Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	343,447	T→G	H188Q (CAT→CAG)	yahK →	broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	343,447	0	T	G	100.0%	109.5 / NA	34	H188Q (CAT→CAG)	yahK	broad specificity NADPH‑dependent aldehyde reductase, Zn‑containing
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (11/23);  total (11/23)

CGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCT  >  NC_000913/343216‑343680                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          cGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGa                                                                                                                                                                                                                                    <  1:795139/239‑1 (MQ=255)            ccTCGCGGCCGCCGGCCGAACCGGGCCATACTCTGGGCGGCTCCTCACAACAGATCGTCGTTCAGGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGcc                                                                                                                                                                                                                         <  2:1150450/238‑1 (MQ=255)                                       atcctctGGGCTGCTACCCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCa                                                                                                                                                                                              <  1:580178/236‑1 (MQ=255)                                       aTACTCTTGGCGTCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCa                                                                                                                                                                                              <  1:804999/239‑1 (MQ=255)                                            ctGGCCGGCTACTCACAACAGATCGTCGTTCATGAGCGTTATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTAc                                                                                                                                                                                         <  1:173713/239‑1 (MQ=255)                                               ggcggcTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGGAAAAAAGTGGGCGTGGTCGGCATCGGGGGGGCTGGGGCGGTGGGGGATGAAGTGGGCCCCGGCGTGGGGGGGCCAGGTGGGGGCAttttca                                                                                                                                                                                        >  2:214815/1‑235 (MQ=255)                                                                         tCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCaa                                                                                                                                                            <  1:183437/239‑1 (MQ=255)                                                                         tCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCa                                                                                                                                                             <  1:811171/238‑1 (MQ=255)                                                                                     tGTTCTGCGTATTCGTCACCCGCAAGGGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGAt                                                                                                                                                                                                                 <  1:585710/174‑1 (MQ=255)                                                                                           gCGTATTCGTCCCCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATg                                                                                                                                          <  1:214815/239‑1 (MQ=255)                                                                                           gCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAg                                                                                                                                                                                                                               >  1:179037/1‑154 (MQ=255)                                                                                           gCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAg                                                                                                                                                                                                                               <  2:179037/154‑1 (MQ=255)                                                                                                   tgcgcccGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGGAGGGATCACCAGGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCGGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt                                                                                                                                  <  1:1013334/235‑1 (MQ=255)                                                                                                   gTCACCCGCAAGAGCGGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAgttgtt                                                                                                                                  <  1:724808/239‑1 (MQ=255)                                                                                                                  agcTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGGCTGGGGCAGGTGGGGGTTAAGATGGGCCCCGCGGAGGGGGGGCCAGTGGGGGGCTTTTTCACTTTTTGGGGAAAAAGACAAACGGGCAAAGCCCCGGGGGGCGGTGTAGATGTTTACACAAGCaaatg                                                                                                                    >  2:4814/1‑235 (MQ=255)                                                                                                                                                                    ccgcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGGCAGATGGGGATTAAGCTGGCCCCCGCGATGGGGGCACATGTGGTGGCATTTACCCATTCTGAGGCAAAACGCGAAGCGGCAAAAACCCCTGGGGCCGATGAAAGTGTTAAATCACGCAAAGCCGATGAGAAGGGGGCGCCTCTTAAGAGTTTCGGTTTTATTTTGAATaaaag                                                                 >  2:871895/1‑235 (MQ=255)                                                                                                                                                                      gcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAg                                                               <  2:153618/239‑1 (MQ=255)                                                                                                                                                                      gcTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAg                                                               <  2:166819/239‑1 (MQ=255)                                                                                                                                                                          cGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATgag                                                                                                             <  1:625469/189‑1 (MQ=255)                                                                                                                                                                          cGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgc                                                           <  1:247870/239‑1 (MQ=255)                                                                                                                                                                           gTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTgcg                                                          <  1:669255/239‑1 (MQ=255)                                                                                                                                                                            tCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATgtgg                                                                                                                                                                                                   <  1:824800/101‑1 (MQ=255)                                                                                                                                                                                gggCAGGCCGGGCCGGGTAAAAAAGTGGGCGTTGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCcaca                                                     <  1:1050373/238‑1 (MQ=255)                                                                                                                                                                                      gccgggccggGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATggggg                                                                                                                                                                                                            <  2:1565/82‑1 (MQ=255)                                                                                                                                                                                      gccgggccggGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGGCTGGGGCAGATGGGGGGTAAGCTTGGCCACGCGGTGGGGGGCCATGGGGGGGCAttttc                                                                                                                                                                                         >  2:824800/1‑99 (MQ=255)                                                                                                                                                                                                gTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGGCAGGTGGGGGTTAAGCTGGCCCACGCGGTGGGGGGACATggggg                                                                                                                                                                                                  >  1:1565/1‑81 (MQ=255)                                                                                                                                                                                                          gggCGTGGTCGGCATC‑GGCGGGCTGGGGCAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCGCAACTGAAGAGGTTCGATTTTCTTTTTAATAAAGTAGCCGCGCCACATAATATCGACGATTTTACCACCtt                             >  2:218068/1‑237 (MQ=255)                                                                                                                                                                                                           ggCGTGGTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCCGGCCCCCGCGGTGGGGGCACATGTGGGGGCATTTACCACTTCTGAGGCCAAACGCGGAGCGGCAAAAGCCCTGGGGGGCGATGAAGTTGTTAACTCACGCAAAGCCGATGAGATGGGGGGCCATCTGAAGAGGTTCGATTTCATTTTTAATACAGGAGCTGCGCCCCCTAAATTCGACGATTTTTCCACCtt                             >  2:280005/1‑236 (MQ=255)                                                                                                                                                                                                                 gTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGcac                                                                                                                                                                                                         >  1:672522/1‑58 (MQ=255)                                                                                                                                                                                                                 gTCGGCATC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGcac                                                                                                                                                                                                         <  2:672522/58‑1 (MQ=255)                                                                                                                                                                                                                        tC‑GGCGGTCTGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATgg              <  1:501322/238‑1 (MQ=255)                                                                                                                                                                                                                                  tGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGcac                                                                                                                                                                                                         <  2:699292/42‑1 (MQ=255)                                                                                                                                                                                                                                  tGGGACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGcac                                                                                                                                                                                                         >  1:699292/1‑42 (MQ=255)                                                                                                                                                                                                                                    ggACAGATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATATCGACGATTTTACCACCTTGCTGAAGCGTGATGGAACCATGACGcg  >  2:329477/1‑238 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          CGGCACCTATAACTCGCCGACGCCGGACGAACCGGGCCATACTCTGGGCGGCTACTCACAACAGATCGTCGTTCATGAGCGATATGTTCTGCGTATTCGTCACCCGCAAGAGCAGCTGGCGGCGGTGGCTCCTTTGTTGTGTGCAGGGATCACCACGTATTCGCCGCTACGTCACTGGCAGGCCGGGCCGGGTAAAAAAGTGGGCGTGGTCGGCATC‑GGCGGTCTGGGACATATGGGGATTAAGCTGGCCCACGCGATGGGGGCACATGTGGTGGCATTTACCACTTCTGAGGCAAAACGCGAAGCGGCAAAAGCCCTGGGGGCCGATGAAGTTGTTAACTCACGCAATGCCGATGAGATGGCGGCTCATCTGAAGAGTTTCGATTTCATTTTGAATACAGTAGCTGCGCCACATAATCTCGACGATTTTACCACCTTGCTGAAGCGTGATGGCACCATGACGCT  >  NC_000913/343216‑343680

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