Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,872,469 C→T 8.2% Y143Y (TAC→TAT)  yeaJ → putative diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,872,4690CT 8.2% 161.8 / 3.6 50Y143Y (TAC→TAT) yeaJputative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base T (2/2);  total (25/25)
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

ACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGT  >  NC_000913/1872335‑1872604
                                                                                                                                       |                                                                                                                                       
aCTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTc                                                                                                                                    >  1:939418/1‑140 (MQ=255)
aCTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTc                                                                                                                                    >  1:516893/1‑140 (MQ=255)
              cccGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTATGGTTCAGGTCTTGGGATgg                                                                                                                      >  2:333250/1‑140 (MQ=255)
              cccGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTATGGTTCAGGTCTTGGGATgg                                                                                                                      >  2:333256/1‑140 (MQ=255)
              cccGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATgg                                                                                                                      >  2:821682/1‑140 (MQ=255)
              cccGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATgg                                                                                                                      >  1:610652/1‑140 (MQ=255)
                 gggCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCg                                                                                                                   >  2:825301/1‑140 (MQ=255)
                       aCGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAgaga                                                                                                             <  1:889236/140‑1 (MQ=255)
                            gCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAAttt                                                                                                        <  1:354249/140‑1 (MQ=255)
                                       aCCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTAtct                                                                                             >  1:477599/1‑140 (MQ=255)
                                          gCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGa                                                                                          >  2:959489/1‑140 (MQ=255)
                                            aTCCACTGACTGCGACACAATTTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAt                                                                                         >  1:335762/1‑140 (MQ=255)
                                            aTCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                                                        >  1:912038/1‑140 (MQ=255)
                                            aTCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                                                        >  1:911969/1‑140 (MQ=255)
                                            aTCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                                                        >  1:621795/1‑140 (MQ=255)
                                            aTCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                                                        >  1:621531/1‑140 (MQ=255)
                                            aTCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                                                        >  1:335779/1‑140 (MQ=255)
                                            aTCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtc                                                                                        >  2:4946/1‑140 (MQ=255)
                                                 ctgactgCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAAt                                                                                   <  2:333974/140‑1 (MQ=255)
                                                 ctgactgCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAAt                                                                                   <  1:685432/140‑1 (MQ=255)
                                                        cGACACAC‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCTCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCt                                                                            <  2:488199/140‑1 (MQ=255)
                                                        cGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCt                                                                            >  1:846135/1‑140 (MQ=255)
                                                        cGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCt                                                                            >  1:154087/1‑140 (MQ=255)
                                                            acaA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGt                                                                        <  2:1053267/140‑1 (MQ=255)
                                                              aa‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGACATTTCACTATTATCTCGATCTCAATGACCGCGATGTCt                                                                      >  2:948912/1‑140 (MQ=255)
                                                                               gcaCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCGCTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGt                                                      <  2:929522/140‑1 (MQ=255)
                                                                                 ccTATTACCCGCTTTTGATCTGGCAGTGGAAGGCAACCGCCACATGGATGCTTAGGGTTCAGGTCTTGGGATGGCCGAAGAGAAACTACACTATTATCTCGACCTCAAGGACCGCTATGTCCAGGTTTAGGAGCCGGTTa                                                    <  2:265849/139‑1 (MQ=255)
                                                                                 aCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACGAGGATGATTACGGTTGAGGTCTTGGGATGGCCGAAGAGTAATTTCACTATTATCTCGATGTCAATGACCGCTATGTCTATTTTTATGAGCCGGttt                                                    <  2:620750/140‑2 (MQ=255)
                                                                                 aCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGTTGATTACGGTTCTGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTa                                                    <  2:1125203/140‑1 (MQ=255)
                                                                                 aCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTa                                                    <  2:282647/140‑1 (MQ=255)
                                                                                 aCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTa                                                    <  2:365751/140‑1 (MQ=255)
                                                                                 aCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTa                                                    <  2:450488/140‑1 (MQ=255)
                                                                                 aCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTa                                                    <  2:499985/140‑1 (MQ=255)
                                                                                       aCCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTg                                              >  1:890506/1‑140 (MQ=255)
                                                                                       aCCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTg                                              >  1:87224/1‑140 (MQ=255)
                                                                                        cccGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGa                                             <  1:678734/140‑1 (MQ=255)
                                                                                               ttGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACttt                                      >  1:755921/1‑140 (MQ=255)
                                                                                                      ggCAGTGGAAGGGAACCGCCACCAGGTTGATTTCGGTTCAGGTATTGGGATGGCGGTAGAGAAATTTCACTATTATCTCGTTCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGa                               <  2:797078/140‑1 (MQ=255)
                                                                                                      ggCAGTGGAAGGCAACCGCCACCAGGATGATTATGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGa                               <  1:333250/140‑1 (MQ=255)
                                                                                                      ggCAGTGGAAGGCAACCGCCACCAGGATGATTATGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGa                               <  1:333256/140‑1 (MQ=255)
                                                                                                      ggCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGa                               <  2:127716/140‑1 (MQ=255)
                                                                                                      ggCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGa                               <  2:925672/140‑1 (MQ=255)
                                                                                                          gTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATaa                           <  2:221787/140‑1 (MQ=255)
                                                                                                          gTGGAAGGCAACCGCCACCAGGATGAGTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATaa                           <  1:882453/140‑1 (MQ=255)
                                                                                                                       gCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTg              <  1:193739/140‑1 (MQ=255)
                                                                                                                         caccaGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCa            <  1:301510/140‑1 (MQ=255)
                                                                                                                               gatgatTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAgg      <  2:1145559/140‑1 (MQ=255)
                                                                                                                                   atTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCTATGAATAACTGGTCATTCCTGCAGTCAGGAAGt  >  2:642503/1‑140 (MQ=255)
                                                                                                                                   atTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGt  >  1:569454/1‑140 (MQ=255)
                                                                                                                                   atTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGAGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGt  >  2:642504/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ACTGCCCATAAATACCCGGGCTTACGCGGCACACTACAAACCGCATCCACTGACTGCGACACAA‑TTGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGT  >  NC_000913/1872335‑1872604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑