Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,193,105	G→A	R81H (CGC→CAC)	yqiK →	protein YqiK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,193,105	0	G	A	100.0%	106.2 / NA	30	R81H (CGC→CAC)	yqiK	protein YqiK
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/11);  total (19/11)

GCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTAAAA  >  NC_000913/3193033‑3193181                                                                         |                                                                             gcgcAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCa                                                                      >  1:532309/1‑81 (MQ=255)   gcAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCAtt                                                                    >  2:532315/1‑81 (MQ=255)                      ttCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGa                                                    >  1:532318/1‑78 (MQ=255)                           gAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCg                                             >  2:532310/1‑80 (MQ=255)                           gAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCg                                             >  1:532315/1‑80 (MQ=255)                           gAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCg                                             >  1:532320/1‑80 (MQ=255)                                  tCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGcgcg                                     <  1:532276/81‑1 (MQ=255)                                  tCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGcgcg                                     <  2:532274/81‑1 (MQ=255)                                   ccccATCAATATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGt                                    >  2:532309/1‑81 (MQ=255)                                          aaTATGAATACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATg                               >  2:532331/1‑79 (MQ=255)                                                aaTACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTagtagt                        <  1:532279/80‑1 (MQ=255)                                                aaTACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTc                       <  1:532277/81‑1 (MQ=255)                                                aaTACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATCTAGTAGTc                       <  1:532278/81‑1 (MQ=255)                                                 aTACTCTGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCg                      <  1:532280/81‑1 (MQ=255)                                                       tGAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGctttct                >  1:532330/1‑81 (MQ=255)                                                        gAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGctttctt               >  1:532336/1‑81 (MQ=255)                                                        gAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGctttctt               >  2:532333/1‑81 (MQ=255)                                                        gAAGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGctttctt               >  1:532322/1‑81 (MQ=255)                                                         aaGCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGctttctt               <  1:532281/80‑1 (MQ=255)                                                           gCTGGAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTtgt            <  2:532281/81‑1 (MQ=255)                                                              ggAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCg         >  1:532327/1‑81 (MQ=255)                                                              ggAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCg         >  2:532322/1‑81 (MQ=255)                                                               gAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGc          <  2:532282/79‑1 (MQ=255)                                                               gAAGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCgg        >  2:532327/1‑81 (MQ=255)                                                                aaGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCggg       >  2:532332/1‑81 (MQ=255)                                                                 aGTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGt      <  2:532283/81‑1 (MQ=255)                                                                  gTCAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTa     <  2:532284/81‑1 (MQ=255)                                                                    cAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTaa    >  2:532328/1‑80 (MQ=255)                                                                    cAGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTaa    >  2:532334/1‑80 (MQ=255)                                                                     aGCCACTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTaaaa  >  1:532326/1‑81 (MQ=255)                                                                         |                                                                             GCGCAATCGTGATGCCGATCTTTCATGAAATAATCCCCATCAATATGAATACTCTGAAGCTGGAAGTCAGCCGCTCAACCATTGATAGCCTGATTACGAAAGATCGTATGCGCGTCGATGTAGTAGTCGCTTTCTTTGTGCGGGTAAAA  >  NC_000913/3193033‑3193181

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