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

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169                                               |                                              cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             <  1:4895249/49‑1 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:1776449/49‑1 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:1032069/1‑50 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:2694661/50‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:3213461/50‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:4990749/50‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:2618912/1‑49 (MQ=255)     aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:4408937/50‑1 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       <  1:1877036/50‑1 (MQ=255)       ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      <  1:5730759/50‑1 (MQ=255)       ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      <  1:1335776/50‑1 (MQ=255)       ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      >  1:1975798/1‑50 (MQ=255)       ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      <  1:3560034/50‑1 (MQ=255)        cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:4504264/50‑1 (MQ=255)        cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:6110261/50‑1 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:49294/1‑50 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:5737175/50‑1 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:377088/50‑1 (MQ=255)           gAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:2430637/1‑48 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:6003339/50‑1 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:6330187/50‑1 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:6198900/50‑1 (MQ=255)             aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 >  1:2430981/1‑49 (MQ=255)             aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                <  1:3870159/50‑1 (MQ=255)             aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                <  1:1019698/50‑1 (MQ=255)              tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               <  1:782281/50‑1 (MQ=255)                tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              <  1:1219594/49‑1 (MQ=255)                 ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:323706/1‑49 (MQ=255)                 ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:3433340/1‑49 (MQ=255)                 ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:5017868/1‑49 (MQ=255)                   gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                          <  1:5782762/50‑1 (MQ=255)                    ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:2560239/1‑48 (MQ=255)                    ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                         <  1:5457571/50‑1 (MQ=255)                     tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        >  1:103428/1‑50 (MQ=255)                     tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        >  1:3949889/1‑50 (MQ=255)                      aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        <  1:3059570/49‑1 (MQ=255)                      aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        <  1:5814233/49‑1 (MQ=255)                        aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:5094125/50‑1 (MQ=255)                        aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:3745602/50‑1 (MQ=255)                        aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:2468256/50‑1 (MQ=255)                          cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:5799902/1‑50 (MQ=255)                          cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   <  1:1450162/50‑1 (MQ=255)                             aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGGCAt                <  1:5539234/50‑1 (MQ=255)                             aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 >  1:1888438/1‑49 (MQ=255)                             aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt                <  1:3271682/50‑1 (MQ=255)                             aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt                <  1:4508608/50‑1 (MQ=255)                               cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              <  1:676349/50‑1 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:5351501/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:900326/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:3848862/50‑1 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:2183297/1‑50 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:3571486/1‑49 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:615063/1‑49 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:4421447/50‑1 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:2279569/50‑1 (MQ=255)                                  aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            >  1:699116/1‑49 (MQ=255)                                  aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            >  1:5796064/1‑49 (MQ=255)                                   cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          <  1:1986071/50‑1 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:5311819/50‑1 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:1882454/1‑49 (MQ=255)                                      ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       <  1:1613883/50‑1 (MQ=255)                                      ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       <  1:2404340/50‑1 (MQ=255)                                        gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg      <  1:1943073/49‑1 (MQ=255)                                         tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    <  1:2925477/50‑1 (MQ=255)                                          aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:2564844/1‑49 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:5719480/49‑1 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  <  1:4244210/50‑1 (MQ=255)                                            gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:4890427/1‑49 (MQ=255)                                            gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:1550814/1‑49 (MQ=255)                                               |                                              CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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