Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,344,156 A→C 5.3% intergenic (‑166/+77) PP_2935 ← / ← ybiT membrane protein/ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,344,1560AC 5.3% 711.4 / 4.1 237intergenic (‑166/+77)PP_2935/ybiTmembrane protein/ABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base A (101/118);  new base C (3/15);  total (104/133)
Fisher's exact test for biased strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 241 total aligned reads displayed.

GGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGAC  >  NC_002947/3344007‑3344305
                                                                                                                                                     |                                                                                                                                                     
gggcggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATa                                                                                                                                                       <  1:5785774/150‑1 (MQ=255)
   cggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATATACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCt                                                                                                                                                    <  1:1076260/150‑1 (MQ=255)
   cggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGCCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCt                                                                                                                                                    <  2:2135811/150‑1 (MQ=255)
        gggCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCg                                                                                                                                               >  2:2750868/1‑150 (MQ=255)
        gggCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCg                                                                                                                                               <  1:2916890/150‑1 (MQ=255)
          gCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCCCGCACCCAAGGCCCCACAAACCGCCAAAACCAGTGCggg                                                                                                                                             >  1:1716508/1‑149 (MQ=255)
              ggTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCt                                                                                                                                         <  1:2047212/150‑1 (MQ=255)
               ttcagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTC‑‑CCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgtg                                                                                                                                      >  1:4908197‑M1/88‑150 (MQ=255)
                gcccATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTACCGAATGATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGAGGTACTCGTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGGCCCtgt                                                                                                                                       <  2:3602727/149‑1 (MQ=255)
                 cccATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgtg                                                                                                                                      <  2:2211775/150‑1 (MQ=255)
                 cagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgtg                                                                                                                                      >  1:358274‑M1/86‑150 (MQ=255)
                  ccATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTgg                                                                                                                                     <  1:143783/150‑1 (MQ=255)
                  ccATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTgg                                                                                                                                     <  1:3443053/150‑1 (MQ=255)
                  agcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTgg                                                                                                                                     >  1:4327691‑M1/85‑150 (MQ=255)
                    atcatGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAAAACCTGTGCGGTCCCTGTGGGa                                                                                                                                   >  2:4183299/1‑150 (MQ=255)
                     tcatGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCCAGGCCCCACAAGCGTCCAAAACCTGTGCGGTCCCGGTGGGAg                                                                                                                                  >  1:6148537/1‑150 (MQ=255)
                       atGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCg                                                                                                                                <  2:5781331/150‑1 (MQ=255)
                        tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                                                               <  2:853429/150‑1 (MQ=255)
                        tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                                                               <  2:6853434/150‑1 (MQ=255)
                        tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCCCGCACCCAAGGCCCCCCAACCCACCACAACCAGATCGGTCCCGGGGGGAGCgg                                                                                                                               >  1:2881079/1‑150 (MQ=255)
                          gCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGc                                                                                                                             <  1:468137/150‑1 (MQ=255)
                           ccGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCCCTCTCCCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCTGTGGGAGCGGGcc                                                                                                                            >  2:6383806/1‑149 (MQ=255)
                            cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                                                           <  1:4225500/150‑1 (MQ=255)
                            cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                                                           <  2:5272521/150‑1 (MQ=255)
                            acgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                                                           >  1:383703‑M1/75‑150 (MQ=255)
                            acgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                                                           >  2:5095458‑M1/75‑150 (MQ=255)
                                 cgcgAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCCCGCGCCCAAGGCCCCCCCACCGGCCCATACCAGTGCCGGCCCCCTGGGGGCGGGCCACCCcg                                                                                                                      >  1:5409973/1‑150 (MQ=255)
                                  gcgAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCACGCCcgc                                                                                                                     >  2:4071011/1‑150 (MQ=255)
                                   cgAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCcgcg                                                                                                                    <  1:2642926/150‑1 (MQ=255)
                                   cgAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCcgcg                                                                                                                    <  1:6109081/150‑1 (MQ=255)
                                    gAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCCAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGa                                                                                                                   >  1:6729778/1‑150 (MQ=255)
                                     aaGGAAACGACAGCAGGCTAACGACTTAGTCGTCGGGCCGTGCCGGTCGGTGGGAGGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGaa                                                                                                                  <  2:781119/150‑1 (MQ=255)
                                       ggAACAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGGCCCTGTGGGAGCGGGCAAGCCCGCGAaga                                                                                                                <  2:4425294/150‑1 (MQ=255)
                                        taaaaGAAAGCAGGCTAACTAATTTTTCGTCGGGCCGTGCGGGTGGGTGGTTTGGGTATTTACTCTTCGCGGGCTTGCCCGCTATCCCAAGGCCCCAAAAGGGTCCAATACGTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAagag                                                                                                               <  1:5995282/149‑1 (MQ=255)
                                        gAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAAAACCAGTGCGGTACCTGTGGGAGCGGGCAAGCCCGCGAagag                                                                                                               >  2:2934687/1‑150 (MQ=255)
                                            aGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGcct                                                                                                           <  2:6277734/150‑2 (MQ=255)
                                             gACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCGTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACACGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCg                                                                                                          <  2:5416101/150‑1 (MQ=255)
                                                agcagGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCGCGCCCAAGGCCCCACAAGCGTCCAAAACCAGTGCGGTCCCGGTGGGAGCGGGCCACCCCGCCAACAGGCCCCGcc                                                                                                       >  2:12574/1‑150 (MQ=255)
                                                agcagGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCCCGCTCCCAACGCCCCCCAAACGTCCCAAAACAGGGCGGTCCCGGGGGGAGCGGGCAAGCCCCCGAAGAGGCCCGGcc                                                                                                       >  1:5818570/1‑150 (MQ=255)
                                                       tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTACCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCAGCCcag                                                                                                >  2:2613066/1‑150 (MQ=255)
                                                        accGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGAGGTCCTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACACGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCgaga                                                                                               <  2:2673188/148‑1 (MQ=255)
                                                         aCGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGAt                                                                                              >  2:6168793/1‑150 (MQ=255)
                                                         aCGAATTATTCGGCGGGCCGTGCCGGTCGGTGGGCTGGGTATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATCCCTGTGGGGTCCCTGTGGGCGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGAt                                                                                              <  1:526980/150‑1 (MQ=255)
                                                          ctaaTTATTCGTCGGGCCGTGCCGGTCGGTGGTATGGGTATGTACTCTTCGCGGGCTTGCCCGCTCTCCAAAGGCCCCACAAGCTTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATg                                                                                             <  2:1410187/148‑1 (MQ=255)
                                                            aaTTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGc                                                                                           <  2:1522360/150‑1 (MQ=255)
                                                              ttattCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGACCGCGAAGAGGCCCGGCCATCCGAGATGGCCa                                                                                         >  2:6807171/1‑150 (MQ=255)
                                                                attCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGCGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAgg                                                                                       >  2:6001269/1‑150 (MQ=255)
                                                                attCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAgg                                                                                       <  2:4494380/150‑1 (MQ=255)
                                                                 ttCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATCCCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATCTCGGATGGCCGGGc                                                                                      <  2:1242704/150‑1 (MQ=38)
                                                                    gtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTTCCCGCGCTCCCAAGGCCCCACAAGCGGCCAAAACCAGTGCGGTGCCGGTGGGAGCCGGCCAGCCCGCCAACAGGCCCCGCCCACCCACGTGGCCCGGcccg                                                                                   >  1:4766186/1‑148 (MQ=255)
                                                                    gtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCAACCAAGGCCCCACAAGCGTCCAATACCAGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGt                                                                                   >  2:1969253/1‑150 (MQ=255)
                                                                           ggggCCGGGCGGGGGGAGGGGCTGGAACTCTCCGGGGGCGGCCCCGCGCCCCCACGCCCCCACAAGCTTCCAACCCCTGGGCGGCCCCTGTGGGTGCGGCCACCCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatca                                                                            <  1:3495258/147‑1 (MQ=255)
                                                                              gcgggtgggtgggtgGGGTATGTACTCTTGGGGGGCTTGCCCGCCCTCCAACGGCCCCAAAAGCGTCCAAAACGTGTGGGGACCCTGTGGGCGGGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatcaatc                                                                         <  1:4172295/147‑1 (MQ=255)
                                                                                 ggtcggtGGGAGGGGGATGCACTCTTCGCGGGCTGGCCCGCTCCCCCACGGCCCCACAAGCGTCCACTCCCTGTGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                      <  1:1059645/150‑1 (MQ=255)
                                                                                 gggcggtGGGAGGGGGTTGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGa                                                                      <  2:2543082/147‑1 (MQ=39)
                                                                                  gtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGACCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGAc                                                                     <  1:6807171/150‑1 (MQ=255)
                                                                                     ggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCGCTCCCCAGGCCCCCCACGCGTCCCATACCTGGGCGGGGCCTGTGGGGGCGGGCCAGCCCGCGGAGGGGCCCGGCCCTCCGAGATGGGGCGGCCGTCCCAGCCATCAGAcccg                                                                  >  2:1767791/1‑148 (MQ=255)
                                                                                      gtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCCGTCCCTGTGGGAGCGGGCCAGCCCGCGAAGAGGCCCGGCCAACCGAGATGGCCCGGCCGTCCCATAAATCAGACCAGc                                                                 >  2:4642164/1‑150 (MQ=255)
                                                                                      gtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCACTCCCACGGCCCCACAAGCTTCCACTACCGGGGCGGTCCGTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGc                                                                 <  1:6201292/150‑1 (MQ=255)
                                                                                      ggagcggggggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCAAAAGCGTCCAATCCCCGTGCGGTCCCTGTGGGTGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGc                                                                 <  2:3180135/143‑1 (MQ=255)
                                                                                             ggggtttTACTTTTGGGGGGTTTGCCCGATCCCCCACGCCCCCAAACGGGCCCATTCCGTGGGGGGGCCGTGGGGTCGGGGGCACCCCCGCTAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCt                                                          <  1:313842/150‑1 (MQ=255)
                                                                                                       tctTCGCGGGCTTGCCCGCTCTCCCAACGCCCCACAAGCGTCCAAAAACGGTGCGGTCCCGGTGGGAGCGGGCAAGCCCCCGAAGAGGCCCGGCCAACCGAGATGGCCCGGCCGTCCCCACAATCAGACCAGCACACCCAGGCTACCCAg                                                >  2:1219040/1‑150 (MQ=255)
                                                                                                       TCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATAACAGTGCGGTCCCAGTGGGAGCGGGCCAGCCCGCGAAGAGGCCCGGCCATCCGAGATTGCCAGGCCGGCCCATCAATCAGACCCGCACACCCAGGCTACCCAg                                                <  2:4773462‑M1/150‑1 (MQ=255)
                                                                                                         TTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAATCGTGCAATACCTGTGCGGTGCCTGTGGGTGCGGGCAAGCCCGCGAATATGCCCGGCCATCCGAGATTGCCAGGCCGTCCCATCAATAATAACAGCACACCCTGGCTAAGCAGGt                                              <  2:4215464‑M1/150‑1 (MQ=255)
                                                                                                          ggGCGGGCTGGCCCGCTCCCCCACGCCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                             >  1:3718802‑M1/3‑150 (MQ=255)
                                                                                                          TCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                             >  2:5134143‑M1/1‑150 (MQ=255)
                                                                                                          TCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                             >  2:2383523‑M1/1‑150 (MQ=255)
                                                                                                          TCGCGGGCTTGCCCGCCCCCCCACGCCCCCACAAGCGCCCAAGCCCGGGGCGGCCCCTGGGGGCGCGGCCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa                                             >  1:1907874‑M1/1‑150 (MQ=255)
                                                                                                            GCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGt                                           <  2:3230329‑M1/150‑1 (MQ=255)
                                                                                                              GGGCTTGCCCGCTCTCCCACGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCa                                         >  2:4316609‑M1/1‑150 (MQ=255)
                                                                                                               GGCTGGCCCGCACGCCCCCGGCCCCACCAGCGTCCAAGACCTGGGGGGCCCGTGTGGGAGGGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCAt                                        >  1:1263936‑M1/1‑150 (MQ=255)
                                                                                                                  TTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAgg                                               >  2:5469626‑M1/1‑140 (MQ=255)
                                                                                                                  TTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAgg                                               <  1:5469626‑M1/140‑1 (MQ=255)
                                                                                                                  TTGCCCGCTCCCCCAAGGCCCCACACGCGTCAAATACGTGTGCGGTCCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCGATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGt                                     >  1:3531182‑M1/1‑150 (MQ=255)
                                                                                                                   TGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCGGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa                                    >  2:5209605‑M1/1‑150 (MQ=255)
                                                                                                                      CCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCCGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGt                                 <  1:2201143‑M1/150‑1 (MQ=255)
                                                                                                                      CCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGt                                 <  2:6652970‑M1/150‑1 (MQ=255)
                                                                                                                      CCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGt                                 <  1:921730‑M1/150‑1 (MQ=255)
                                                                                                                        GCTCCCCCAAGGCCCCACAAGCGCCCAACCCCTGTGCGGTCCCTGTGGGCGCGGGCACGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTAc                               >  1:1106606‑M1/1‑150 (MQ=255)
                                                                                                                           CTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGc                            >  2:5938496‑M1/1‑150 (MQ=255)
                                                                                                                              ccAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCCGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGa                         >  1:2390320/1‑150 (MQ=255)
                                                                                                                               cgaGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                        <  1:5042785/148‑1 (MQ=255)
                                                                                                                               cAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                        >  1:5933005/1‑150 (MQ=255)
                                                                                                                               cAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGaa                        >  1:5048610/1‑150 (MQ=255)
                                                                                                                                aaGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAg                       >  2:3489490/1‑150 (MQ=255)
                                                                                                                                aaGGCCCCAAAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAg                       <  1:4424646/150‑1 (MQ=255)
                                                                                                                                 aGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGt                      <  1:6771581/150‑1 (MQ=255)
                                                                                                                                  ggCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTc                     <  2:2538447/150‑1 (MQ=255)
                                                                                                                                    ccccACACGCGTCAAATACCTGGGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTcca                   <  1:4071011/150‑1 (MQ=255)
                                                                                                                                    cccaACATGCGTCCAATACCTGTGGGTTCCCTGTGGGAGGGGGAAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTcca                   <  1:2966196/150‑1 (MQ=255)
                                                                                                                                     cccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccac                  >  2:6302117/1‑150 (MQ=255)
                                                                                                                                     cccACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccac                  <  2:389937/150‑1 (MQ=255)
                                                                                                                                      ccACACGCGTCCAATACCTGTGCGGTCCCGGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacc                 <  2:5818860/150‑1 (MQ=255)
                                                                                                                                       cacaAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTccacca                <  2:5808582/150‑1 (MQ=255)
                                                                                                                                         caAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGCCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCat                                                                              <  2:5590375/86‑1 (MQ=255)
                                                                                                                                          aaGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCcaccacac             <  2:2084307/150‑1 (MQ=255)
                                                                                                                                              gTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATc         >  2:3800221/1‑150 (MQ=255)
                                                                                                                                                  aaTCCCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGcc     <  2:1951548/150‑1 (MQ=255)
                                                                                                                                                  aaTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGcc     >  1:4018050/1‑150 (MQ=255)
                                                                                                                                                     aCCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGAc  <  1:6638939/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTAGGTACCGCTGAAGTCCACCACACCATCGGCCGAC  >  NC_002947/3344007‑3344305

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