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

TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438077‑4438169                                                |                                              tCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             >  1:1878306/1‑50 (MQ=255)   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:4105369/1‑50 (MQ=255)      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:3819536/50‑1 (MQ=255)       aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       >  1:4065918/1‑50 (MQ=255)        ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      >  1:144335/1‑50 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:3178049/50‑1 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:1102109/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:2089542/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:1930084/1‑50 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:4286861/1‑50 (MQ=255)           tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                   <  1:1503037/50‑1 (MQ=255)           tGAATTTCCGTAAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                   <  1:2674905/50‑1 (MQ=255)            gAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:1073104/1‑48 (MQ=255)             aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                  <  1:2784136/49‑1 (MQ=255)              aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 >  1:302477/1‑49 (MQ=255)              aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                <  1:1274553/50‑1 (MQ=255)                ttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               >  1:1234814/1‑49 (MQ=255)                  ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:3406924/1‑49 (MQ=255)                   cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:1569714/1‑50 (MQ=255)                   cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           <  1:2569384/50‑1 (MQ=255)                   cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:2378719/1‑50 (MQ=255)                     ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                         <  1:4204218/50‑1 (MQ=255)                       aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       >  1:4215372/1‑50 (MQ=255)                       aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       >  1:180745/1‑50 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      >  1:1217731/1‑50 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      >  1:1145982/1‑50 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:1839430/50‑1 (MQ=255)                         aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     >  1:1812173/1‑50 (MQ=255)                           cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:1077687/1‑50 (MQ=255)                           cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:458251/1‑50 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:4125458/1‑49 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:846978/50‑1 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 <  1:629207/50‑1 (MQ=255)                              aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt                <  1:3853703/50‑1 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               <  1:1283302/50‑1 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:3097145/1‑50 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:3261817/1‑50 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:4078218/1‑50 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:3455044/1‑50 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:1331241/1‑50 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:1572429/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:822815/1‑50 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:1886912/50‑1 (MQ=255)                                  gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:1251911/1‑49 (MQ=255)                                   aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            >  1:2502467/1‑49 (MQ=255)                                   aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            >  1:2465745/1‑49 (MQ=255)                                    cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          <  1:2302333/50‑1 (MQ=255)                                    cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          <  1:2429345/50‑1 (MQ=255)                                       ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:1863031/1‑50 (MQ=255)                                       ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       <  1:310852/50‑1 (MQ=255)                                        tGTATGCCCGGGATATTCTACAGAAAAAACCGTCTGACATACCGCAAAg       >  1:4074773/1‑49 (MQ=255)                                        tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:1000575/1‑49 (MQ=255)                                         gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg      <  1:4211794/49‑1 (MQ=255)                                           aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:3866775/1‑49 (MQ=255)                                            tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:1498815/49‑1 (MQ=255)                                            tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:1677301/49‑1 (MQ=255)                                             gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:619932/1‑49 (MQ=255)                                                |                                              TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438077‑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: ‑