Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 484,203 C→G 20.0% G80A (GGC→GCC)  yeaG ← protein kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947484,2030CG20.0% 55.9 / 9.0 25G80A (GGC→GCC) yeaGprotein kinase
Reads supporting (aligned to +/- strand):  ref base C (10/10);  new base G (3/2);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCATCAGCTGCTTGAGTTTTTCCGCAAGCGAGGACTTGCCGCCGCCTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTCGGCCGCCGATGCGTAGGCACTGCGATCC  >  NC_002947/484058‑484352
                                                                                                                                                 |                                                                                                                                                     
ttCCATCAGCTGCTTGAGTTTTTCCGAAAGCGAGGACTTGGCGCCTCGTACCGTGCGCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGCCCATGCCCTCTTccatcccatg                                                                                                                                                   <  2:529669/150‑1 (MQ=255)
   cATCAGCTGCTTGAGTTTTTCCGCAAGCGAGGACTTGCCGCCGCCTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGaa                                                                                                                                                <  2:212491/150‑1 (MQ=255)
             ttGAGTTTTTCCGCAAGCGAGGACTTGCCGCCGCCTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCGCTCTTCCATGGCATGGAAga                                                                                                                                              <  2:363600/142‑2 (MQ=255)
             ttGAGTTTTTCCGCAAGCGAGGACTTGCCGCCGCCTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGGCATGGAAga                                                                                                                                              >  1:363600/1‑141 (MQ=255)
                                        ccgccgccTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGTGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCttgttg                                                                                                           <  1:290363/150‑1 (MQ=255)
                                          gccgtcTACCTGGCTCAGTAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGACGTACGACACGATCTGGTCGATGCACTCTTCCATTCCATGGTAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGga                                                                                                         <  2:125151/150‑1 (MQ=255)
                                          gccgccTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGga                                                                                                         >  2:373268/1‑150 (MQ=255)
                                             gccTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGagaa                                                                                                      <  1:373268/150‑1 (MQ=255)
                                             gccTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCCTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGagaa                                                                                                      >  1:98629/1‑150 (MQ=255)
                                                        ccAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACa                                                                                           >  1:350349/1‑150 (MQ=255)
                                                            cagATAGAGGATCTGTTTCTTCTCCTCCAGGCCTTTGTCCGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCt                                                                                       <  1:502099/150‑1 (MQ=255)
                                                                 agagGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGt                                                                                  >  2:265431/1‑150 (MQ=255)
                                                                                 ctcttcCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGAt                                                                  <  1:265431/150‑1 (MQ=255)
                                                                                  tcttcCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATc                                                                 >  1:572021/1‑150 (MQ=255)
                                                                                            cTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATgc                                                                                                 >  1:291492/1‑108 (MQ=255)
                                                                                            cTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATgc                                                                                                 <  2:291492/108‑1 (MQ=255)
                                                                                                                   gACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGt                                <  2:350349/150‑1 (MQ=255)
                                                                                                                     cacGATCTGGTCGATGCACTCTTCCATGGCATGGaa                                                                                                                                                <  1:248870/36‑1 (MQ=255)
                                                                                                                     cacGATCTGGTCGATGCACTCTTCCATGGCATGGaa                                                                                                                                                >  2:248870/1‑36 (MQ=255)
                                                                                                                     cacGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTc                              >  1:223470/1‑150 (MQ=255)
                                                                                                                      acGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTg                                                                                <  1:23878/99‑1 (MQ=255)
                                                                                                                      acGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTc                              <  1:555028/149‑1 (MQ=255)
                                                                                                                      acGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTc                              >  2:555028/1‑149 (MQ=255)
                                                                                                                      acGATCTGGTCGAGGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTg                                                                                >  2:23878/1‑99 (MQ=255)
                                                                                                                                     cACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTCGGCCGCCGATGCGTAg              >  2:393695/1‑150 (MQ=255)
                                                                                                                                      cgactTCCATGGCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTCGGCCGCCGATGCGTAgg             >  1:13514/4‑150 (MQ=255)
                                                                                                                                              atgccatgGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTCGGCCGCCGATGCGTAGGCACTGCGa     <  2:572021/150‑1 (MQ=255)
                                                                                                                                                 ccatgGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTCGGCCGCCGATGCGTAGGCACTGCGATcc  >  2:506546/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTCCATCAGCTGCTTGAGTTTTTCCGCAAGCGAGGACTTGCCGCCGCCTACCGGGCCCAGCAGATAGAGGATCTGTTTCTTCTCTTCCAGGCCTTGGGCGGCGTGGCGGAAGTACGACACGATCTGGTCGATGCACTCTTCCATGCCATGGAAGTCGGCAAAGGCCGGATAGCGGCGGATAACCTTGTTGGAGAAGATGCGCGACAGCCTGGAGTTGGTAGAGGTGTCGATCAGCTCTGGCTCACCGATGGCCAGCAACAGCCGTTCGGCCGCCGATGCGTAGGCACTGCGATCC  >  NC_002947/484058‑484352

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