Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,038,809	A→G	A216A (GCA→GCG)	appC →	cytochrome bd‑II oxidase, subunit I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,038,809	0	A	G	100.0%	129.1 / NA	38	A216A (GCA→GCG)	appC	cytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/16);  total (22/16)

TTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCG  >  W3110S.gb/1038753‑1038873                                                         |                                                                 tttATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGc                                                    >  1:4401597/1‑71 (MQ=255) tttATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGc                                                    >  1:405943/1‑71 (MQ=255) tttATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGc                                                    >  1:3481579/1‑71 (MQ=255) tttATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGc                                                    >  1:2541815/1‑71 (MQ=255) tttATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGc                                                    >  1:2426849/1‑71 (MQ=255)  ttattaTGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTg                                                     >  1:4227386/1‑69 (MQ=255)   tattatGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCa                                                  >  1:4363999/1‑71 (MQ=255)               cAGCGCCTGGTATTTACTGCGCGGACGGGAGCGAAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCt                                      >  1:4513933/1‑71 (MQ=255)                   gcCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:1048688/71‑1 (MQ=255)                   gcCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:1837410/71‑1 (MQ=255)                    ccTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:357138/70‑1 (MQ=255)                        gTATTTACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTAGGTAc                               <  1:3230204/69‑1 (MQ=255)                           tttACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGt                                             >  1:3128262/1‑52 (MQ=255)                           tttACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:3481351/63‑1 (MQ=255)                           tttACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:2930331/63‑1 (MQ=255)                           tttACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:210532/63‑1 (MQ=255)                           tttACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:567467/63‑1 (MQ=255)                           tttACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTg                           >  1:754015/1‑70 (MQ=255)                             tACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCtt                                     >  1:1791224/1‑58 (MQ=255)                             tACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  >  1:3236479/1‑61 (MQ=255)                             tACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGc                         >  1:1804067/1‑70 (MQ=255)                             tACTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGc                         >  1:1825898/1‑70 (MQ=255)                              aCTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGa                       >  1:2826444/1‑71 (MQ=255)                              aCTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGa                       >  1:2446282/1‑71 (MQ=255)                               cTGCGCGGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGAt                      <  1:2744776/71‑1 (MQ=255)                                  cgcgGCCGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTc                                    <  1:3575367/54‑1 (MQ=255)                                  cgcgGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTAt                   >  1:1497463/1‑71 (MQ=255)                                   gcgGACGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATc                  <  1:3998976/71‑1 (MQ=255)                                       aCGGGAGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTc                                       <  1:587993/46‑1 (MQ=255)                                            aGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCgg                                  <  1:3160902/46‑1 (MQ=255)                                            aGCGCAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTg         >  1:1136286/1‑71 (MQ=255)                                             gcgcAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGAt                      <  1:4176860/57‑1 (MQ=255)                                             gcgcAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTg         >  1:1333431/1‑70 (MQ=255)                                             gcgcAATGTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGc        <  1:1628871/71‑1 (MQ=255)                                                   tGTCGCGTTACGCTCGTTTGCCATCGGTTCCGGCTTCGGTAc                               <  1:4091086/42‑1 (MQ=38)                                                    gTCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCg  >  1:430500/1‑70 (MQ=255)                                                     tCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCg  >  1:3223336/1‑69 (MQ=255)                                                     tCGCGTTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCg  >  1:985039/1‑69 (MQ=255)                                                         |                                                                 TTTATTATGGCGATCAGCGCCTGGTATTTACTGCGCGGACGGGAGCGCAATGTCGCATTACGCTCGTTTGCCATCGGTTCCGTCTTCGGTACTCTGGCGATTATCGGTACCCTGCAACTCG  >  W3110S.gb/1038753‑1038873

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