Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,606,903 A→G 26.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,9030AG26.7% 94.4 / 7.1 58S267P (TCC→CCC) rsmC16S rRNA m(2)G1207 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (9/32);  new base G (14/3);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 5.16e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTC  >  NC_000913/4606665‑4607132
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
gaTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCCTCGTGGAACGGCGGGTTGga                                                                                                                                                                                                                                                >  2:69892/1‑239 (MQ=255)
  aatttAACCTTTCTTCCCCGGGGGGGGTCTAAGGGCGGGATCCCCCTTGAAGCGGCCTGTTGGCGCGATCCCTTCGTGGAAGCCAAATGTCTCATCCAGCCCGGCCGGATAAGGCAGGAAGGCGGTCGCTAAAAGTCGCAGCTGGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                               <  1:195090/235‑1 (MQ=255)
          cTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTTCAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCCGGCTGGTTTTCCTCCCAACGGGGAACGGGGGCGTGGGgatgtgcctt                                                                                                                                                                                                                                      >  1:898359/1‑233 (MQ=255)
             tcttcGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCGCCATCCCGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGgatgagtcaatca                                                                                                                                                                                                                                   >  2:59990/1‑231 (MQ=255)
              cttcGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGTTGGGGAGTAT‑Catattca                                                                                                                                                                                                                                 >  2:608157/1‑236 (MQ=255)
                cccgctgtggggTCATAATGGCGCGATCCCCCTTGAAGCCCCCTGTTTGCCCGACCACTTCGGGGAAGCCAAATGCCTCATCCAGCCCGTCCGGATAAGGCAGGAAGGCGTCCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCaaaa                                                                                                                                                                                                                               <  1:967400/230‑1 (MQ=255)
                             tttAATGGCGCGTTACCCCTTTAAGCGCCCGTTTTGCGGGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCTTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCAcc                                                                                                                                                                                                                 <  1:98782/236‑1 (MQ=255)
                             taTAATGGGGCGATACTCCTTGAAGCGCCCTGTTTGCGCGGTCCCTTCGGGGAAGCCAAAGGTCTAACCCAGCACGTCCGGATAAGGCAGGAAGGCGGTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCAcc                                                                                                                                                                                                                 <  1:370407/237‑1 (MQ=255)
                              aTAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCAACGGGGAACGGCGGGGTGGGGGAGATACTTATCaaaaaaacattcacccc                                                                                                                                                                                                                  >  1:168104/1‑226 (MQ=255)
                                     cgcgATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGAGTTTGCGCCCCATCCAGGCTGGTTTGCATCCCCTCGGGGAACGCCGGGGTGGGGGTGAT‑CAGATCaaaaacacacttcccccccggaaag                                                                                                                                                                                                          >  1:987714/1‑218 (MQ=255)
                                           accccTTGAAGCCCCCTGTTTGCGCGATCGCTTCGGGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTg                                                                                                                                                                                                 <  2:999406/236‑1 (MQ=255)
                                            ccccTTGAAGCGCCCTGTTTGCGCGATCACTTCGGGGAAGCCAAATGTCTCATCCAGCCCGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGc                                                                                                                                                                                                <  2:898080/237‑1 (MQ=255)
                                                 tggaccgcctcttttGCGCCATCCCTTCGTGGAACCCAAATTTCTCATCCAGCACGCCCGGATATGGATGGAAGGCGTTCGCGACAATTCCCAGCTCGCCGGCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCGCCTCGGAAAAG‑ACGTTGC‑TGGCa                                                                                                                                                                                          <  1:377598/227‑1 (MQ=255)
                                                   aaGCGCCCTGTTTGCGGGTTCACTTCGGGGAACGCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATGCGCAGCTCGCCGCCGCTATTAAGATGACGCCCCGCACCGCGAATCAGGGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCaaa                                                                                                                                                                                        <  1:220485/239‑1 (MQ=255)
                                                        cccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCAACCCATCGTGGAACGGCGGGTTGGAGATGAT‑CACATCACAAC‑GG‑CCTTTCACCCCGGGAAAG‑ACGGTGT‑TGGCcaaaaac                                                                                                                                                                                    >  2:953450/1‑234 (MQ=255)
                                                         ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCCCCGCGAATCAGCGTTTGCCCCGCATCCAGGCTGGTTTGCATCCCATCGTGGACCGGCGGGTGTGGGATGAT‑TATTTCAACAA‑GACCCTTCCCCCCCGGaaaaaaggtgccgggaaaaacat                                                                                                                                                                                      >  1:1016778/1‑219 (MQ=255)
                                                         ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGAT‑CCTATCCAAAA‑GA‑CCTTTCCCCCCGGGAAAG‑ACGTTGC‑TGGCAAAAAATTc                                                                                                                                                                                  >  2:254971/1‑239 (MQ=255)
                                                         ccTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGCTGGAGATGAT‑CATTTCAAAAC‑GA‑CACTTTACCCCGGGAAAA‑ACGATGC‑TGGCAAAGACCTc                                                                                                                                                                                  >  2:495453/1‑239 (MQ=255)
                                                           tttttgggggttaccttgggggcaGCCAAAGGTCTCATCCAGCACTTCCGGATAAGGCGGGAAGGCGTTCGCTCCGATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑c                                                                                                                                                                               <  1:528439/216‑1 (MQ=255)
                                                           tttttGGGCGATTCCTTTGGGAAGGCAAAATGTTTCAGCCAGCACGCCCGGAAAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑c                                                                                                                                                                               <  1:648191/237‑1 (MQ=255)
                                                                     aTCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATTAT‑CATTTCACAAC‑GA‑CCTTTTCCCCCGGAAAAA‑ACGTTGC‑TGGGAAAGACTTAC‑CCTTCC‑acac                                                                                                                                                                     >  2:388093/1‑238 (MQ=255)
                                                                       cACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                               >  2:416620/1‑169 (MQ=255)
                                                                       cACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgag                                                                                                                                                                                                                                               <  1:416620/169‑1 (MQ=255)
                                                                           tggTGGAAGCCAAATTTCCAATCCACGCCGTCCGGATAAGGAAGGAAGGGGTTCGCTGCAATTCGCAGTCCCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGc                                                                                                                                                              <  2:748183/237‑1 (MQ=255)
                                                                              tGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGCTTGCATCCCATCGTGGAACGGCGGGGTGGGGCGGTT‑ACTATCAAAAC‑CA‑CCCTTCACCCCCGACACA‑ACGCTGCGTGGCAAA‑ACCTCACCCCTCA‑AAACCCTTGGcgc                                                                                                                                                            >  1:956637/1‑237 (MQ=255)
                                                                                      aaaTGTCTCAGCCAGCACGGCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTg                                                                                                                                                   <  2:701680/239‑1 (MQ=255)
                                                                                       aaTGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGGTGAT‑CATACCAAAAA‑CA‑CCCTTCTACTCGGAAAAG‑AAGTTGA‑TGGGAAAGACTTCA‑ACCTCA‑ACAACggcgaccgccccagtagg                                                                                                                                                  >  1:844909/1‑222 (MQ=255)
                                                                                             tcATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGTGTGGGCGTGAA‑CCTATTAAAAA‑AACCCTTCCACCCCGAAAAAG‑A‑GTGTC‑GCGCACAAACCTCA‑CCCTCC‑ACACCCATGGCGGCaaaggggggcgggcg                                                                                                                                            >  1:864438/1‑226 (MQ=255)
                                                                                                tCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCt                                                                                                                                                                                                                <  2:980437/172‑1 (MQ=255)
                                                                                                    gcccGTCCGGATAAGGCAGGAATGCGTTCGCTACAATTCCGATCCAGCCGCCGCTATTAAGATGCCCCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCt                                                                                                                                     <  2:671176/236‑1 (MQ=255)
                                                                                                     cACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGAT‑TCTATCAAAAC‑GC‑CCTTTCACCTCGGGAAAG‑ACGTTGC‑TGGCAAAAACTTCA‑CCTTCA‑AAACCCTTGGGCGCAAAGGATGCGCGGCTGGGTTCt                                                                                                                                     >  2:114369/1‑238 (MQ=255)
                                                                                                         ttcggaaagggcgaggaGGCGGTCGCAACAATTCGTACGCCGCCGCCGCTATTAAGATGACGTACCGCCCCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCAATCGAGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGGTGC‑TGGCAAAGACTTCA‑CCTGCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgc                                                                                                                                <  1:608157/223‑1 (MQ=255)
                                                                                                          ccGGATAAGGCAGGAAGGCTTTCCCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTAccgcc                                                                                                                               <  2:168104/239‑1 (MQ=255)
                                                                                                                 aGGCAGGAAGGGGTTCGCTACAATTCGAAGCCCCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCa‑g                                                                                                                       <  2:924349/239‑1 (MQ=255)
                                                                                                                  ggCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGAT‑CGATTCAAAACACA‑CTTTCCACCT‑GGAAAAGAAGGTTGC‑T‑GCGAAAACTTCC‑CCTTCC‑ACACCACTGGCCGCCAGAGTGGGGCgggtggccatcaccccccgggcaag                                                                                                                       >  1:306598/1‑221 (MQ=255)
                                                                                                                       gAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGGTGC‑TGGCAAAGACTTCA‑CCTTCA‑AAACCGTTGGCCGCAAGTGTTGCGCGGCTGGGTTTTACCGCCCGGGCC‑GAGACa                                                                                                                  >  2:207327/1‑238 (MQ=255)
                                                                                                                          ggCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGAGAGGAT‑CCTCTCAAAAA‑GA‑CCCTTCAACCCCGAAAAA‑AAGTTGG‑TGGCAAAGACTTCA‑ACTCCA‑AAACCCTTGTCCGCACAGGATTCGCGGCGGGGTTCTCCCCCCcgggggacgaccaccc                                                                                                               >  1:668236/1‑227 (MQ=255)
                                                                                                                            cggtCGCTACAATTCGCAGCCCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGca                                                                                                             <  2:440127/235‑1 (MQ=255)
                                                                                                                            cGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTc                                                                                                                                                                                  <  2:713659/172‑1 (MQ=255)
                                                                                                                            cGTTCCCTAAAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGca                                                                                                             <  2:956637/238‑1 (MQ=255)
                                                                                                                              ttCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGCAAAG‑ACGTTGC‑TGGGAAAGACTTCC‑ACCTCA‑ACACCATTGGCCGCAAGGGTTTGCCGGCTGGCCTCTACCCCCGGGGCA‑GAGACaacccgaaa                                                                                                          >  1:826954/1‑231 (MQ=255)
                                                                                                                                 gCTACAATTCGCCGCTGCCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGt                                                                                                      <  2:826954/239‑1 (MQ=255)
                                                                                                                                                                  tcccgccccgcccCCCGAAGCACGGTTGGCGCCGCAGCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATAACAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGCCTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAgg                                                                      <  1:114369/235‑1 (MQ=255)
                                                                                                                                                                    aCGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGAGATGAT‑CATATCAAAAA‑GA‑CCTTTCACCCCGGAAAAG‑ACGTTGC‑TGGCCAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGGGTTGCGCGGCGGGCTTCTACCGCCCGGGCA‑AAGAAACCGCCC‑AAAGGGAGACGAAGTTTTGGCGAATGGCGGGCGAAGGCaa                                                                   >  1:225453/1‑239 (MQ=255)
                                                                                                                                                                                 cgAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGAT‑CATATCCAAAC‑GA‑CCTTTCACCTCCGGAAAG‑ACGTTGC‑TTGCAAAGAATTCA‑CCTTCACCCCCCTTTGCCCGCAAGGGTT‑TGCGGCTGGCTTCTAACCACCGGGCA‑GGGGCAACGCAC‑AAGGGGGAGCGAGATTTTGGCGAACAGCGGGCaaaaaacaccgaaaaaaaccc                                                      >  1:596915/1‑221 (MQ=255)
                                                                                                                                                                                  gAACCAGGGTTGGCGCCGCATCCGGGCTGGTTTGCACCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGcc                                                      <  2:987714/238‑1 (MQ=255)
                                                                                                                                                                                         ggTTTGCGCCTCAACCGGGCTGGGTTGCCTCCCATCGTGGACCGGGGGGTGGGGGATGGA‑CATATCAAAAC‑GG‑CCTTTCCCCTGGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTCCA‑ACACCGTGGGCCGCAAGTGTTGCGCGGCTGGCTTCTCCCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCcac                                              <  1:897351/238‑1 (MQ=255)
                                                                                                                                                                                         cGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAA‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑AAACCGTTGGCCGCAAGGGTTGCGCGGCTGGCTTCTACCGCCCGGGCA‑GAGACATCGCAC‑AAGGGGGGACGAATTTTCGGCGAATGGCGCGCCAAGGCAACTTAAAgaaacccccgccccc                                              >  2:756928/1‑225 (MQ=255)
                                                                                                                                                                                              tgccccgctCCGGGGTGGTTTGCATCCAAGCGTGGAAGGGGGGGTTGGGGATGAT‑AAAATCAAAAC‑GA‑CCTTTCCCCTCGGAAAAG‑ACGTGGC‑GGGCAAAGCCTTCC‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATTGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                         <  1:388093/231‑1 (MQ=255)
                                                                                                                                                                                              gcgcCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑TG‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCg                                         <  2:306598/239‑1 (MQ=255)
                                                                                                                                                                                                   gcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGGGGGTTGGAGATGAT‑CAGATCAAAAC‑GA‑CCTTGCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTGC‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATc                                    <  2:919470/239‑1 (MQ=255)
                                                                                                                                                                                                     aaccAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATc                                    <  2:5702/235‑1 (MQ=255)
                                                                                                                                                                                                     aTCCGGGCTGGTTTGCACCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCa                                  <  2:596915/239‑1 (MQ=255)
                                                                                                                                                                                                            cttgTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATACCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACtt                            <  2:646286/235‑1 (MQ=255)
                                                                                                                                                                                                               gTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑AAGGTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGGTGCGCGGCTGGCTTCTCCCGCCGGGGCAGGAGAC‑TCGCCCAAAGGGGAG‑CGAAATTTCGGCGACAGGGCCGCAAACGGAACCGAAAACAGCCCCGCCCCCCAccccgacaccacccccttac                         >  1:602453/1‑221 (MQ=255)
                                                                                                                                                                                                                 ttGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACCCCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCt                       <  2:579132/238‑1 (MQ=255)
                                                                                                                                                                                                                        ccATCGTGGAACGGGGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGtgtg               <  2:1016778/239‑1 (MQ=255)
                                                                                                                                                                                                                         catgggggAAGGGGGGGTTGGGAATGAA‑CATATCAAAAC‑GG‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCC‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc              <  2:429674/233‑1 (MQ=255)
                                                                                                                                                                                                                         cATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGc              <  2:225453/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cggGTTGGAGATGAT‑CATAGCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAACACTTCA‑CCTACA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTGTCGTGTGCGGAGTTAACGTc  <  2:237373/238‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCa‑g                                                                                                                       <  2:25323/120‑1 (MQ=255)
                                                                                                                                                                                                                                        ggTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCa‑g                                                                                                                       >  1:25323/1‑120 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
GATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTCTCATCCAGCACGTCCGGATAAGGCAGGAAGGCGTTCGCTACAATTCGCAGCTCGCCGCCGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGAT‑CATATCAAAAC‑GA‑CCTTTCACCTCGGAAAAG‑ACGTTGC‑TGGCAAAGACTTCA‑CCTTCA‑ACACCGTTGGCCGCAAGTGTTGCGCGGCTGGCTTCTACCGCCGGGGCA‑GAGACATCGCAC‑AAGGTGAGACGAATTTTCGGCGAATGGCGCGCAAAGGCAACTGAAAGCACGCCCGCGCCACAGCCGACATCCAGCACTTTACCTTTCGTGTGCGGAGTTAACGTC  >  NC_000913/4606665‑4607132

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: ‑