Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,252,810 T→G 14.3% intergenic (+22/+85) 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,8100TG14.3% 75.0 / 5.4 35intergenic (+22/+85)yeiI/nupXputative kinase/nucleoside permease
Reads supporting (aligned to +/- strand):  ref base T (18/12);  new base G (2/3);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCG  >  NC_000913/2252666‑2252953
                                                                                                                                                |                                                                                                                                               
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                                                                      >  3:78771/1‑148 (MQ=255)
                                                                        aacCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAgg                                                                      >  5:311955/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)
                                                                                                                                                |                                                                                                                                               
TACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCG  >  NC_000913/2252666‑2252953

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: