Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,252,814 C→A 14.3% intergenic (+26/+81) yeiI → / ← nupX putative kinase/nucleoside permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,252,8140CA14.3% 74.3 / 3.7 34intergenic (+26/+81)yeiI/nupXputative kinase/nucleoside permease
Reads supporting (aligned to +/- strand):  ref base C (17/12);  new base A (2/3);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCCGC  >  NC_000913/2252666‑2252960
                                                                                                                                                    |                                                                                                                                                  
tACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAc                                                                                                                                                    >  3:2372/1‑149 (MQ=255)
           cgACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATcc                                                                                                                                         <  4:2372/149‑1 (MQ=255)
            gACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTCGCATCCg                                                                                                                                        >  2:315750/1‑149 (MQ=255)
              cAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGc                                                                                                                                      <  3:264221/149‑1 (MQ=255)
               aGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTCGCATCCGGca                                                                                                                                     <  1:315750/149‑2 (MQ=255)
                gcgcCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTCGCATCCGGcaa                                                                                                                                    <  8:273227/149‑3 (MQ=38)
                   cccGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTAc                                                                                                                                  >  7:362665/1‑148 (MQ=255)
                    ccGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACa                                                                                                                                 >  1:189150/1‑148 (MQ=255)
                                      cgcgGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGc                                                                                                                <  8:122619/147‑1 (MQ=255)
                                         ggCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGt                                                                                                            >  8:85729/1‑148 (MQ=255)
                                               ttCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATg                                                                                                       >  8:303864/1‑147 (MQ=255)
                                                 cgcgcgCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAg                                                                                                    <  8:362665/148‑1 (MQ=255)
                                                                gCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTgcg                                                                                    <  7:303864/149‑1 (MQ=255)
                                                                  ttAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCa                                                                                  >  6:38034/1‑149 (MQ=255)
                                                                     aacaacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTTACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATc                                                                                <  4:141435/148‑1 (MQ=255)
                                                                        aacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCt                                                                                        >  6:329292/1‑137 (MQ=255)
                                                                        aacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCt                                                                                        <  5:329292/137‑1 (MQ=255)
                                                                        aacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAgg                                                                             >  5:311955/1‑148 (MQ=255)
                                                                        aacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAgg                                                                             >  3:78771/1‑148 (MQ=255)
                                                                            ctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCaa                                                                                                              >  4:61580/1‑111 (MQ=255)
                                                                            ctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCaa                                                                                                              <  3:61580/111‑1 (MQ=255)
                                                                               ccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGaca                                                                     >  5:40991/1‑149 (MQ=255)
                                                                                      tGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGc                                                              >  4:138539/1‑149 (MQ=255)
                                                                                               cgcgCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAAt                                                     <  2:189150/149‑1 (MQ=255)
                                                                                                               ccAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCa                                                                           <  2:292573/111‑1 (MQ=255)
                                                                                                                  aTGGAGTAACAATGTTCACCGATGCTACGCTTAACTCCTCTTATACCGCCTACAAATCAGCCAAATACAGCAAGAAAATCAACGCATGATAAGCTTACTCATCAGGCAATGACACTCAAGCTAAACCAATAAAGACACCGGAAATggac                                  >  4:146836/1‑147 (MQ=255)
                                                                                                                      tgtAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCAc                              <  2:37338/149‑1 (MQ=255)
                                                                                                                            aaTGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCATGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGCCACTCAAGCGAAACCAATTAAGCACCCGGCAATGGTCGCACtcatca                        >  7:348591/1‑149 (MQ=39)
                                                                                                                                   ccGGATAAGACGCGTTAAGCGTCGCATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTgg                  >  3:88695/1‑148 (MQ=39)
                                                                                                                                      gATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAg              >  4:370291/1‑149 (MQ=255)
                                                                                                                                      gATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAg              >  1:19565/1‑149 (MQ=255)
                                                                                                                                       aTGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGc             <  6:40991/149‑1 (MQ=255)
                                                                                                                                         gCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGt           >  4:165443/1‑149 (MQ=255)
                                                                                                                                          agACGCGTTAAGCGTCGCATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTc          <  4:88695/148‑1 (MQ=38)
                                                                                                                                           gACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcg         >  6:20738/1‑149 (MQ=255)
                                                                                                                                                  accgcgTCGTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCGTGCGCATCAGGCACTGACACTCAAGCTATACCAATAAAGAACCGGGCAATGGTCGCACTCATCAAGTTGTCAAGCGTcgccgccgc  >  6:296035/3‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCCGC  >  NC_000913/2252666‑2252960

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