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

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169                                               |                                              cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             <  1:934158/49‑1 (MQ=255) cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             >  1:1449883/1‑49 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:867592/50‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:292453/50‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:1532677/1‑49 (MQ=255)    aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         <  1:2705420/50‑1 (MQ=255)    aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         <  1:2219480/50‑1 (MQ=255)     aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        >  1:2367996/1‑50 (MQ=255)     aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        >  1:1058779/1‑50 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       <  1:345113/50‑1 (MQ=255)        cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:2068881/50‑1 (MQ=255)             aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                <  1:679383/50‑1 (MQ=255)              tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               <  1:218980/50‑1 (MQ=255)                tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              <  1:1237590/49‑1 (MQ=255)                 ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             <  1:2233150/49‑1 (MQ=255)                      aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       >  1:1625229/1‑50 (MQ=255)                         aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                    <  1:434703/50‑1 (MQ=255)                          cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:959067/1‑50 (MQ=255)                             aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt                <  1:645209/50‑1 (MQ=255)                              aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:1902604/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:728422/1‑50 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:1834684/1‑49 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:1728629/50‑1 (MQ=255)                                   cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          <  1:1437042/50‑1 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:450013/1‑49 (MQ=255)                                     aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:688938/1‑48 (MQ=255)                                          aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:1858618/1‑49 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  <  1:1311259/50‑1 (MQ=255)                                            gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:351920/1‑49 (MQ=255)                                               |                                              CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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