Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 14,173 T→C 19.6% A2A (GCT→GCC dnaJ → chaperone Hsp40, co‑chaperone with DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb14,1730TC19.6% 117.3 / 25.8 56A2A (GCT→GCCdnaJchaperone Hsp40, co‑chaperone with DnaK
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base C (9/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 9.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGC  >  W3110S.gb/14112‑14238
                                                             |                                                                 
ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                                          <  1:3448436/71‑1 (MQ=255)
ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCCAAGCAAGat                                                          >  1:794248/1‑71 (MQ=255)
  cGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatta                                                        <  1:1102303/71‑1 (MQ=255)
     aGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAttatta                                                     >  1:1681196/1‑71 (MQ=255)
         gAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACgag                                                 <  1:621422/71‑1 (MQ=255)
           aTTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCCAAGCAAGATTATTACGAGAt                                               >  1:1535975/1‑71 (MQ=255)
            tttCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatt                                                         <  1:1977082/60‑1 (MQ=255)
              tCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAttt                                             >  1:3776451/1‑70 (MQ=255)
              tCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAttt                                             <  1:2345396/70‑1 (MQ=255)
                   tcCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCCAAGCAAGATTATTACGAGATTTTAGGCg                                       >  1:1546995/1‑71 (MQ=255)
                          gTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAttt                                             <  1:3964040/58‑1 (MQ=255)
                          gTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAttt                                             <  1:2731251/58‑1 (MQ=255)
                          gTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAttt                                             <  1:2198692/58‑1 (MQ=255)
                             cATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                                          <  1:2562901/42‑1 (MQ=255)
                             cATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa                              >  1:182097/1‑70 (MQ=255)
                             cATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa                              >  1:3895327/1‑70 (MQ=255)
                               ttCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAg                           >  1:1372735/1‑71 (MQ=255)
                                tCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                                          >  1:2386719/1‑39 (MQ=255)
                                tCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                                          >  1:3875258/1‑39 (MQ=255)
                                 cATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                           >  1:1260598/1‑53 (MQ=255)
                                 cATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                           >  1:2043603/1‑53 (MQ=255)
                                 cATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                           >  1:2345826/1‑53 (MQ=255)
                                 cATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                           >  1:899531/1‑53 (MQ=255)
                                 cATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                         >  1:685020/1‑55 (MQ=255)
                                 cATCTAGGGGCAATTTAAAAAAGATGGCCAAGCAAGatta                                                        >  1:303842/1‑40 (MQ=255)
                                  aTCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCgg                        <  1:3772927/71‑1 (MQ=255)
                                   tCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCgg                        <  1:2230792/70‑1 (MQ=255)
                                     tAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACgaga                                                >  1:2378490/1‑44 (MQ=255)
                                     tAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACgaga                                                >  1:537717/1‑44 (MQ=255)
                                     tAGGGGCAATTTAAAAAAGATGGCCAAGCAAGattat                                                       >  1:826376/1‑37 (MQ=255)
                                           cAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTg               <  1:2532025/71‑1 (MQ=255)
                                              tttAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaa                                <  1:2266422/51‑1 (MQ=255)
                                              tttAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaa                                <  1:2241061/51‑1 (MQ=255)
                                              tttAAAAAAGATGGCCAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCgg                        <  1:2152643/59‑1 (MQ=255)
                                               ttAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           >  1:315393/1‑71 (MQ=255)
                                               ttAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           >  1:348139/1‑71 (MQ=255)
                                               ttAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           >  1:2234669/1‑71 (MQ=255)
                                               ttAAAAAAGATGGCCAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa                                 >  1:2781868/1‑49 (MQ=255)
                                                tAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                           >  1:3521206/1‑38 (MQ=255)
                                                tAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                           >  1:2362470/1‑38 (MQ=255)
                                                tAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa                              >  1:2615203/1‑51 (MQ=255)
                                                tAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           >  1:1398469/1‑70 (MQ=255)
                                                tAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           <  1:3783522/70‑1 (MQ=255)
                                                tAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           >  1:1536474/1‑70 (MQ=255)
                                                tAAAAAAGATGGCCAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           >  1:3728996/1‑70 (MQ=255)
                                                tAAAAAAGATGGCCAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt           >  1:3858625/1‑70 (MQ=255)
                                                   aaaaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa       <  1:3198778/71‑1 (MQ=255)
                                                   aaaaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa       <  1:4033047/71‑1 (MQ=255)
                                                   aaaaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa       <  1:195660/71‑1 (MQ=255)
                                                     aaGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAg                           <  1:1241630/49‑1 (MQ=255)
                                                     aaGATGGCCAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa       >  1:811724/1‑69 (MQ=255)
                                                      aGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAg    <  1:3601413/71‑1 (MQ=255)
                                                        aTGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa       <  1:929106/66‑1 (MQ=255)
                                                        aTGGCCAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGc  <  1:1434114/71‑1 (MQ=255)
                                                         tGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGa       <  1:1743550/65‑1 (MQ=255)
                                                          ggCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaaa     >  1:551222/1‑66 (MQ=255)
                                                             |                                                                 
GGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGC  >  W3110S.gb/14112‑14238

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: