Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,252,821 T→C 13.7% intergenic (+33/+74) 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,8210TC13.7% 73.3 / 2.6 36intergenic (+33/+74)yeiI/nupXputative kinase/nucleoside permease
Reads supporting (aligned to +/- strand):  ref base T (19/12);  new base C (2/3);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCCGCCAGCGC  >  NC_000913/2252677‑2252966
                                                                                                                                                |                                                                                                                                                 
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                                                                                   >  3:78771/1‑148 (MQ=255)
                                                             aacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAgg                                                                                   >  5:311955/1‑148 (MQ=255)
                                                                 ctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCaa                                                                                                                    <  3:61580/111‑1 (MQ=255)
                                                                 ctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCaa                                                                                                                    >  4:61580/1‑111 (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)
                                                                                                                                             tCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCCGCCAgcgc  >  1:138552/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCCGCCAGCGC  >  NC_000913/2252677‑2252966

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: