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

AAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATT  >  NC_000913/4438079‑4438168                                              |                                             aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                            >  1:1075489/1‑48 (MQ=255)  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                          >  1:799921/1‑49 (MQ=255)  aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                          <  1:445988/49‑1 (MQ=255)       cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                     >  1:1659382/1‑49 (MQ=255)        gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                     >  1:1800075/1‑48 (MQ=255)        gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                    <  1:1917806/49‑1 (MQ=255)         tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                   <  1:1964188/49‑1 (MQ=255)           aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                  <  1:831246/48‑1 (MQ=255)           aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                 <  1:1746898/49‑1 (MQ=255)             tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                >  1:1346736/1‑48 (MQ=255)               tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                              <  1:1693045/48‑1 (MQ=255)                  gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                          >  1:587531/1‑49 (MQ=255)                   ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                         <  1:1078593/49‑1 (MQ=255)                    tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                          >  1:62078/1‑47 (MQ=255)                     aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                        <  1:1772615/48‑1 (MQ=255)                      gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                       <  1:22094/48‑1 (MQ=255)                       aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                     <  1:1402266/49‑1 (MQ=255)                       aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                     <  1:1394865/49‑1 (MQ=255)                        aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                    >  1:1731788/1‑49 (MQ=255)                         cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                   <  1:327377/49‑1 (MQ=255)                          aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                  <  1:1060020/49‑1 (MQ=255)                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                >  1:1247358/1‑48 (MQ=255)                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt               <  1:1812721/49‑1 (MQ=255)                             aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt               <  1:1374414/49‑1 (MQ=255)                              cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa              <  1:1307139/49‑1 (MQ=255)                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa         <  1:46818/49‑1 (MQ=255)                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa          >  1:1088628/1‑48 (MQ=255)                                   gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa          >  1:1610455/1‑48 (MQ=255)                                    aTTGTATACCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa         >  1:1608795/1‑48 (MQ=39)                                     ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa         >  1:536708/1‑47 (MQ=255)                                       gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg      <  1:1876066/48‑1 (MQ=255)                                           gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt  >  1:528787/1‑48 (MQ=255)                                           gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt  >  1:766197/1‑48 (MQ=255)                                              |                                             AAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATT  >  NC_000913/4438079‑4438168

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: ‑