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

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169                                               |                                              cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             <  1:6033688/49‑1 (MQ=255) cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             >  1:4593328/1‑49 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:345521/49‑1 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:4960578/49‑1 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:5469924/49‑1 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:2790650/1‑50 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:797761/50‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:6357258/50‑1 (MQ=255)    aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         >  1:5784952/1‑50 (MQ=255)     aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:1679738/50‑1 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       >  1:3082213/1‑50 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       >  1:3275277/1‑50 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       >  1:4947572/1‑50 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       >  1:1808248/1‑50 (MQ=255)       ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      <  1:4551158/50‑1 (MQ=255)       ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      >  1:137035/1‑50 (MQ=255)        cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:3975615/50‑1 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:4905441/50‑1 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:1405016/1‑50 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:695767/50‑1 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:282129/50‑1 (MQ=255)          tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:2053396/1‑49 (MQ=255)          tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:986985/1‑49 (MQ=255)          tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                   <  1:5640667/50‑1 (MQ=255)          tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                   <  1:1173423/50‑1 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:4149454/50‑1 (MQ=255)                 ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             <  1:1949946/49‑1 (MQ=255)                 ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:413583/1‑49 (MQ=255)                  cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           <  1:3908578/50‑1 (MQ=255)                  cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaa                            <  1:1300588/49‑1 (MQ=255)                  cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaa                            <  1:1943864/49‑1 (MQ=255)                  cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaa                            <  1:2476182/49‑1 (MQ=255)                   gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:5090908/1‑49 (MQ=255)                   gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                          <  1:3308569/50‑1 (MQ=255)                   gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                          <  1:6702005/50‑1 (MQ=255)                      aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       >  1:4694386/1‑50 (MQ=255)                       gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:1715386/50‑1 (MQ=255)                       gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:5096179/50‑1 (MQ=255)                        aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:2392878/50‑1 (MQ=255)                         aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                    <  1:1216261/50‑1 (MQ=255)                          cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   <  1:616783/50‑1 (MQ=255)                          cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:4112690/1‑50 (MQ=255)                           aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:2501964/1‑49 (MQ=255)                           aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:2107187/1‑49 (MQ=255)                           aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:4413773/50‑1 (MQ=255)                           aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:4821264/50‑1 (MQ=255)                              aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               <  1:3778160/50‑1 (MQ=255)                               cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              >  1:1918145/1‑50 (MQ=255)                               cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              >  1:5111361/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:2602249/50‑1 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:5020283/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:2974661/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:3977955/50‑1 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:4857987/50‑1 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:441868/50‑1 (MQ=255)                                  aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            >  1:6714874/1‑49 (MQ=255)                                   cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa           >  1:98313/1‑49 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:3619369/50‑1 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:4680838/50‑1 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:4836039/1‑49 (MQ=255)                                     aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa        <  1:1209644/50‑1 (MQ=255)                                       tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:5972350/1‑49 (MQ=255)                                       tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:956373/1‑49 (MQ=255)                                        gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg      <  1:4939447/49‑1 (MQ=255)                                         tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa     <  1:5196627/49‑1 (MQ=255)                                         tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:5606044/1‑50 (MQ=255)                                         tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:2598417/1‑50 (MQ=255)                                          aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:117818/1‑49 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:1106538/49‑1 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  <  1:4932963/50‑1 (MQ=255)                                               |                                              CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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