Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 15,286 C→T 100% 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,2860CT82.6% 49.3 / 5.8 23D373D (GAC→GATdnaJchaperone Hsp40, co‑chaperone with DnaK
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (9/10);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGAT  >  W3110S.gb/15219‑15355
                                                                   |                                                                     
cAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTgacgac                                                                       >  1:1261027/1‑68 (MQ=255)
cAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTgacgac                                                                       >  1:139438/1‑68 (MQ=255)
cAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCt                                                                     >  1:272177/1‑70 (MQ=255)
cAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCt                                                                     >  1:1499077/1‑70 (MQ=255)
cAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCt                                                                     >  1:2252194/1‑70 (MQ=255)
    cGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGAccc                                                                >  1:832039/1‑71 (MQ=255)
    cGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGAccc                                                                >  1:1568400/1‑71 (MQ=255)
         aGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGa                                                                   <  1:1726249/63‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                             <  1:1184662/46‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                             <  1:822395/46‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                             <  1:3505091/46‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                             <  1:344157/46‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                             <  1:3092769/46‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                             <  1:3048064/46‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCt                                                             <  1:2202629/46‑1 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTgg                                    >  1:2104770/1‑71 (MQ=255)
                                agCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGATCTACCCGCt                                                             <  1:2095734/45‑1 (MQ=38)
                                  cttcttTGATGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTaa                                                           >  1:1251809/1‑46 (MQ=255)
                                     cttTGATGGTGTGAAGAAGTTTTTTGACGACCTGa                                                                   >  1:2864016/1‑35 (MQ=255)
                                     cttTGATGGTGTGAAGAAGTTTTTTGACGACCTGa                                                                   >  1:3085021/1‑35 (MQ=255)
                                           tGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTAAc                                                          >  1:257748/1‑38 (MQ=255)
                                           tGGTGTGAAGAAGTTTTTTGACGATCTGACCCGCTAAc                                                          >  1:2081577/1‑38 (MQ=255)
                                                       gTTTTTTGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGt             >  1:1517157/1‑71 (MQ=255)
                                                          ttttGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGc           <  1:3296909/70‑1 (MQ=255)
                                                          ttttGACGATCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCg          >  1:1740048/1‑71 (MQ=255)
                                                                 gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGa   <  1:318608/71‑1 (MQ=255)
                                                                  acCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGat  <  1:2257816/71‑1 (MQ=255)
                                                                   |                                                                     
CAACCGGCGAGCACAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGAT  >  W3110S.gb/15219‑15355

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: