Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 604,473 C→G 7.0% P101R (CCA→CGA) ‡ nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947604,4730CG 7.0% 1214.2 / 56.5 364P101R (CCA→CGA) ‡nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base C (172/166);  new base G (11/14);  total (183/181)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 374 total aligned reads displayed.

TCGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGAC  >  NC_002947/604326‑604621
                                                                                                                                                   |                                                                                                                                                    
tCGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAg                                                                                                                                                    >  2:4831725/1‑150 (MQ=255)
  gAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATAGACAAAGCGCTGGTGCCTTGCAGGACCATCGCACTGGAAGAGCTCGACCCAGTc                                                                                                                                                  >  1:7680683/1‑150 (MQ=255)
   aaGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCg                                                                                                                                                 >  2:6498395/1‑150 (MQ=255)
    aGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGa                                                                                                                                                <  1:7322230/150‑1 (MQ=255)
    aGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGa                                                                                                                                                <  1:6443491/150‑1 (MQ=255)
    aGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGa                                                                                                                                                >  1:5101699/1‑150 (MQ=255)
     gcgcAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAg                                                                                                                                               <  1:586952/150‑1 (MQ=255)
     gcgcAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAg                                                                                                                                               <  2:6092052/150‑1 (MQ=255)
       gcAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgct                                                                                                                                             <  1:6890150/150‑1 (MQ=255)
         aaTTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctct                                                                                                                                           <  1:5517580/150‑3 (MQ=255)
           ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttc                                                                                                                                         <  1:4353726/150‑5 (MQ=255)
           ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCAAgctggct                                                                                                                                         <  2:3352855/150‑1 (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)
                 gTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGGCGAGCTGGCTGtgctg                                                                                                                                   >  2:7577278/1‑150 (MQ=255)
                    gATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                                <  1:6016981/150‑1 (MQ=255)
                    gATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                                <  1:7567371/150‑1 (MQ=255)
                    gATAACGATTTCACCGATGTCGACGGTGCCTATGTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                                <  2:7039885/150‑1 (MQ=255)
                      tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                              >  2:3988249/1‑150 (MQ=255)
                      tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                              >  1:8251040/1‑150 (MQ=255)
                      tAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCt                                                                                                                              >  1:2328419/1‑150 (MQ=255)
                            tttCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCAc                                                                                                                        <  2:8207557/150‑1 (MQ=255)
                            tttCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCAc                                                                                                                        >  1:3388123/1‑150 (MQ=255)
                            tttCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCAc                                                                                                                        <  2:7603432/150‑1 (MQ=255)
                              tCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTCGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCt                                                                                                                      <  2:2614521/150‑1 (MQ=255)
                                aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGTGCCTGTCCACCTgg                                                                                                                    >  2:2203839/1‑150 (MQ=255)
                                  cGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGa                                                                                                                  <  1:5578629/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                                                                                                           >  1:1183738/1‑150 (MQ=255)
                                         gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATc                                                                                                           <  1:2134095/150‑1 (MQ=255)
                                         gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATc                                                                                                           <  1:3819399/150‑1 (MQ=255)
                                            ggTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                        <  1:3137002/150‑1 (MQ=255)
                                            ggTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                        <  2:8004900/150‑1 (MQ=255)
                                            ggTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                        >  2:3504157/1‑150 (MQ=255)
                                            ggTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                        <  1:1452148/150‑1 (MQ=255)
                                             gTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                       <  2:1931007/150‑1 (MQ=255)
                                             gTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                       <  2:4265039/150‑1 (MQ=255)
                                             gTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                       <  2:150110/150‑1 (MQ=255)
                                                 cTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGt                                                                                                   >  1:5239853/1‑150 (MQ=255)
                                                      tCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACg                                                                                              <  2:6977767/150‑1 (MQ=255)
                                                       ccGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGt                                                                                             >  2:5559633/1‑150 (MQ=255)
                                                         gcgAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGc                                                                                           >  1:8457264/1‑150 (MQ=255)
                                                          cgAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGcc                                                                                          >  2:4163937/1‑150 (MQ=255)
                                                          cgAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGcc                                                                                          >  2:4572185/1‑150 (MQ=255)
                                                                ccTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACcg                                                                                    >  2:6356810/1‑150 (MQ=255)
                                                                    cATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtc                                                                                <  1:5036099/150‑1 (MQ=255)
                                                                    cATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtc                                                                                <  1:1871436/150‑1 (MQ=255)
                                                                    cATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtc                                                                                <  1:1500112/150‑1 (MQ=255)
                                                                     aTGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtcg                                                                               >  2:7868316/1‑150 (MQ=255)
                                                                        gggTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGa                                                                            <  2:526309/150‑1 (MQ=255)
                                                                         ggTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGAt                                                                           <  2:4471333/150‑1 (MQ=255)
                                                                             ccGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAAc                                                                       >  2:3754346/1‑150 (MQ=255)
                                                                              cGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACg                                                                      >  2:1459526/1‑150 (MQ=255)
                                                                              cGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACg                                                                      <  1:1014447/150‑1 (MQ=255)
                                                                                     cAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgtgt                                                               <  2:6578775/150‑1 (MQ=255)
                                                                                     cAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgtgt                                                               >  1:782120/1‑150 (MQ=255)
                                                                                      aaGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                              <  1:3710844/150‑1 (MQ=255)
                                                                                      aaGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                              >  2:5808595/1‑150 (MQ=255)
                                                                                      aaGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                              >  2:4630318/1‑150 (MQ=255)
                                                                                         ggCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGaa                                                           >  2:3200978/1‑150 (MQ=255)
                                                                                         ggCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGaa                                                           >  1:4600406/1‑150 (MQ=255)
                                                                                           cGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACt                                                         >  1:359681/1‑150 (MQ=255)
                                                                                           cGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACt                                                         <  2:1720454/150‑1 (MQ=255)
                                                                                             aaaTCGACAAAGCTCTGGTGCCTTGCATGACCAGCGCACTGTAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGTTCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                       <  2:8634859/150‑1 (MQ=255)
                                                                                             aaaTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                       >  2:1561494/1‑150 (MQ=255)
                                                                                               aTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGcc                                                     >  1:6217519/1‑150 (MQ=255)
                                                                                                 cGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCaa                                                   <  1:8818494/150‑1 (MQ=255)
                                                                                                   aCAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGAGGAACTGGCCAAgg                                                 >  2:7100290/1‑150 (MQ=255)
                                                                                                       aGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCtt                                             >  1:2302204/1‑150 (MQ=255)
                                                                                                        gcgcTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTc                                            <  2:2944960/150‑1 (MQ=255)
                                                                                                           cTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgc                                         <  1:3076979/150‑1 (MQ=255)
                                                                                                            tGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcg                                        >  2:3698119/1‑150 (MQ=255)
                                                                                                             ggTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcgc                                       >  1:2664631/1‑150 (MQ=255)
                                                                                                             ggTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcgc                                       <  2:7172658/150‑1 (MQ=255)
                                                                                                                   ttGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGa                                 >  2:4306870/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)
                                                                                                                                     ggAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa               <  2:5239853/150‑1 (MQ=255)
                                                                                                                                      gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAAc              >  2:4973759/1‑150 (MQ=255)
                                                                                                                                       gccagCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGTCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACg             >  1:4864557/4‑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          >  2:4288805/1‑150 (MQ=255)
                                                                                                                                          agCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg          >  1:8118319/1‑150 (MQ=255)
                                                                                                                                            cTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTa        <  2:2754301/150‑1 (MQ=255)
                                                                                                                                             tCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTAc       >  2:910856/1‑150 (MQ=255)
                                                                                                                                             tCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTAc       >  1:6614828/1‑150 (MQ=255)
                                                                                                                                               gACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTg     >  1:4455850/1‑150 (MQ=255)
                                                                                                                                                aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTCCTgg    >  1:8596104/1‑150 (MQ=255)
                                                                                                                                                aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTgg    <  2:2899249/150‑1 (MQ=255)
                                                                                                                                                aCCCAGTCGAGCTAGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTgg    <  2:1242362/150‑1 (MQ=255)
                                                                                                                                                  ccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGAc  >  1:6071769/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGGAC  >  NC_002947/604326‑604621

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: