Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	651,373	G→A	R162C (CGT→TGT)	citC ←	citrate lyase synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	651,373	0	G	A	100.0%	108.4 / NA	31	R162C (CGT→TGT)	citC	citrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/15);  total (16/15)

GGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGC  >  NC_000913/651295‑651449                                                                               |                                                                             ggAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGt                                                                            >  2:123656/1‑81 (MQ=255) ggAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGt                                                                            <  1:123614/81‑1 (MQ=255)        tGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCAtt                                                                     <  1:123615/81‑1 (MQ=255)        tGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCAtt                                                                     <  1:123616/81‑1 (MQ=255)         gaagaaTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTc                                                                    >  1:123650/1‑81 (MQ=255)          aagaaTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCg                                                                   <  2:123617/81‑1 (MQ=255)            gaaTCTTCTTTGACTAAAAACAGATGCAACCGGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTa                                                                 <  2:123618/81‑1 (MQ=255)              aTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTaaa                                                               <  1:123619/81‑1 (MQ=255)                     ttGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTg                                                        >  2:123654/1‑81 (MQ=255)                       gACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGc                                                      >  2:123646/1‑81 (MQ=255)                              aaCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTtcatc                                               <  2:123620/81‑1 (MQ=255)                                   aTGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAAt                                          <  1:123621/81‑1 (MQ=255)                                      cAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATg                                         >  2:123661/1‑79 (MQ=255)                                       aaCCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAg                                      >  1:123658/1‑81 (MQ=255)                                          cAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCa                                   <  2:123622/81‑1 (MQ=255)                                            gTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAAt                                 <  1:123623/81‑1 (MQ=255)                                               gCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCtt                              >  1:123669/1‑81 (MQ=255)                                                cACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTg                             >  1:123653/1‑81 (MQ=255)                                                   tgtgCCGCAGCCCGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTc                          >  2:123673/1‑81 (MQ=255)                                                    gtgCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTcc                         >  1:123655/1‑81 (MQ=255)                                                       ccGCAGCCTGTTGAATCAGATAACAGTGACCGTTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTg                      <  2:123624/81‑1 (MQ=255)                                                          cAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGAt                   <  1:123625/81‑1 (MQ=255)                                                              cTGTTGAATCAGATAACAGTGACCATTCGTAGAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACg               <  2:123627/81‑1 (MQ=255)                                                              cTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACg               <  2:123626/81‑1 (MQ=255)                                                               tGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACg               >  1:123662/1‑80 (MQ=255)                                                               tGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACg               >  2:123668/1‑80 (MQ=255)                                                                  tGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAt           >  1:123667/1‑81 (MQ=255)                                                                   gAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAtt          >  2:123665/1‑81 (MQ=255)                                                                     aTCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAtttt        >  1:123666/1‑81 (MQ=255)                                                                      tCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAAttttt       >  1:123664/1‑81 (MQ=255)                                                                           tAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGc  <  1:123628/81‑1 (MQ=255)                                                                               |                                                                             GGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGC  >  NC_000913/651295‑651449

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