Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,438,124	T→C	W140R (TGG→CGG)	ytfI →	protein YtfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,438,124	0	T	C	100.0%	78.1 / NA	21	W140R (TGG→CGG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (6/15);  total (6/15)

AAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATT  >  NC_000913/4438079‑4438168                                              |                                             aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                           <  1:2202292/49‑1 (MQ=255) aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                          <  1:1908941/50‑1 (MQ=255)  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                         <  1:401915/50‑1 (MQ=255)  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                         <  1:1970662/50‑1 (MQ=255)      ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                     >  1:319351/1‑50 (MQ=255)      ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                     <  1:2188700/50‑1 (MQ=255)           aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                <  1:2067344/50‑1 (MQ=255)             tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                              <  1:1423913/50‑1 (MQ=255)               tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                             >  1:1753197/1‑49 (MQ=255)                ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                            <  1:1542576/49‑1 (MQ=255)                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                          >  1:1541806/1‑49 (MQ=255)                    tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                       >  1:453075/1‑50 (MQ=255)                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                      >  1:1041119/1‑50 (MQ=255)                        aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                   <  1:2114972/50‑1 (MQ=255)                           aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                <  1:1103107/50‑1 (MQ=255)                              cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc             <  1:1720043/50‑1 (MQ=255)                               aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc            <  1:1892494/50‑1 (MQ=255)                                gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg           <  1:925318/50‑1 (MQ=255)                                     ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg      <  1:316743/50‑1 (MQ=255)                                      tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg      >  1:1206593/1‑49 (MQ=255)                                          tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt  <  1:1414427/49‑1 (MQ=255)                                              |                                             AAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATT  >  NC_000913/4438079‑4438168

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