Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 4,611,582 A→G 21.5% S267P (TCC→CCC)  rsmC ← 16S rRNA (guanine(1207)‑N(2))‑methyltransferase RsmC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077794,611,5820AG21.5% 124.4 / 4.9 56S267P (TCC→CCC) rsmC16S rRNA (guanine(1207)‑N(2))‑methyltransferase RsmC
Reads supporting (aligned to +/- strand):  ref base A (7/35);  new base G (11/2);  total (19/37)
Fisher's exact test for biased strand distribution p-value = 1.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTC  >  NC_007779/4611345‑4611811
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
aTAATTCGCCTTTCTGCGCCTGGCGGGTCATAATGGCTCGATACACCTTGAAGCGCCGTGTTTGCGCGATCGCTTCGTGGAGGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCGCCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                         <  2:35324/239‑1 (MQ=255)
 tAATTACCCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCCCTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGga                                                                                                                                                                                                                                          <  1:1101009/237‑1 (MQ=255)
 tAAATAACCTTGTGTCGCCTGGCGGTTCATAATGGCGCGATACACCTTGAAGCCCCCTTTTTGCGCGATCACTTCGTGGAAGCCAAATGTATCATCCAGCACGTCCGGATAAGGCAGGAAGGCCTTCGCTTCAATTCGCAGCACGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                         <  2:377352/238‑1 (MQ=255)
      aaCCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATc                                                                                                                                                                                                                                   <  1:1057065/239‑1 (MQ=255)
        ccTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCat                                                                                                                                                                                                                                 <  1:286874/239‑1 (MQ=255)
             cttcGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCaa                                                                                                                                                                                                                          <  1:409450/239‑1 (MQ=255)
                            cATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAgatgat                                                                                                                                                                                                                                    <  1:426598/216‑1 (MQ=255)
                                    cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACACTTCGCAGCTCGCCGCCGCTAATAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATACCATCGTGGAACGGCGGGGTGGGGATTGTCATATTAAAACGACCTTTCACCTCGGgaaaa                                                                                                                                                                                                  >  1:720421/1‑238 (MQ=255)
                                    cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAg                                                                                                                                                                                                  >  1:1424191/1‑239 (MQ=255)
                                    cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGATCATATCAAAACCACCCTTCAcccagggcaaa                                                                                                                                                                                                  >  1:1443341/1‑230 (MQ=255)
                                              cTTGAAGCGCCCTGTTTGCGAGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgaga                                                                                                                                                                                                                                        <  1:891636/194‑1 (MQ=255)
                                                tGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGGTGAACATTTCAAAACGACCTTTCCCCTCGGAAAAGACGTTGCTGGcc                                                                                                                                                                                      >  1:507969/1‑238 (MQ=255)
                                                 gAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCaa                                                                                                                                                                                     <  1:362222/239‑1 (MQ=255)
                                                   aGCGCCGTGACTGCGCCATAAATTCGTGGAAGTTAACTGTCACCTCCAGCACGTACACATCAGGCAGGAAGGACTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCAGCGTGGAACGGGGGGGTGGAGAGGATCACAACAACAAGACCCATACACACGGGAAAAACGTTGCCGGcggcg                                                                                                                                                                                   >  2:1369819/1‑235 (MQ=255)
                                                   aGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                         >  1:1570823/1‑188 (MQ=255)
                                                   aGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                         <  2:1570823/188‑1 (MQ=255)
                                                   aGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGGAGGGCATAACAAAAAGACCTTTCTCCCCCGCAAAAACAATGGTGGCGAAg                                                                                                                                                                                   >  2:1176780/1‑239 (MQ=255)
                                                      gcCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCCCCTCGGAAAAGACGTTGCTGGCAAAGACt                                                                                                                                                                                <  1:1369819/239‑1 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCAAAAAGAACCTTCACCTCGGGAAAAAAGGTGCTGGCAAAAACtt                                                                                                                                                                               >  1:1717517/1‑238 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                              <  1:888192/239‑1 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                              <  1:368136/239‑1 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                              <  2:556356/239‑1 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTAGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGTGAGGATTAAATCAAAACAACCCTTTCACTCGGGAAAAAAATTTCTGGCAAAAAATTc                                                                                                                                                                              >  1:1042102/1‑239 (MQ=255)
                                                           ttttggcggATCACTGCGGTGAAGCCAAATGTCTCATCCAGCACTTCCGGATAAGGCAGGAGGGCGTTCGCTACAAGTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCAcc                                                                                                                                                                           <  2:417969/231‑1 (MQ=255)
                                                           gTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCAcc                                                                                                                                                                           <  2:1424191/239‑1 (MQ=255)
                                                            tttGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                              <  1:29929/235‑1 (MQ=255)
                                                                    aTCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAacac                                                                                                                                                                   >  1:1032765/1‑238 (MQ=255)
                                                                            gTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCg                                                                                                                                                          <  1:623199/239‑1 (MQ=255)
                                                                                 aGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCTGGAAGGCGTTCGCTACAATTCGCAGCCCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGAGGAGCACATTCAAAAGAACCTTTCCCCCGGAAAAGACGGTGCTGGGAAACAATTTACCCTTAACACCGTTGGCCGccacg                                                                                                                                                     >  2:286550/1‑235 (MQ=255)
                                                                                   ccAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAgtgt                                                                                                                                                   >  1:1193541/1‑239 (MQ=255)
                                                                                       aTGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGGGGGGATGATCATATTAAAACAACCCTTCTCCCCGGCAAAAACGTTGTTGGGAAAAACTTCACCTCCAACACCCATGGCCGccacagctgt                                                                                                                                                >  2:222515/1‑229 (MQ=255)
                                                                                          tctcATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                              <  1:1354046/205‑1 (MQ=255)
                                                                                           ctcATCCAGCACGTCCGGATAAGGCAGGGAGGCGTTCGCGCACATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCggc                                                                                                                                           <  2:1490609/239‑1 (MQ=255)
                                                                                                    cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGGTGGACACATCAAAAAAACCCTTTTCCCCCGAAAAAAAGGTTCTGGGCAAAACCTCTCCTTCAACACACTCGGCCCCCAGGGTTGggcgcgccggctcc                                                                                                                                   >  2:1391609/1‑228 (MQ=255)
                                                                                                    cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCat                                                                                                                                                                                                                                 <  1:731918/147‑1 (MQ=255)
                                                                                                    cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGATCat                                                                                                                                                                                                                                 >  2:731918/1‑147 (MQ=255)
                                                                                                    cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGAGGAGCATATCCAAAACACCCTTTCCCCCGGCAAAAAAGGTGCTGGCCAAAAACTCATCGTCAACACCACTGACCGCGAGAGTTgcggggcgggactct                                                                                                                                   >  2:22065/1‑233 (MQ=255)
                                                                                                     aCGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCCCCGCGAAGCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCt                                                                                                                                   <  1:1498350/237‑1 (MQ=255)
                                                                                                     aCGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGTTCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCt                                                                                                                                   <  1:1176780/237‑1 (MQ=255)
                                                                                                         cctgagaaGGCGGGAAGGCGTTCGCTACAATCGGCAGCTTCGCGCCGCTATTAATATGACGCACCGCCCCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCGCCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgcc                                                                                                                             <  2:976682/233‑1 (MQ=255)
                                                                                                                aGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCCCCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                      <  2:615080/239‑1 (MQ=255)
                                                                                                                 aggaaggaGGCGTTCCGCGCCAATGCGCAGCTCACCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                      <  2:972994/232‑1 (MQ=255)
                                                                                                                       aaGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACa                                                                                                                 <  1:23169/237‑1 (MQ=255)
                                                                                                                       aaGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACAt                                                                                                                <  2:1032765/238‑1 (MQ=255)
                                                                                                                           cGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCCAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCCCCCGGGCAGAGACATCGccc                                                                                                           >  1:1813139/1‑237 (MQ=255)
                                                                                                                            gTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGcac                                                                                                           <  1:776188/238‑1 (MQ=255)
                                                                                                                                             aGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGt                                                                                                      >  1:191893/1‑227 (MQ=255)
                                                                                                                                             aGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGt                                                                                                      <  2:191893/227‑1 (MQ=255)
                                                                                                                                                             tAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCa                                                                          <  2:304488/239‑1 (MQ=255)
                                                                                                                                                                       caccgcaccgCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgcc                                                                                                                             <  1:834556/178‑1 (MQ=255)
                                                                                                                                                                           gcaccgcGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATCTCAAAAAGCCCTTTCCCCTCCGGCAAAACGGGGCCGGCAAAAACCTCACCTTTAACACCGCTTGCCGCAAGGGTTGCGGGGCTGGCCTCTACCGCCGGGGGCGCGACACCCCACCAAGGGGGACGAATTTTCTGCGAACAGGGCGCAAAggacaccaacag                                                            >  2:404699/1‑229 (MQ=255)
                                                                                                                                                                                         cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACa                                                                                                                 >  1:693535/1‑172 (MQ=255)
                                                                                                                                                                                         cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACa                                                                                                                 <  2:693535/172‑1 (MQ=255)
                                                                                                                                                                                              gggcggcATCCAGGCTGGTTTGCATCCCATCGTGGAAGGGGGGGTTGGGGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                         <  2:1717517/234‑1 (MQ=255)
                                                                                                                                                                                               cgccgcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAACCGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAGGACTTCACCTTCAACACCGTTGGCCGCAAGTGTGGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                         <  2:1416496/238‑1 (MQ=255)
                                                                                                                                                                                                   gcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATc                                    <  1:55226/239‑1 (MQ=255)
                                                                                                                                                                                                                 ttGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCt                       <  1:543964/238‑1 (MQ=255)
                                                                                                                                                                                                                     aTCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCg                   <  1:486208/238‑1 (MQ=255)
                                                                                                                                                                                                                        ccATCGTGGAAGGGGGGGTTGGGGATGATAAAATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCCCCTTCAACACGGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg               <  2:1813139/239‑1 (MQ=255)
                                                                                                                                                                                                                        ccATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg               <  1:1450220/239‑1 (MQ=255)
                                                                                                                                                                                                                         cATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc              <  1:843797/239‑1 (MQ=255)
                                                                                                                                                                                                                                 aaCGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTa       <  1:1806873/238‑1 (MQ=255)
                                                                                                                                                                                                                                              agaTGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACTTc  <  1:409631/230‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTC  >  NC_007779/4611345‑4611811

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