Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,452,160	T→C	T43A (ACT→GCT)	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,160	0	T	C	100.0%	53.1 / NA	16	T43A (ACT→GCT)	rplD	50S ribosomal subunit protein L4
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/8);  total (8/8)

AAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTAT  >  NC_000913/3452032‑3452292                                                                                                                                 |                                                                                                                                     aaaGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAg                                                                                                                           <  1:404764/140‑1 (MQ=255)       cACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAg                                                                                                                     <  1:35265/140‑1 (MQ=255)          gCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTg                                                                                                                  >  1:55088/1‑140 (MQ=255)                    cGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGcaac                                                                                                        >  1:323811/1‑140 (MQ=255)                        aGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGcaacaaca                                                                                                    <  1:274801/140‑1 (MQ=255)                                     gATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCa                                                                                       >  1:336714/1‑140 (MQ=255)                                               gAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGtt                                                                             <  2:72033/140‑1 (MQ=255)                                                    cgcACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCTCGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAAt                                                                        >  2:17918/1‑140 (MQ=255)                                                      cACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCa                                                                      <  2:323811/140‑1 (MQ=255)                                                                         gcgcCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCgg                                                   >  2:285679/1‑140 (MQ=255)                                                                             cACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAAc                                               <  2:176612/140‑1 (MQ=255)                                                                               cGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAg                                             <  2:335933/140‑1 (MQ=255)                                                                                          cGGAACCAGTTACTTCAGCATGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAgc                                        <  1:339629/134‑1 (MQ=255)                                                                                          cGGAACCAGTTACTTCAGCATGAGTCTTCTGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAgc                                        >  2:339629/1‑134 (MQ=255)                                                                                                                      ctGAGCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTg      >  1:211588/1‑140 (MQ=255)                                                                                                                          gCACGAGCACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTAt  >  2:197737/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                     AAAGGTCACGCCACCAGAACGCCAGATCGGGCTCTTGATAGAACCAGAACGCGCACGGCCGGTGCCTTTCTGGCGCCACGGTTTTTTACCGGAACCAGTTACTTCAGCACGAGTCTTCTGAGCACGAGTACCCTGACGAGCACCAGCTGCATAAGCAACAACAACCTGGTGAACCAGCGCTTCGTTGAAATCACGACCGAAGGTAGTTTCGGAAACAGTCAGCGCGCTCTGCGCGTCTTTCAATACTAATTCCATTGCTAT  >  NC_000913/3452032‑3452292

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