Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	616,523	G→A	A789A (GCG→GCA)	entF →	apo‑serine activating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	616,523	0	G	A	100.0%	76.8 / NA	22	A789A (GCG→GCA)	entF	apo‑serine activating enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/12);  total (10/12)

CGCGGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTC  >  NC_000913/616461‑616598                                                               |                                                                            cgcgGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTgccg                                                           >  1:127992/1‑81 (MQ=255) cgcgGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTgcc                                                            <  2:127955/80‑1 (MQ=255)   cgGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTgccgcc                                                         <  2:127957/81‑1 (MQ=255)   cgGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTgccgcc                                                         <  2:127956/81‑1 (MQ=255)       agcagTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGgt                                                     <  2:127958/81‑1 (MQ=255)         cagTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGgtgt                                                   <  1:127959/81‑1 (MQ=255)              tgCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCg                                               <  1:127960/80‑1 (MQ=255)                  gATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGa                                          >  1:128003/1‑81 (MQ=255)                       gTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATctct                                     <  1:127961/81‑1 (MQ=255)                          aTCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTAtc                                  >  1:128012/1‑81 (MQ=255)                             cGGTATGGAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTAtctc                                <  2:127962/80‑1 (MQ=255)                                   ggAATACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCa                         >  2:128004/1‑81 (MQ=255)                                     aaTACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCAtt                       >  1:128007/1‑81 (MQ=255)                                     aaTACGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCAtt                       >  1:128015/1‑81 (MQ=255)                                        aCGGGTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCaa                    >  1:128002/1‑81 (MQ=255)                                            gTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTgg                >  2:128001/1‑81 (MQ=255)                                            gTCTGCGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTgg                >  1:128013/1‑81 (MQ=255)                                                 cGTATTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAg           <  1:127963/81‑1 (MQ=255)                                                    aTTCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGCTCTCTATCTCACTGGCCTTCTACTGgcgcc            >  2:128016/1‑76 (MQ=255)                                                     ttCTTGATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCt       <  1:127964/81‑1 (MQ=255)                                                          gATGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTAtc    <  2:127965/79‑1 (MQ=255)                                                           aTGCAATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTAtctc  <  1:127966/80‑1 (MQ=255)                                                               |                                                                            CGCGGCAGCAGTGTGCCGATTGGTTATCCGGTATGGAATACGGGTCTGCGTATTCTTGATGCGATGATGCATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTC  >  NC_000913/616461‑616598

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