Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	343,357	T→C	C158C (TGT→TGC)	yahK →	aldehyde reductase, NADPH‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	343,357	0	T	C	100.0%	94.1 / NA	29	C158C (TGT→TGC)	yahK	aldehyde reductase, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (11/18);  total (11/18)

CCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATG  >  NC_000913/343232‑343474                                                                                                                              |                                                                                                                      ccGACGCCGTCCGAACCGGGCCATACACTGGGAGGCTACTCCCAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGa                                                                                                                <  2:240006/133‑1 (MQ=255) ccGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGa                                                                                                                >  2:393493/1‑133 (MQ=255) ccGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGa                                                                                                                >  1:353738/1‑133 (MQ=255) ccGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCAACCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGa                                                                                                                <  2:317784/133‑1 (MQ=255) ccGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGAGCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGa                                                                                                                >  2:306344/1‑133 (MQ=255)                  gggCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTcgccg                                                                                               <  1:393493/133‑1 (MQ=255)                    gCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCt                                                                                             <  2:331482/133‑1 (MQ=255)                    gCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCt                                                                                             <  1:83000/133‑1 (MQ=255)                             tGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTg                                                                                    >  2:65065/1‑133 (MQ=255)                                gcggcTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCa                                                                                 >  1:423578/1‑133 (MQ=255)                                    cTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAggcc                                                                             >  2:63884/1‑133 (MQ=255)                                    cTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAgg                                                                               >  1:424488/1‑131 (MQ=255)                                    cTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAgg                                                                               <  2:424488/131‑1 (MQ=255)                                         cacaACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAggccgggcc                                                                        <  1:310931/133‑1 (MQ=255)                                                 aTCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTaaaa                                                                >  2:382035/1‑133 (MQ=255)                                                    gtcgtTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAg                                                             <  1:368630/133‑1 (MQ=255)                                                         tCATGAGCGATATGTGCTGCGTTTTCGTCATCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGc                                                        <  1:183258/133‑1 (MQ=255)                                                           aTGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGt                                                      >  1:490615/1‑133 (MQ=255)                                                              aGCGATATGTTCTGCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGt                                                   <  2:423578/133‑1 (MQ=255)                                                                cGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCg                                                 >  1:398851/1‑133 (MQ=255)                                                                   tatGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCa                                              >  1:367164/1‑133 (MQ=255)                                                                       ttCTGCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATcgg                                          <  2:85808/133‑1 (MQ=255)                                                                                      aCCCGCCAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                           <  1:306344/133‑1 (MQ=255)                                                                                      aCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                           <  2:5750/133‑1 (MQ=255)                                                                                      aCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACatat                           <  1:382035/133‑1 (MQ=255)                                                                                           cAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGACCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGa                      <  1:444519/133‑1 (MQ=255)                                                                                                              tGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGAt   <  1:65065/133‑1 (MQ=255)                                                                                                              tGGCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCACCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGAt   <  2:353738/133‑1 (MQ=255)                                                                                                               ggCTCCTTTGTTGTGCGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATg  <  2:367164/133‑1 (MQ=255)                                                                                                                              |                                                                                                                      CCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATCGGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATG  >  NC_000913/343232‑343474

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