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

ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTT  >  NC_000913/4606665‑4607125
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
aTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑Ggag                                                                                                                                                                                                                                           <  1:717120/239‑1 (MQ=255)
           ttcttcGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCCTCGTGGAACGCCGGGTC‑GGGGGTGATCaaatt                                                                                                                                                                                                                                >  1:796078/1‑238 (MQ=255)
               tccccgcgggggTTATAATGGCGCGATACCCCGTGAAGCGCCCTTTTTGCGCGATCCTCTCGTGGAATCCAAATGTCTCACCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑aaaa                                                                                                                                                                                                                           <  2:903560/230‑1 (MQ=255)
                  ccTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCCCATCCAGGCTGGTTTGCATCCCACCGGGGAACGGCGGGGG‑GGGGAGGATAATATA‑aaaaaa                                                                                                                                                                                                                         >  2:619537/1‑236 (MQ=255)
                                   gcgcGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCAGATC‑AAAACAACCTTT‑CACCC‑CGG‑aaaa                                                                                                                                                                                                    >  1:781463/1‑239 (MQ=255)
                                    cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCCCCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGG‑GGGGAGGGGCATATC‑AAAAAAACCTTT‑TCCCC‑CGG‑gaaaa                                                                                                                                                                                                   >  2:879021/1‑238 (MQ=255)
                                         tACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGGGGGTT‑GGGGGTGGTACTATC‑AAAAAGACCTTT‑CAACC‑CGG‑aaaaaagggt                                                                                                                                                                                              >  2:259217/1‑233 (MQ=255)
                                                ttaaGCGCCCTTTTGGGGCGATCCGTTCGGGGAGGCCAAATTTCTCCTCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCa                                                                                                                                                                                       <  2:767885/237‑1 (MQ=255)
                                                tGAAGCGCCCTGTTTGCGCGACCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGAGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCa                                                                                                                                                                                       <  1:818942/239‑1 (MQ=255)
                                                   aGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGGTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGTT‑GGAGATGCTCATATC‑AAAAAGACCTTT‑AACCT‑CGG‑CAAAAACGTTTCTGGCACAg                                                                                                                                                                                    >  1:771385/1‑239 (MQ=255)
                                                        ccggTTTGCGCGTTCACTTCGTGGAAGCCAAATGTCTCACCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTc                                                                                                                                                                               <  1:768416/236‑1 (MQ=255)
                                                        ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGAT‑TGGGATGATCACATC‑ACAACGACCTTTCCACCT‑GGG‑AAAA‑ACGTTGCCGGCaaaagcct                                                                                                                                                                                >  1:103708/1‑233 (MQ=255)
                                                           gTTTGCGCGATCACGTCGTGGAAGCCAAATGTCTCATCCAGCCCGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTGGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCAcc                                                                                                                                                                            <  1:19099/239‑1 (MQ=255)
                                                             ttGCGCGATCACTTCGTGGAAGCCAAAGTGCTCATCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAAGTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCtt                                                                                                                                                                          <  1:136656/239‑1 (MQ=255)
                                                                        cTTCGTTGAACCCAAATTTCTCACCCAGCCCGTCCGGATAAGGCCGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTg                                                                                                                                                              <  2:781463/239‑1 (MQ=255)
                                                                          tCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGc                                                                                                                                                            <  1:59518/239‑1 (MQ=255)
                                                                             tGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGc                                                                                                                                                         <  1:325418/239‑1 (MQ=255)
                                                                              ggAAGCCAAATGTCACATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGc                                                                                                                                                         <  1:848066/238‑1 (MQ=255)
                                                                                 aGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGCTGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCACCGTGGAACGGCGGGCT‑GG‑GTGGGGGATTTCAAAAACGCCCCAT‑CACCT‑CGG‑GAAAACCGCTGTTGGCGACAACTTCACCCTC‑AACACCCTTGGCCGCaaag                                                                                                                                                     >  2:714456/1‑237 (MQ=255)
                                                                                   ccAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGGGATGGTCATATC‑AAAACGCCCTTT‑CTCCC‑CCG‑AAAAAACGTTGTTGGCGAAAACTTTACCCTC‑AACACCGTTGGCCGccacggt                                                                                                                                                   >  1:791776/1‑233 (MQ=255)
                                                                                       aTGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTgcg                                                                                                                                               <  1:619537/239‑1 (MQ=255)
                                                                                           cttctCCAGCACGGCCGGAAAAGGCAGGAAGGCGTTCGGTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCgg                                                                                                                                            <  1:750849/234‑1 (MQ=255)
                                                                                           ctcATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCCCTAAAATTCGCAGCTCGGCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCgg                                                                                                                                            <  1:744512/238‑1 (MQ=255)
                                                                                                    cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGT‑GGGCATGATCAAATC‑ACAACAACCTTT‑CACCT‑CGGCAAAACAC‑TGGCTGGCAAAGACTTCACCTTC‑AAAACCCATCACCCCAAGGGTTGCCCggcggggcttca                                                                                                                                  >  2:155404/1‑233 (MQ=255)
                                                                                                         ccGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGGGATGATCATATC‑TAAAAGACCTTT‑CCCCT‑CGG‑AAAAAACGTTGCGGGCCAAAACTTTACCCTCAAACACC‑TTGGCCGCAAGCGGTGGGGGGCTGGG‑TTCTAccgcc                                                                                                                            >  2:92618/1‑239 (MQ=255)
                                                                                                         acGGATAAGGCAGGAAGGCGTTCGCTACAATCCCCAGCCCCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTAccgc                                                                                                                             <  1:816279/237‑1 (MQ=255)
                                                                                                               aaGGCAGGAAGGCTGTCGCTACAATTCGCAGCTCGCCGCCGCTATGAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCG‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCa                                                                                                                      <  1:664736/239‑1 (MQ=255)
                                                                                                                aGGGAGGAAGGCGTTCGCTACAATTCGCACGTCGCCGCCGCTATTAAGATGCCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCag                                                                                                                     <  1:496301/239‑1 (MQ=255)
                                                                                                                     ggAAGGCGTTCGTTACAATTCGCAGCTCGGCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑Ac                                                                                                                <  1:455271/238‑1 (MQ=255)
                                                                                                                       aaGGCGTTCGGTTCAATTCGCAGGCCGCCGCCGATATTAAGATGACGCACCGCACCGCGAATCCGCGTTTGCACCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACAt                                                                                                              <  1:833453/238‑1 (MQ=255)
                                                                                                                           cGTTCGCTACAATTCCCAGTCCGCCGCCGGTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGca                                                                                                          <  1:795044/238‑1 (MQ=255)
                                                                                                                                                  ccgccgcTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑Ga                                                                                <  1:607860/239‑1 (MQ=255)
                                                                                                                                                      cgcTATTAAGGTGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATgg                                                                            <  1:714456/239‑1 (MQ=255)
                                                                                                                                                      cgcTATTAAGATGGCGCACCGGACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATgg                                                                            <  1:22154/239‑1 (MQ=255)
                                                                                                                                                         tATTAAGATGACGCCCCGCACCGGGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTc                                                                                      <  1:53156/227‑1 (MQ=255)
                                                                                                                                                            tAAGCTGACCCCCCGCACCGTGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑a                                                                    <  1:273453/239‑1 (MQ=255)
                                                                                                                                                                 tgccgcaccgcccCGGGAATCAGCGTTTGCGCCGCCTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGGCGTTGCTGGCAAAGACTTCACCGTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGc                                                               <  1:801521/236‑1 (MQ=255)
                                                                                                                                                                   aCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑GACCT‑CGG‑AAAAGACGGTTCTGGCAAAGACTTCACCTTC‑CAAACCGATGGGCGCAAGAGGTGCGGGGCTGACGTTTCACCCCCGGGCCAAAGCACATC‑C‑CAAGGGGAGAGG‑GATTTTTGGCAGAAGGGGGCG‑C‑aaaaca                                                               >  2:742273/1‑235 (MQ=255)
                                                                                                                                                                     gcaccgcaccGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGGT‑TGAGATTAACATATC‑AAAACGAACCTT‑CCCCTCCGG‑AAAAACCGTTCCGGGCAAA‑ACTTCAACCTC‑AACACCCTTGGCCGCCAAGGTTGCGCGGCCGGC‑CTCTAACGCCGGGGCCAAG‑ACATTCCCCAAGGGGGGCCGAAACTTTCGG‑‑GCACGGCGCG‑CAAAAGGCAAc                                                            >  2:200774/1‑239 (MQ=255)
                                                                                                                                                                          gcaccgaGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGGGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAg                                                      <  1:667935/239‑1 (MQ=255)
                                                                                                                                                                                 gAATCCGCGTTTGCGCCGCATCCAGGCTTGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGcc                                                <  1:852588/238‑1 (MQ=255)
                                                                                                                                                                                   aTCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGGGATGATCATATC‑AAAAAGACCTTT‑CACCT‑CGG‑AAAAGACGGTGCTGGGAAAAACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGGTGGC‑TTCTTCCGCCCGGGCAAAG‑GAATCGCCACAGGGGAGACG‑GATTTTTGGCGGAAGGGCG‑G‑C‑AAAAGAAACAGaaaaaaaacccgc                                             >  2:852588/1‑229 (MQ=255)
                                                                                                                                                                                        cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCa                                                                                                                                                                                       <  2:706304/103‑1 (MQ=255)
                                                                                                                                                                                        cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCa                                                                                                                                                                                       >  1:706304/1‑103 (MQ=255)
                                                                                                                                                                                         gTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTgcg                                                                                                                                               <  1:717386/141‑1 (MQ=255)
                                                                                                                                                                                            tGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGT‑GGGGATGATCATATC‑AAAACGACCTTT‑CACCC‑CGG‑GAAAAACGTTGCTGGCAAAGACTTTACCTTC‑AACACCCTTGGGCGCAAGGGTTGGGCGGCCGGC‑TTTTACCCCCGGGGCGGAG‑ACCTCGCACAAGGGGAGGCG‑AAATTTCTGC‑GAAAGGGCGG‑C‑AAAAGAAACAGAAAACAAGCCCGCCCCCCAAcc                                    >  2:904944/1‑239 (MQ=255)
                                                                                                                                                                                             ggcccgcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CCCCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                   <  1:897806/236‑1 (MQ=255)
                                                                                                                                                                                             gcgcCGGATCCAGGGTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATTATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                   <  1:803377/239‑1 (MQ=255)
                                                                                                                                                                                             gcgcCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGGGATGATCATATC‑AAAACGACCTTT‑CACCC‑CGG‑GAAAGACGTTTCTGGCAAAGAATTCACCTTC‑CACACCGTTGGCCGCAAGTGTTGCGCGGCGGGC‑TTCTAACCCCGGGGACGAG‑ACATCCCACAACGAGGGACG‑GATTTTTGGC‑GAAGAGCGCGAC‑AAAGACAACTAAAAACCCGccccgcccccagcg                                    >  2:312557/1‑228 (MQ=255)
                                                                                                                                                                                                   gaTCCAGGGTGGTTTGCACCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCCCCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATcc                             <  1:92618/238‑1 (MQ=255)
                                                                                                                                                                                                   cATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATcc                             <  1:606822/239‑1 (MQ=255)
                                                                                                                                                                                                    aTCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCa                            <  1:146287/239‑1 (MQ=255)
                                                                                                                                                                                                              gTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTa                    >  1:95039/1‑237 (MQ=255)
                                                                                                                                                                                                              gTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGGGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACCAGGTGAGACG‑CATTTTCGGC‑GAATGGCGCC‑C‑CAAGGCAAATGAAAGCACGCCCGCGCCCCAGCCGACATcctggctctttt                   >  1:4251/1‑228 (MQ=255)
                                                                                                                                                                                                                      cccATGGGGGAACGGCGGGTT‑GGAGATGATCATGTC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCCCCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATcc                             <  1:625293/220‑1 (MQ=255)
                                                                                                                                                                                                                       cctgtgtgtAACGGGGGGTT‑GGGTAGGACAATATC‑AAAACGTCCTTT‑CCCCT‑CGG‑AAAAGAAGTTTCTGGCAAAGGCTTCCCCCTC‑AACCCCGTTGGCCGCAAGTGTTGGGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCAAAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg         <  2:157097/234‑1 (MQ=255)
                                                                                                                                                                                                                        cATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CCCCG‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc        <  1:856072/239‑1 (MQ=255)
                                                                                                                                                                                                                        cATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑GAAAGACGTTGCTGGCCAAGACTTCACCTTC‑AACACCCTTGGCCGCAAGGGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCCGAG‑ACATCGCACAAGGTGGGACG‑AAGTTTCGTC‑GAATGGGGCG‑C‑CAAGGCAACTGCAAAACCCCCCCCCCCACCCGCCACATCCACCACCTTACCATTctcgtgg        >  2:642747/1‑233 (MQ=255)
                                                                                                                                                                                                                         catgggggaCGGGGGGTTGGGGGATGAACATATC‑AAAACTCCCTTC‑CGCCT‑CGG‑AAAAGTCGTGGCGGGAAAGGACTGCCCCTCC‑ACCACCGTTGGCCGAAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTGCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc        <  2:796078/231‑1 (MQ=255)
                                                                                                                                                                                                                            gTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CCCCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAg    <  1:21866/239‑1 (MQ=255)
                                                                                                                                                                                                                               gAACGGGGGGTT‑GGAAATGAACAAATC‑AAAACGCCCTTT‑CCACT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTGGGCCGCAAGTGTGGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGtt  <  1:839849/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
ATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTT‑GGAGATGATCATATC‑AAAACGACCTTT‑CACCT‑CGG‑AAAAGACGTTGCTGGCAAAGACTTCACCTTC‑AACACCGTTGGCCGCAAGTGTTGCGCGGCTGGC‑TTCTACCGCCGGGGCAGAG‑ACATCGCACAAGGTGAGACG‑AATTTTCGGC‑GAATGGCGCG‑C‑AAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTT  >  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: ‑