Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,799,883	T→C	K39K (AAA→AAG)	rpmI ←	50S ribosomal subunit protein L35

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,799,883	0	T	C	100.0%	122.4 / NA	35	K39K (AAA→AAG)	rpmI	50S ribosomal subunit protein L35
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/17);  total (18/17)

CTCTCCTGTATCTATATTCTAATTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGTTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAGCGCTTAGCAGCACCGCGTACGGTCTTAATTTTT  >  NC_000913/1799755‑1799994                                                                                                                                 |                                                                                                                ctctCCTGTATCTATATTCTAATTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCtt                                                                                                       >  2:1096192/1‑139 (MQ=255)              atatTCTAATTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAAt                                                                                          >  1:970600/1‑139 (MQ=255)                atTCTAATTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAAtgt                                                                                        >  2:901765/1‑139 (MQ=255)                   cTAATTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCttttt                                                                                                    <  1:1080484/124‑1 (MQ=255)                   cTAATTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCttttt                                                                                                    >  2:1080484/1‑124 (MQ=255)                     aaTTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGc                                                                                   >  1:898449/1‑139 (MQ=255)                      aTTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCa                                                                                  >  2:1113025/1‑139 (MQ=255)                      aTTAAAAAGTTAAAAACGTTAACGGCTGATGCGTACGGCGGGCAGGCGGTTACCAGGCCCAGATCGCCGTTGGAAACCAGGGCTTTCGGACGCAGGTGACAGTTACGCTTGGTCGCTTTATTGTTCAGGATGTGACGCa                                                                                  >  2:468822/1‑139 (MQ=255)                       ttAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCttttt                                                                                                    >  1:414228/1‑120 (MQ=255)                       ttAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCttttt                                                                                                    <  2:414228/120‑1 (MQ=255)                                   aaaCGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCgtgctt                                                                     <  2:898449/139‑1 (MQ=255)                                          aaCGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTa                                                              <  1:1113025/139‑1 (MQ=255)                                                 tATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAaccac                                                       >  2:280473/1‑139 (MQ=255)                                                       tACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGATGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAAtg                                                                                         <  1:702399/99‑1 (MQ=255)                                                       tACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGATGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAAtg                                                                                         >  2:702399/1‑99 (MQ=255)                                                         cGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCg                                               >  2:902205/1‑139 (MQ=255)                                                         cGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCg                                               >  1:70166/1‑139 (MQ=255)                                                            caggcaCGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGtt                                            >  1:566886/1‑139 (MQ=255)                                                               gcaCGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGc                                                                                   <  2:720148/97‑1 (MQ=255)                                                               gcaCGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGc                                                                                   >  1:720148/1‑97 (MQ=255)                                                               gcaCGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCgtgctt                                                                     <  1:873977/111‑1 (MQ=255)                                                               gcaCGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCgtgctt                                                                     >  2:873977/1‑111 (MQ=255)                                                                    cgATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGaa                                    <  2:70166/139‑1 (MQ=255)                                                                       ttACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAgcg                                 <  1:280473/139‑1 (MQ=255)                                                                       ttACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAgcg                                 <  1:902205/139‑1 (MQ=255)                                                                       ttACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAgcg                                 <  1:901765/139‑1 (MQ=255)                                                                       ttACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAgcg                                 <  2:566886/139‑1 (MQ=255)                                                                               cccAGATCTCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGAc                                                                                     >  1:723687/1‑79 (MQ=255)                                                                               cccAGATCGCCTTTGTAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGAc                                                                                     <  2:723687/79‑1 (MQ=255)                                                                                        ccTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCgtgc                                                                       >  1:118447/1‑84 (MQ=255)                                                                                        ccTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCgtgc                                                                       <  2:118447/84‑1 (MQ=255)                                                                                        ccTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAGCGCTTAGCAGCACCGCGTa                <  2:1094065/139‑1 (MQ=255)                                                                                         cTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAGCGCTTAGCAGCACCGCGTAc               <  1:1096192/139‑1 (MQ=255)                                                                                                   cATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAGCGCTTAgcagca                       >  1:994323/1‑121 (MQ=255)                                                                                                   cATGGCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAGCGCTTAgcagca                       <  2:994323/121‑1 (MQ=255)                                                                                                      ggCTTTCGGACGCAGGTGACGTTTACGCTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAGCGCTTAGCAGCACCGCGTACGGTCTTAAttttt  >  1:649484/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                CTCTCCTGTATCTATATTCTAATTAAAAAGTTAAAAACGTTAACGGCTTATGCGTACGGCAGGCACGCGATTACCAGGCCCAGATCGCCTTTGGAAACCATGGCTTTCGGACGCAGGTGACGTTTACGTTTGGTCGCTTTTTTGGTCAGAATGTGACGCAGGTTAGCGTGCTTGTGCTTAAAACCACCTTTACCGGTTTTTTTGAAGCGCTTAGCAGCACCGCGTACGGTCTTAATTTTT  >  NC_000913/1799755‑1799994

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