Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,209,713	G→A	intergenic (+25/‑61)	rrsE → / → gltV	16S ribosomal RNA/tRNA‑Glu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,209,713	0	G	A	100.0%	43.7 / NA	23	intergenic (+25/‑61)	rrsE/gltV	16S ribosomal RNA/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/11);  total (12/11)

TTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCG  >  NC_000913/4209596‑4209830                                                                                                                      |                                                                                                                      ttCGGGAGGGCGCTTACCACTTAGTTATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGa                                                                                                 <  1:279763/140‑1 (MQ=2)   cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGaaa                                                                                               >  1:714498/1‑140 (MQ=17)    gggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAg                                                                                              >  1:719305/1‑140 (MQ=17)         ggCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGaaa                                                                                         >  1:108359/1‑140 (MQ=17)                       tgtgATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAGCCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTg                                                                           >  2:681740/1‑140 (MQ=17)                            ttCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAg                                                                      <  2:719305/140‑1 (MQ=18)                            ttCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAg                                                                      >  1:123361/1‑140 (MQ=17)                             tCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCGCCTTGCCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGc                                                                     <  2:12261/140‑1 (MQ=12)                              cATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCATCTTGCGAAGCa                                                                    <  2:631849/140‑1 (MQ=16)                                              tCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCt                                                    <  1:558414/140‑1 (MQ=17)                                                   aCAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTc                                               <  2:108598/140‑1 (MQ=17)                                                      aGGTAACCGTAGGGGAGCCTGCGGTTGGATCGCCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                                            <  1:681740/140‑1 (MQ=12)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                                            <  2:374301/140‑1 (MQ=17)                                                      aGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTag                                            <  1:426022/140‑1 (MQ=17)                                                       ggTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTaga                                           >  2:152735/1‑140 (MQ=17)                                                                 ggggAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGac                                 >  1:641111/1‑140 (MQ=17)                                                                       ccTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGcc                           >  1:580456/1‑140 (MQ=17)                                                                                  aTCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCa                     >  1:658547/1‑135 (MQ=17)                                                                                  aTCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCa                     <  2:658547/135‑1 (MQ=17)                                                                                            taccttaAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAgggg      >  1:541524/1‑140 (MQ=17)                                                                                            taccttaAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAgggg      >  1:88720/1‑140 (MQ=18)                                                                                            taccttaAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAgggg      >  2:651297/1‑140 (MQ=18)                                                                                                ttaAAGAAGCGTTCTTTGCAGTACTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCg  <  2:580456/140‑1 (MQ=17)                                                                                                                      |                                                                                                                      TTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTTCTTTGCAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTCTTGCGAAGCAGACTGATACGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCG  >  NC_000913/4209596‑4209830

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