Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,864,768 T→C 50.0% intergenic (+34/+14) 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,7680TC50.0% ‑4.1 / 39.8 28intergenic (+34/+14)yeaD/yeaED‑hexose‑6‑phosphate epimerase‑like protein/aldo‑keto reductase, methylglyoxal to acetol, NADPH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base C (7/7);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

TCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGCTTTT  >  NC_000913/1864663‑1864914
                                                                                                         |                                                                                                                                                  
tCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                                                                                  <  1:695986/108‑1 (MQ=255)
tCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATcc                                                                                                                                                  >  2:695986/1‑108 (MQ=255)
                     aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCg                                                                                                                                                 <  4:31009/88‑1 (MQ=38)
                     aaaGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCg                                                                                                                                                 >  3:31009/1‑88 (MQ=38)
                       agagAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                                <  1:576536/87‑1 (MQ=38)
                       agagAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGATGCATTGCACCCgg                                                                                                                                                >  2:576536/1‑87 (MQ=38)
                                                                                         tgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTa                                                                                      <  1:830231/79‑1 (MQ=255)
                                                                                         tgCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTa                                                                                      >  2:830231/1‑79 (MQ=255)
                                                                                           ccGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTttgtt                            <  4:90456/135‑1 (MQ=255)
                                                                                           ccGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTttgtt                            >  3:90456/1‑135 (MQ=255)
                                                                                            cggatGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCAt                                                                              <  4:390358/80‑1 (MQ=38)
                                                                                            cggatGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTttgt                             >  3:1114215/5‑133 (MQ=255)
                                                                                            cggatGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTttgt                             <  4:1114215/129‑1 (MQ=255)
                                                                                            cGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc                                                                  >  3:1104241/1‑96 (MQ=255)
                                                                                            cGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGc                                                                  <  4:1104241/96‑1 (MQ=255)
                                                                                            cGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCAt                                                                              >  3:390358/1‑84 (MQ=38)
                                                                                            cGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                                                                     >  2:798509/1‑93 (MQ=38)
                                                                                            cGGATGCATTGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTc                                                                     <  1:798509/93‑1 (MQ=38)
                                                                                                  cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGcac                                          <  1:386163/114‑1 (MQ=255)
                                                                                                  cAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGcac                                          >  2:386163/1‑114 (MQ=255)
                                                                                                    ttGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATg                                                                                           <  3:289691/62‑1 (MQ=255)
                                                                                                    ttGCACCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATg                                                                                           >  4:289691/2‑63 (MQ=255)
                                                                                                     tGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                                                                                >  3:42638/1‑73 (MQ=255)
                                                                                                     tGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGc                                                                                <  4:42638/73‑1 (MQ=255)
                                                                                                     tGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGCt     >  2:538433/1‑148 (MQ=255)
                                                                                                      gCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTga                                                               >  3:206708/1‑89 (MQ=255)
                                                                                                      gCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTga                                                               <  4:206708/89‑1 (MQ=255)
                                                                                                        aTCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGCtttt  <  1:538433/148‑1 (MQ=255)
                                                                                                         |                                                                                                                                                  
TCAGAAACGCAAAAAGTGACCAAAGAGAAACCTGCACATCTGGCGCAATCCATTCGCGTTGCGAAACGTTAATTTACGTTAATGTTGTGTGCCGGGTGCAATGCATCCGGCACACAACATCACACCATATCCAGCGCAGTTTTTCCTTTTGGTGCCGGATATGCCTTATCCAGCATAGCTAATTCCGCTGAAGAAAGTTCGACCTCAAGCACAGCCGCATTTTGTTGGACATGGGCAATCGTGGCCGCTTTT  >  NC_000913/1864663‑1864914

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