Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	197,387	G→C	5.2%	S281T (AGT→ACT)	yaeL →	zinc metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	197,387	0	G	C	5.2%	233.3 / 2.4	76	S281T (AGT→ACT)	yaeL	zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base C (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATT  >  W3110S.gb/197321‑197457                                                                   |                                                                       ccTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTccc                                                                    <  1:2658221/71‑1 (MQ=255)  cTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCt                                                                   >  1:5470426/1‑71 (MQ=255)    ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:7466262/1‑71 (MQ=255)    ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:7457029/1‑71 (MQ=255)    ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:6567932/1‑71 (MQ=255)    ttGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 >  1:4188778/1‑71 (MQ=255)     tgtgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTg                                                                 <  1:945900/70‑1 (MQ=255)     tgtgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                >  1:3317047/1‑71 (MQ=255)       tgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCt                                                              <  1:3834946/71‑1 (MQ=255)       tgATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCt                                                              <  1:1385580/71‑1 (MQ=255)         aTGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCttt                                                            <  1:4209517/71‑1 (MQ=255)          tGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                >  1:1515477/1‑66 (MQ=255)          tGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCttt                                                            >  1:3452303/1‑70 (MQ=255)           gCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                          >  1:5793100/1‑71 (MQ=255)             ttGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                        >  1:7486774/1‑71 (MQ=255)               gTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                        <  1:4798441/69‑1 (MQ=255)                tCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTccc                                                                    <  1:2919299/56‑1 (MQ=255)                tCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCtt                                                             <  1:2054973/63‑1 (MQ=255)                  cGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAAt                                                   >  1:7789048/1‑71 (MQ=255)                  cGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAAt                                                   <  1:3041545/71‑1 (MQ=255)                     gATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATccc                                                >  1:3979165/1‑71 (MQ=255)                      aTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCg                                               <  1:3207197/71‑1 (MQ=255)                      aTAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGGCCCTTGTCTTTGACATTAATCCCg                                               <  1:3427132/71‑1 (MQ=255)                          cccGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             >  1:2190408/1‑69 (MQ=255)                             gggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTT‑TCTTTGACATTAATCCCGGAGAGt                                         >  1:1452639/1‑69 (MQ=255)                              ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             <  1:4899522/65‑1 (MQ=255)                              ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:5833893/71‑1 (MQ=255)                              ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:3343685/71‑1 (MQ=255)                              ggTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTa                                        >  1:2436857/1‑70 (MQ=255)                               gTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaaa                                      <  1:1821578/71‑1 (MQ=255)                                tAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAAc                                     >  1:3602359/1‑71 (MQ=255)                                 aaaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTccc                                                                    >  1:7649139/1‑39 (MQ=255)                                 aaaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             <  1:4636894/62‑1 (MQ=255)                                  aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                          >  1:1416368/1‑48 (MQ=255)                                  aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGAc                                                         <  1:532508/49‑1 (MQ=255)                                  aaTCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGAc                                                         <  1:523108/49‑1 (MQ=255)                                     ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTg                                                           >  1:1695216/1‑44 (MQ=255)                                     ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTaa                                                    >  1:204824/1‑51 (MQ=255)                                     ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:3307857/1‑70 (MQ=255)                                     ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:897626/1‑70 (MQ=255)                                     ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:3053792/1‑70 (MQ=255)                                     ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:1993208/1‑70 (MQ=255)                                     ccTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCggg                                 >  1:237401/1‑70 (MQ=255)                                      cTTAGCGTTAGAAATCGAAAGGCAGTGGAGTCCCTTGTCTTTGACa                                                        >  1:7258646/1‑46 (MQ=255)                                      cTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACa                                                        >  1:4291081/1‑46 (MQ=255)                                      cTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTa                               <  1:1676505/71‑1 (MQ=255)                                        tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                <  1:5600713/36‑1 (MQ=255)                                        tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGt                                                                <  1:3673409/36‑1 (MQ=255)                                        tAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTa                                        >  1:5032656/1‑60 (MQ=255)                                         aGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGga                                             >  1:5177872/1‑54 (MQ=255)                                          gCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGa                                                          >  1:678515/1‑40 (MQ=255)                                            gTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTa                         <  1:7501580/71‑1 (MQ=255)                                             ttAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTa                         >  1:6033077/1‑70 (MQ=255)                                                 aaaTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGc                     >  1:1378508/1‑70 (MQ=255)                                                 aaaTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCg                    <  1:3119412/71‑1 (MQ=255)                                                   aTCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATgg                           <  1:3044962/62‑1 (MQ=255)                                                      gAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCTATTgg               <  1:1684330/71‑1 (MQ=255)                                                       aaaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGt              <  1:5635083/71‑1 (MQ=255)                                                       aaaGACAAGGGACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTg                >  1:7653797/1‑69 (MQ=25)                                                       aaaGACAAGGGACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTg                >  1:4534769/1‑69 (MQ=25)                                                        aaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTTAAGCGATTGGtt             >  1:6626766/1‑71 (MQ=255)                                                        aaGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGtt             <  1:3489967/71‑1 (MQ=255)                                                         aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:6244903/44‑1 (MQ=255)                                                         aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTaa                                       <  1:1576187/44‑1 (MQ=255)                                                         aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGttt            <  1:2736806/71‑1 (MQ=255)                                                         aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGttt            >  1:4424518/1‑71 (MQ=255)                                                         aGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGt              <  1:3372853/69‑1 (MQ=255)                                                           gCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTg          <  1:3276904/71‑1 (MQ=255)                                                            cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGt                          >  1:7321514/1‑54 (MQ=255)                                                            cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt         <  1:5553728/71‑1 (MQ=255)                                                            cAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt         <  1:4922697/71‑1 (MQ=255)                                                             aGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGt         <  1:6787180/70‑1 (MQ=255)                                                               gggAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg      <  1:285673/71‑1 (MQ=255)                                                               gggAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCgg      >  1:4511659/1‑71 (MQ=255)                                                               gggACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTaa                        <  1:6114528/53‑1 (MQ=38)                                                               gggACTCCCCTGCCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTaa                        <  1:1824813/53‑1 (MQ=38)                                                                   gTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTAtt  <  1:890886/71‑1 (MQ=255)                                                                   |                                                                       CCTTTGTGATGCTTGTCCGGGATAACCCGGGTAAATCCTTAGCGTTAGAAATCGAAAGGCAGGGGAGTCCCTTGTCTTTGACATTAATCCCGGAGAGTAAACCGGGTAATGGTAAAGCGATTGGTTTTGTCGGTATT  >  W3110S.gb/197321‑197457

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