Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,520,695 C→G 5.1% A231A (GCC→GCG)  PP_4853 → DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,520,6950CG 5.1% 283.5 / 2.9 90A231A (GCC→GCG) PP_4853DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base C (46/39);  new base G (3/2);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTACCTGCAGGTGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTAGTGCTGGC  >  NC_002947/5520547‑5520844
                                                                                                                                                    |                                                                                                                                                     
tgTACCTGCAGGTGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCg                                                                                                                                                      >  2:2352738/1‑150 (MQ=255)
      tGCAGGTGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCAcc                                                                                                                                                >  1:679030/1‑150 (MQ=255)
         aGGTGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGcc                                                                                                                                             >  1:275886/1‑150 (MQ=255)
         aGGTGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGcc                                                                                                                                             <  2:1373759/150‑1 (MQ=255)
           gTGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTg                                                                                                                                           >  2:1506310/1‑150 (MQ=255)
            tGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTgg                                                                                                                                          <  1:2193275/150‑1 (MQ=255)
             gcgggtggCTATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTggg                                                                                                                                         <  1:683223/147‑1 (MQ=255)
             gCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGGTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGTTCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTggg                                                                                                                                         <  2:1054165/150‑1 (MQ=255)
             gCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGTCGCTGGTGTGGATGCAGGCCGTGCACCGCCTggg                                                                                                                                         <  1:1067973/150‑1 (MQ=255)
             gCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTggg                                                                                                                                         <  1:2352738/150‑1 (MQ=255)
               tggtggCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGc                                                                                                                                       >  1:1483342/1‑150 (MQ=255)
                      cATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCg                                                                                                                                >  2:2338621/1‑150 (MQ=255)
                         cgtcgtGTTGCTGCCGCTGTTCCTGCTGTCAGGCAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCGGTGCACCGCCTGGGGCCCAGCCGcat                                                                                                                             <  2:2084301/150‑2 (MQ=255)
                               gtTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACt                                                                                                                       <  2:1329800/150‑1 (MQ=255)
                                ttGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTg                                                                                                                      >  1:906390/1‑150 (MQ=255)
                                  gctgcCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTGCTGTATTCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGtt                                                                                                                    <  2:74553/150‑1 (MQ=255)
                                   ctgcCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGttc                                                                                                                   >  1:1343058/1‑150 (MQ=255)
                                        gCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCaa                                                                                                              <  2:1483342/150‑1 (MQ=255)
                                               ctgctgTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTtgctg                                                                                                       <  1:884583/150‑1 (MQ=255)
                                                tgctgTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTtgctgc                                                                                                      <  1:1334399/150‑1 (MQ=255)
                                                     tCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTg                                                                                                 >  1:1839209/1‑150 (MQ=255)
                                                       aGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGATGGTGGGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGt                                                                                               >  1:1620136/1‑150 (MQ=255)
                                                       aGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCGGTATGCGTGGGTACTGGCGGCGAGGATCGCGCCGCTGGTGTGGATGCAGGCCGTCCACCGTCGGGGGCCCAGCCGAACCACACGTTTCTTCAATTTGCTTCCGCTggg                                                                                               >  2:1004549/1‑149 (MQ=255)
                                                         aCAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTg                                                                                                 >  2:919781/1‑146 (MQ=255)
                                                         aCAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTg                                                                                                 <  1:919781/146‑1 (MQ=255)
                                                           aaGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCAcc                                                                                           >  1:1841998/1‑150 (MQ=255)
                                                             gACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACcgc                                                                                         >  1:449193/1‑150 (MQ=255)
                                                                 ggCCTGAACAGCCACAACATCGGCCTGTTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTgg                                                                                                <  1:1934861/139‑1 (MQ=255)
                                                                 ggCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTgg                                                                                                >  2:1934861/1‑139 (MQ=255)
                                                                  gCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACcaca                                                                                                                         <  1:2271418/113‑1 (MQ=255)
                                                                  gCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACcaca                                                                                                                         >  2:2271418/1‑113 (MQ=255)
                                                                      gAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGAGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATcgc                                                                                >  2:2392922/1‑150 (MQ=255)
                                                                        aCAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGc                                                                                                                                                   <  1:1550668/81‑1 (MQ=255)
                                                                        aCAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGc                                                                                                                                                   >  2:1550668/1‑81 (MQ=255)
                                                                           gCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGcac                                                                                                                             <  2:1222960/100‑1 (MQ=255)
                                                                           gCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGcac                                                                                                                             >  1:1222960/1‑100 (MQ=255)
                                                                             cacaACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGc                                                                                                   >  2:2369843/1‑124 (MQ=255)
                                                                             cacaACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGc                                                                                                   <  1:2369843/124‑1 (MQ=255)
                                                                             cacaACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCgtg                                                                         >  2:862395/1‑150 (MQ=255)
                                                                              acaacaTCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGcc                                                                                                                                             <  2:2119624/81‑1 (MQ=255)
                                                                              acaacaTCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGcc                                                                                                                                             >  1:2119624/1‑81 (MQ=255)
                                                                               caacaTCGGCCTGGTGCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCgtggt                                                                       <  2:1676480/150‑1 (MQ=255)
                                                                               caacaTCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCgtggt                                                                       >  1:217516/1‑150 (MQ=255)
                                                                                aacaTCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCgtggtg                                                                      <  2:2013481/150‑1 (MQ=255)
                                                                                 acaTCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGc                                                                     <  1:862395/150‑1 (MQ=255)
                                                                                   aTCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTc                                                                   >  2:1308495/1‑150 (MQ=255)
                                                                                       gCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGTCGTGGTGCTcgacg                                                               >  2:1484788/1‑150 (MQ=255)
                                                                                       gCCTGGTCCTGTATGCCTGCGTACTGGCTTCGAGGATCGCGCCGCTGGTGTGGATGCAGTCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGCTCGCCGCCGTGGTGCTcgccg                                                               >  2:490914/1‑150 (MQ=255)
                                                                                         cTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGTCGTGGTGCTCGACGaa                                                             <  1:1484788/150‑1 (MQ=255)
                                                                                         cTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGaa                                                             <  2:449193/150‑1 (MQ=255)
                                                                                         cTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCAACGCGCTGATCGCCGCCTTGGTGCTCTTCGaa                                                             <  2:1354244/150‑1 (MQ=255)
                                                                                          tGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGaac                                                            >  2:2168208/1‑150 (MQ=255)
                                                                                             tCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGaacaac                                                         >  1:1952846/1‑150 (MQ=255)
                                                                                             tCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGaacaac                                                         <  2:966055/150‑1 (MQ=255)
                                                                                              ccTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCTCCACTTTGTTCTTCAACCTGCTGTCGCTTGTCACCGCGCTGGTCGCCGCCGTGGTGCTCGACGAACAACt                                                        <  2:275886/150‑1 (MQ=255)
                                                                                              ccTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACt                                                        <  1:2066860/150‑1 (MQ=255)
                                                                                                tatATGCCTGCGTACTGGCTTCGATGATAGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTATTCAACTTGCTGCCGCTGGTCAAAGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTgg                                                      >  2:1698743/3‑150 (MQ=255)
                                                                                                tGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCgt                                                                          <  2:168199/130‑1 (MQ=255)
                                                                                                tGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCgt                                                                          >  1:168199/1‑130 (MQ=255)
                                                                                                    tGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGaacaac                                                         <  1:754079/143‑1 (MQ=255)
                                                                                                    tGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGaacaac                                                         >  2:754079/1‑143 (MQ=255)
                                                                                                      ccTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGtt                                                >  1:1951778/1‑150 (MQ=255)
                                                                                                        tGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGACCAACTGGCCAGTTac                                              >  2:405578/1‑150 (MQ=255)
                                                                                                           ggaCTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTTGCCAGTTaccac                                           <  2:549186/148‑1 (MQ=255)
                                                                                                              cTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTg                                        <  1:2392922/150‑1 (MQ=255)
                                                                                                               tGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTc                                                                   <  1:258289/122‑1 (MQ=255)
                                                                                                               tGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTc                                                                   >  2:258289/1‑122 (MQ=255)
                                                                                                               tGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTgg                                       >  2:433297/1‑150 (MQ=255)
                                                                                                                   ttCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTcgg                                   >  1:342015/1‑150 (MQ=255)
                                                                                                                   ttCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTcgg                                   >  2:1134400/1‑150 (MQ=255)
                                                                                                                   ttCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTcgg                                   >  1:1639009/1‑150 (MQ=255)
                                                                                                                      gatgatCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCGCCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATcgcc                                                                               >  2:2178274/1‑103 (MQ=255)
                                                                                                                      gatgatCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCGCCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATcgcc                                                                               <  1:2178274/103‑1 (MQ=255)
                                                                                                                       atgatCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCggcggc                               >  1:1532708/1‑150 (MQ=255)
                                                                                                                         gatCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCtt                             >  2:1858208/1‑150 (MQ=255)
                                                                                                                         gatCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCtt                             <  2:2085363/150‑1 (MQ=255)
                                                                                                                          atCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTg                            <  1:2168208/150‑1 (MQ=255)
                                                                                                                           tCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGc                           >  2:1725625/1‑150 (MQ=255)
                                                                                                                            cgcgCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCt                          <  2:1951778/150‑1 (MQ=255)
                                                                                                                             gcgcTGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTg                         <  1:484711/150‑1 (MQ=255)
                                                                                                                                 cgcTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCt                     <  1:1725625/150‑1 (MQ=255)
                                                                                                                                  gcTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCtt                    >  2:1162264/1‑150 (MQ=255)
                                                                                                                                           gATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCTTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGt           >  1:1668569/1‑150 (MQ=255)
                                                                                                                                            aTGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTa          >  1:664737/1‑150 (MQ=255)
                                                                                                                                             cccAGGCGGTGCACGGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTAg         >  2:1532373/3‑150 (MQ=255)
                                                                                                                                              ccAGGCGGTGCACGGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTacc                                             <  2:702642/112‑1 (MQ=255)
                                                                                                                                              ccAGGCGGTGCACGGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTacc                                             >  1:702642/2‑113 (MQ=255)
                                                                                                                                              ccAGGCGGTGCACGGCCTGGGGCCCAGCCGCACCACACTGTTCTTAAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTCACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTAGt        >  2:499810/2‑150 (MQ=255)
                                                                                                                                                aGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTAGTGc      >  1:58354/1‑150 (MQ=255)
                                                                                                                                                 ggCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTAGTGCt     <  2:1841998/150‑1 (MQ=255)
                                                                                                                                                    cGTGCACCGCCTTGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTAGTGCTGGc  <  1:1308495/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTACCTGCAGGTGCTGGTGGCCATCGTCGTGTTGCTGCCGCTGTTCCTGCTGTCAGACAAGACCGGCCTGAACAGCCACAACATCGGCCTGGTCCTGTATGCCTGCGTACTGGCTTCGATGATCGCGCCGCTGGTGTGGATGCAGGCCGTGCACCGCCTGGGGCCCAGCCGCACCACACTGTTCTTCAACTTGCTGCCGCTGGTCACCGCGCTGATCGCCGCCGTGGTGCTCGACGAACAACTGGCCAGTTACCACCTGGTCGGCGGCTTGCTGACCTTGGCCGGCGTAGTGCTGGC  >  NC_002947/5520547‑5520844

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑