Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 604,474 A→G 6.7% P101P (CCA→CCG) ‡ nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947604,4740AG 6.7% 1181.8 / 45.7 374P101P (CCA→CCG) ‡nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base A (181/168);  new base G (11/14);  total (192/182)
Fisher's exact test for biased strand distribution p-value = 5.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 383 total aligned reads displayed.

TCGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAA  >  NC_002947/604326‑604623
                                                                                                                                                    |                                                                                                                                                     
tCGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAg                                                                                                                                                      >  2:4831725/1‑150 (MQ=255)
   aaGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCg                                                                                                                                                   >  1:5033667/1‑150 (MQ=255)
    aGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGa                                                                                                                                                  <  1:5046154/150‑1 (MQ=255)
    aGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGa                                                                                                                                                  <  1:7322230/150‑1 (MQ=255)
    aGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGa                                                                                                                                                  <  1:6443491/150‑1 (MQ=255)
     gcgcAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAg                                                                                                                                                 <  2:6092052/150‑1 (MQ=255)
       gcAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgct                                                                                                                                               >  1:4481493/1‑150 (MQ=255)
         aaTTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctct                                                                                                                                             <  1:5517580/150‑3 (MQ=255)
           ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttc                                                                                                                                           <  2:4556067/150‑5 (MQ=255)
           ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttc                                                                                                                                           <  1:4353726/150‑5 (MQ=255)
           ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgctggct                                                                                                                                           >  2:7944862/1‑150 (MQ=255)
           ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCAAgctggct                                                                                                                                           <  2:3352855/150‑1 (MQ=255)
            tCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttcc                                                                                                                                          <  1:7656154/150‑6 (MQ=255)
            tCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGAGCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttcc                                                                                                                                          <  2:5149903/150‑6 (MQ=255)
             ccGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttcca                                                                                                                                         <  1:3186656/150‑7 (MQ=255)
             ccGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgctcttcca                                                                                                                                         <  2:2697973/150‑7 (MQ=255)
               gggTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGtgc                                                                                                                                       >  2:1353693/1‑150 (MQ=255)
                  tGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGtgctgc                                                                                                                                    <  2:4226072/150‑1 (MQ=255)
                    gATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                                  <  1:7567371/150‑1 (MQ=255)
                      tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                                >  2:3988249/1‑150 (MQ=255)
                      tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                                >  1:2328419/1‑150 (MQ=255)
                      tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                                >  1:8251040/1‑150 (MQ=255)
                        aCGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGt                                                                                                                              <  2:7396862/150‑1 (MQ=255)
                         cGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTc                                                                                                                             >  1:1494160/1‑150 (MQ=255)
                            tttCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCAc                                                                                                                          <  2:7603432/150‑1 (MQ=255)
                            tttCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCAc                                                                                                                          <  2:8207553/150‑1 (MQ=255)
                            tttCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCAc                                                                                                                          >  1:3388123/1‑150 (MQ=255)
                              tCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTCGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCt                                                                                                                        <  2:2614521/150‑1 (MQ=255)
                               cACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTg                                                                                                                       >  2:7165842/1‑150 (MQ=255)
                                aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGTGCCTGTCCACCTgg                                                                                                                      >  2:2203839/1‑150 (MQ=255)
                                  cGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGa                                                                                                                    <  1:7118529/150‑1 (MQ=255)
                                       tCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGa                                                                                                               <  1:1610124/150‑1 (MQ=255)
                                       tCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGa                                                                                                               <  1:2482434/150‑1 (MQ=255)
                                         gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATc                                                                                                             <  1:4271420/150‑1 (MQ=255)
                                         gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATc                                                                                                             <  2:7084007/150‑1 (MQ=255)
                                            ggTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                          <  1:3137002/150‑1 (MQ=255)
                                            ggTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                          >  2:3504157/1‑150 (MQ=255)
                                             gTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                         <  2:1931007/150‑1 (MQ=255)
                                             gTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                         <  2:150110/150‑1 (MQ=255)
                                               gCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAcgc                                                                                                       >  2:4407610/1‑150 (MQ=255)
                                                 cTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGt                                                                                                     <  1:6220479/150‑1 (MQ=255)
                                                       ccGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGt                                                                                               <  1:8280663/150‑1 (MQ=255)
                                                          cgAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGcc                                                                                            >  2:4163937/1‑150 (MQ=255)
                                                              aTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTAGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCATACGCGTCGACGTGCCTTAc                                                                                        <  2:2835597/150‑1 (MQ=255)
                                                                ccTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACcg                                                                                      >  2:6356810/1‑150 (MQ=255)
                                                                    cATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtc                                                                                  <  1:1871436/150‑1 (MQ=255)
                                                                     aTGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtcg                                                                                 >  2:7868316/1‑150 (MQ=255)
                                                                        gggTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGa                                                                              <  2:526309/150‑1 (MQ=255)
                                                                          gTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATc                                                                            >  2:1625239/1‑150 (MQ=255)
                                                                           tCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCa                                                                           <  1:398186/150‑1 (MQ=255)
                                                                             ccGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAAc                                                                         >  2:3754346/1‑150 (MQ=255)
                                                                              cGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCATTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACg                                                                        >  2:3641836/1‑150 (MQ=255)
                                                                              cGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACg                                                                        <  1:1014447/150‑1 (MQ=255)
                                                                              cGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACg                                                                        >  2:1459526/1‑150 (MQ=255)
                                                                                gCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTGGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGaa                                                                      <  1:2110253/150‑1 (MQ=255)
                                                                                 cAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAg                                                                     <  1:5652144/150‑1 (MQ=255)
                                                                                     cAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgtgt                                                                 <  2:6578775/150‑1 (MQ=255)
                                                                                     cAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgtgt                                                                 >  1:782120/1‑150 (MQ=255)
                                                                                      aaGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                                >  2:4630318/1‑150 (MQ=255)
                                                                                         ggCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGaa                                                             >  2:1291765/1‑150 (MQ=255)
                                                                                         ggCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGaa                                                             >  2:3201136/1‑150 (MQ=255)
                                                                                         ggCGAAATCGACAAAGCGCTGGTGCCTCGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGaa                                                             >  1:5981743/1‑150 (MQ=255)
                                                                                           cGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACt                                                           >  1:359681/1‑150 (MQ=255)
                                                                                           cGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACt                                                           <  2:1720709/150‑1 (MQ=255)
                                                                                             aaaTCGACAAAGCTCTGGTGCCTTGCATGACCAGCGCACTGTAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGTTCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                         <  2:8634859/150‑1 (MQ=255)
                                                                                                 cGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCaa                                                     <  1:8818494/150‑1 (MQ=255)
                                                                                                  gACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAg                                                    >  1:1253824/1‑150 (MQ=255)
                                                                                                   aCAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAgg                                                   >  2:7033019/1‑150 (MQ=255)
                                                                                                   aCAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAgg                                                   >  1:3720436/1‑150 (MQ=255)
                                                                                                      aaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCt                                                >  2:6514415/1‑150 (MQ=255)
                                                                                                       aGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCtt                                               >  1:2302112/1‑150 (MQ=255)
                                                                                                        gcgcTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTc                                              <  2:2944960/150‑1 (MQ=255)
                                                                                                           cTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgc                                           <  1:3076979/150‑1 (MQ=255)
                                                                                                             ggTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcgc                                         >  1:2664631/1‑150 (MQ=255)
                                                                                                               tGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCa                                       <  1:4670965/150‑1 (MQ=255)
                                                                                                                   ttGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGa                                   >  2:4306866/1‑150 (MQ=255)
                                                                                                                   ttGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTAGGCGCCACTGa                                   >  2:5565633/1‑150 (MQ=255)
                                                                                                                          aCCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCcac                            >  1:6055007/1‑150 (MQ=255)
                                                                                                                            cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCcacac                          >  2:3143643/1‑149 (MQ=255)
                                                                                                                            cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                          >  2:8473413/1‑150 (MQ=255)
                                                                                                                            cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                          >  1:6747803/1‑150 (MQ=255)
                                                                                                                            cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                          <  2:6803013/150‑1 (MQ=255)
                                                                                                                                gCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTc                      <  2:1092094/150‑1 (MQ=255)
                                                                                                                                 cACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCg                     >  2:582028/1‑150 (MQ=255)
                                                                                                                                     ggAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa                 <  2:5239853/150‑1 (MQ=255)
                                                                                                                                     ggAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa                 <  1:5131062/150‑1 (MQ=255)
                                                                                                                                      gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAAc                >  2:4973759/1‑150 (MQ=255)
                                                                                                                                        ccagCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACgg              >  2:544869/3‑150 (MQ=255)
                                                                                                                                         gagCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGc             <  2:2051199/150‑1 (MQ=255)
                                                                                                                                          agCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg            >  1:8716405/1‑150 (MQ=255)
                                                                                                                                          agCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg            >  1:8118319/1‑150 (MQ=255)
                                                                                                                                          agCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg            >  2:4288805/1‑150 (MQ=255)
                                                                                                                                            cTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTa          >  1:5283962/1‑150 (MQ=255)
                                                                                                                                             tCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTAc         >  2:910856/1‑150 (MQ=255)
                                                                                                                                             tCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTAc         >  1:6614828/1‑150 (MQ=255)
                                                                                                                                                aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTCCTgg      >  1:8596104/1‑150 (MQ=255)
                                                                                                                                                aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTgg      <  2:2899249/150‑1 (MQ=255)
                                                                                                                                                  ccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGAc    >  1:6071769/1‑150 (MQ=255)
                                                                                                                                                    aGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACaa  >  2:6944501/1‑150 (MQ=255)
                                                                                                                                                    aGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACaa  >  1:1176999/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGACAA  >  NC_002947/604326‑604623

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