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

CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169                                               |                                              cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             >  1:390208/1‑49 (MQ=255) cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             <  1:6573813/49‑1 (MQ=255) cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             <  1:198940/49‑1 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:3399132/49‑1 (MQ=255)  aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:1288852/49‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:653050/50‑1 (MQ=255)   aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:6431394/50‑1 (MQ=255)    aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         <  1:1013764/50‑1 (MQ=255)    aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         >  1:2506699/1‑50 (MQ=255)    aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                         <  1:168888/50‑1 (MQ=255)     aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:6656821/50‑1 (MQ=255)     aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:4043515/50‑1 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       <  1:3468109/50‑1 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       <  1:2019455/50‑1 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       <  1:1983680/50‑1 (MQ=255)      aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                       >  1:5264032/1‑50 (MQ=255)        cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:2316573/50‑1 (MQ=255)        cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:5227707/50‑1 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:5432044/50‑1 (MQ=255)         gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:6289937/1‑50 (MQ=255)          tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                   <  1:2553592/50‑1 (MQ=255)            aaTTTCCGTTAGAGCAAAACAGACGATGGTATGCCCGGGATATTCTACAg                                 <  1:2212405/50‑1 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:5284636/1‑47 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:1180826/1‑47 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:4343112/1‑47 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                  <  1:2295682/49‑1 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                  <  1:2801772/49‑1 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:1765571/50‑1 (MQ=255)            aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:3826638/50‑1 (MQ=255)              tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAgc                                >  1:2701520/1‑48 (MQ=255)              tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                                >  1:892628/1‑49 (MQ=255)               ttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               >  1:4547855/1‑49 (MQ=255)                tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              >  1:1419036/1‑49 (MQ=255)                tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              <  1:2276584/49‑1 (MQ=255)                tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              <  1:3453075/49‑1 (MQ=255)                  cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaa                            <  1:5178825/49‑1 (MQ=255)                   gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                          <  1:5190156/50‑1 (MQ=255)                   gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                          <  1:3582500/50‑1 (MQ=255)                    ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                         <  1:4614400/50‑1 (MQ=255)                      aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       >  1:3434205/1‑50 (MQ=255)                       gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:5251394/50‑1 (MQ=255)                       gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      >  1:1160500/1‑50 (MQ=255)                       gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:3273440/50‑1 (MQ=255)                       gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:2132885/50‑1 (MQ=255)                        aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:3422899/50‑1 (MQ=255)                        aaCAAAACAGACGATGGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:1226871/50‑1 (MQ=255)                         aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                    <  1:744437/50‑1 (MQ=255)                          cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   <  1:1461903/50‑1 (MQ=255)                           aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:6660018/1‑49 (MQ=255)                           aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:3690619/50‑1 (MQ=255)                           aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:181180/50‑1 (MQ=255)                            aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 >  1:2446908/1‑50 (MQ=255)                             aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 >  1:3125957/1‑49 (MQ=255)                             aaCAGACGATTGTATGCACGGGATATTCTACAGAAAAAACAGTCTGACa                 >  1:380157/1‑49 (MQ=39)                              aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               >  1:289053/1‑50 (MQ=255)                               cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              <  1:1922407/50‑1 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:1965570/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:1310764/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:6363759/50‑1 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:6709915/50‑1 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:5883072/1‑50 (MQ=255)                                aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:2668075/50‑1 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:6061604/1‑49 (MQ=255)                                 gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:5714641/50‑1 (MQ=255)                                  aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            >  1:3236523/1‑49 (MQ=255)                                  aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa           <  1:3435836/50‑1 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:6786677/50‑1 (MQ=255)                                    gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:3864034/50‑1 (MQ=255)                                     aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa        <  1:718969/50‑1 (MQ=255)                                      ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:9527/1‑50 (MQ=255)                                        gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg      <  1:4174469/49‑1 (MQ=255)                                        gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg      <  1:2703657/49‑1 (MQ=255)                                        gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg      <  1:828997/49‑1 (MQ=255)                                        gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa     >  1:3870927/1‑50 (MQ=255)                                        gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa     >  1:7067496/1‑50 (MQ=255)                                         tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGaa    <  1:782646/50‑2 (MQ=255)                                         tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:1350552/1‑50 (MQ=255)                                         tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt    >  1:1672779/1‑50 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:5224190/49‑1 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:3143801/49‑1 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:4335529/49‑1 (MQ=255)                                           tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  <  1:1531146/50‑1 (MQ=255)                                            gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:6591008/1‑49 (MQ=255)                                               |                                              CAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438078‑4438169

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