Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,110,393 A→G 8.2% N84N (AAT→AAC hslO ← heat shock protein Hsp33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,110,3930AG 8.2% 159.7 / 5.3 61N84N (AAT→AAChslOheat shock protein Hsp33
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (3/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CCAGCGTTTTCAGGTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGATATCACCATC  >  W3110S.gb/4110323‑4110461
                                                                      |                                                                    
ccAGCGTTTTCAGGTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTa                                                                      <  1:291441/71‑1 (MQ=255)
   gCGTTTTCAGGTCGGCATTTTCTGGAATTTCGCCCTGCAGGCGCGCCACACCGCGCATCTGCTGGTTATTg                                                                   >  1:3183848/1‑71 (MQ=255)
    cGTTTTCAGGTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGt                                                                  <  1:2603007/71‑1 (MQ=255)
    cGTTTTCAGGTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGt                                                                  <  1:3423771/71‑1 (MQ=255)
         tCAGGTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCg                                                             >  1:3742450/1‑71 (MQ=255)
          cAGGTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGt                                                            >  1:971602/1‑71 (MQ=255)
            ggTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTta                                                          <  1:3549897/71‑1 (MQ=255)
              tCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTtaat                                                        <  1:2103188/71‑1 (MQ=255)
                ggCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGt                                                                  <  1:420919/59‑1 (MQ=255)
                ggCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGt                                                                  <  1:673281/59‑1 (MQ=255)
                 gCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAAc                                                     <  1:3172567/71‑1 (MQ=255)
                    ttttCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGt                                                                  <  1:1908721/55‑1 (MQ=255)
                       tCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTGTTGTTACCGTTAATAACCGCCa                                                <  1:2465860/70‑1 (MQ=255)
                       tCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTGTTGTTACCGTTAATAACCGCCa                                                <  1:2223475/70‑1 (MQ=255)
                        cTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTAc                                                               >  1:1234037/1‑54 (MQ=255)
                         tGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTg                                                                   >  1:1674008/1‑49 (MQ=255)
                           gAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGAtt                                            <  1:1256376/70‑1 (MQ=255)
                           gAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGAtt                                            >  1:2598377/1‑70 (MQ=255)
                           gAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTc                                           <  1:194016/71‑1 (MQ=255)
                             aTTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCAt                                         <  1:1738090/71‑1 (MQ=255)
                              tttCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCa                                                <  1:392242/63‑1 (MQ=255)
                                tCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCg                                       >  1:2416808/1‑70 (MQ=255)
                                tCGCCCTGCACGCGCGCCACACCGCGCATCTGCTCGTTATTGTTACCGTTAATAACCGCCAGATTCATCg                                       >  1:491604/1‑70 (MQ=255)
                                 cGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCgg                                      <  1:3187614/70‑1 (MQ=255)
                                 cGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCgg                                      <  1:3593501/70‑1 (MQ=255)
                                      tGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCa                                                <  1:67685/55‑1 (MQ=255)
                                      tGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTc                                >  1:3756651/1‑71 (MQ=255)
                                       gCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTc                                >  1:3172146/1‑70 (MQ=255)
                                        cACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCgg                                      >  1:3354273/1‑63 (MQ=255)
                                        cACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGc                              >  1:3194326/1‑71 (MQ=255)
                                        cACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGc                              >  1:2341562/1‑71 (MQ=255)
                                            cgcgcCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTg                          <  1:1210474/71‑1 (MQ=255)
                                               gcCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAg                       >  1:3749010/1‑71 (MQ=255)
                                                 cacaCCGCGCATCTGCTGGTTATTGTTACCGTtaat                                                        <  1:2555523/36‑1 (MQ=255)
                                                  acacCGCGCATCTGCTGGTTATTGTTACCGTTAATAAcc                                                    >  1:3224451/1‑39 (MQ=255)
                                                  acacCGCGCATCTGCTGGTTATTGTTACCGTTAATAAcc                                                    >  1:1425700/1‑39 (MQ=255)
                                                  acacCGCGCATCTGCTGGTTATTGTTACCGTTAATAAcc                                                    >  1:1619994/1‑39 (MQ=255)
                                                     ccGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCt                           <  1:1507799/61‑1 (MQ=255)
                                                        cgcATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGt              <  1:2251674/71‑1 (MQ=255)
                                                        cgcATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGt              <  1:3348794/71‑1 (MQ=255)
                                                         gcATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTg             <  1:3207979/71‑1 (MQ=255)
                                                            tCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGAcc                                   <  1:678208/46‑1 (MQ=255)
                                                              tgctgGTTGTTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCTTGCAGCTGTACGGTGatat         >  1:524141/1‑70 (MQ=255)
                                                              tgctgGTTGTTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCTTGCAGCTGTACGGTGatat         >  1:3784802/1‑70 (MQ=255)
                                                              tgctgGTTGTTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCTTGCAGCTGTACGGTGat           >  1:83205/1‑68 (MQ=255)
                                                              tgctgGTTATTGTTACCGTTAATAACCGCCAGATTCAt                                         >  1:691275/1‑38 (MQ=255)
                                                              tgctgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGatat         >  1:2875273/1‑70 (MQ=255)
                                                              tgctgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGatat         >  1:1627629/1‑70 (MQ=255)
                                                              tgctgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGatat         >  1:1589868/1‑70 (MQ=255)
                                                              tgctgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGatat         >  1:594360/1‑70 (MQ=255)
                                                              tgctgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGatat         >  1:3281347/1‑70 (MQ=255)
                                                              tgctgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGatat         >  1:105678/1‑70 (MQ=255)
                                                                 tgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGt              <  1:541547/62‑1 (MQ=255)
                                                                 tgGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGt              <  1:58522/62‑1 (MQ=255)
                                                                   gTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTa                  <  1:3534170/56‑1 (MQ=255)
                                                                   gTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTa                  <  1:2199067/56‑1 (MQ=255)
                                                                    ttATTGTTACCGTTAATAACCGCCAGATTCATCgg                                      <  1:1709877/35‑1 (MQ=255)
                                                                    ttATTGTTACCGTTAATAACCGCCAGATTCATCgg                                      <  1:3517045/35‑1 (MQ=255)
                                                                    ttATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGt              <  1:1849868/59‑1 (MQ=255)
                                                                    ttATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGt              <  1:1239895/59‑1 (MQ=255)
                                                                     tATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGATATCACCATc  >  1:1840409/1‑70 (MQ=255)
                                                                     tATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGATATCACCATc  >  1:2755970/1‑70 (MQ=255)
                                                                      aTTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGatat         >  1:2922042/1‑62 (MQ=255)
                                                                      |                                                                    
CCAGCGTTTTCAGGTCGGCATTTTCTGGAATTTCGCCCTGCACGCGCGCCACACCGCGCATCTGCTGGTTATTGTTACCGTTAATAACCGCCAGATTCATCGGACCGTCGCCCTGCAGCTGTACGGTGATATCACCATC  >  W3110S.gb/4110323‑4110461

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