Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	649,372	A→G	F67S (TTC→TCC)	citF ←	citrate lyase, citrate‑ACP transferase (alpha) subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	649,372	0	A	G	100.0%	220.9 / NA	61	F67S (TTC→TCC)	citF	citrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/30);  total (31/30)

GAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCA  >  NC_000913/649238‑649497                                                                                                                                       |                                                                                                                              gAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATgggagga                                                                                                                           <  1:916208/139‑1 (MQ=255)     gtgtTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAAcc                                                                                                                       >  2:13554/1‑139 (MQ=255)              cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                              <  1:13554/139‑1 (MQ=255)                   gcgcATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTTGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTg                                                                                                         <  1:676297/139‑1 (MQ=255)                     gcATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                       >  1:552331/1‑139 (MQ=255)                              tCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGc                                                                                              >  2:732394/1‑139 (MQ=255)                                         gCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                   <  2:552331/139‑1 (MQ=255)                                         gCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGct                                                                                   <  1:34974/139‑1 (MQ=255)                                          cTGGACGCCAGGGTCAGGTTTTTAAAGCTCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGctt                                                                                  >  1:886553/1‑139 (MQ=255)                                           tGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGGAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            <  1:58290/112‑1 (MQ=255)                                           tGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGGAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            >  2:58290/1‑112 (MQ=255)                                                    gggTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGTAAACCGTCATg                                                                                                                 <  2:830367/98‑1 (MQ=255)                                                    gggTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 >  1:830367/1‑98 (MQ=255)                                                          ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  2:886553/139‑1 (MQ=255)                                                          ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:732394/139‑1 (MQ=255)                                                          ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:1084855/139‑1 (MQ=255)                                                          ggTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:1372423/139‑1 (MQ=255)                                                              tttAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGt                                                              >  1:175268/1‑139 (MQ=255)                                                                     cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 <  1:519794/81‑1 (MQ=255)                                                                     cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                 >  2:519794/1‑81 (MQ=255)                                                                       cATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            <  1:506801/84‑1 (MQ=255)                                                                       cATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            >  2:506801/1‑84 (MQ=255)                                                                        aTCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            <  2:772985/83‑1 (MQ=255)                                                                        aTCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                            >  1:772985/1‑83 (MQ=255)                                                                        aTCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  2:1124149/1‑125 (MQ=255)                                                                        aTCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:1124149/125‑1 (MQ=255)                                                                        aTCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgc                                                    >  2:1092496/1‑139 (MQ=255)                                                                         tCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcg                                                   >  2:1057988/1‑139 (MQ=255)                                                                         tCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCATCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcg                                                   >  2:865680/1‑139 (MQ=255)                                                                           ttCGCGATGACGTCCATCACCATATTGACGGTCAGATCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 <  1:865680/139‑1 (MQ=255)                                                                                 aTGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                       <  2:58064/79‑1 (MQ=255)                                                                                 aTGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTa                                                                                                       >  1:58064/1‑79 (MQ=255)                                                                                  tGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTc                                          <  1:1057988/139‑1 (MQ=255)                                                                                   gACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCa                                         <  1:1092496/139‑1 (MQ=255)                                                                                       tCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCt                                     >  1:236170/1‑139 (MQ=255)                                                                                       tCCATCACCAAATTGACGGTGAGGTCACCGCCACGGAAAGCGTGCTGGGCGGACACAGTCATTGCCTCCTATAAACCAGCGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGcca                                     >  2:1387608/1‑138 (MQ=255)                                                                                            cACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGaaa                                                                                                                         <  1:1252357/50‑1 (MQ=255)                                                                                            cACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGaaa                                                                                                                         >  2:1252357/1‑50 (MQ=255)                                                                                              ccATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                               <  2:1050884/58‑1 (MQ=255)                                                                                              ccATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                               >  1:1050884/1‑58 (MQ=255)                                                                                                 tatTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  2:952270/1‑100 (MQ=255)                                                                                                 tatTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCACTTCTTCCAGGTTGGCGca                                                                  <  1:952270/100‑1 (MQ=255)                                                                                                   ttGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCg                                                                                     >  1:50320/1‑79 (MQ=255)                                                                                                   ttGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCg                                                                                     <  2:50320/79‑1 (MQ=255)                                                                                                      aCGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  <  1:198820/95‑1 (MQ=255)                                                                                                      aCGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                  >  2:198820/1‑95 (MQ=255)                                                                                                         gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCg                                                         >  1:1080380/1‑101 (MQ=255)                                                                                                         gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCg                                                         <  2:1080380/101‑1 (MQ=255)                                                                                                         gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcg                                                     >  1:1335511/1‑105 (MQ=255)                                                                                                         gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcg                                                     <  2:1335511/105‑1 (MQ=255)                                                                                                         gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCg                   >  2:248504/1‑139 (MQ=255)                                                                                                            aggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCg                                                                                     <  2:489947/70‑1 (MQ=255)                                                                                                            aggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCg                                                                                     >  1:489947/1‑70 (MQ=255)                                                                                                              gtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGt                                  <  2:471577/119‑1 (MQ=255)                                                                                                              gtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGt                                  >  1:471577/1‑119 (MQ=255)                                                                                                              gtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTc              >  2:923611/1‑139 (MQ=255)                                                                                                                  ccGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGa          >  2:1089893/1‑139 (MQ=255)                                                                                                                        cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    <  1:838132/123‑1 (MQ=255)                                                                                                                        cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGc                    >  2:838132/1‑123 (MQ=255)                                                                                                                          gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa  <  2:175268/139‑1 (MQ=255)                                                                                                                               cGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 >  2:123515/1‑87 (MQ=255)                                                                                                                               cGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATcgcgcg                                                 <  1:123515/87‑1 (MQ=255)                                                                                                                                       |                                                                                                                              GAATGTGTTCTACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCA  >  NC_000913/649238‑649497

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