Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,866,193 T→C 5.6% intergenic (‑202/+28) mtlD ← / ← mtlA mannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding/fused mannitol‑specific PTS enzyme IIABC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,866,1930TC 5.6% 254.6 / 5.8 89intergenic (‑202/+28)mtlD/mtlAmannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding/fused mannitol‑specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base T (39/45);  new base C (3/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCAAATTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTG  >  AP009048/3866123‑3866263
                                                                      |                                                                      
gCAAATTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCt                                                                        >  1:4660938/1‑71 (MQ=255)
gCAAATTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCt                                                                        >  1:633914/1‑71 (MQ=255)
  aaaTTCAACAAATTGTGATGCATCCTTTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtc                                                                      <  1:1771011/71‑1 (MQ=255)
  aaaTTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtc                                                                      >  1:1315432/1‑71 (MQ=255)
  aaaTTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCCtc                                                                      >  1:3381390/1‑71 (MQ=255)
   aaTTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtct                                                                     >  1:3831922/1‑71 (MQ=255)
   aaTTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtct                                                                     >  1:3277383/1‑71 (MQ=255)
    aTTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtctc                                                                    >  1:797523/1‑71 (MQ=255)
    aTTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtct                                                                     >  1:4457822/1‑70 (MQ=255)
     ttCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtctc                                                                    >  1:3778873/1‑70 (MQ=255)
        aacaaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAt                                                                 <  1:3588521/70‑1 (MQ=255)
         acaaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCtt                                                                       <  1:4868854/63‑1 (MQ=255)
          caaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                              <  1:1200245/71‑1 (MQ=255)
          caaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCACACCTTCTCCAtgtg                                                              >  1:1312947/1‑71 (MQ=255)
             aTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTtct                                                                     >  1:1122711/1‑61 (MQ=255)
             aTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                              >  1:838799/1‑68 (MQ=255)
             aTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                              >  1:4292489/1‑68 (MQ=255)
              ttGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGgaga                                                          <  1:3645101/71‑1 (MQ=255)
                gtgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                         <  1:1332377/70‑1 (MQ=255)
                 tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                       <  1:2966907/71‑1 (MQ=255)
                 tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                       >  1:3778970/1‑71 (MQ=255)
                 tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAgg                                                        <  1:2653601/70‑1 (MQ=255)
                 tgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAgg                                                        <  1:1952358/70‑1 (MQ=255)
                   aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                      <  1:2345324/70‑1 (MQ=255)
                   aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                      <  1:5082437/70‑1 (MQ=255)
                   aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTg                                                     >  1:2532348/1‑71 (MQ=255)
                      cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTg                                                     >  1:1215358/1‑68 (MQ=255)
                      cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGgga                                                  <  1:3630135/71‑1 (MQ=255)
                        tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgt                                                               <  1:1160235/56‑1 (MQ=255)
                        tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggatt                                                <  1:2979297/71‑1 (MQ=255)
                          cTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                               <  1:144184/70‑1 (MQ=255)
                           tGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                               <  1:4112772/69‑1 (MQ=255)
                                gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTgg                                                             <  1:3275491/50‑1 (MQ=255)
                                gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                        <  1:4312243/71‑1 (MQ=255)
                                gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                        <  1:3174406/71‑1 (MQ=255)
                                gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                        <  1:5236795/71‑1 (MQ=255)
                                 ccTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                        <  1:827125/70‑1 (MQ=255)
                                   tGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                         <  1:3803544/51‑1 (MQ=255)
                                   tGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttacttac                                     >  1:680712/1‑71 (MQ=255)
                                   tGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAtta                                          >  1:2467193/1‑66 (MQ=255)
                                    gATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTTCATGTGGAGAGGGTGGGATTGGAttacttac                                     >  1:4207061/1‑70 (MQ=255)
                                    gATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAt                                                                 <  1:973503/42‑1 (MQ=255)
                                    gATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGgaga                                                          >  1:1454356/1‑49 (MQ=255)
                                     aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACg                                    <  1:3487071/70‑1 (MQ=255)
                                     aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                   >  1:1068446/1‑71 (MQ=255)
                                     aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                   >  1:2330691/1‑71 (MQ=255)
                                      tAAGCGTAGCGCATCAGGCAATTAACCACACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGAc                                  >  1:4852193/1‑71 (MQ=255)
                                        aGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGgaga                                                          >  1:5347440/1‑45 (MQ=255)
                                        aGCGTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGgagaa                                                         <  1:1879473/46‑2 (MQ=21)
                                         gCGTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGGAGAagg                                                       <  1:347914/47‑1 (MQ=18)
                                          cGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATTTGGAGAGGGTGggattgga                                             >  1:3617350/1‑56 (MQ=255)
                                          cGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGc                              <  1:4021157/71‑1 (MQ=255)
                                           gTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                        <  1:2238873/60‑1 (MQ=255)
                                           gTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGGAGAAGGTg                                                     >  1:2077999/1‑47 (MQ=17)
                                              gcgcATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                               >  1:1271706/1‑50 (MQ=255)
                                              gcgcATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcag                           >  1:4445413/1‑70 (MQ=255)
                                               cgcATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcagc                          <  1:3026443/70‑1 (MQ=255)
                                                 cATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                       >  1:4387983/1‑71 (MQ=255)
                                                 cATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                       >  1:216850/1‑71 (MQ=255)
                                                  aTCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                      <  1:2471884/39‑1 (MQ=255)
                                                  aTCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGtt                      >  1:4544558/1‑71 (MQ=255)
                                                  aTCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                       <  1:375539/70‑1 (MQ=255)
                                                     aGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTcc                    <  1:3100337/70‑1 (MQ=255)
                                                     aGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTc                     >  1:2038205/1‑69 (MQ=255)
                                                      ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcagca                         >  1:4469543/1‑64 (MQ=255)
                                                      ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                  <  1:1964857/71‑1 (MQ=255)
                                                       gCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGATCTGCCAGCAGTTCCAg                  >  1:2510022/1‑70 (MQ=255)
                                                       gCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                  >  1:4504366/1‑70 (MQ=255)
                                                        cAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                <  1:2791445/71‑1 (MQ=255)
                                                         aaTTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcagc                          <  1:101380/60‑1 (MQ=255)
                                                            tAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttactta                                      >  1:1445527/1‑45 (MQ=255)
                                                            tAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACtt             <  1:4414191/70‑1 (MQ=255)
                                                             aaCCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCa           >  1:141234/1‑71 (MQ=255)
                                                              aCCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCa           <  1:591826/70‑1 (MQ=255)
                                                               tcccACCCTCTCCACATGGAGAAGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc         >  1:1312819/2‑71 (MQ=25)
                                                               ccccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc         <  1:2218806/71‑1 (MQ=255)
                                                               ccccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc         >  1:2937224/1‑71 (MQ=255)
                                                                cccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATc         >  1:3572972/1‑70 (MQ=255)
                                                                 ccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGt                       <  1:4366777/55‑1 (MQ=255)
                                                                 ccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTc            <  1:1991114/66‑1 (MQ=255)
                                                                  cACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCAc      <  1:3502840/71‑1 (MQ=255)
                                                                   aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCAc      >  1:3682374/1‑70 (MQ=255)
                                                                   aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACg     >  1:1283496/1‑71 (MQ=255)
                                                                    ccTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCcag                           <  1:2719237/48‑1 (MQ=255)
                                                                    ccTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTcc                    >  1:1005395/1‑55 (MQ=255)
                                                                    ccTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGc    <  1:1727949/71‑1 (MQ=255)
                                                                     cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGc                 <  1:2420796/57‑1 (MQ=255)
                                                                     cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACt              <  1:1174665/60‑1 (MQ=255)
                                                                     cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt   <  1:4943847/71‑1 (MQ=255)
                                                                     cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt   <  1:4450701/71‑1 (MQ=255)
                                                                     cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt   <  1:1681338/71‑1 (MQ=255)
                                                                     cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt   >  1:1874980/1‑71 (MQ=255)
                                                                      ttCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGAc                                  <  1:1412/39‑1 (MQ=255)
                                                                      ttCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtg  >  1:4335628/1‑71 (MQ=255)
                                                                      |                                                                      
GCAAATTCAACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTG  >  AP009048/3866123‑3866263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: