Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,887,033	G→A	G69D (GGC→GAC)	mnmE →	tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,887,033	0	G	A	100.0%	147.6 / NA	40	G69D (GGC→GAC)	mnmE	tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/16);  total (24/16)

CATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCT  >  NC_000913/3886887‑3887177                                                                                                                                                   |                                                                                                                                                 cATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGAccc                                                                                                                                               >  1:563736/1‑150 (MQ=255)    ccTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAAc                                                                                                                                           <  2:603106/151‑1 (MQ=255)       gcgcATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACt                                                                                                                                          <  1:287465/149‑1 (MQ=255)       gcgcATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTc                                                                                                                                         <  2:769684/150‑1 (MQ=255)               ccGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCg                                                                                                                                 <  1:485222/150‑1 (MQ=255)                   cTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGa                                                                                                                             >  2:418340/1‑150 (MQ=255)                     tCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGa                                                                                                                             >  2:293298/1‑148 (MQ=255)                      cAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGa                                                                                                                             >  2:175148/1‑147 (MQ=255)                        aaGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGAtgt                                                                                                                       >  1:550244/1‑151 (MQ=255)                           cccGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCt                                                                                                                    >  2:454759/1‑151 (MQ=255)                           cccGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCt                                                                                                                    <  1:245964/151‑1 (MQ=255)                                    ttGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGc                                                                                                            >  1:649638/1‑150 (MQ=255)                                           aaCCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCAt                                                                                                    >  1:46482/1‑151 (MQ=255)                                            aCCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATg                                                                                                   <  1:353531/151‑1 (MQ=255)                                               gTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCg                                                                                                <  1:407335/151‑1 (MQ=255)                                                 gCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGt                                                                                              >  2:243677/1‑151 (MQ=255)                                                        aaaaTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCgggg                                                                                        >  1:500616/1‑148 (MQ=255)                                                            tGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATcc                                                                                    >  2:114707/1‑150 (MQ=255)                                                                     cgcgCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACctgct                                                                          <  2:724548/151‑1 (MQ=255)                                                                         cTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGtt                                                                       <  2:586565/150‑1 (MQ=255)                                                                          tACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTa                                                                      >  1:211349/1‑150 (MQ=255)                                                                           aCGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTaa                                                                     >  2:660424/1‑150 (MQ=255)                                                                            cGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTaaaa                                                                   <  2:645403/151‑1 (MQ=255)                                                                            cGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCCCGACCTGCTGTTaaa                                                                    >  1:276376/1‑150 (MQ=255)                                                                                 aTTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCa                                                               >  2:472983/1‑150 (MQ=255)                                                                                 aTTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCAt                                                              >  2:595146/1‑151 (MQ=255)                                                                                    aTCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTc                                                            <  2:74103/150‑1 (MQ=255)                                                                                           gTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGATTCCGTTCACCGGCGAAGATGTGCTGGAACTGCCAGCTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGAccct                                                     <  1:499119/150‑3 (MQ=255)                                                                                              tAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTCGACCTGCTGTTAAAACGCATTCTGACCAATccc                                                 >  1:653286/1‑151 (MQ=255)                                                                                               aaaGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCg                                                <  1:400280/151‑1 (MQ=255)                                                                                                aaGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTccc                                                 >  2:755340/1‑149 (MQ=255)                                                                                                   aCGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGAACTGCTGTTAAAACGAATTCTGACCATTCCCGGcc                                             >  1:491603/1‑150 (MQ=255)                                                                                                                   gCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGGAGATGTGCTGGAACTGCAAGGTCATGGCGGGTCCGAGATCCTCGAACTGCTGTTAAAACGCACTCTGATCCTTCCCGGCCTGCGGATTGCTAGccc                             >  1:624023/1‑150 (MQ=255)                                                                                                                          cAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAg                      >  2:470588/1‑150 (MQ=255)                                                                                                                            ggggATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGttt                   <  1:188025/151‑1 (MQ=255)                                                                                                                                     gcTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAAcg           >  2:331298/1‑150 (MQ=255)                                                                                                                                      cTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAAcgc          >  2:336027/1‑150 (MQ=255)                                                                                                                                        aTGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAAcgcgcg       <  1:114675/151‑1 (MQ=255)                                                                                                                                              cccTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCt  <  2:20626/150‑1 (MQ=255)                                                                                                                                                cTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCAATCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGGTTCt  >  1:440226/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCT  >  NC_000913/3886887‑3887177

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