Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,934,313 G→A 19.0% A53V (GCA→GTA)  eda ← multifunctional 2‑keto‑3‑deoxygluconate 6‑phosphate aldolase, 2‑keto‑4‑hydroxyglutarate aldolase, and oxaloacetate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,934,3130GA19.0% 44.2 / 5.2 21A53V (GCA→GTA) edamultifunctional 2‑keto‑3‑deoxygluconate 6‑phosphate aldolase, 2‑keto‑4‑hydroxyglutarate aldolase, and oxaloacetate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (12/5);  new base A (2/2);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TCAGCACCGTACCGGCACCCACAATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCACCCCACCAGCAACCAACGCTTTTGCCATCGG  >  W3110S.gb/1934245‑1934380
                                                                    |                                                                   
tCAGCACCGTACCGGCACCCACAATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTg                                                                     >  1:2965862/1‑69 (MQ=255)
tCAGCACCGTACCGGCACCCACAATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTg                                                                     >  1:2585457/1‑69 (MQ=255)
      ccatACCGGCACCCACAATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTGCACACt                                                               >  1:544488/4‑69 (MQ=255)
      ccatACCGGCACCCACAATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTGCACACt                                                               >  1:548736/4‑69 (MQ=255)
                    acaATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTACACACTCGGTACGCAGAGTc                                                 >  1:2023977/1‑69 (MQ=255)
                    acaATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTACACACTCGGTACGCAGAGTc                                                 >  1:2945220/1‑69 (MQ=255)
                       aTCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACt                                              <  1:3143028/69‑1 (MQ=255)
                       aTCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACt                                              <  1:2187108/69‑1 (MQ=255)
                              cAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGaa                                       >  1:2412813/1‑69 (MQ=255)
                                   aCTTCTTTGGCGATAGCACGGATAGCGTCAACTACACACTCGGTACGCAGAGTCACTTCCAGAACgcgc                                  <  1:1378088/69‑1 (MQ=255)
                                      tcttTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAAcgcg                                   <  1:1034330/65‑1 (MQ=255)
                                       cttTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAg                                         >  1:471357/1‑58 (MQ=255)
                                       cttTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAg                                         >  1:832065/1‑58 (MQ=255)
                                       cttTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAg                                         <  1:2327105/58‑1 (MQ=255)
                                       cttTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAg                                         <  1:1330771/58‑1 (MQ=255)
                                       cttTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCAccc                              >  1:2818514/1‑69 (MQ=255)
                                           ggCGATAGCACGGATAGCGTCAACTACACACTCGGTACGCAGAGTCACTTCCAGAACgcgc                                  <  1:2729379/61‑1 (MQ=255)
                                                 aGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCACCCCACCAGc                       >  1:2572504/1‑66 (MQ=255)
                                                     cGGATAGCGTCAACTGCACACTCGGTACGCAGAGt                                                  >  1:1222546/1‑35 (MQ=255)
                                                     cGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCACCCCACCAGcaaccaac                >  1:2589541/1‑69 (MQ=255)
                                                     cGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCACCCCACCAGcaaccaac                >  1:2554941/1‑69 (MQ=255)
                                                           gCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCACCCCACCAGCAACCAACGCtttt          >  1:2653720/1‑69 (MQ=255)
                                                                   tGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCACCCCACCAGCAACCAACGCTTTTGCCATCgg  >  1:3721017/1‑69 (MQ=255)
                                                                    |                                                                   
TCAGCACCGTACCGGCACCCACAATCGCTTCAGGCACTTCTTTGGCGATAGCACGGATAGCGTCAACTGCACACTCGGTACGCAGAGTCACTTCCAGAACGCGCACCCCACCAGCAACCAACGCTTTTGCCATCGG  >  W3110S.gb/1934245‑1934380

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: