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%	167.3 / NA	48	F67S (TTC→TCC)	citF	citrate lyase, citrate‑ACP transferase (alpha) subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/23);  total (25/23)

TACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATT  >  NC_000913/649248‑649507                                                                                                                             |                                                                                                                                        tACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATg                                                                                                                           <  2:222010/139‑1 (MQ=255)    cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                        <  2:32391/139‑1 (MQ=255)    cAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                        <  1:543523/139‑1 (MQ=255)      gcggcgCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                                      <  2:964660/139‑1 (MQ=255)      gcggcgCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTc                                                                                                                      <  1:813768/139‑1 (MQ=255)                         cAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGcgacga                                                                                                   >  1:321177/1‑139 (MQ=255)                          aGGGAGCTGGACGCCAGGGTCAGGGTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGaa                                                                                                  >  2:891794/1‑139 (MQ=255)                                  ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                                         >  2:415362/1‑108 (MQ=255)                                  ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGcc                                                                                                                         <  1:415362/108‑1 (MQ=255)                                  ggACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGcttct                                                                                          >  2:1079134/1‑139 (MQ=255)                                          gggTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGtt                                                                                  >  2:645858/1‑139 (MQ=255)                                                           cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGcga                                                                                                      >  2:737544/1‑102 (MQ=255)                                                           cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGcga                                                                                                      <  1:737544/102‑1 (MQ=255)                                                           cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAgg                                                                                    <  2:1049004/120‑1 (MQ=255)                                                           cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAgg                                                                                    >  1:1049004/1‑120 (MQ=255)                                                           cccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGAt                                                                 >  1:1048945/1‑139 (MQ=255)                                                            ccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                        >  1:1040272/1‑83 (MQ=255)                                                            ccATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCg                                                                                                                        <  2:1040272/83‑1 (MQ=255)                                                                       aTGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCttttt                                                     <  2:631470/139‑1 (MQ=255)                                                                         gACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCa                                                   <  1:510660/139‑1 (MQ=255)                                                                          aCGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTaa                                                                                                                >  1:185367/1‑77 (MQ=255)                                                                          aCGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTaa                                                                                                                <  2:185367/77‑1 (MQ=255)                                                                           cGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                            >  2:1270016/1‑112 (MQ=255)                                                                           cGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                            <  1:1270016/112‑1 (MQ=255)                                                                           cGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGc                                                 <  1:645858/139‑1 (MQ=255)                                                                             tCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCt                                               <  1:450880/139‑1 (MQ=255)                                                                             tCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCt                                               <  2:194034/139‑1 (MQ=255)                                                                                    ccATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGt                                        >  2:749532/1‑139 (MQ=255)                                                                                       tatTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                            <  1:308554/100‑1 (MQ=255)                                                                                       tatTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                            >  2:308554/1‑100 (MQ=255)                                                                                       tatTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTg                                     >  2:82653/1‑139 (MQ=255)                                                                                            aCGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                            >  1:477946/1‑95 (MQ=255)                                                                                            aCGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGca                                                                            <  2:477946/95‑1 (MQ=255)                                                                                            aCGGTCAGGTCACCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTg                                >  1:674143/1‑139 (MQ=255)                                                                                               gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGTGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCg                             >  2:691293/1‑139 (MQ=255)                                                                                               gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCg                             >  2:290683/1‑139 (MQ=255)                                                                                               gtcaggtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTACCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCg                             >  2:804105/1‑139 (MQ=255)                                                                                                    gtcaCCGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGATCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTc                        >  1:562897/1‑139 (MQ=255)                                                                                                         cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGaa                   >  1:169215/1‑139 (MQ=255)                                                                                                         cGCCACGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGaa                   <  1:1079134/139‑1 (MQ=255)                                                                                                              cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGaa                                                                                                  >  1:479784/1‑55 (MQ=255)                                                                                                              cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGaa                                                                                                  <  2:479784/55‑1 (MQ=255)                                                                                                              cGGAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCgg              >  1:784585/1‑139 (MQ=255)                                                                                                                gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGTGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa            <  1:691293/139‑1 (MQ=255)                                                                                                                gAAAGCGTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa            <  2:784585/139‑1 (MQ=255)                                                                                                                      gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa            <  1:1134897/133‑1 (MQ=255)                                                                                                                      gTGATGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa            >  2:1134897/1‑133 (MQ=255)                                                                                                                         aTGGGAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCAt   >  2:207535/1‑139 (MQ=255)                                                                                                                          tGGGAGGAAAACGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTATCAGGTTGGCGCAGAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAATCGGCAAGAACGCAtt  <  2:321177/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                        TACCAGCGGCGCATGGCAATCACTCAGGGAGCTGGACGCCAGGGTCAGGTTTTTAAAGCCCATCTTCGCGATGACGTCCATCACCATATTGACGGTCAGGTCACCGCCACGGAAAGCGTGATGGAAGGAAACCGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATT  >  NC_000913/649248‑649507

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