Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,162,833	C→T	Q146Q (CAG→CAA)	trmA ←	tRNA m(5)U54 methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,162,833	0	C	T	100.0%	110.9 / NA	33	Q146Q (CAG→CAA)	trmA	tRNA m(5)U54 methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/8);  total (25/8)

CGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGG  >  NC_000913/4162766‑4162908                                                                    |                                                                            cGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTg                                                                >  2:758636/1‑81 (MQ=255)      tAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATcg                                                           >  1:758628/1‑81 (MQ=255)          ccTGATCCAGCTCGATTTTGGTTTTAGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATcgcg                                                         >  2:758626/1‑79 (MQ=255)             gATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAggg                                                    >  2:758627/1‑81 (MQ=255)               tCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGcc                                                  >  1:758599/1‑81 (MQ=255)                ccAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCt                                                 >  2:758637/1‑81 (MQ=255)                 cAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTc                                                <  2:758600/81‑1 (MQ=255)                  aGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTcc                                               >  1:758630/1‑81 (MQ=255)                   gCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCg                                              >  2:758638/1‑81 (MQ=255)                     tCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCt                                            >  1:758639/1‑81 (MQ=255)                      cGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCtt                                           >  2:758625/1‑81 (MQ=255)                       gATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTc                                          >  1:758646/1‑81 (MQ=255)                        aTTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTcc                                         >  1:758648/1‑81 (MQ=255)                        aTTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTcc                                         >  1:758601/1‑81 (MQ=255)                         ttttGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCt                                        <  1:758602/81‑1 (MQ=255)                           ttGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTg                                       <  2:758603/80‑1 (MQ=255)                              gTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGc                                   <  1:758604/81‑1 (MQ=255)                               ttttCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGcc                                  <  2:758605/81‑1 (MQ=255)                                 ttCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCa                                 <  1:758606/80‑1 (MQ=255)                                  tCGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACt                               >  2:758644/1‑81 (MQ=255)                                   cGTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACt                               >  2:758642/1‑80 (MQ=255)                                    gTTGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACtca                             >  2:758634/1‑81 (MQ=255)                                      tGCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACtcat                            >  2:758631/1‑80 (MQ=255)                                       gCCCGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACtcatca                          >  1:758643/1‑81 (MQ=255)                                          cGACCAATCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATcc                       <  1:758607/81‑1 (MQ=255)                                                aTCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGcttc                 >  1:758645/1‑81 (MQ=255)                                                aTCAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGcttc                 <  1:758608/81‑1 (MQ=255)                                                  cAGATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGcttct                >  1:758649/1‑80 (MQ=255)                                                    gATGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTAt             >  2:758640/1‑81 (MQ=255)                                                     aTGCACATTCAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATg            >  2:758651/1‑81 (MQ=255)                                                        cacaTTCAGATTTTGCGCGCGCAGTGCATCGCGCAGAGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGta         >  1:758656/1‑81 (MQ=255)                                                              cAGATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAg   >  1:758647/1‑81 (MQ=255)                                                                gATTTTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAgg  >  1:758653/1‑80 (MQ=255)                                                                    |                                                                            CGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACGCCACTCATCATCCAGCTTCTTATGGTATAGCAGG  >  NC_000913/4162766‑4162908

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: ‑