Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,502,074 T→C 28.8% E258G (GAA→GGA)  katG ← catalase/hydroperoxidase HPI(I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,502,0740TC28.8% 82.8 / 46.3 59E258G (GAA→GGA) katGcatalase/hydroperoxidase HPI(I)
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base C (7/10);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

CATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTCGCCGCT  >  W3110S.gb/3502006‑3502144
                                                                    |                                                                      
cATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttct                                                                      <  1:1035058/71‑1 (MQ=255)
cATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttct                                                                      >  1:4245014/1‑71 (MQ=255)
  tttGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttcttc                                                                    >  1:4149348/1‑71 (MQ=255)
     gATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgt                                                                  <  1:3117681/70‑1 (MQ=255)
     gATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgt                                                                  <  1:3443471/70‑1 (MQ=255)
     gATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgt                                                                  <  1:3897427/70‑1 (MQ=255)
     gATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTtcgt                                                                  <  1:396550/70‑1 (MQ=255)
         tCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGtt                                                              <  1:557118/70‑1 (MQ=255)
         tCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTc                                                             <  1:2200269/71‑1 (MQ=255)
         tCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGtt                                                              <  1:3552703/70‑1 (MQ=255)
           ggACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCa                                                            >  1:2604683/1‑70 (MQ=255)
           ggACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTTGTCGTTCAt                                                           >  1:2275828/1‑71 (MQ=255)
             aCCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGc                                                         >  1:2651880/1‑71 (MQ=255)
              ccGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTttctt                                                                     >  1:887559/1‑58 (MQ=255)
               cGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGccc                                                       <  1:280246/71‑1 (MQ=255)
                        gggTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGcc                                               >  1:1855380/1‑70 (MQ=255)
                         ggTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCg                                              <  1:895705/70‑1 (MQ=255)
                         ggTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCg                                              <  1:1025109/70‑1 (MQ=255)
                            tttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAgg                                          >  1:524449/1‑71 (MQ=255)
                            tttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAgg                                          >  1:807900/1‑71 (MQ=255)
                             ttACCCAGCGTATGACCACCGGCAATCAGCGCCACGGTTTCTTCGTCGTTCa                                                            >  1:186017/1‑52 (MQ=255)
                             ttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCa                                                            >  1:4848375/1‑52 (MQ=255)
                             ttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAgct                                         >  1:992846/1‑69 (MQ=255)
                             ttACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAgg                                          <  1:336715/70‑1 (MQ=255)
                                cccAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcgc                                      <  1:1780338/71‑1 (MQ=255)
                                   aGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCa                                                            >  1:4919149/1‑46 (MQ=255)
                                   aGCGTATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcgcgcg                                   <  1:2419573/71‑1 (MQ=255)
                                    gCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGccc                                                       <  1:2256298/50‑1 (MQ=255)
                                     cGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa                                 >  1:1662719/1‑71 (MQ=255)
                                      gTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa                                 >  1:2520136/1‑70 (MQ=255)
                                      gTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa                                 >  1:1392579/1‑70 (MQ=255)
                                       tATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGtt                                                              <  1:626983/40‑1 (MQ=255)
                                       tATGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGtt                                                              <  1:2583735/40‑1 (MQ=255)
                                         tGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAgc                             >  1:4721288/1‑71 (MQ=255)
                                         tGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAgc                             >  1:3170768/1‑71 (MQ=255)
                                         tGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAg                              <  1:1204028/70‑1 (MQ=255)
                                         tGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAg                              <  1:3284656/70‑1 (MQ=255)
                                         tGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAg                              <  1:3526618/70‑1 (MQ=255)
                                         tGACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGTc                                        >  1:965603/1‑60 (MQ=255)
                                          gACCACCCGCAATCAGCGCCACGGTTCCTTCGTCttt                                                              <  1:201982/37‑1 (MQ=37)
                                          gACCACCCGCAATCAGCGCCACGGTTCCTTCGTCGtt                                                              <  1:3631382/37‑1 (MQ=255)
                                            ccaccCGCAATCAGCGCCACGGTTCCTTCGTCGTTCa                                                            >  1:4427999/1‑37 (MQ=255)
                                            ccaccCGCAATCAGCGCCACGGTTCCTTCGTCGTTCa                                                            >  1:2220315/1‑37 (MQ=255)
                                            ccaccCGCAATCAGCGCCACGGTTCCTTCGTCGTTCa                                                            >  1:2015283/1‑37 (MQ=255)
                                                ccGCAATCAGCGCCACGGTTCCTTCGTCGTTCATGCCCa                                                      <  1:4459707/39‑1 (MQ=255)
                                                      tCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCg                                              <  1:2307776/41‑1 (MQ=255)
                                                      tCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCg                                              <  1:4837301/41‑1 (MQ=255)
                                                      tCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCg                                              <  1:1442970/41‑1 (MQ=255)
                                                      tCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGa                                             <  1:2461428/42‑1 (MQ=255)
                                                      tCAGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGaa                <  1:5077409/71‑1 (MQ=255)
                                                        aGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGaa                                            >  1:4098791/1‑41 (MQ=255)
                                                        aGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAg                                           >  1:3377251/1‑42 (MQ=255)
                                                        aGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAg                                           >  1:312253/1‑42 (MQ=255)
                                                        aGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcgcg                                     >  1:2929432/1‑48 (MQ=255)
                                                        aGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTcgcg                                     >  1:237870/1‑48 (MQ=255)
                                                        aGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAg              >  1:2816007/1‑71 (MQ=255)
                                                        aGCGCCACGGTTCCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGaaa               >  1:2514573/1‑70 (MQ=255)
                                                             cACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGtt         >  1:4645824/1‑71 (MQ=255)
                                                             cACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGtt         >  1:3060699/1‑71 (MQ=255)
                                                              aCGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAgct                            <  1:1443291/51‑1 (MQ=255)
                                                              aCGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAgct                            <  1:2563603/51‑1 (MQ=255)
                                                              aCGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAgct                            <  1:1082161/51‑1 (MQ=255)
                                                              aCGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTc        >  1:410715/1‑71 (MQ=255)
                                                                  tttCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa                                 <  1:3959507/42‑1 (MQ=255)
                                                                  tttCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGa                                 <  1:2207442/42‑1 (MQ=255)
                                                                  tttCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCgg           <  1:4526560/64‑1 (MQ=255)
                                                                  tttCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCgg           <  1:1813559/64‑1 (MQ=255)
                                                                  tttCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCgg           <  1:513396/64‑1 (MQ=255)
                                                                  tttCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTcgccg    >  1:2987917/1‑71 (MQ=255)
                                                                   ttcttcGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAg                              >  1:5106252/1‑44 (MQ=255)
                                                                   ttcttcGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAg                              >  1:933993/1‑44 (MQ=255)
                                                                    tcttcGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTCGCCGCt  >  1:1966203/1‑71 (MQ=255)
                                                                    |                                                                      
CATTTGATGTCGGACCGGCACCGTGGGTTTTACCCAGCGTATGACCACCCGCAATCAGCGCCACGGTTTCTTCGTCGTTCATGCCCATGTTGCCGAAGGTCGCGCGGATAGCTGCTGCCGCAGAAAGCGGTTCGCCGCT  >  W3110S.gb/3502006‑3502144

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