Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,452,218	G→A	F23F (TTC→TTT)	rplD ←	50S ribosomal subunit protein L4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,452,218	0	G	A	100.0%	32.6 / NA	13	F23F (TTC→TTT)	rplD	50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/5);  total (8/5)

CTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTATCTCCTTACGCCTTCACAGCTGGTTTAACGATCAGGTCGCTACCGGTTGCACCCGGGACAGCACCTTTAACCAGCAGCAGGTTGCGCTCAGCGTCAACGCGT  >  NC_000913/3452009‑3452393                                                                                                                                                                                                                  |                                                                                                                                                                                ntgtgGTCCTGCGAACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGct                                                                                                                                        >  1:273732/2‑251 (MQ=255)                           ggTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGGCTTCTGATCACGAGTGCCCGGACGAGCACCATCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTc                                                                                                                                                <  1:389461/217‑1 (MQ=255)                           ggTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGCGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATACGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAgcg                                                                                                                                            >  1:66061/1‑221 (MQ=255)                                                                               cGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAACTCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTATCTCCTTACGCCTTCACAGCTGGTTTAACGATCAGGTCGCTACCgg                                                          >  1:68037/1‑251 (MQ=255)                                                                                        ccTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAAGTCCATTGCTATCTcc                                                                                                   <  1:120060/201‑1 (MQ=255)                                                                                        ccTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTAAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCCTAAGCAACAACAAACTGGTGACCCAGCGCTTCGTTAAAATCAc                                                                                                                                                                          >  1:70624/1‑130 (MQ=255)                                                                                                                           tACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGAGCTTAGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTAt                                                                                                       >  1:74057/1‑162 (MQ=255)                                                                                                                                        gTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTATCTCCTTACGCCTTCACAGCTGGTTTAACGATCAGGTCGCTACCGGTTGCACCCGGGACAGCACCTTTAACCAGCAGCAGGTTGCGCTCAGCGTCAACGCGt  <  1:290100/250‑1 (MQ=255)                                                                                                                                                  nnnnnnnnnncTNNNNNGCANCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCGTTCAATACTAATTCCATTGCTATCTCCTTACGCCTTca                                                                                        <  2:273732/144‑1 (MQ=255)                                                                                                                                                       ntACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATc                                                                                                                                                                            >  1:42984/2‑65 (MQ=255)                                                                                                                                                                       cAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTATCTCCTTACGCCTTCACAGCTGGTTTAACGATCAg                                                                     >  1:207669/1‑152 (MQ=255)                                                                                                                                                                                    aacaacaacCTGGTGAACCAGCGCTTCGTTAAAATCACGCCCGAAGGTAGTTTAGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTATCTCCTTACGCCTTCACCGCTGGTTTacc                                                                           >  1:301481/1‑131 (MQ=255)                                                                                                                                                                                                gTGAACCAGCGCTTCGTTAAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTAGCTCCTTACGCCTTCACAGCTGGTTTAACGATCAGGTCGCTACCGGTTg                                                       <  1:156136/141‑1 (MQ=255)                                                                                                                                                                                                                  |                                                                                                                                                                                CTGTGGTCCTGCGGACGAGCAGCAAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTATCTCCTTACGCCTTCACAGCTGGTTTAACGATCAGGTCGCTACCGGTTGCACCCGGGACAGCACCTTTAACCAGCAGCAGGTTGCGCTCAGCGTCAACGCGT  >  NC_000913/3452009‑3452393

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