Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	648,446	T→C	T376A (ACC→GCC)	citF ←	citrate lyase alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	648,446	0	T	C	100.0%	138.9 / NA	39	T376A (ACC→GCC)	citF	citrate lyase alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/25);  total (14/25)

TACGACCGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGATTGCGCAGCATG  >  NC_000913/648320‑648578                                                                                                                               |                                                                                                                                     tACGACCGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAgc                                                                                                                         <  1:602594/140‑1 (MQ=255)  aCGACCGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAgcg                                                                                                                        <  1:581225/140‑1 (MQ=255)   cGACCGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAgcgc                                                                                                                       >  2:495567/1‑140 (MQ=255)     aCCGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAgcgcgc                                                                                                                     <  1:1081035/140‑1 (MQ=255)     aCCGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAgcgcgc                                                                                                                     <  1:932715/140‑1 (MQ=255)      ccGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCt                                                                                                                    <  2:514358/140‑1 (MQ=255)      ccGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCt                                                                                                                    >  2:309240/1‑140 (MQ=255)      ccGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCt                                                                                                                    <  2:140164/140‑1 (MQ=255)      ccGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCt                                                                                                                    <  1:495567/140‑1 (MQ=255)       cgcgTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCt                                                                                                                    <  1:1011534/139‑1 (MQ=255)       cgcgTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCt                                                                                                                    >  2:1011534/1‑139 (MQ=255)         cgTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAg                                                                                                                 >  1:1021515/1‑140 (MQ=255)         cgTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAg                                                                                                                 <  1:948840/140‑1 (MQ=255)                   gcgcGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACg                                                                                                       >  2:898186/1‑140 (MQ=255)                    cgcgACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGt                                                                                                      >  2:1060732/1‑140 (MQ=255)                         cgaTGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCgaga                                                                                                 <  1:863889/140‑1 (MQ=255)                           atgatgGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGGCGAGACg                                                                                               <  2:713219/140‑1 (MQ=255)                              atgGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATc                                                                                            <  2:158735/140‑1 (MQ=255)                                      cgcAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATg                                                                                    >  2:552656/1‑140 (MQ=255)                                      cgcAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATg                                                                                    <  1:552656/140‑1 (MQ=255)                                      cgcAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATg                                                                                    >  2:522007/1‑140 (MQ=255)                                        cAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCGGTACCACCACGTCGAGACGATCAACCGATgcg                                                                                  <  2:236949/140‑1 (MQ=255)                                              cAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTc                                                                            <  1:492962/140‑1 (MQ=255)                                                  cgcgGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAAc                                                                        <  1:101294/140‑1 (MQ=255)                                                     ggTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAAcccc                                                                     <  1:1060732/140‑1 (MQ=255)                                                     ggTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAAcccc                                                                     <  2:868096/140‑1 (MQ=255)                                                                  ccaccGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTg                                                        <  2:579741/140‑1 (MQ=255)                                                                      cGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTg                                                    <  1:898186/140‑1 (MQ=255)                                                                          aGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGgcgc                                                <  2:123560/140‑1 (MQ=255)                                                                          aGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGgcgc                                                <  1:536899/140‑1 (MQ=255)                                                                                          ccGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGAtt                                <  1:427897/140‑1 (MQ=255)                                                                                          ccGTCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAACTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGAtt                                <  2:1021515/140‑1 (MQ=255)                                                                                             tCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTggg                             >  2:454808/1‑140 (MQ=255)                                                                                             tCAGAGCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTggg                             >  1:652628/1‑140 (MQ=255)                                                                                                gagCCGGTCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGtt                          >  2:755664/1‑140 (MQ=255)                                                                                                       tCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGcc                   >  1:414343/1‑140 (MQ=255)                                                                                                       tCAGCGCGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGcc                   >  1:126614/1‑140 (MQ=255)                                                                                                            gcGTTAACGTTGAACTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGa              <  2:609460/139‑1 (MQ=255)                                                                                                                        aaCTGGGCGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGATTGCGCAGCATg  >  1:576557/1‑140 (MQ=255)                                                                                                                               |                                                                                                                                     TACGACCGCGTACCAGCGGCGCGACGATGATGGAAAGCGCAGAGGCAATCGCGGTATCGCAGTGACCACCGGAAGCACCACGCAGTACGCCGTCAGAGCCGGTCAGCACGTTAACGTTGAACTGGGTGTCAATTTCCAGCGCGCTCAGTACCACCACGTCGAGACGATCAACCGATGCGCCTTTCGAACCCCAGTTAGCGTACTGGTTGGCGCTGATTTCGATGTGATTGGGGTTACGGGCCAGCGATTGCGCAGCATG  >  NC_000913/648320‑648578

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