Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,758 G→A 34.8% intergenic (+24/+24) yeaD → / ← yeaE D‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,864,7580GA34.8% 19.2 / 20.4 23intergenic (+24/+24)yeaD/yeaED‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (6/9);  new base A (4/4);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

TATGCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATC  >  NC_000913/1864611‑1864902
                                                                                                                                                   |                                                                                                                                                
tatGCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGt                                                                                                                                                 <  4:580305/149‑1 (MQ=255)
 atGCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTg                                                                                                                                                <  3:208612/149‑1 (MQ=255)
     cGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAAt                                                                                                                                            <  1:602417/149‑1 (MQ=255)
     cGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAAt                                                                                                                                            <  4:1416077/149‑1 (MQ=255)
          gatgGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCacc                                                                                                                                       >  2:534547/1‑147 (MQ=255)
             ggCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                    <  1:534547/149‑1 (MQ=255)
                    aaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                    >  1:685262/1‑142 (MQ=255)
                    aaaCATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                    <  2:685262/142‑1 (MQ=255)
                                   tAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                                                                      >  3:1062323/1‑125 (MQ=255)
                                   tAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                                                                      <  4:1062323/125‑1 (MQ=255)
                                    aGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAAt                                                                                                                                            >  2:538528/1‑118 (MQ=255)
                                    aGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAAt                                                                                                                                            <  1:538528/118‑1 (MQ=255)
                                        aCGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                    >  3:1208293/1‑122 (MQ=255)
                                        aCGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                    <  4:1208293/122‑1 (MQ=255)
                                                      aGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGt                                                                                                                                                 <  1:38695/95‑1 (MQ=255)
                                                      aGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGt                                                                                                                                                 >  2:38695/1‑95 (MQ=255)
                                                                                                      cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGc                                                                                                                                          <  2:673609/54‑1 (MQ=255)
                                                                                                      cATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGc                                                                                                                                          >  1:673609/1‑54 (MQ=255)
                                                                                                                                             tgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAg                            >  3:38017/1‑125 (MQ=255)
                                                                                                                                             tgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAg                            <  4:38017/125‑1 (MQ=255)
                                                                                                                                             tgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGaaa                                              <  3:725914/107‑1 (MQ=39)
                                                                                                                                             tgCCGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGaaa                                              >  4:725914/1‑107 (MQ=38)
                                                                                                                                               ccGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATc  >  4:615105/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TATGCCGGATGATGGCTACAAAACATTTGTTTGTGTAGAAACGGCTTACGCTTCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATC  >  NC_000913/1864611‑1864902

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