Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 496,170 A→G 45.9% E609E (GAA→GAG htpG → molecular chaperone HSP90 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb496,1700AG45.9% 3.5 / 49.9 37E609E (GAA→GAGhtpGmolecular chaperone HSP90 family
Reads supporting (aligned to +/- strand):  ref base A (8/12);  new base G (10/7);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

TCAGCGAGTGGGTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTAATGCCGGAT  >  W3110S.gb/496109‑496231
                                                             |                                                             
tCAGCGAGTGGGTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGa                                                        <  1:1319403/69‑1 (MQ=255)
    cGAGTGGGTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGAACCt                                                    <  1:436620/69‑1 (MQ=255)
      aGTGGGTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGAACCTGt                                                  <  1:3446632/69‑1 (MQ=255)
      aGTGGGTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGAACCTGt                                                  <  1:739962/69‑1 (MQ=255)
         gggTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTa                                               <  1:1746918/69‑1 (MQ=255)
          ggTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTAt                                              >  1:1129735/1‑69 (MQ=255)
          ggTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTAt                                              >  1:795753/1‑69 (MQ=255)
            tAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGAACCTg                                                   >  1:2651062/1‑62 (MQ=255)
             aGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATtcg                                           >  1:3610063/1‑69 (MQ=255)
             aGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATtcg                                           >  1:1364814/1‑69 (MQ=255)
                      gctgGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGa                                   >  1:899909/1‑68 (MQ=255)
                      gctgGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGa                                   >  1:2126015/1‑68 (MQ=255)
                         ggATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGAAcc                                <  1:1055616/68‑1 (MQ=255)
                         ggATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCa                               >  1:3222319/1‑69 (MQ=255)
                             cAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGAACCAgctg                           >  1:847545/1‑69 (MQ=255)
                             cAGGCGCTGCTGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGAACCAgctg                           >  1:878092/1‑69 (MQ=255)
                                  gctgctGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATg                                    <  1:2351629/55‑1 (MQ=255)
                                  gctgctGGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATg                                    <  1:1090640/55‑1 (MQ=255)
                                   ctgctgGCAGAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATg                                    <  1:652569/54‑1 (MQ=255)
                                        ggCAGAACGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTg                <  1:3443583/69‑1 (MQ=255)
                                           aGAACGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGttt                    >  1:1934054/1‑62 (MQ=255)
                                           aGAACGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGttt                    >  1:3436925/1‑62 (MQ=255)
                                            gAACGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGt            <  1:2618250/69‑1 (MQ=255)
                                            gAACGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGt            <  1:926043/69‑1 (MQ=255)
                                            gAACGCGGCACCCTGGAGGATCCGAAGCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTg                <  1:3544946/65‑1 (MQ=255)
                                             aaCGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGt            <  1:1791419/68‑1 (MQ=255)
                                              aCGCGGCACGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTaa          >  1:2196488/1‑69 (MQ=255)
                                              aCGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTaa          >  1:2431733/1‑69 (MQ=255)
                                               cgcgGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATg                                    >  1:860704/1‑42 (MQ=255)
                                               cgcgGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTAAt         >  1:3307261/1‑69 (MQ=255)
                                               cgcgGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTAAt         <  1:445213/69‑1 (MQ=255)
                                               cgcgGCACGCTGGAAGATCCGAACCTGTTTAATCGTCGTATg                                    <  1:3008469/42‑1 (MQ=255)
                                                    cACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTa           >  1:1015730/1‑62 (MQ=255)
                                                     aCGCTGGAGGATCCGAACCTGTTTATTCGTCGTATGa                                   >  1:3089319/1‑37 (MQ=255)
                                                      cGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTAATGCCGGAt  <  1:245283/69‑1 (MQ=255)
                                                      cGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTAATGCCGGAt  <  1:1668190/69‑1 (MQ=255)
                                                      cGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTAATGCCGGAt  <  1:1140377/69‑1 (MQ=255)
                                                             |                                                             
TCAGCGAGTGGGTAGAACTGCTGCTGGATCAGGCGCTGCTGGCAGAACGCGGCACGCTGGAAGATCCGAACCTGTTTATTCGTCGTATGAACCAGCTGCTGGTTTCCTGATGTAATGCCGGAT  >  W3110S.gb/496109‑496231

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: