Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,632	G→A	intergenic (‑231/+37)	leuQ ← / ← rsmC	tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,632	0	G	A	100.0%	92.9 / NA	30	intergenic (‑231/+37)	leuQ/rsmC	tRNA‑Leu/16S rRNA m(2)G1207 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/17);  total (13/17)

ATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTC  >  NC_000913/4606513‑4606756                                                                                                                        |                                                                                                                             aTCAGGAACCTAGCGGGCAGGTTTAGTGCGCAGGCGGGGGCTGATGCCAGATGCGACGCTGGGGCGTCTTATCTGGCCTACGAAGGGCTAGCGTGCAGGTTTTGTAGGTAGGATAAGGAATTCACGCCGCATc                                                                                                                 <  2:167833/133‑1 (MQ=21)                  cAGGTTTGGTATGCATGCGGGGGCCGATCCCGCATGCGACGCTGGCGCGTCTTATCGGGCTTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCg                                                                                                <  2:159273/133‑1 (MQ=38)                     gTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCgaga                                                                                             <  1:460149/133‑1 (MQ=255)                     gTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTAACGCCGCATCCGACTCGGTATTCGGCgaga                                                                                             >  1:155770/1‑133 (MQ=255)                               catgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCt                                                                                   >  1:429717/1‑133 (MQ=255)                               catgcGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATACACGCCGCCTCCAACATGGTATTCGGCGAGATAATTAACCt                                                                                   >  2:201792/1‑133 (MQ=37)                                  gcGGGGGCAGATGCCAGATGCGACGCTTGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTttc                                                                                >  2:457623/1‑133 (MQ=255)                                           gatgccagattcGACCCTGGCGCGTCTTCTCTAGCCTACGAGGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTgg                                                                       <  2:155770/133‑1 (MQ=38)                                           gatgccagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTgg                                                                       <  2:429290/133‑1 (MQ=255)                                                cagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTgg                                                                       >  1:411088/1‑128 (MQ=255)                                                cagatgcGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTgg                                                                       <  2:411088/128‑1 (MQ=255)                                                            tGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGcgcg                                                      <  1:72636/133‑1 (MQ=255)                                                                  gTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATacac                                                <  2:429717/133‑1 (MQ=255)                                                                             ccTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGcc                                     >  2:326643/1‑133 (MQ=255)                                                                             ccTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGcc                                     >  1:345200/1‑133 (MQ=255)                                                                                 cGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATATTTGCGCGATACACCTTGATGCGCCCTGt                                 >  2:229855/1‑133 (MQ=255)                                                                                     gggCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTgc                             <  1:450526/133‑1 (MQ=255)                                                                                      ggCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTgcg                            >  2:387510/1‑133 (MQ=255)                                                                                          aaCGTGCAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGAt                        <  1:201792/133‑1 (MQ=255)                                                                                          aaCGTGCAGGTTTTATAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGAt                        <  2:502981/133‑1 (MQ=255)                                                                                              tGCAGGTTATGTAGGACGAATACGGCATTCTCGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACt                    <  1:138507/133‑1 (MQ=255)                                                                                                cAGGTTTTGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACAACTTGAAGCGCCATGTTTGCGCTATCACTTc                  >  1:197500/1‑133 (MQ=255)                                                                                                       tGTAGGTCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGGAGCGCCCGGTTTGCGCGATCATTTCGTGGAAg           >  2:138507/1‑133 (MQ=255)                                                                                                             tCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAAt     >  2:25036/1‑133 (MQ=255)                                                                                                             tCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAAt     >  2:398984/1‑133 (MQ=255)                                                                                                             tCGGATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAAt     <  1:280868/133‑1 (MQ=255)                                                                                                             tCGGATAAGGCATTCACGCCGCATCCGACAAGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAAt     <  1:229855/133‑1 (MQ=255)                                                                                                               ggATAAGGCATTCACGCCGCATACGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGt   <  1:212230/133‑1 (MQ=255)                                                                                                                gATAAGGCATTCACGTCGCATCCGACACTGTATTCGGCGAGTTAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtc  <  2:197500/133‑1 (MQ=255)                                                                                                                gATAAGGCATTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGtc  <  1:326643/133‑1 (MQ=255)                                                                                                                        |                                                                                                                             ATCAGGAAACTAGCGTGCAGGTTTGGTATGCATGCGGGGGCAGATGCCAGATGCGACGCTGGCGCGTCTTATCTGGCCTACGAAGGGCTAACGTGCAGGTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACACGGTATTCGGCGAGATAATTAACCTTTCTTCGCCTGGCGGGTCATAATGGCGCGATACACCTTGAAGCGCCCTGTTTGCGCGATCACTTCGTGGAAGCCAAATGTC  >  NC_000913/4606513‑4606756

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