Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,903 A→G 25.2% 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,9030AG25.2% 137.8 / 9.2 70S267P (TCC→CCC) rsmC16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (8/42);  new base G (15/4);  total (24/46)
Fisher's exact test for biased strand distribution p-value = 1.81e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTA  >  NC_000913/4606666‑4607127
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
aTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                             <  1:717120/239‑1 (MQ=255)
     taccgttttttcgtcgtgggggTTATAAGGGTGCGATACCGCTTGCAGCCCCCTGTTTGCGGCAACACTTGGGGGAAGCAAAAGGTCTCATCCCCCACGTCCGGATAAGGCAGGAAGTGTGTCGCTCCAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAgatgat                                                                                                                                                                                                                                        <  2:94737/220‑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)
                                                   aGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGGTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAAAGGCGGGGGGGGGGTGATCACATTAAAACAACCCTTCCCCCCCGaaaaaaaagtggagggcaaa                                                                                                                                                                                      >  2:856072/1‑223 (MQ=255)
                                                        cctgtttggtggtgctcttgggggAGGCAAAATGTCTCACCCGGCACGCGCGGATAAGGCGGGAGGGCGTTCGCTACAATTCGCAGCGCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                                 <  2:227090/217‑1 (MQ=255)
                                                        ccggTTTGCGCGTTCACTTCGTGGAAGCCAAATGTCTCACCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                                 <  1:768416/236‑1 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCCTCGAGGAACGGCGGGGTGGGGGGAAGGCTATCAAAAAAACCCTTCACCCCGGAAAAGaacgtgttgggaaaaacata                                                                                                                                                                                 >  2:795044/1‑220 (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)
                                                                                                   ccccccccggaatggggcggAAGGCCTTCGCACCAATTCCGCGCCCGCCGCCGCTATTAAGATGGCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGCCCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTCCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCt                                                                                                                                    <  2:841037/221‑1 (MQ=255)
                                                                                                    cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGCATGATCAAATCACAACAACCTTTCACCTCGGCAAAACACTGGCTGGCAAAGACTTCACCTTCAAAACCCATCACCCCAAGGGTTGCCCggcggggcttca                                                                                                                                    >  2:155404/1‑233 (MQ=255)
                                                                                                        gccgggatggggcggAAGGTGTTCGCCAAAATTCGCCACCCGCCCCCCGTTTTAAGGTGCCGCACCGCACCGAGAATCAGCGTTGGCAGCGCATCCAGGCATGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATAATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgc                                                                                                                               <  2:202123/226‑1 (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)
                                                                                                                                                     ccggttttaagggtgcgacccgccccgcGAAACAGGGTTGGCGCCACATCCAGGCTGGTTTGCGTCCCAACGGTGAACGGCGGGTTGGAGATGATCATCTCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATg                                                                               <  2:791776/220‑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)
                                                                                                                                                                   aCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCAAAACGACCTTTCCCCTCCGGAAAGACGTTGCTGGCAAAAACTTTACCTTCACAACCCTTTGCCGCAAGTGTTGCGCGGCGGGGCTTTACCGCCCGGGGGAAGACACCCCCCAAAGGGGGCCGAAATTTTGGGGAATgggcgggtctccttt                                                                 >  2:53156/1‑226 (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)
                                                                                                                                                                                        ataaaggacagctaccAGGCTGGTTGGCATCCCATGGTGGAACGGGGGGTTGGGGAGGATCAAAACAAAACGACCTTCCCCCCGGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTGGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCca                                           <  2:4251/224‑1 (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)
                                                                                                                                                                                                                                aaCGGCGGGTTGGGGATGATCATGTCAAAACGAACTTTCCCCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGGGTTGCGCGGCTGGCTTCTACCCCCGGGGCGGAGACAACCCACAAGGGGGGACGAAGTTTCTGGGCAAGGCGCGCAAAAGCAAACGAAAGAACGCCCCGCCCCACGCCCCCCACCAACACTTTAccccggcccttataaaaaa  >  2:625293/1‑221 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            
ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCAGAGACATCGCACAAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTA  >  NC_000913/4606666‑4607127

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