Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,248,780 C→T R280Q (CGG→CAG)  yeiE ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,248,7800CT97.1% 109.6 / ‑3.7 34R280Q (CGG→CAG) yeiEputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (17/16);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

AAGCAGCCAATCTCTTTTCTGGCGAATGCATGCTGAAGGCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCCGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGTGCCTGCCTGCAATTGATC  >  NC_000913/2248649‑2248889
                                                                                                                                   |                                                                                                             
aaGCAGCCAATCTCTTTTCTGGCGAATGCATGCTGAAGGCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAgcgc                                                                                                        >  1:137162/1‑139 (MQ=255)
 agcagcCAATCTCTTTTCTGGCGAATGCATGCTGAAGGCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGcgcg                                                                                                       >  2:427569/1‑139 (MQ=255)
                     gCGAATGCATGCTGAAGGCGCTGCGTCACGAAATCGACACAGCACCATCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCACTAGTCCAGAAAGCGCCGTAGCGCGTTGGAAAGGTGTTTTTGCCg                                                                                   <  1:412194/139‑1 (MQ=255)
                     gCGAATGCATGCTGAAGGCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCg                                                                                   <  2:137162/139‑1 (MQ=255)
                     gCGAATGCATGCTGAAGGCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCg                                                                                   >  1:521711/1‑139 (MQ=255)
                                     ggCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCcac                                                                   <  2:9667/139‑1 (MQ=255)
                                     ggCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCcac                                                                   >  1:9667/1‑139 (MQ=255)
                                       cgcTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACaa                                                                 <  1:427569/139‑1 (MQ=255)
                                           gCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTa                                                             <  2:321271/139‑1 (MQ=255)
                                                               cACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATc                                                       <  1:131344/125‑1 (MQ=255)
                                                               cACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATc                                                       >  2:131344/1‑125 (MQ=255)
                                                                       tGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCgatg                                                                                >  2:221585/1‑92 (MQ=255)
                                                                       tGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCgatg                                                                                <  1:221585/92‑1 (MQ=255)
                                                                              gTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACt                          >  2:458332/1‑139 (MQ=255)
                                                                                   gCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCg                                                                                   <  1:256070/77‑1 (MQ=255)
                                                                                   gCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCg                                                                                   >  2:256070/1‑77 (MQ=255)
                                                                                   gCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAAtg                     >  1:164604/1‑139 (MQ=255)
                                                                                   gCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAAtg                     >  2:325498/1‑139 (MQ=255)
                                                                                    cAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGc                                                                                     >  1:86269/1‑74 (MQ=255)
                                                                                    cAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGc                                                                                     >  1:391210/1‑74 (MQ=255)
                                                                                    cAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGc                                                                                     <  2:86269/74‑1 (MQ=255)
                                                                                    cAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGc                                                                                     <  2:391210/74‑1 (MQ=255)
                                                                                    cAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCcac                                                                   <  1:261254/92‑1 (MQ=255)
                                                                                    cAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCcac                                                                   >  2:261254/1‑92 (MQ=255)
                                                                                      aCTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCg                                                                                   >  1:545565/1‑74 (MQ=255)
                                                                                      aCTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCg                                                                                   <  2:545565/74‑1 (MQ=255)
                                                                                       cTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGtgcc                 >  2:505286/1‑139 (MQ=255)
                                                                                           aCGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGtgcctgcc             <  1:325498/139‑1 (MQ=255)
                                                                                              cgGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGTGCctgcctgc          <  2:511652/139‑1 (MQ=255)
                                                                                                gCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGtgc                  <  2:408448/129‑1 (MQ=255)
                                                                                                gCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGtgc                  >  1:408448/1‑129 (MQ=255)
                                                                                                      ttGCGGGATCGCAATAGTCCAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGTGCCTGCCTGCAATTGATc  <  2:164604/139‑1 (MQ=255)
                                                                                                                         cAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAAtg                     <  2:260498/101‑1 (MQ=255)
                                                                                                                         cAGAAAGCGCTGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAAtg                     >  1:260498/1‑101 (MQ=255)
                                                                                                                                   |                                                                                                             
AAGCAGCCAATCTCTTTTCTGGCGAATGCATGCTGAAGGCGCTGCGTCACGAAATCGACACAGCACCAGCATGTTCTTGTACAGCAACTTAACGCGGCACATTTGCGGGATCGCAATAGTCCAGAAAGCGCCGTAGCGCGTTGGAAAGGTGTTTTTGCCGATGATGTATACGCCACAACGTACGCATCAGGCGCGGCAGAGGGACCGCAACTTCACTTAATGTGCCTGCCTGCAATTGATC  >  NC_000913/2248649‑2248889

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