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

TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAT  >  NC_000913/4438077‑4438167                                                |                                            tCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                           >  1:4319723/1‑50 (MQ=255)  cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                           >  1:489517/1‑49 (MQ=255)   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                         >  1:7245954/1‑50 (MQ=255)   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                         >  1:5730954/1‑50 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                        <  1:6632094/50‑1 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                        <  1:5496217/50‑1 (MQ=255)     aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                       <  1:467143/50‑1 (MQ=255)     aaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGata                                       <  1:2649616/50‑1 (MQ=255)      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                      >  1:2476855/1‑50 (MQ=255)       aCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATAtt                                     >  1:3618493/1‑50 (MQ=255)        ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                    >  1:1756177/1‑50 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                   <  1:2164670/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                  <  1:987604/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                  <  1:6321527/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                  <  1:3709039/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                  >  1:200877/1‑50 (MQ=255)             aaTTTCCGTTAGAACAAAACAGGCGATTGTATGCCCGGGATATTCTACAg                               <  1:1829509/50‑1 (MQ=255)             aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                  >  1:4474496/1‑47 (MQ=255)             aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                <  1:6826426/49‑1 (MQ=255)             aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                               <  1:3204483/50‑1 (MQ=255)               tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                              >  1:7308346/1‑49 (MQ=255)               tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGa                              >  1:5770991/1‑49 (MQ=255)                ttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                             >  1:3099991/1‑49 (MQ=255)                 tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                            <  1:5798237/49‑1 (MQ=255)                 tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                            <  1:4494826/49‑1 (MQ=255)                  ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                           >  1:2873027/1‑49 (MQ=255)                   cGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                         >  1:3411286/1‑50 (MQ=255)                    gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAAc                        <  1:1417092/50‑1 (MQ=255)                     ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                         >  1:2871822/1‑48 (MQ=255)                     ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                       <  1:3904786/50‑1 (MQ=255)                      tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                      >  1:1103083/1‑50 (MQ=255)                       aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                     >  1:1571965/1‑50 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                    <  1:265375/50‑1 (MQ=255)                         aaCAAAACAGACGGTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                   <  1:1408726/50‑1 (MQ=255)                         aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                   <  1:5492001/50‑1 (MQ=255)                         aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                   <  1:1942098/50‑1 (MQ=255)                          aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                   >  1:3136435/1‑49 (MQ=255)                          aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                  <  1:3826943/50‑1 (MQ=255)                          aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                  <  1:5918620/50‑1 (MQ=255)                           cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                 <  1:6262704/50‑1 (MQ=255)                           cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                 <  1:4780841/50‑1 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                <  1:2988032/50‑1 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa               >  1:5528003/1‑50 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa               <  1:5725575/50‑1 (MQ=255)                              aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa               >  1:6570012/1‑49 (MQ=255)                              aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt              <  1:7108940/50‑1 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa             <  1:4812503/50‑1 (MQ=255)                                cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc            >  1:7102521/1‑50 (MQ=255)                                cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc            >  1:3529831/1‑50 (MQ=255)                                cAGACGATTGTATGCCCGGGATAGTCTACAGAAAAAACAGTCTGACATAc            <  1:5989230/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc           <  1:5770584/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc           >  1:6899332/1‑50 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc           <  1:4371762/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc           <  1:6562412/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc           <  1:3732522/50‑1 (MQ=255)                                  gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc           >  1:3000624/1‑49 (MQ=255)                                  gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc           >  1:398666/1‑49 (MQ=255)                                  gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg          <  1:7373155/50‑1 (MQ=255)                                   aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg          >  1:3788737/1‑49 (MQ=255)                                   aCGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg          >  1:4301633/1‑49 (MQ=255)                                     gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa       <  1:5762432/50‑1 (MQ=255)                                     gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa        >  1:3924538/1‑49 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa      <  1:732714/50‑1 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa        >  1:1906209/1‑48 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa        >  1:7197129/1‑48 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa        >  1:2940955/1‑48 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa        >  1:750018/1‑48 (MQ=255)                                       ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg     <  1:3892005/50‑1 (MQ=255)                                        tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg     >  1:2928068/1‑49 (MQ=255)                                         gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg    <  1:2444721/49‑1 (MQ=255)                                         gTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAgg    <  1:2653812/49‑1 (MQ=255)                                           aTGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAt  >  1:5996626/1‑49 (MQ=255)                                                |                                            TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAT  >  NC_000913/4438077‑4438167

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