Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,874	C→T	noncoding (67/1542 nt)	rrsC →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,874	0	C	T	100.0%	51.8 / NA	27	noncoding (67/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/7);  total (20/7)

GCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCC  >  NC_000913/3941657‑3942119                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      gcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACg                                                                                                                                                                                                                      <  2:442624/251‑1 (MQ=14)  cGGCACTGCTCTTTAACTATCTATCAGACAATCTGCGAGGGCACTCGAAGATCCGGATTCTTAAAGTCGCAAGACGAAAAATGAAAACCAAATCTCACGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCACACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAAAGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACg                                                                                                                                                                                                                      >  2:582760/1‑250 (MQ=11)  cGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGCAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGACGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACg                                                                                                                                                                                                                      >  1:346342/1‑250 (MQ=14)  cGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACg                                                                                                                                                                                                                      >  1:141302/1‑250 (MQ=14)                                                                                                                aatgcATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGAGTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCTCATAACGTCGCAAGATCCAAGAGCTGGACCTTc                                                                                                        <  2:386783/246‑1 (MQ=11)                                                                                                                aaTTCATTACGAAGTTTAATTCTTTGAGCAGCAAACTGTTAAATTGAAGAGTGTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGATGAGTGGCGGACGGGTGAGTAATGTCTGGGAATCTGCCGGATGGAGGGGGATAACAGCTGGAAACGGGAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGCGGCCCTTCg                                                                                                       <  2:276554/251‑1 (MQ=11)                                                                                                                   tCATTACGAAGTTTAATTCTTTGAGCATCACACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAAAAAATGCAAGTCGAATGGTACCAGGCAACAGCTTGCTGTTTCGCTGACGAGTGGCGGTAGGGTGCGTAATGTCTCGGCAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCACGAACAAGGAGGAGGACCTTCgcgc                                                                                                    >  1:274501/1‑248 (MQ=11)                                                                                                                         cGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGAGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGGTAACTACTGGAAACGGAAGCTAATACCGCATAACGTCGCAAGGCCAAAGAGGGGGACCTTCGGGCCTCTTg                                                                                              >  1:375411/1‑251 (MQ=9)                                                                                                                                  aTTCTTTGAGCATCAAACTTTTAAATTGAAGAGATTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGAAGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAAAGTCGCTAGACAAAAGAGGGGGACCTTCGGGCCTCTTGcactccgat                                                                                     >  2:381153/1‑243 (MQ=11)                                                                                                                                    tCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTTAACGCTAACGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGATGATTCGCTGACGAGTGGCGGAAGGGGGAGTAATGTCTAGGAAACTGCCTGATGGAGGAGGATAAACACTGGAAACGGTAGCTAATACCGCATAACGTCGGAAGACCAAAGAGAGGGCCCATCCGGCCTCTTGCCATCAGAtgt                                                                                   >  1:342639/1‑251 (MQ=11)                                                                                                                                               cAAACTTTTACATTGAAGAGTTTGCTCATGGCTCAGATTGAACGCTGGCGGCAGGACTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGCCTGGTGGAGGGGGATAACTACTGGAAACGGACGCTAATACCGCCTAACGTCGCAAGAACAAAGAGGGGGGCCTTCGGGCCTCTTGCCATCAGGTGTGCCCAGATggg                                                                        >  2:421139/1‑251 (MQ=11)                                                                                                                                               cAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGCACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGAACAAAGCGGGGGACCTTCg                                                                                                       >  2:98379/1‑220 (MQ=9)                                                                                                                                               cAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCg                                                                                                       <  1:98379/220‑1 (MQ=14)                                                                                                                                                aaaCTTTTAAATTGAAGAGATTGCTCATGGCTCAGATTGAACGCTGGCGACAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTATCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCATGAGGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGACCCTTGCCATCTGCTGTGcaccaatgggg                                                                       >  2:356449/1‑241 (MQ=11)                                                                                                                                                 aaCTTTTAAATTGAAGCGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGAt                                                                      >  1:497835/1‑251 (MQ=18)                                                                                                                                                        taaTTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACAAATGCAAGTCGCATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGACACTGCCTGATGGAGGGGGATAAATACTGGAAACGGTAGCTAATACCGCATAACGTCGAAAGACCACAGAGGGGGACCTTAGGGCCTCTTGCCATCAGTTGTGCCCCGATGGGGTTAGCt                                                                 >  2:495719/2‑249 (MQ=11)                                                                                                                                                        aaaTTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCt                                                                 >  1:57400/1‑249 (MQ=18)                                                                                                                                                                gagTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTAATGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGGCCTTCGGGCCTCTTGCCATCAGATGTGCCCAGCTGGGATTAGCTAGTAGGTggg                                                       >  1:357825/1‑251 (MQ=255)                                                                                                                                                                 agTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGGGAGAAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACGACTGGAAACGGTAGCTAATACCGCATAACGTGGCAAGACCAAAGAGGGGGACCTTCGGGCCTTTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGGAGGTgggg                                                      >  2:421528/1‑251 (MQ=11)                                                                                                                                                                  gTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGt                                                     <  1:232560/251‑1 (MQ=37)                                                                                                                                                                       aTCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGCATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACTGGTTAGTAATGTCTGGGAAACTGCCTGCTGGCTGGGGCCAACTACTGGAAACGGTAGCTAATACCGCATTACGTCTCACGACCAAATAGGGGGACCTTCGGGTCTCTTTTCACCCACTGTGCCCCGATGGGATTAGCTCGTAGGTGCGGTAACgg                                                >  2:8546/1‑251 (MQ=11)                                                                                                                                                                         cATGGCTAAGATTGAACGCTGGCGGCAGGCCTAACACATGCACGTCGAATGGTAACAGCAAACAGCTTGCTGTTTCGCTGACGCGTGGCGCACGGGTGAGTAATGTCCGGGAAACTGCCTGCTGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGGAGGCGGGGTAACGGCt                                              >  2:524055/1‑251 (MQ=11)                                                                                                                                                                                 aGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAgg                                      <  1:102284/251‑1 (MQ=38)                                                                                                                                                                                  gATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGc                                     >  1:426165/1‑251 (MQ=38)                                                                                                                                                                                             ggcggcAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCGGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCGCTTGCCATCAGAGGTGCCCAGATGGGATTAGCGAGTAGGTGGGGTAACGGCTCACCGAGGCGACGAGCCCTa                          <  2:346342/251‑1 (MQ=11)                                                                                                                                                                                                                  aaGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAAATGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACAGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGCTGTGCCCAGATGGGATTAGCTGGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGAcc      >  1:354060/1‑250 (MQ=255)                                                                                                                                                                                                                     tCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGAGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGcc  >  1:51575/1‑251 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                      GCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCC  >  NC_000913/3941657‑3942119

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