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%	99.7 / NA	29	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 (14/15);  total (14/15)

GTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCAGTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACTGTGACTC  >  NC_000913/3030683‑3030833                                                                            |                                                                            gTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATg                                                                        >  1:554406/1‑81 (MQ=255) gTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATg                                                                        <  2:554357/81‑1 (MQ=255)  tCCAGTTGTACGCTATCGGACTGGAAAGTAAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATgg                                                                       <  2:554358/81‑1 (MQ=255)       ttGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGc                                                                  <  2:554359/81‑1 (MQ=255)         gTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCa                                                                >  1:554410/1‑81 (MQ=255)             gCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCg                                                             >  2:554411/1‑80 (MQ=255)                 tCGGATTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGAt                                                          <  2:554360/79‑1 (MQ=255)                 tCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTc                                                        <  2:554361/81‑1 (MQ=255)                   ggACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCg                                                       <  2:554362/80‑1 (MQ=255)                           aaGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGc                                              >  1:554403/1‑81 (MQ=255)                             gTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCt                                            >  2:554400/1‑81 (MQ=255)                              tCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCtt                                           <  2:554363/81‑1 (MQ=255)                                agagCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCttt                                          <  2:554366/80‑1 (MQ=255)                                agagCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCttt                                          <  2:554364/80‑1 (MQ=255)                                agagCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTc                                         <  2:554365/81‑1 (MQ=255)                                  agCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTccc                                       >  2:554402/1‑81 (MQ=255)                                             gACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGa                            <  2:554368/81‑1 (MQ=255)                                             gACAAGCCGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGa                            <  2:554367/81‑1 (MQ=255)                                                   ccGTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGCCCTTTCCCTACAACGTGGATAcgcg                      >  1:554425/1‑81 (MQ=255)                                                     gTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGt                    <  1:554369/81‑1 (MQ=255)                                                     gTAACGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGt                    >  1:554407/1‑81 (MQ=255)                                                         cGCAAGGGGCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGa                 >  1:554413/1‑80 (MQ=255)                                                               gggCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCAc          >  1:554424/1‑81 (MQ=255)                                                                ggCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACt         >  1:554418/1‑81 (MQ=255)                                                                 gCATTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACtg        >  1:554419/1‑81 (MQ=255)                                                                   aTTGTTGCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACtgt       >  1:554435/1‑80 (MQ=255)                                                                     tgttgCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACTGTGAc    >  1:554422/1‑81 (MQ=255)                                                                      gttgCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACTGTGAc    <  2:554370/80‑1 (MQ=255)                                                                       ttgCACTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACTGTGACTc  <  2:554371/81‑1 (MQ=255)                                                                            |                                                                            GTCCAGTTGTACGCTATCGGACTGGAAAGTCAGAGCATTAACCGGACAAGCCGTAACGCAAGGGGCATTGTTGCAGTGATGGCAGGCCACCGGATTCGCAGCCTGGCCTTTCCCTACAACGTGGATACGCGGTCGAAAGTCACTGTGACTC  >  NC_000913/3030683‑3030833

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