Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,953,658	G→T	A73S (GCA→TCA)	cmoB →	tRNA U34 carboxymethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,953,658	0	G	T	100.0%	102.6 / NA	30	A73S (GCA→TCA)	cmoB	tRNA U34 carboxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/13);  total (17/13)

TTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGG  >  NC_000913/1953583‑1953734                                                                            |                                                                             tttCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAggg                                                                         <  1:333468/81‑1 (MQ=255)       ccTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCaa                                                                      <  2:333469/78‑1 (MQ=255)        cTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAAtt                                                                   >  2:333510/1‑80 (MQ=255)        cTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAAtt                                                                   >  2:333509/1‑80 (MQ=255)         tGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTa                                                                  >  1:333501/1‑80 (MQ=255)            aaTTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACg                                                              >  2:333506/1‑81 (MQ=255)                 aaCCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCAtt                                                          <  1:333471/80‑1 (MQ=255)                 aaCCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAAGGCAAATTAAACGCAtt                                                          <  2:333470/80‑1 (MQ=255)                    cGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGaaa                                                      <  1:333472/81‑1 (MQ=255)                         cGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCt                                                  >  1:333514/1‑80 (MQ=255)                          gTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTg                                                 <  1:333473/80‑1 (MQ=255)                             ttGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATgc                                             >  1:333518/1‑81 (MQ=255)                              tGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATgcg                                            >  1:333507/1‑81 (MQ=255)                               gATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATgc                                             <  2:333474/79‑1 (MQ=255)                                     ttGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCt                                      <  2:333475/80‑1 (MQ=255)                                     ttGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCt                                      <  1:333476/80‑1 (MQ=255)                                      tGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGa                                    >  2:333521/1‑81 (MQ=255)                                           aGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCg                               >  1:333504/1‑81 (MQ=255)                                            gCGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCg                               <  2:333477/80‑1 (MQ=255)                                             cGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTg                             >  2:333512/1‑81 (MQ=255)                                             cGTTACGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTg                             >  2:333517/1‑81 (MQ=255)                                                  cGGCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGgcgc                         <  2:333478/80‑1 (MQ=255)                                                   ggCAGAAAGCGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCa                        >  1:333508/1‑80 (MQ=255)                                                            cGAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGtt              <  1:333479/81‑1 (MQ=255)                                                             gAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTtc             >  1:333522/1‑81 (MQ=255)                                                             gAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTtc             >  1:333516/1‑81 (MQ=255)                                                             gAAGAGCCACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTtc             >  2:333528/1‑81 (MQ=255)                                                                    cACTGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGt      <  2:333480/81‑1 (MQ=255)                                                                       tGAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATg   >  1:333526/1‑81 (MQ=255)                                                                        gAGCTCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATgg  >  2:333513/1‑81 (MQ=255)                                                                            |                                                                             TTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTGAAACGCTGATGCGCAACCTGATGCCGTGGCGCAAAGGGCCGTTCTCACTGTATGG  >  NC_000913/1953583‑1953734

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