Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	649,472	G→A	Q34* (CAA→TAA)	citF ←	citrate lyase alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	649,472	0	G	A	100.0%	133.2 / NA	38	Q34* (CAA→TAA)	citF	citrate lyase alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/12);  total (26/12)

AAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCGTTCTTGTCGTT  >  NC_000913/649396‑649546                                                                             |                                                                           aaaCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACttt                                                                        >  2:123235/1‑81 (MQ=255) aaaCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACttt                                                                        >  1:123239/1‑81 (MQ=255)  aaCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACttt                                                                        <  2:123200/80‑1 (MQ=255)  aaCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACttt                                                                        <  1:123198/80‑1 (MQ=255)  aaCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTg                                                                       <  2:123199/81‑1 (MQ=255)   aCCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTgg                                                                      >  2:123248/1‑81 (MQ=255)    ccAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGc                                                                     >  2:123251/1‑81 (MQ=255)     cAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGc                                                                     >  1:123243/1‑80 (MQ=255)       gagCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAg                                                                  >  1:123250/1‑81 (MQ=255)        agCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGt                                                                 >  1:123242/1‑81 (MQ=255)           gacgaATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCt                                                              >  1:123230/1‑81 (MQ=255)             cgaATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTgg                                                            <  2:123201/81‑1 (MQ=255)              gaATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTgg                                                            >  1:123257/1‑80 (MQ=255)              gaATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGa                                                           <  1:123202/81‑1 (MQ=255)                  cGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGc                                                       >  2:123232/1‑81 (MQ=255)                         tcCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGc                                                    >  2:123244/1‑77 (MQ=255)                         tcCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGc                                                    >  1:123249/1‑77 (MQ=255)                          ccAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGc                                                    >  1:123236/1‑76 (MQ=255)                                tGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCAtt                                         >  1:123233/1‑81 (MQ=255)                                 ggCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCAtt                                         >  1:123245/1‑80 (MQ=255)                                  gcgcACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCa                                       <  1:123203/81‑1 (MQ=255)                                      acaGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAgcgc                                   >  2:123255/1‑81 (MQ=255)                                       caGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAgcgc                                   >  1:123240/1‑80 (MQ=255)                                               gATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTcc                          <  1:123204/81‑1 (MQ=255)                                               gATCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTc                           >  1:123263/1‑80 (MQ=255)                                                aTCGCGCGCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCa                         <  1:123205/81‑1 (MQ=255)                                                  cgcgcgCTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAgg                       >  1:123282/1‑81 (MQ=255)                                                  cgcgcgCTTTTTCAGCCTGGTAGGATTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGCCGATTCCAgg                       >  2:123247/1‑81 (MQ=255)                                                      cgcTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGc                   >  1:123270/1‑81 (MQ=255)                                                       gcTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGc                   >  1:123281/1‑80 (MQ=255)                                                        cTTTTTCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTa                 <  2:123206/81‑1 (MQ=255)                                                           tttCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTAcc               <  2:123207/80‑1 (MQ=255)                                                            ttCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCg             >  2:123259/1‑81 (MQ=255)                                                             tCAGCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCg             <  2:123208/80‑1 (MQ=255)                                                                gCCTGGTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCGTTCt         <  2:123209/81‑1 (MQ=255)                                                                    ggTAGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCGTTCTTGTc     >  1:123277/1‑81 (MQ=255)                                                                       aGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCGTTCTTGTCGtt  >  1:123258/1‑81 (MQ=255)                                                                       aGGTTTACTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCGTTCTTGTCGtt  >  2:123262/1‑81 (MQ=255)                                                                             |                                                                           AAACCAGAGCGACGAATCGCTTCTTCCAGGTTGGCGCACAGTTTGCGATCGCGCGCTTTTTCAGCCTGGTAGGTTTGCTTTGGCGAGTTCTGGAAAGCGGCAAGATCGCATTCAGCGCGACGATTCCAGGCCGCTACCCGTTCTTGTCGTT  >  NC_000913/649396‑649546

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