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%	59.0 / 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)

CACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGC  >  NC_000913/3941630‑3942118                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAAc                                                                                                                                                                                                                                                 >  1:188527/1‑251 (MQ=14)     ccgctccgctGAGAAAAAGCAAAGCGGCACTGCGCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATGCTTAACGTCGCAAGACGAAAAATGAAGACCAAGTCTCAAGAGTGAACACGGAATTCATTACGAAGTTTAATTCTTTGAGCATCAATCTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGa                                                                                                                                                                                                                                             <  1:81182/251‑1 (MQ=11)                 gaAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTg                                                                                                                                                                                                                                 <  2:223623/251‑1 (MQ=18)                    aaagcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGttt                                                                                                                                                                                                                              >  1:389696/1‑251 (MQ=21)                     aaacaaagcGGCACTGCTCTTTAACAAATTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTc                                                                                                                                                                                                                             <  2:430553/248‑1 (MQ=18)                      agcaaagcGGCACTGCTCTTTAACAATTTATCAGACAATCTGTATGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGAGAACACGTAATTCATTACGAAGTTTAATTCTTTGCGCATCAAACTTTTAAATTACAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGCATGGTAACAGGCAACAACTTGCTGTTTCg                                                                                                                                                                                                                            >  1:230989/1‑251 (MQ=11)                               gcactgcactTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATGCTTAACGACGCAAGACGAAAAATGAATGCCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAGACTTTTAAATTGAAGAGGTTGAGCATGGCTCAGGTTGAACGGTGGCGGCAGGCCTAAGACATGCAAGTCGAATGGTAGCAGGAAACAGCTTGCTGTTTCGCTGACGAGt                                                                                                                                                                                                                   <  1:178578/251‑1 (MQ=11)                                                                                                                                                      aaGTTTAATTCTTTGAGCATCAAACTTTTACATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACCTGCAAGTCGAATGGTAACAGGACACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTATGGGAAACTGCCTGATGGAGGGGGGAAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGTGACCTTCGGGCCTCTTGcc                                                                                           >  1:380249/1‑251 (MQ=11)                                                                                                                                                                     aGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGCTGGGccc                                                                              >  1:384748/1‑249 (MQ=17)                                                                                                                                                                       cATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGGAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGTGCCTCTTGCCATCAGATGTGccc                                                                              >  2:254912/1‑247 (MQ=18)                                                                                                                                                                        aTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTATGGGAACCTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATg                                                                         >  1:296821/1‑251 (MQ=18)                                                                                                                                                                        aTCAAAATTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCCTCTGATTTGCCCCGATg                                                                         >  2:204186/1‑251 (MQ=255)                                                                                                                                                                           aaaCTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAAAGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCAAGATggg                                                                       >  1:114389/1‑250 (MQ=18)                                                                                                                                                                             aCTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTGCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGAtt                                                                    <  2:355693/251‑1 (MQ=18)                                                                                                                                                                                   aaaTTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATTGGATTAGCt                                                                >  2:80253/1‑249 (MQ=18)                                                                                                                                                                                     aTTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGGAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTCCTGGCACCGGTAGCTAATACCGCATAACGTCGCACGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAgatcgtt                                                            >  2:32681/1‑245 (MQ=11)                                                                                                                                                                                       tGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGGAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAgg                                                          >  2:266913/1‑251 (MQ=25)                                                                                                                                                                                        gAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGCAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGAAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGGCCTTCGGGCCTCTTGCCCTCAGATGTGCCCCGATGGGATTAGCTAGTAggg                                                         >  2:375468/1‑250 (MQ=255)                                                                                                                                                                                        gAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGCCGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAAATACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAggg                                                         >  1:35985/1‑250 (MQ=255)                                                                                                                                                                                             gTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGt                                                    <  1:72855/251‑1 (MQ=37)                                                                                                                                                                                                    cATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGACGAATGGTAACAAAAACCAGCTTGCTGTTTCGCTGACGAGTGGCGGAAGGGTGAGTAATGTCTGGGATACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTCGCTAATACCGCATAAAGGCGCAAGACCCAAAG‑GGGGGACCTTCGGGCCTATTGCCATCAGATGTGCCCAGATGGGTTTAGCTAGTAGGTGGGGTAACGGcc                                             >  2:428816/1‑250 (MQ=11)                                                                                                                                                                                                     aTGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGCCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGGTGGGCCCCGATGGGATTAGCTAGTAGGGGGGGTAACGGcac                                            >  1:28949/1‑249 (MQ=11)                                                                                                                                                                                                                    cGCTGGCGGCAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCCGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATc                             >  1:385762/1‑251 (MQ=37)                                                                                                                                                                                                                          gggcAGGCCTAACACATGCAAGGCGAATGGTAACAGGAAACAGCGTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCAGAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATTTGCCCAGATGGTATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGTCGATTCCTAGc                       <  2:170889/250‑1 (MQ=14)                                                                                                                                                                                                                           ggcAGGCCTAACACATGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCt                      >  1:287720/1‑251 (MQ=38)                                                                                                                                                                                                                                          tGCAAGTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGCGAGGAATGTCTGGGAAACCGCCTGATGGAGGGGGATAACTACCGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGCGGGGGCCCTTCGGGCCTCCTGCCCTCAGCTGTGCCCCGATGGGATTCGCTAGTAGGGGGGGTAACGGCCCACCCAGGCGACGATCCCTCGCCGGCCCGAGAGGacg        >  2:374589/1‑248 (MQ=11)                                                                                                                                                                                                                                               gTCGAATGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc  >  1:366893/1‑251 (MQ=38)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAACAGCTTGCTGTTTCGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACC‑AAAGAGGGGGACCTTCGGGCCTCTTGCCATCAGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGC  >  NC_000913/3941630‑3942118

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