Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,030,758	G→C	H58Q (CAC→CAG)	ygfT ←	fused putative oxidoreductase, Fe‑S subunit and nucleotide‑binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,030,758	0	G	C	100.0%	119.3 / NA	35	H58Q (CAC→CAG)	ygfT	fused putative oxidoreductase, Fe‑S subunit and nucleotide‑binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (19/16);  total (19/16)

CTCGTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCAGTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACTGTGACT  >  NC_000913/3030680‑3030832                                                                               |                                                                           cTCGTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTg                                                                          <  2:497214/81‑1 (MQ=255)    gTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGa                                                                         <  1:497216/79‑1 (MQ=255)       cAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATgg                                                                      >  2:497242/1‑79 (MQ=255)            gTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGcc                                                                <  2:497216/80‑1 (MQ=255)             tACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCAc                                                              <  1:497218/81‑1 (MQ=255)              aCGCTATCGGGCTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCAcc                                                             <  1:497219/81‑1 (MQ=255)              aCGCTAACGGACTGGAAAGACAGGGCATAAACAGGACAAGCCGTAAAGCAAGGGGAAATGCTGCACTGATGGCAGGCCAcc                                                             >  2:497280/1‑81 (MQ=255)                     cGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCg                                                      >  1:497251/1‑81 (MQ=255)                       gACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCa                                                    <  1:497220/81‑1 (MQ=255)                        aCTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAg                                                   >  2:497261/1‑81 (MQ=255)                            gAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGcc                                                 <  2:497220/79‑1 (MQ=255)                            gAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTg                                               >  1:497271/1‑81 (MQ=255)                             aaaGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTgg                                              >  1:497260/1‑81 (MQ=255)                                 tCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCtt                                          >  1:497275/1‑81 (MQ=255)                                  cAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCttt                                         >  1:497256/1‑81 (MQ=255)                                   agagCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCttt                                         <  2:497221/80‑1 (MQ=255)                                      gCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTccc                                      >  1:497261/1‑80 (MQ=255)                                      gCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTccc                                      >  2:497256/1‑80 (MQ=255)                                          tAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACaa                                 >  2:497269/1‑81 (MQ=255)                                            aCCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACg                               <  2:497222/81‑1 (MQ=255)                                               ggACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTg                             >  2:497278/1‑80 (MQ=255)                                                gACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGa                           >  1:497258/1‑81 (MQ=255)                                                gACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGa                           <  1:497224/81‑1 (MQ=255)                                                     gCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATAcgc                      <  1:497225/81‑1 (MQ=255)                                                     gCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATAcgc                      <  1:497226/81‑1 (MQ=255)                                                        gTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGt                   >  1:497269/1‑81 (MQ=255)                                                         tAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTc                  <  2:497226/81‑1 (MQ=255)                                                         tAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTc                  <  2:497228/81‑1 (MQ=255)                                                         tAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTc                  <  1:497228/81‑1 (MQ=255)                                                                 ggggCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCa          >  1:497280/1‑81 (MQ=255)                                                                   ggCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACt        >  2:497266/1‑81 (MQ=255)                                                                   ggCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACt        >  1:497272/1‑81 (MQ=255)                                                                      aTTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACtgtg     <  2:497229/81‑1 (MQ=255)                                                                      aTTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACtgt      >  2:497264/1‑80 (MQ=255)                                                                         gttgCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACATCGTGGATACGCGGTCGAAAGTCACTGTGACt  >  2:497265/1‑81 (MQ=255)                                                                               |                                                                           CTCGTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCAGTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACTGTGACT  >  NC_000913/3030680‑3030832

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