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	96.5%	95.2 / ‑4.0	29	R162C (CGT→TGT)	citC	citrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (13/15);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGC  >  NC_000913/651295‑651449                                                                               |                                                                             ggAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGt                                                                            <  2:135276/81‑1 (MQ=255)        tGAAGAATCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCAtt                                                                     <  1:135277/81‑1 (MQ=255)          aagaaTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCg                                                                   >  1:135315/1‑81 (MQ=255)             aaTCTTCTTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTa                                                                 >  1:135316/1‑80 (MQ=255)               tcttctTTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAg                                                              >  2:135332/1‑81 (MQ=255)                  tcttTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGa                                                           >  1:135322/1‑81 (MQ=255)                   cttTGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAAGGGa                                                           <  2:135278/80‑1 (MQ=255)                      tGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGAtt                                                         <  2:135281/79‑1 (MQ=255)                      tGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGAt                                                          <  2:135279/78‑1 (MQ=255)                      tGACTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTg                                                        <  1:135280/80‑1 (MQ=255)                        aCTAAAAACAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCg                                                     >  1:135319/1‑81 (MQ=255)                            aaaaCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTtca                                                 <  2:135283/81‑1 (MQ=255)                            aaaaCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTtc                                                  <  2:135282/80‑1 (MQ=255)                              aaCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTtcat                                                >  1:135338/1‑80 (MQ=255)                               aCAGATGCAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTtcatca                                              >  2:135323/1‑81 (MQ=255)                                      cAACCAGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATg                                         >  2:135327/1‑79 (MQ=255)                                         ccAGTCGCACTGTGCCGCAGCCTGTTGAATCGGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGcc                                    <  2:135284/81‑1 (MQ=255)                                           aGTCGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCaa                                  >  2:135325/1‑81 (MQ=255)                                              cGCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCt                               >  1:135337/1‑81 (MQ=255)                                               gCACTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCt                               <  1:135285/80‑1 (MQ=255)                                                 aCTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGt                            <  1:135286/81‑1 (MQ=255)                                                 aCTGTGCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGt                            >  2:135329/1‑81 (MQ=255)                                                     tgCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTccc                        >  1:135326/1‑81 (MQ=255)                                                     tgCCGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTcc                         <  1:135287/80‑1 (MQ=255)                                                       ccGCAGCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCt                       <  2:135288/80‑1 (MQ=255)                                                            gCCTGTTGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGa                 <  2:135289/81‑1 (MQ=255)                                                                 ttGAATCAGATAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGaa             <  2:135290/80‑1 (MQ=255)                                                                           tAACAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGc  <  2:135291/81‑1 (MQ=255)                                                                            aaCAGTGACCATTCGTAAAGGGATTGGCGTTCATCACAATGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGc  >  1:135328/1‑80 (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: ‑