Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,972,184 C→T 6.4% intergenic (+62/+57) PP_1762 → / ← yciK hypothetical protein/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,972,1840CT 6.4% 639.7 / 21.1 209intergenic (+62/+57)PP_1762/yciKhypothetical protein/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (95/98);  new base T (3/10);  total (100/109)
Fisher's exact test for biased strand distribution p-value = 8.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 211 total aligned reads displayed.

GACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTT  >  NC_002947/1972036‑1972333
                                                                                                                                                    |                                                                                                                                                     
gACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCa                                                                                                                                                      >  2:6783063/1‑150 (MQ=255)
 aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                     <  1:1346863/150‑1 (MQ=255)
 aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                     <  2:2354199/150‑1 (MQ=255)
 aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                     <  1:6539352/150‑1 (MQ=255)
 aCTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCaa                                                                                                                                                     <  1:5807391/150‑1 (MQ=255)
   tGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCACCCACAGGATCACCGCCACCCACAAgg                                                                                                                                                   >  2:3733056/1‑150 (MQ=255)
           cctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGcc                                                                                                                                           <  2:5687432/150‑1 (MQ=255)
            gaagagatcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                          <  2:914779‑M1/61‑1 (MQ=255)
            ctgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGccc                                                                                                                                          <  2:6624855/150‑1 (MQ=255)
             tgttcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCCCGCCCCCCCGAACACCCCAACCCACCACGATTGCGcccc                                                                                                                                         >  1:5142152/1‑149 (MQ=255)
               gtcctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgc                                                                                                                                       <  2:4105993/149‑1 (MQ=255)
                 cctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgcct                                                                                                                                     <  2:3115867/150‑1 (MQ=255)
                 cctgGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgcct                                                                                                                                     <  2:88379/150‑1 (MQ=255)
                  atcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCcctgcctg                                                                                                                                    >  1:5479448‑M1/84‑150 (MQ=255)
                   tcttgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggctggcgttcaacgcccaataaacctgcttCGTCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                                                   >  1:4992050‑M1/83‑150 (MQ=255)
                   tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                                                   >  2:3845789‑M1/83‑150 (MQ=255)
                   tcatgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgt                                                                                                                                   >  1:614395‑M1/83‑150 (MQ=255)
                    ggAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                                                  <  2:5661805/150‑1 (MQ=255)
                    ggAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                                                  <  1:4267518/150‑1 (MQ=255)
                    ggAGCGCTGCATGCAGCTGGAAGGGCGCCAACAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                                                  <  2:2577301/150‑1 (MQ=255)
                    catgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCtgtg                                                                                                                                  >  1:5547137‑M1/82‑150 (MQ=255)
                     gAGCGCTGCATGCAGCTGGAAGGGCGCCATCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                                 <  1:5507682/150‑1 (MQ=255)
                     atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGTGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                                 >  1:711173‑M1/81‑150 (MQ=255)
                     atgccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTgg                                                                                                                                 >  2:2712571‑M1/81‑150 (MQ=255)
                      aGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTggg                                                                                                                                <  1:2300466/150‑1 (MQ=255)
                       gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                               <  1:6017149/150‑1 (MQ=255)
                       gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                               <  2:3209552/150‑1 (MQ=255)
                       gcgcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                               <  2:1258774/150‑1 (MQ=255)
                       gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                               >  1:2834479‑M1/79‑150 (MQ=255)
                       gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                               >  1:1052230‑M1/79‑150 (MQ=255)
                       gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                               >  2:5141857‑M1/79‑150 (MQ=255)
                       gccggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGa                                                                                                                               >  2:4483373‑M1/79‑150 (MQ=255)
                         gcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                                                             <  2:3742970/150‑1 (MQ=255)
                         gcTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGc                                                                                                                             <  1:2460108/150‑1 (MQ=255)
                         cggtgtacctgtacctggtggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTTCGCCCTGCCTGTGGGAGc                                                                                                                             >  1:3453500‑M1/77‑150 (MQ=255)
                          ggtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgcatCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCg                                                                                                                            >  2:1438299‑M1/76‑150 (MQ=255)
                          cTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCg                                                                                                                            <  1:351825/150‑1 (MQ=255)
                           tgcatgcaGCTGGATGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCTCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCgg                                                                                                                           <  2:1021685/150‑1 (MQ=17)
                           tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGGGTGCCCCCCCCCCCACGATCCCCCCCACCCACAAAGCTTTGGCCCAGCCGGGGGGTGCgg                                                                                                                           >  1:5387892/1‑150 (MQ=255)
                           tgcatgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                                           <  2:5449462/150‑1 (MQ=255)
                           gtgtacctgtacctgatggggccggacagcacaggggtgaacgggcaggcgttcaacgcccaataaacctgtacCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCgg                                                                                                                           >  1:1857941‑M1/75‑150 (MQ=255)
                               tgcaGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCaa                                                                                                                       <  2:6167459/150‑1 (MQ=255)
                                  aGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCGTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGCTTGCGCCCTGCCTGGGGGAGCGGGCCAGcc                                                                                                                    >  2:4694314/1‑150 (MQ=255)
                                     tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGTGCGGGCAAGCCcgc                                                                                                                 <  2:4831363/150‑1 (MQ=14)
                                     tGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGTTTGCGCCCTGCCTGTGGGAGCGGGCAAGCCcgc                                                                                                                 <  2:541871/150‑1 (MQ=255)
                                       gAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCGGTGGGAGCGGGCAAGCCCGCGa                                                                                                               >  2:5505917/1‑150 (MQ=255)
                                                    aaGGCCTTATGCGGCCGGCGTTGGCTGCGGGCTGATCGATTTCCTGTACCGGCCTCTGCGCGGGCTGGCCGGCTCCCACAGGATCACCGCAACCCTCACGGATTGCGCCCTGCCTGTGGGCGCGGGCAAGCCCGCGAAGAGGCCGATGCa                                                                                                  <  2:3305760/150‑1 (MQ=38)
                                                           tATGCGGCGGGCTTTGGCTGGGGGCTTAGCGCTTTCGTGTCCCGCCCTCTTCGGGGGCTTGCCCGCACCCACTGGTCCACCCCAACCCAAATGGTTTGCGCCCTCCCTGTGGGCGGGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTAt                                                                                           <  2:213067/150‑1 (MQ=255)
                                                           caggcggcgGGCGTGGGCTGCGGGCTGATCGCTGGCCTGTACCGGCCTCTCCGCGGGCTTGCCCGCCCCCACCGGATCACCGCACCCCTCACGGATTGCGCCCGGCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTAt                                                                                           <  2:6585890/147‑1 (MQ=37)
                                                                  gcgggGTGGGCGGGGGGTTTATCGCTTTCGTTCCCCGGCCTCGGCGGGGGCGGGCCCGCCCCCCCCGTATCCCCCAAACCCTCTGGGTTTGCGCCCGGCGGGTGGGCGGGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGt                                                                                    <  2:2852241/149‑1 (MQ=255)
                                                                             ggcgggcggAGCGCGGTCCTGTACCGGCCTCCGCGGGGGCTGCCCCGCCCCCACCGGACCCCCGCAACCCTCAGGGAGTGCGCCCTGCCTGTGGGCGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGc                                                                         <  2:5190002/149‑1 (MQ=255)
                                                                                       tCGCTTTCCTGTACCGGCCTCTTCGCGGGCGTGCCCCCGCCCACACGCTCACCCCAACCCACCAAGCTTTCCCCCCGCCCGTGGGGGCCGGCCAGCCCCCCCAAAGGCCCAGGCCGGTGTATTTGGCGGTGTACCCCAGCCCGGGCAccc                                                               >  2:4050535/1‑150 (MQ=255)
                                                                                         gCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACCGGATCACCGCAACCCTCAAGGATTGCGCCCGGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCt                                                             <  2:6237079/150‑1 (MQ=255)
                                                                                          cTTTCCTGTACCGGCCTCTTCGCGGGGGTGGCCGCTGCCAAAGGTTTACAGCAAACCAGAAGGATTTGGCCCTTGCGGTTGGTGGGGGGAAGCCCGCGTAGAGGGCGATGGAGGTTTTTTTGGGGTTTTAAGGCTGGCCGTTAAacccgg                                                            >  2:2932904/1‑148 (MQ=255)
                                                                                                               GCGGGCAAGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggt                                       <  1:5970778‑M1/150‑1 (MQ=255)
                                                                                                                   GCTTGCCCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag                                   <  2:2079125‑M1/150‑1 (MQ=255)
                                                                                                                   GCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacag                                   <  2:2471619‑M1/150‑1 (MQ=255)
                                                                                                                    CTTGCCCGCTCCCACCGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCGGTGGGAGCGGGCCAGCCCGCGAAGAGGCCGATGCAGGTGGATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacagg                                  <  2:5090110‑M1/150‑1 (MQ=255)
                                                                                                                     TTGCCCGCTCCCCCCGGATCACCCAACCCCGCACGGATTGCGCCCTGCGTGTGGGCGCGGGAAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                                 >  1:4566917‑M1/1‑150 (MQ=255)
                                                                                                                     TTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATcaggtacaggt                                 <  2:2607427‑M1/150‑1 (MQ=255)
                                                                                                                       GCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAATCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCaggtacaggtac                               >  2:1604682‑M1/1‑150 (MQ=255)
                                                                                                                       GCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCaggtacaggtac                               <  1:952985‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                             <  1:4548312‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGCAGGGCGCAATCCTTGAGGGTTGCGGTGATCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                             <  2:653562‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                             <  1:6138373‑M1/150‑1 (MQ=255)
                                                                                                                         CCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTacac                             <  1:1334719‑M1/150‑1 (MQ=255)
                                                                                                                          CGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACAcc                            <  1:4171083‑M1/150‑1 (MQ=255)
                                                                                                                           GCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCg                           <  1:6767599‑M1/150‑1 (MQ=255)
                                                                                                                           GCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCg                           <  2:5356895‑M1/150‑1 (MQ=255)
                                                                                                                             TCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGc                         <  1:327765‑M1/150‑1 (MQ=255)
                                                                                                                              cccccctGATCACCCCAACCCTCACGTTTCGCCCCCGGCCTGTGGGCGCGGGCACGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCa                        >  1:2083006‑M1/8‑150 (MQ=255)
                                                                                                                               CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt                       <  1:1620449‑M1/150‑1 (MQ=255)
                                                                                                                               CCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCAt                       <  2:556590‑M1/150‑1 (MQ=255)
                                                                                                                                CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg                      <  1:4970523‑M1/150‑1 (MQ=255)
                                                                                                                                CACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATg                      <  1:3778188‑M1/150‑1 (MQ=255)
                                                                                                                                    GGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtct                  <  2:4630003‑M1/150‑1 (MQ=255)
                                                                                                                                    GGATCACCGCAACCCCCCAGGAGGGCGCCCTGCCTGGGGGTGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtct                  >  2:5972797‑M1/1‑150 (MQ=255)
                                                                                                                                     GATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGAtctc                 <  2:1680774‑M1/150‑1 (MQ=255)
                                                                                                                                        gggcGCAATCCTTGAGGGTTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTc              >  2:2182506/4‑150 (MQ=38)
                                                                                                                                         aCCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCg             >  1:3257617/1‑150 (MQ=255)
                                                                                                                                           cGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCggg           >  2:251942/1‑150 (MQ=255)
                                                                                                                                             cAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCa         >  2:3595314/1‑150 (MQ=255)
                                                                                                                                                 ccTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCgg     >  2:2403475/1‑150 (MQ=255)
                                                                                                                                                 ccTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCgg     >  1:3608567/1‑150 (MQ=255)
                                                                                                                                                  cTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg    >  1:3990504/1‑150 (MQ=255)
                                                                                                                                                  cTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg    >  1:857613/1‑150 (MQ=255)
                                                                                                                                                  cTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCggg    >  1:6571983/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:6388729/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  1:6624872/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:1745208/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  1:2513261/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:5390978/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:510607/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:4968285/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:1104616/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:4732160/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  2:3355746/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgt   >  1:6294309/1‑150 (MQ=255)
                                                                                                                                                   tCAAGGATTGCGCCCAGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTAATTGGGCGTTGAACGCCTCCCCGTTCACCCCTGAGCTGTCCGCCCCCATCACGTACAGGAACACCGGCAACATCTCTTCGGGCAGCggg  >  2:2324770/1‑149 (MQ=255)
                                                                                                                                                    cAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGgtt  >  2:1857297/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACTGCATGTTCCTGTTCCTGGAGCGCTGCATGCAGCTGGAAGGGCGCCAGCAAGGCCGTATGCGGCCGGCGTTGGCTGCGGGCTGATCGCTTTCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGATCACCGCAACCCTCAAGGATTGCGCCCTGCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCGATGCAGGTTTATTGGGCGTTGAACGCCTGCCCGTTCACCCCTGTGCTGTCCGGCCCCATCAGGTACAGGTACACCGGCATGATCTCTTCGGGCAGCGGGTT  >  NC_002947/1972036‑1972333

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