Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	649,456	T→C	E39G (GAA→GGA)	citF ←	citrate lyase alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	649,456	0	T	C	100.0%	168.8 / NA	47	E39G (GAA→GGA)	citF	citrate lyase alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/22);  total (25/22)

CGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTA  >  NC_000913/649380‑649531                                                                             |                                                                            cGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGc                                                                         <  1:138903/81‑1 (MQ=255)   tCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCt                                                                       <  1:138904/81‑1 (MQ=255)       gCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTgg                                                                     >  2:138949/1‑79 (MQ=255)       gCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTgg                                                                     >  1:138966/1‑79 (MQ=255)        ccGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTa                                                                   <  2:138906/80‑1 (MQ=255)        ccGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAg                                                                  <  1:138905/81‑1 (MQ=255)        ccGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAg                                                                  >  2:138962/1‑81 (MQ=255)          gTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGt                                                                <  1:138907/81‑1 (MQ=255)          gTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGt                                                                >  1:138970/1‑81 (MQ=255)             cTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGttt                                                              <  2:138908/80‑1 (MQ=255)             cTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTg                                                             <  2:138909/81‑1 (MQ=255)             cTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTg                                                             >  2:138959/1‑81 (MQ=255)               gTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCt                                                           >  2:138974/1‑81 (MQ=255)                tAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGc                                                            <  2:138911/79‑1 (MQ=255)                tAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCtt                                                          <  1:138912/81‑1 (MQ=255)                tAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCt                                                           <  2:138910/80‑1 (MQ=255)                 aaaCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCttt                                                         >  2:138980/1‑81 (MQ=255)                   aCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTgg                                                       >  2:138976/1‑81 (MQ=255)                     cAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCg                                                     <  1:138913/81‑1 (MQ=255)                        agCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGt                                                  >  1:138982/1‑81 (MQ=255)                          cgacgaATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTc                                                >  1:138984/1‑81 (MQ=255)                                aTCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGaaa                                          <  2:138914/81‑1 (MQ=255)                                    cttcttCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCgg                                      >  1:138965/1‑81 (MQ=255)                                    cttcttCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCgg                                      >  1:138985/1‑81 (MQ=255)                                       cttcCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGc                                     <  2:138915/79‑1 (MQ=255)                                       cttcCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCa                                    <  2:138916/80‑1 (MQ=255)                                             ggTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGc                             <  2:138917/81‑1 (MQ=255)                                             ggTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGc                             <  2:138918/81‑1 (MQ=255)                                              gTTGGCGCACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCa                            >  1:138967/1‑81 (MQ=255)                                                  gcgcACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCAtt                          >  1:138951/1‑79 (MQ=255)                                                   cgcACAGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAg                       >  1:138977/1‑81 (MQ=255)                                                     cacaGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAgcg                     >  1:138978/1‑81 (MQ=255)                                                     cacaGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAgcg                     >  2:138975/1‑81 (MQ=255)                                                       caGTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAgcgc                    >  2:138971/1‑80 (MQ=255)                                                         gTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGcgac                 >  2:138979/1‑81 (MQ=255)                                                         gTTTGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGcgac                 <  1:138919/81‑1 (MQ=255)                                                           ttGCGATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGcgacga               >  1:138999/1‑81 (MQ=255)                                                               gATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTc            >  2:138981/1‑80 (MQ=255)                                                               gATCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTc            >  1:138992/1‑80 (MQ=255)                                                                 tCGCGCGCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAg         <  2:138920/81‑1 (MQ=255)                                                                  cgcgcgCTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAgg        >  2:138969/1‑81 (MQ=255)                                                                    cgcgcTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGcc      <  2:138921/81‑1 (MQ=255)                                                                     gcgcTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACAATTCCAGGCCg     <  1:138922/81‑1 (MQ=255)                                                                      cgcTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGc    >  2:138973/1‑81 (MQ=255)                                                                       gcTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGc    >  2:138972/1‑80 (MQ=255)                                                                        cTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTa  <  2:138924/81‑1 (MQ=255)                                                                        cTTTTCCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTa  <  1:138923/81‑1 (MQ=255)                                                                             |                                                                            CGTCATGCCGTCCTGTAAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTA  >  NC_000913/649380‑649531

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