Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,599,304 G→A 100% M57I (ATG→ATA bcp → thiol peroxidase, thioredoxin‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,599,3040GA81.8% 71.0 / 11.6 33M57I (ATG→ATAbcpthiol peroxidase, thioredoxin‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/6);  new base A (13/14);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 5.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTC  >  W3110S.gb/2599239‑2599367
                                                                 |                                                               
tATTTCTACCCGAAAGCCATGACCCCCGGCTGTACGGTACAGGCCTGCGGCTTACGCGATAACATAGAt                                                              <  1:3605286/69‑1 (MQ=255)
  tttCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATg                                                             <  1:1373817/68‑1 (MQ=255)
  tttCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATg                                                             <  1:1602099/68‑1 (MQ=255)
  tttCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATg                                                             <  1:3154653/68‑1 (MQ=255)
    tCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAg                                                           >  1:232196/1‑68 (MQ=255)
           gAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                                   <  1:1960983/69‑1 (MQ=255)
           gAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                                   <  1:1984604/69‑1 (MQ=255)
           gAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGaaaaa                                                   <  1:1999461/69‑1 (MQ=255)
                     aCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGtt                                         >  1:121031/1‑69 (MQ=255)
                     aCCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGtt                                         <  1:579078/69‑1 (MQ=255)
                       ccccGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCggg                                             <  1:708170/63‑1 (MQ=255)
                              tGTACCGTACAGGCCTGCGGCTTACGCGATAACATg                                                                 <  1:3656038/36‑1 (MQ=255)
                                  ccGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGt                             <  1:3584104/68‑1 (MQ=255)
                                  ccGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGt                             >  1:483274/1‑68 (MQ=255)
                                  ccGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCAGGCGTTGATGTGCTGGGt                             >  1:2193474/1‑68 (MQ=255)
                                   cGTACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTg                                        <  1:732430/56‑1 (MQ=255)
                                     tACAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATc                          >  1:2463961/1‑68 (MQ=255)
                                       cAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                             <  1:698592/47‑1 (MQ=255)
                                       cAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                             <  1:1339997/47‑1 (MQ=255)
                                       cAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGtt                                         <  1:54479/51‑1 (MQ=255)
                                       cAGGCCTGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGc                       >  1:3575980/1‑69 (MQ=255)
                                       cAGGCCTCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCggg                                             <  1:2890762/46‑1 (MQ=255)
                                          gCCTGCGGCTTACGCGATAACATAGATGAGTTGaaa                                                     >  1:879055/1‑36 (MQ=255)
                                          gCCTGCGGCTTACGCGATAACATAGATGAGTTGaaa                                                     >  1:894274/1‑36 (MQ=255)
                                          gCCTGCGGCTTACGCGATAACATAGATGAGTTGaaa                                                     >  1:3433315/1‑36 (MQ=255)
                                          gCCTGCGGCTTACGCGATAACATAGATGAGTTGaaa                                                     >  1:1946027/1‑36 (MQ=255)
                                             tGCGGCTTACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGAt                 <  1:1747324/69‑1 (MQ=255)
                                                   ttACGCGATAACATAGATGAGTTGAAAAAAGCGGGCGtt                                         <  1:2866010/39‑1 (MQ=255)
                                                      cgcgATAACATGGATGAGTTGAACAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGa                  <  1:1487146/59‑1 (MQ=255)
                                                      cgcgATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGa                  <  1:2711484/59‑1 (MQ=255)
                                                      cgcgATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGa         >  1:2112459/1‑68 (MQ=255)
                                                        cgATAACATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGaaa       <  1:3124498/68‑1 (MQ=255)
                                                            aaCATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAActc  >  1:2454091/1‑69 (MQ=255)
                                                             aCATAGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTAt                           >  1:23883/1‑43 (MQ=255)
                                                                 |                                                               
TATTTCTACCCGAAAGCCATGACCCCCGGCTGTACCGTACAGGCCTGCGGCTTACGCGATAACATGGATGAGTTGAAAAAAGCGGGCGTTGATGTGCTGGGTATCAGCACCGATAAACCCGAAAAACTC  >  W3110S.gb/2599239‑2599367

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: