Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,471,692	G→T	A608E (GCG→GAG)	fusA ←	elongation factor G

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,471,692	0	G	T	100.0%	127.6 / NA	35	A608E (GCG→GAG)	fusA	elongation factor G
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/21);  total (14/21)

ACCGATAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCGCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCAGAGGAGTCAACGTCATGGT  >  NC_000913/3471616‑3471767                                                                             |                                                                            aCCGATAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCttt                                                                         <  1:654400/81‑1 (MQ=255) aCCGATAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCtt                                                                          <  2:654399/80‑1 (MQ=255)  ccGATAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTc                                                                        <  2:654401/81‑1 (MQ=255)   cGATAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCt                                                                       <  2:654402/81‑1 (MQ=255)    gATAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCtctttctt                                                                      <  1:654403/81‑1 (MQ=255)     aTAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCtctttctt                                                                      >  1:654443/1‑80 (MQ=255)     aTAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTa                                                                     <  1:654404/81‑1 (MQ=255)     aTAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTa                                                                     <  1:654405/81‑1 (MQ=255)         cGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGc                                                                 <  2:654406/81‑1 (MQ=255)                ggTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCcttctt                                                           >  2:654433/1‑80 (MQ=255)                ggTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCttt                                                          <  2:654407/81‑1 (MQ=255)                 gtgtTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTAGTTTCTCTTTCTTAAAGCCTTCTTTa                                                         <  2:654408/81‑1 (MQ=255)                    ttCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAg                                                      <  2:654410/81‑1 (MQ=255)                    ttCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAg                                                      <  2:654409/81‑1 (MQ=255)                        ctTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGa                                                  >  1:654464/1‑81 (MQ=255)                           ccGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATa                                                <  1:654411/80‑1 (MQ=255)                                gTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAag                                            >  1:654452/1‑79 (MQ=255)                                 tttCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAagcag                                         >  1:654449/1‑81 (MQ=255)                                    cTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAagcagc                                        >  1:654446/1‑79 (MQ=255)                                    cTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGcc                                       <  1:654412/80‑1 (MQ=255)                                       cTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGtt                                   <  1:654414/81‑1 (MQ=255)                                       cTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGtt                                   <  2:654413/81‑1 (MQ=255)                                            aCCTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAAc                              <  1:654415/81‑1 (MQ=255)                                             ccTTCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACg                             >  2:654465/1‑81 (MQ=255)                                                tCATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGc                            >  2:654458/1‑79 (MQ=255)                                                 cATGATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGcc                           <  2:654416/79‑1 (MQ=255)                                                 cATGATCGGCTCAAGCAGAACTGGTTTCTCATTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCa                          <  2:654417/80‑1 (MQ=255)                                                    gATCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGtt                       >  1:654444/1‑80 (MQ=255)                                                      tCGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGtt                       >  2:654454/1‑78 (MQ=255)                                                       cGGCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCa                     >  2:654451/1‑79 (MQ=255)                                                        ggCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCagag                  >  2:654450/1‑81 (MQ=255)                                                         gCTCAAGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCagag                  >  1:654455/1‑80 (MQ=255)                                                             aaGCAGAACTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCAGAGGAGTc             <  2:654418/81‑1 (MQ=255)                                                             aaGCAGAACTGGGTTCTCTTTCTTAAGGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCAGAGGAGt              <  2:654419/80‑1 (MQ=255)                                                                   aaCTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCAGAGGAGTCAACGTc       <  2:654420/81‑1 (MQ=255)                                                                    aCTGGTTTCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCAGAGGAGTCAACGTCa      >  1:654453/1‑81 (MQ=255)                                                                          ttCTCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCAGAGGAGTCAACGTCATGGt  >  1:654469/1‑79 (MQ=255)                                                                             |                                                                            ACCGATAACGTCACCGGTGTTCTCTTCCGGAGTTTCTACTTCAACCTTCATGATCGGCTCAAGCAGAACTGGTTTCGCTTTCTTAAAGCCTTCTTTAAAGGCGATAGAAGCAGCCAGTTTAAACGCCAGTTCAGAGGAGTCAACGTCATGGT  >  NC_000913/3471616‑3471767

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