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%	216.8 / NA	55	W140R (TGG→CGG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/28);  total (27/28)

TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATT  >  NC_000913/4438077‑4438168                                                |                                             tCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                            >  1:2100437/1‑50 (MQ=255)  cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                            <  1:4368158/49‑1 (MQ=255)  cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                            <  1:1918950/49‑1 (MQ=255)   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                          >  1:4415709/1‑50 (MQ=255)   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                          >  1:1466613/1‑50 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                         <  1:2157571/50‑1 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                          >  1:3591161/1‑49 (MQ=255)      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                       >  1:377120/1‑50 (MQ=255)      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                       <  1:1944237/50‑1 (MQ=255)      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                       <  1:644470/50‑1 (MQ=255)       aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                      <  1:706550/50‑1 (MQ=255)       aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                      >  1:632179/1‑50 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                    <  1:971915/50‑1 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                    <  1:563954/50‑1 (MQ=255)           tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                   >  1:2273915/1‑49 (MQ=255)              aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                >  1:283111/1‑49 (MQ=255)                ttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                              >  1:2440825/1‑49 (MQ=255)                 tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                             <  1:3021972/49‑1 (MQ=255)                  ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                            <  1:2644845/49‑1 (MQ=255)                  ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                            >  1:3050323/1‑49 (MQ=255)                   cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                          <  1:4902893/50‑1 (MQ=255)                    gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                         <  1:4861081/50‑1 (MQ=255)                      tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                       >  1:2375605/1‑50 (MQ=255)                      tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                       >  1:3476026/1‑50 (MQ=255)                       aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                       <  1:1021809/49‑1 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                      <  1:2438915/49‑1 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                      <  1:913941/49‑1 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                      <  1:4351562/49‑1 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                     <  1:3473762/50‑1 (MQ=255)                         aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAGAAACAGTCt                    <  1:1147446/50‑1 (MQ=255)                         aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                    <  1:2723104/50‑1 (MQ=255)                          aCAAACCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                   <  1:2868601/50‑1 (MQ=255)                          aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                   <  1:3552944/50‑1 (MQ=255)                          aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                   <  1:4632990/50‑1 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                  >  1:1435412/1‑49 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                  >  1:4056719/1‑49 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAGAACAGTCTGACa                <  1:3580132/50‑1 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                >  1:3669098/1‑50 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                >  1:4293620/1‑50 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                >  1:1878619/1‑50 (MQ=255)                              aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                >  1:3726388/1‑49 (MQ=255)                              aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                >  1:2907806/1‑49 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa              <  1:385754/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc            >  1:3152393/1‑50 (MQ=255)                                  gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg           <  1:1493345/50‑1 (MQ=255)                                   aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg           >  1:3452366/1‑49 (MQ=255)                                     gATGGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa        <  1:4902863/50‑1 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa         >  1:3001344/1‑48 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa         >  1:1672235/1‑48 (MQ=255)                                       ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg      <  1:1105801/50‑1 (MQ=255)                                        tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg      >  1:1074137/1‑49 (MQ=255)                                          tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa    <  1:2553646/49‑1 (MQ=255)                                           aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt   >  1:4873004/1‑49 (MQ=255)                                           aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt   >  1:4634488/1‑49 (MQ=255)                                            tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt  >  1:3742797/1‑49 (MQ=255)                                                |                                             TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATT  >  NC_000913/4438077‑4438168

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