Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,902 A→G 21.7% S267P (TCC→CCC)  rsmC ← 16S rRNA m(2)G1207 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,9020AG21.7% 126.5 / 5.8 60S267P (TCC→CCC) rsmC16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (8/37);  new base G (11/3);  total (20/40)
Fisher's exact test for biased strand distribution p-value = 5.92e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTT  >  NC_000913/4606665‑4607125
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
aTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                           <  1:717120/239‑1 (MQ=255)
           ttcttcGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCCTCGTGGAACGCCGGGTCGGGGGTGATCaaatt                                                                                                                                                                                                                                >  1:796078/1‑238 (MQ=255)
               tccccgcgggggTTATAATGGCGCGATACCCCGTGAAGCGCCCTTTTTGCGCGATCCTCTCGTGGAATCCAAATGTCTCACCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCaaaa                                                                                                                                                                                                                           <  2:903560/230‑1 (MQ=255)
                  ccTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCCCATCCAGGCTGGTTTGCATCCCACCGGGGAACGGCGGGGGGGGGAGGATAATATAaaaaaa                                                                                                                                                                                                                         >  2:619537/1‑236 (MQ=255)
                                   gcgcGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCAGATCAAAACAACCTTTCACCCCGGaaaa                                                                                                                                                                                                    >  1:781463/1‑239 (MQ=255)
                                    cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCCCCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGGGGGGAGGGGCATATCAAAAAAACCTTTTCCCCCGGgaaaa                                                                                                                                                                                                   >  2:879021/1‑238 (MQ=255)
                                         tACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGGGGGTTGGGGGTGGTACTATCAAAAAGACCTTTCAACCCGGaaaaaagggt                                                                                                                                                                                              >  2:259217/1‑233 (MQ=255)
                                                ttaaGCGCCCTTTTGGGGCGATCCGTTCGGGGAGGCCAAATTTCTCCTCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCa                                                                                                                                                                                       <  2:767885/237‑1 (MQ=255)
                                                tGAAGCGCCCTGTTTGCGCGACCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGAGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCa                                                                                                                                                                                       <  1:818942/239‑1 (MQ=255)
                                                   aGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGGTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGTTGGAGATGCTCATATCAAAAAGACCTTTAACCTCGGCAAAAACGTTTCTGGCACAg                                                                                                                                                                                    >  1:771385/1‑239 (MQ=255)
                                                        ccggTTTGCGCGTTCACTTCGTGGAAGCCAAATGTCTCACCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                               <  1:768416/236‑1 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGATTGGGATGATCACATCACAACGACCTTTCCACCTGGGAAAAACGTTGCCGGCaaaagcct                                                                                                                                                                                >  1:103708/1‑233 (MQ=255)
                                                           gTTTGCGCGATCACGTCGTGGAAGCCAAATGTCTCATCCAGCCCGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTGGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCAcc                                                                                                                                                                            <  1:19099/239‑1 (MQ=255)
                                                             ttGCGCGATCACTTCGTGGAAGCCAAAGTGCTCATCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAAGTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCtt                                                                                                                                                                          <  1:136656/239‑1 (MQ=255)
                                                                        cTTCGTTGAACCCAAATTTCTCACCCAGCCCGTCCGGATAAGGCCGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTg                                                                                                                                                              <  2:781463/239‑1 (MQ=255)
                                                                          tCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGc                                                                                                                                                            <  1:59518/239‑1 (MQ=255)
                                                                             tGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGc                                                                                                                                                         <  1:325418/239‑1 (MQ=255)
                                                                              ggAAGCCAAATGTCACATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGc                                                                                                                                                         <  1:848066/238‑1 (MQ=255)
                                                                                 aGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGCTGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCACCGTGGAACGGCGGGCTGGGTGGGGGATTTCAAAAACGCCCCATCACCTCGGGAAAACCGCTGTTGGCGACAACTTCACCCTCAACACCCTTGGCCGCaaag                                                                                                                                                     >  2:714456/1‑237 (MQ=255)
                                                                                   ccAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGGTCATATCAAAACGCCCTTTCTCCCCCGAAAAAACGTTGTTGGCGAAAACTTTACCCTCAACACCGTTGGCCGccacggt                                                                                                                                                   >  1:791776/1‑233 (MQ=255)
                                                                                       aTGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTgcg                                                                                                                                               <  1:619537/239‑1 (MQ=255)
                                                                                           cttctCCAGCACGGCCGGAAAAGGCAGGAAGGCGTTCGGTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCgg                                                                                                                                            <  1:750849/234‑1 (MQ=255)
                                                                                           ctcATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCCCTAAAATTCGCAGCTCGGCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCgg                                                                                                                                            <  1:744512/238‑1 (MQ=255)
                                                                                                    cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGCATGATCAAATCACAACAACCTTTCACCTCGGCAAAACACTGGCTGGCAAAGACTTCACCTTCAAAACCCATCACCCCAAGGGTTGCCCggcggggcttca                                                                                                                                  >  2:155404/1‑233 (MQ=255)
                                                                                                         ccGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCTAAAAGACCTTTCCCCTCGGAAAAAACGTTGCGGGCCAAAACTTTACCCTCAAACACCTTGGCCGCAAGCGGTGGGGGGCTGGGTTCTAccgcc                                                                                                                            >  2:92618/1‑239 (MQ=255)
                                                                                                         acGGATAAGGCAGGAAGGCGTTCGCTACAATCCCCAGCCCCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgc                                                                                                                             <  1:816279/237‑1 (MQ=255)
                                                                                                               aaGGCAGGAAGGCTGTCGCTACAATTCGCAGCTCGCCGCCGCTATGAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCGCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCa                                                                                                                      <  1:664736/239‑1 (MQ=255)
                                                                                                                aGGGAGGAAGGCGTTCGCTACAATTCGCACGTCGCCGCCGCTATTAAGATGCCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCag                                                                                                                     <  1:496301/239‑1 (MQ=255)
                                                                                                                     ggAAGGCGTTCGTTACAATTCGCAGCTCGGCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGAc                                                                                                                <  1:455271/238‑1 (MQ=255)
                                                                                                                       aaGGCGTTCGGTTCAATTCGCAGGCCGCCGCCGATATTAAGATGACGCACCGCACCGCGAATCCGCGTTTGCACCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACAt                                                                                                              <  1:833453/238‑1 (MQ=255)
                                                                                                                           cGTTCGCTACAATTCCCAGTCCGCCGCCGGTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGca                                                                                                          <  1:795044/238‑1 (MQ=255)
                                                                                                                                                  ccgccgcTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGa                                                                                <  1:607860/239‑1 (MQ=255)
                                                                                                                                                      cgcTATTAAGGTGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATgg                                                                            <  1:714456/239‑1 (MQ=255)
                                                                                                                                                      cgcTATTAAGATGGCGCACCGGACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATgg                                                                            <  1:22154/239‑1 (MQ=255)
                                                                                                                                                         tATTAAGATGACGCCCCGCACCGGGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTc                                                                                      <  1:53156/227‑1 (MQ=255)
                                                                                                                                                            tAAGCTGACCCCCCGCACCGTGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCa                                                                    <  1:273453/239‑1 (MQ=255)
                                                                                                                                                                 tgccgcaccgcccCGGGAATCAGCGTTTGCGCCGCCTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGGCGTTGCTGGCAAAGACTTCACCGTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGc                                                               <  1:801521/236‑1 (MQ=255)
                                                                                                                                                                   aCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTGACCTCGGAAAAGACGGTTCTGGCAAAGACTTCACCTTCCAAACCGATGGGCGCAAGAGGTGCGGGGCTGACGTTTCACCCCCGGGCCAAAGCACATCCCAAGGGGAGAGGGATTTTTGGCAGAAGGGGGCGCaaaaca                                                               >  2:742273/1‑235 (MQ=255)
                                                                                                                                                                     gcaccgcaccGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGGTTGAGATTAACATATCAAAACGAACCTTCCCCTCCGGAAAAACCGTTCCGGGCAAAACTTCAACCTCAACACCCTTGGCCGCCAAGGTTGCGCGGCCGGCCTCTAACGCCGGGGCCAAGACATTCCCCAAGGGGGGCCGAAACTTTCGGGCACGGCGCGCAAAAGGCAAc                                                            >  2:200774/1‑239 (MQ=255)
                                                                                                                                                                          gcaccgaGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAg                                                      <  1:667935/239‑1 (MQ=255)
                                                                                                                                                                                 gAATCCGCGTTTGCGCCGCATCCAGGCTTGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGcc                                                <  1:852588/238‑1 (MQ=255)
                                                                                                                                                                                   aTCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCAAAAAGACCTTTCACCTCGGAAAAGACGGTGCTGGGAAAAACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGGTGGCTTCTTCCGCCCGGGCAAAGGAATCGCCACAGGGGAGACGGATTTTTGGCGGAAGGGCGGCAAAAGAAACAGaaaaaaaacccgc                                             >  2:852588/1‑229 (MQ=255)
                                                                                                                                                                                        cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCa                                                                                                                                                                                       <  2:706304/103‑1 (MQ=255)
                                                                                                                                                                                        cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCa                                                                                                                                                                                       >  1:706304/1‑103 (MQ=255)
                                                                                                                                                                                         gTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTgcg                                                                                                                                               <  1:717386/141‑1 (MQ=255)
                                                                                                                                                                                            tGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGATCATATCAAAACGACCTTTCACCCCGGGAAAAACGTTGCTGGCAAAGACTTTACCTTCAACACCCTTGGGCGCAAGGGTTGGGCGGCCGGCTTTTACCCCCGGGGCGGAGACCTCGCACAAGGGGAGGCGAAATTTCTGCGAAAGGGCGGCAAAAGAAACAGAAAACAAGCCCGCCCCCCAAcc                                    >  2:904944/1‑239 (MQ=255)
                                                                                                                                                                                             ggcccgcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCCCCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                   <  1:897806/236‑1 (MQ=255)
                                                                                                                                                                                             gcgcCGGATCCAGGGTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATTATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                   <  1:803377/239‑1 (MQ=255)
                                                                                                                                                                                             gcgcCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCAAAACGACCTTTCACCCCGGGAAAGACGTTTCTGGCAAAGAATTCACCTTCCACACCGTTGGCCGCAAGTGTTGCGCGGCGGGCTTCTAACCCCGGGGACGAGACATCCCACAACGAGGGACGGATTTTTGGCGAAGAGCGCGACAAAGACAACTAAAAACCCGccccgcccccagcg                                    >  2:312557/1‑228 (MQ=255)
                                                                                                                                                                                                   gaTCCAGGGTGGTTTGCACCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCCCCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATcc                             <  1:92618/238‑1 (MQ=255)
                                                                                                                                                                                                   cATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATcc                             <  1:606822/239‑1 (MQ=255)
                                                                                                                                                                                                    aTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCa                            <  1:146287/239‑1 (MQ=255)
                                                                                                                                                                                                              gTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTa                    >  1:95039/1‑237 (MQ=255)
                                                                                                                                                                                                              gTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGGGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACCAGGTGAGACGCATTTTCGGCGAATGGCGCCCCAAGGCAAATGAAAGCACGCCCGCGCCCCAGCCGACATcctggctctttt                   >  1:4251/1‑228 (MQ=255)
                                                                                                                                                                                                                      cccATGGGGGAACGGCGGGTTGGAGATGATCATGTCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCCCCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATcc                             <  1:625293/220‑1 (MQ=255)
                                                                                                                                                                                                                       cctgtgtgtAACGGGGGGTTGGGTAGGACAATATCAAAACGTCCTTTCCCCTCGGAAAAGAAGTTTCTGGCAAAGGCTTCCCCCTCAACCCCGTTGGCCGCAAGTGTTGGGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCAAAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg         <  2:157097/234‑1 (MQ=255)
                                                                                                                                                                                                                        cATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCCCCGCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc        <  1:856072/239‑1 (MQ=255)
                                                                                                                                                                                                                        cATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGGAAAGACGTTGCTGGCCAAGACTTCACCTTCAACACCCTTGGCCGCAAGGGTTGCGCGGCTGGCTTCTACCGCCGGGGCCGAGACATCGCACAAGGTGGGACGAAGTTTCGTCGAATGGGGCGCCAAGGCAACTGCAAAACCCCCCCCCCCACCCGCCACATCCACCACCTTACCATTctcgtgg        >  2:642747/1‑233 (MQ=255)
                                                                                                                                                                                                                         catgggggaCGGGGGGTTGGGGGATGAACATATCAAAACTCCCTTCCGCCTCGGAAAAGTCGTGGCGGGAAAGGACTGCCCCTCCACCACCGTTGGCCGAAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTGCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc        <  2:796078/231‑1 (MQ=255)
                                                                                                                                                                                                                            gTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCCCCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAg    <  1:21866/239‑1 (MQ=255)
                                                                                                                                                                                                                               gAACGGGGGGTTGGAAATGAACAAATCAAAACGCCCTTTCCACTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTGGGCCGCAAGTGTGGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGtt  <  1:839849/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTT  >  NC_000913/4606665‑4607125

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