Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	372,966	A→G	T16A (ACC→GCC)	mhpF →	acetaldehyde‑CoA dehydrogenase II, NAD‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	372,966	0	A	G	100.0%	133.6 / NA	38	T16A (ACC→GCC)	mhpF	acetaldehyde‑CoA dehydrogenase II, NAD‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/20);  total (18/20)

TGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372830‑373099                                                                                                                                         |                                                                                                                                      tGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGcc                                                                                                                                     >  1:695739/1‑139 (MQ=255)     tGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATc                                                                                                                                 >  2:227826/1‑139 (MQ=255)         ggCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtgat                                                                                                                             >  1:740603/1‑139 (MQ=255)             aTCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGAtt                                                                                                                         >  1:623093/1‑139 (MQ=255)                     gAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAAtttt                                                                                                                 <  1:227826/139‑1 (MQ=255)                           cATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCa                                                                                                           <  2:623093/139‑1 (MQ=255)                                tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGTAGTCTGTCATGAGTAAGCGTAAATTCGCCATTATCGGTTCTGGCAACATTTGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  2:511216/139‑1 (MQ=255)                                tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:749678/139‑1 (MQ=255)                                   aGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAg                                                                                                    >  1:611647/1‑138 (MQ=255)                                   aGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAg                                                                                                    <  2:611647/138‑1 (MQ=255)                                       aTAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:604311/132‑1 (MQ=255)                                       aTAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  2:604311/1‑132 (MQ=255)                                         aGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTgg                                                                                             >  1:620109/1‑139 (MQ=255)                                               aGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATgg                                                                                       >  1:94737/1‑139 (MQ=255)                                                              ttCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGTGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCa                                                                        <  2:236722/139‑1 (MQ=255)                                                                       aGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTc                                                               >  2:253182/1‑139 (MQ=255)                                                                                              gAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGcc                                        <  1:252225/139‑1 (MQ=255)                                                                                               aGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCa                                       <  2:159209/139‑1 (MQ=255)                                                                                                 tAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGa                                     >  2:147887/1‑139 (MQ=255)                                                                                                    gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGTCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGt                                  >  2:141918/1‑139 (MQ=255)                                                                                                    gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTc                                                           <  1:199743/114‑1 (MQ=255)                                                                                                    gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTc                                                           >  2:199743/1‑114 (MQ=255)                                                                                                    gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGt                                  >  1:475508/1‑139 (MQ=255)                                                                                                       tAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATg                               <  2:31507/139‑1 (MQ=255)                                                                                                       tAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATg                               <  2:498711/139‑1 (MQ=255)                                                                                                              gCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCg                        >  2:52678/1‑139 (MQ=255)                                                                                                                  ttATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGccac                    <  2:620109/139‑1 (MQ=255)                                                                                                                      cGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCaccacc                <  1:147887/139‑1 (MQ=255)                                                                                                                          tCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATg            >  2:269365/1‑139 (MQ=255)                                                                                                                          tCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATg            >  2:553494/1‑139 (MQ=255)                                                                                                                          tCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATg            <  1:269365/139‑1 (MQ=255)                                                                                                                            tGGCAACATTTGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGaa          <  2:740603/139‑1 (MQ=255)                                                                                                                            tGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGccac                    >  2:95557/1‑129 (MQ=255)                                                                                                                            tGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGccac                    <  1:95557/129‑1 (MQ=255)                                                                                                                                   aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcg                                              >  1:486170/1‑96 (MQ=255)                                                                                                                                   aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcg                                              <  2:486170/96‑1 (MQ=255)                                                                                                                                   aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGa   <  1:553494/139‑1 (MQ=255)                                                                                                                                    ttGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAt  <  2:695739/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      TGAATGCGGGCGATCGTTTTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372830‑373099

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