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

TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438077‑4438169                                                |                                              tCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             >  1:1165284/1‑50 (MQ=255)  cAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCgg                                             <  1:955755/49‑1 (MQ=255)   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCggg                                            <  1:5085764/49‑1 (MQ=255)   aaaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           <  1:3696672/50‑1 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:3365659/50‑1 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:3886843/50‑1 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGat                                          <  1:4604529/50‑1 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGa                                           >  1:358967/1‑49 (MQ=255)    aaaaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATACCCGGGa                                           >  1:2878117/1‑49 (MQ=38)      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:1168884/50‑1 (MQ=255)      aaCCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGatat                                        <  1:2003937/50‑1 (MQ=255)        ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      >  1:3405647/1‑50 (MQ=255)        ccGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      <  1:2093288/50‑1 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTc                                      >  1:1447812/1‑49 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:266805/50‑1 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:4812813/50‑1 (MQ=255)         cGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCt                                     <  1:4917017/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:212131/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:4108434/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:3886949/1‑50 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:929390/50‑1 (MQ=255)          gTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:966930/50‑1 (MQ=255)          gTGAATCTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    <  1:5575838/50‑1 (MQ=255)           tGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTAc                                   <  1:839496/50‑1 (MQ=255)            gAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:3212088/1‑48 (MQ=255)            gAATTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTa                                    >  1:3196369/1‑48 (MQ=255)             aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                  <  1:5018938/49‑1 (MQ=255)             aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACa                                  <  1:856546/49‑1 (MQ=255)             aaTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 <  1:1950332/50‑1 (MQ=255)              aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 >  1:1109143/1‑49 (MQ=255)              aTTTCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAg                                 >  1:2348604/1‑49 (MQ=255)               tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               <  1:454743/50‑1 (MQ=255)               tttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               <  1:375558/50‑1 (MQ=255)                ttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               >  1:88978/1‑49 (MQ=255)                ttCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaa                               >  1:5286084/1‑49 (MQ=255)                 tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              <  1:1662726/49‑1 (MQ=255)                 tCCGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaa                              >  1:1817906/1‑49 (MQ=255)                  ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             <  1:134942/49‑1 (MQ=255)                  ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:1669091/1‑49 (MQ=255)                  ccGTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaa                             >  1:3034431/1‑49 (MQ=255)                    gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:2632254/1‑49 (MQ=255)                    gTTAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGaaaaaa                           >  1:2271996/1‑49 (MQ=255)                     ttAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACa                         <  1:1766557/50‑1 (MQ=255)                      tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        >  1:2856605/1‑50 (MQ=255)                      tAGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        >  1:4707347/1‑50 (MQ=255)                       aGAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAg                        <  1:3754821/49‑1 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGt                       <  1:5254433/49‑1 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:153816/50‑1 (MQ=255)                        gAACAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTc                      <  1:1651065/50‑1 (MQ=255)                         aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     <  1:955650/50‑1 (MQ=255)                         aaCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     >  1:1923577/1‑50 (MQ=255)                          aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCt                     >  1:1892847/1‑49 (MQ=255)                          aCAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTg                    <  1:3245139/50‑1 (MQ=255)                           cAAAACAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:1843485/1‑50 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:4093634/1‑49 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGa                   >  1:1930732/1‑49 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:1176388/50‑1 (MQ=255)                            aaaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGAc                  <  1:3073929/50‑1 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 <  1:5234980/50‑1 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 >  1:2516429/1‑50 (MQ=255)                             aaaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACa                 <  1:3641876/50‑1 (MQ=255)                              aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt                <  1:4247239/50‑1 (MQ=255)                              aaCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACAt                <  1:3950149/50‑1 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               <  1:2372398/50‑1 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               <  1:135188/50‑1 (MQ=255)                               aCAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATa               <  1:2979636/50‑1 (MQ=255)                                cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              <  1:3675121/50‑1 (MQ=255)                                cAGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAc              <  1:982342/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:1696076/50‑1 (MQ=255)                                 aGACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             <  1:677859/50‑1 (MQ=255)                                  gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATAcc             >  1:5252901/1‑49 (MQ=255)                                  gACGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCg            <  1:973377/50‑1 (MQ=255)                                    cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          <  1:2604288/50‑1 (MQ=255)                                    cGATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGa           >  1:818230/1‑49 (MQ=255)                                     gATTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaa         <  1:5559147/50‑1 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaaaa        <  1:2209564/50‑1 (MQ=255)                                      aTTGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGaa          >  1:705624/1‑48 (MQ=255)                                       ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:4685112/1‑50 (MQ=255)                                       ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       <  1:3658777/50‑1 (MQ=255)                                       ttGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:5096557/1‑50 (MQ=255)                                        tGTATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAg       >  1:1031787/1‑49 (MQ=255)                                          tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa     <  1:2222409/49‑1 (MQ=255)                                          tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa     <  1:914824/49‑1 (MQ=255)                                          tATGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGa     <  1:4295419/49‑1 (MQ=255)                                            tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   <  1:5550522/49‑1 (MQ=255)                                            tGCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGAtt   >  1:19730/1‑49 (MQ=255)                                             gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:78411/1‑49 (MQ=255)                                             gCCCGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTa  >  1:2979468/1‑49 (MQ=255)                                                |                                              TCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGGGATATTCTACAGAAAAAACAGTCTGACATACCGAAAAGGATTA  >  NC_000913/4438077‑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: ‑