Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,244,017 G→A 46.9% Q86Q (CAG→CAA exuT → hexuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,244,0170GA46.9% 2.6 / 42.2 32Q86Q (CAG→CAAexuThexuronate transporter
Reads supporting (aligned to +/- strand):  ref base G (12/5);  new base A (8/7);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTG  >  W3110S.gb/3243950‑3244086
                                                                   |                                                                     
tGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAt                                                                    >  1:394978/1‑71 (MQ=255)
    gcgTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCt                                                                <  1:1724248/71‑1 (MQ=255)
       tAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAg                                                                       <  1:1397450/61‑1 (MQ=255)
        aaCACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTatat                                                            >  1:1499628/1‑71 (MQ=255)
              gtgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTc                                                                  >  1:436179/1‑59 (MQ=255)
              gtgGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTc                                                                  >  1:1585095/1‑59 (MQ=255)
                 gcggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAg                                                    <  1:1695266/70‑1 (MQ=255)
                  cggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGcc                                                  >  1:1247437/1‑71 (MQ=255)
                  cggcAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGcc                                                  >  1:1268325/1‑71 (MQ=255)
                        cTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTct                                            >  1:125195/1‑71 (MQ=255)
                         tGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTAt                                                                    <  1:24801/46‑1 (MQ=255)
                          gcgcCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATatcatc                                                        >  1:1057756/1‑57 (MQ=255)
                               aaCTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCtt                                     <  1:299207/71‑1 (MQ=255)
                                    tGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACg                                >  1:1382984/1‑71 (MQ=255)
                                    tGATGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACg                                >  1:459229/1‑71 (MQ=255)
                                       tGGAAGAGTTAAACATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTc                             <  1:513943/71‑1 (MQ=255)
                                                 aaaCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTct                                            >  1:787290/1‑46 (MQ=255)
                                                 aaaCATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTtata                                  >  1:1133496/1‑56 (MQ=255)
                                                 aaaCATCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:1126477/1‑71 (MQ=255)
                                                  aaCATCTCCACCCAACAGTATTCCTATATCATCGCa                                                     >  1:1576487/1‑36 (MQ=255)
                                                  aaCATCTCCACCCAACAGTATTCCTATATCATCGCa                                                     >  1:395049/1‑36 (MQ=255)
                                                   aCATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCgg                   >  1:423131/1‑69 (MQ=255)
                                                    cATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAg                <  1:69305/71‑1 (MQ=255)
                                                     aTCTCCACCCAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGc               <  1:1231689/71‑1 (MQ=255)
                                                      tctcCACCCAACAATATTCCTATATCATCGCAGCCTa                                                <  1:506732/37‑1 (MQ=255)
                                                      tctcCACCCAACAATATTCCTATATCATCGCAGCCTa                                                <  1:855836/37‑1 (MQ=255)
                                                            cccAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTAt        >  1:1313283/1‑71 (MQ=255)
                                                             ccAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATg       <  1:1450592/71‑1 (MQ=255)
                                                             ccAACAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATg       <  1:698483/71‑1 (MQ=255)
                                                              caacaaTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATg       >  1:1292393/1‑70 (MQ=255)
                                                              cAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGt      <  1:1763400/71‑1 (MQ=255)
                                                               aaCAATATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTa     >  1:1673777/1‑71 (MQ=255)
                                                                  aGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCaa                       >  1:1073606/1‑50 (MQ=255)
                                                                  aGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCTGTAGCAGGTTATGTACTg  >  1:102010/1‑71 (MQ=255)
                                                                   |                                                                     
TGACGCGTAACACTGTGGCGGCAGCTGCGCCAACTCTGATGGAAGAGTTAAACATCTCCACCCAACAGTATTCCTATATCATCGCAGCCTATTCTGCTGCTTATACGGTCATGCAACCGGTAGCAGGTTATGTACTG  >  W3110S.gb/3243950‑3244086

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: