Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,125,971 A→G 50.0% V236A (GTT→GCT)  glcD ← glycolate oxidase subunit, FAD‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,125,9710AG50.0% ‑4.0 / 34.7 26V236A (GTT→GCT) glcDglycolate oxidase subunit, FAD‑linked
Reads supporting (aligned to +/- strand):  ref base A (4/9);  new base G (3/10);  total (7/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

ATTGGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCACGCC  >  W3110S.gb/3125907‑3126026
                                                                |                                                       
aTTGGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCcg                                                     <  1:3372284/69‑1 (MQ=255)
aTTGGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCcg                                                     <  1:3199738/69‑1 (MQ=255)
   ggCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAGCCCgcgc                                                  <  1:1237638/69‑1 (MQ=255)
   ggCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAGCCCgcgc                                                  >  1:302936/1‑69 (MQ=255)
   ggCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCgcgc                                                  >  1:1734925/1‑69 (MQ=255)
   ggCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCgcgc                                                  >  1:2205274/1‑69 (MQ=255)
   ggCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCgcgc                                                  >  1:3127879/1‑69 (MQ=255)
        tgatgTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACgg                                             <  1:924249/69‑1 (MQ=255)
            gTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGggcgg                                         >  1:1551811/1‑69 (MQ=255)
             tCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGggcggc                                        <  1:1435851/69‑1 (MQ=255)
             tCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGggcggc                                        <  1:3519742/69‑1 (MQ=255)
                        aaGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGcagcag                             <  1:95356/69‑1 (MQ=255)
                        aaGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGcagcag                             <  1:707328/69‑1 (MQ=255)
                        aaGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGcagcag                             <  1:541962/69‑1 (MQ=255)
                           tCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGttt                          <  1:2949212/69‑1 (MQ=255)
                           tCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGttt                          >  1:1203981/1‑69 (MQ=255)
                            ccGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGtttt                         <  1:3393518/69‑1 (MQ=255)
                               gCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGtttt                         >  1:3330023/1‑66 (MQ=255)
                                 tttttCTACCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGt                    <  1:291780/69‑1 (MQ=255)
                                        aCCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTc              <  1:1612893/68‑1 (MQ=255)
                                        aCCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTc              <  1:1292580/68‑1 (MQ=255)
                                        aCCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTc              <  1:1288383/68‑1 (MQ=255)
                                        aCCGAGTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTc              <  1:741478/68‑1 (MQ=255)
                                             gTCAAAGCTGGCTAACAGAGCCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCggtggt        <  1:483299/69‑1 (MQ=255)
                                                aaaGCTGGCTAACAGAGCCCGCGCCACGGGCGGATTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCa      <  1:1134117/68‑1 (MQ=255)
                                                   gCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCACGcc  <  1:889069/69‑1 (MQ=255)
                                                                |                                                       
ATTGGCGATGATGTCACCAACCGCAAGTCCGGCTTTTTCTACCGAGTCAAAGCTGGCTAACAGAACCCGCGCCACGGGCGGCTTCGGCAGCAGTTTTACCGTCACTTCGGTGGTCACGCC  >  W3110S.gb/3125907‑3126026

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: