Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,474,425	T→G	K43T (AAA→ACA)	rpsL ←	30S ribosomal subunit protein S12

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,474,425	0	T	G	100.0%	86.0 / NA	26	K43T (AAA→ACA)	rpsL	30S ribosomal subunit protein S12
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/13);  total (13/13)

CCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATA  >  NC_000913/3474183‑3474663                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               ccTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGAAGAACAAGGCATACTTTACGCAGCGCGGCGTTCGGTTTTGTAGGAGTg                                                                                                                                                                                                                                        >  1:40771/1‑251 (MQ=255)  attAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGCCGAACACGGCATACTTTACGCAGCGAGGGGTTCGGTTTTGTAGGAGTgg                                                                                                                                                                                                                                       >  2:10216/2‑251 (MQ=255)      ggACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACTCCGAAGCAGTCAAGCGCACCACGTCCGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTGCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGgtagt                                                                                                                                                                                                                                   <  2:31585/251‑1 (MQ=255)      ggACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGgtagt                                                                                                                                                                                                                                   <  2:198061/251‑1 (MQ=255)        acttttcccGCCATACTGGGAACGAGCCTGCTTCCGGTCTTTACCGCCGGAGCAGTCAAGCGCCCCACGTACGGTGTGGTAACGAACGCCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGtat                                                                                                                                                                                                                                 <  2:33135/243‑1 (MQ=255)                                                  aaCGCCGCAGCAGTCCAGCGCCCAACGTACGGTGAGGTAACGAACACCCGGGAGGTCTTTAAAACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAACCCGTTAGTCAGCCGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAATACATACGCCACGTTTTTGCGGGAATGCTTcc                                                                                                                                                                                        >  2:233506/1‑250 (MQ=255)                                                     gCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAgcg                                                                                                                                                                                    >  2:85580/1‑251 (MQ=255)                                                               gcAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACCCGACCGCCACGGATCAGGATCACGTAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAATTCACGTCGAAACCGTGAGTCAGACGAACTCGTCATACTTTACGCATCGCGGAGTTCGTTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGtt                                                                                                                                                                          <  1:20346/250‑1 (MQ=255)                                                                   ccccaccacGTCCGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGCCCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCtt                                                                                                                                                                      <  2:46947/248‑1 (MQ=255)                                                                                    tgGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGc                                                                                                                                                     >  1:145934/1‑251 (MQ=255)                                                                                             acacCCGGGGGGTCTTTAACCCGACCGCCCCGGATCATGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGTAAGTCACGTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACTCATCGCGTAGTTCGGTTTTGTATGAGTGGTAGTAGATACGCGAGTACATACGCTACGTTTTTTCGGGCAGGCTTCCAGCGCAGGCGTGTTGCTTTTCGCAACTTTGCTAGCACGTGGttt                                                                                                                                            <  2:246851/251‑1 (MQ=255)                                                                                                             tAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATAAACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGtt                                                                                                                             >  1:118014/1‑250 (MQ=255)                                                                                                                  cgcccGCCACGGATAAGTATCCCGGAGTGCTCCTGTAGGTGGTGACCTTCACCACCGATGTATGAAGTCGCTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTGCGGTTTTGTAGGAGTGGTATTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACTGTGCTTTTCGCAACTTTGCTAGCACGTGGTTTTCGTACCAGCTGTTTAACTGt                                                                                                                       <  2:150127/248‑1 (MQ=255)                                                                                                                    aCCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTg                                                                                                                     >  2:210161/1‑251 (MQ=255)                                                                                                                                 aggatcaCGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCt                                                                                                        >  1:144159/1‑251 (MQ=255)                                                                                                                                   gatcaCGGAGTGATCCTGCAGGTTGTGACCTTCACCACCGCTGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTcc                                                                                                      >  2:16552/1‑251 (MQ=255)                                                                                                                                                      aGGTTTTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTc                                                                                   <  2:200595/251‑1 (MQ=255)                                                                                                                                                                                                ccGTTAGTCAGACGAACCCGGCATACTTTACGCCGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGcc                                         >  1:202543/1‑251 (MQ=255)                                                                                                                                                                                                     aGTCAGACGAACACGGCATTCTTTACGCAGCGCGTAGTTCGGGTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCGTCCAGCGCAAGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTCCTTCGTAAACACGTAATAAAACGGCCTCACACAATATGAGTACGCCGAAt                                     <  2:1660/250‑1 (MQ=255)                                                                                                                                                                                                        cAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTCGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTa                                 <  1:143387/251‑1 (MQ=255)                                                                                                                                                                                                         aGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAAtttt                                  >  1:90636/1‑249 (MQ=255)                                                                                                                                                                                                              aaCACGGCATTCTTTACGCAGCGCTGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATTAGGACGCCGAATTTTAGGGCGa                           <  2:67045/251‑1 (MQ=255)                                                                                                                                                                                                              aaCACGGCATACTTTCCGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCg                            >  2:176051/1‑250 (MQ=255)                                                                                                                                                                                                                  cGGCATACTTTACGCAGCGCGGAGTTCGGGTTTGGAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGGTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTGTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGcc                       <  1:223506/251‑1 (MQ=255)                                                                                                                                                                                                                       tACTTTACGCAGCGCGGAGTTCGGTTTTGTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAACTAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGaaa                   >  2:215981/1‑250 (MQ=255)                                                                                                                                                                                                                                       gAGTTCGGTTTTGTAGGTGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAAtata  <  1:186596/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               CCTTAGGACGCTTCACGCCATACTTGGAACGAGCCTGCTTACGGTCTTTAACGCCGGAGCAGTCAAGCGCACCACGTACGGTGTGGTAACGAACACCCGGGAGGTCTTTAACACGACCGCCACGGATCAGGATCACGGAGTGCTCCTGCAGGTTGTGACCTTCACCACCGATGTAGGAAGTCACTTCGAAACCGTTAGTCAGACGAACACGGCATACTTTACGCAGCGCGGAGTTCGGTTTTTTAGGAGTGGTAGTATATACACGAGTACATACGCCACGTTTTTGCGGGCATGCTTCCAGCGCAGGCACGTTGCTTTTCGCAACTTTGCGAGCACGTGGTTTGCGTACCAGCTGGTTAACTGTTGCCATTAAATAGCTCCTGGTTTTAGCTTTTGCTTCGTAAACACGTAATAAAACGTCCTCACACAATATGAGGACGCCGAATTTTAGGGCGATGCCGAAAAGGTGTCAAGAAATATA  >  NC_000913/3474183‑3474663

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