Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 15,286 C→T 61.5% D373D (GAC→GAT dnaJ → chaperone Hsp40, co‑chaperone with DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb15,2860CT61.5% 14.4 / 26.6 26D373D (GAC→GATdnaJchaperone Hsp40, co‑chaperone with DnaK
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base T (6/10);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

CCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTG  >  W3110S.gb/15222‑15345
                                                                |                                                           
ccGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTgt                                                         <  1:500701/69‑2 (MQ=255)
ccGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGa                                                         <  1:863896/69‑1 (MQ=255)
ccGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGa                                                         <  1:243956/69‑1 (MQ=255)
ccGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGa                                                         >  1:2729180/1‑69 (MQ=255)
      aGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                   <  1:255543/69‑1 (MQ=255)
       gCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTa                                                  <  1:1128336/69‑1 (MQ=255)
       gCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTa                                                  <  1:3091172/69‑1 (MQ=255)
         acaacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTAAc                                                <  1:2117863/69‑1 (MQ=255)
             caGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTa                                                  >  1:2043749/1‑63 (MQ=255)
             caGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTa                                                  >  1:2954278/1‑63 (MQ=255)
                   gcgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTAACCTCCCCaaaa                                      <  1:1078504/69‑1 (MQ=255)
                   gcgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTAACCTCCCCaaaa                                      <  1:1074059/69‑1 (MQ=255)
                    cgcTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTAACCTCCCCaaaa                                      <  1:1295203/68‑1 (MQ=255)
                         aaaGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGc                                >  1:3678471/1‑69 (MQ=255)
                                  cttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCa                       >  1:271551/1‑69 (MQ=255)
                                  cttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCa                       >  1:2759324/1‑69 (MQ=255)
                                          gtgtGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTg                 >  1:1979791/1‑67 (MQ=255)
                                               aagaagTTTTTTGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTaaa           <  1:726632/68‑1 (MQ=255)
                                                agaagTTTTTTGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTaaa           <  1:701753/67‑1 (MQ=255)
                                                   agTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAg      <  1:2035974/69‑1 (MQ=255)
                                                    gTTTTTTGACGATCTGACCCGCTAACCTCCCCGAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAgg     >  1:2068755/1‑69 (MQ=255)
                                                    gTTTTTTGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAgg     >  1:359754/1‑69 (MQ=255)
                                                       ttttGACGATCTGACCCGCTAACCTCCCCAAAAGCCTg                                 >  1:1210033/1‑38 (MQ=255)
                                                       ttttGACGATCTGACCCGCTAACCTCCCCAAAAGCCTg                                 >  1:705538/1‑38 (MQ=255)
                                                       ttttGACGATCTGACCCGCTAACCTCCCCAAAAGCCTg                                 >  1:2109636/1‑38 (MQ=255)
                                                       ttttGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTg  >  1:1998841/1‑69 (MQ=255)
                                                                |                                                           
CCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTG  >  W3110S.gb/15222‑15345

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: