Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,873,346 C→T 32.2% L416L (CTC→CTT yeaI → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,873,3460CT32.2% 37.5 / 22.7 31L416L (CTC→CTTyeaIpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (4/6);  total (12/19)
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.70e-01

AGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTT  >  W3110S.gb/1873278‑1873411
                                                                    |                                                                 
aGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTTaa                                                                 <  1:20590/71‑1 (MQ=255)
aGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCaa                                                                 <  1:1397283/71‑1 (MQ=255)
aGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCaa                                                                 <  1:1571452/71‑1 (MQ=255)
      cGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGc                                                           >  1:965411/1‑71 (MQ=255)
      cGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGc                                                           >  1:447299/1‑71 (MQ=255)
      cGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGc                                                           >  1:389247/1‑71 (MQ=255)
       gCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGcc                                                          <  1:1512982/71‑1 (MQ=255)
          ccGAAGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTTAATAGTGCCCAg                                                       >  1:729117/1‑71 (MQ=255)
               gatATTTTAGCGCGACTGGAGGGTGTGGTGTTTGGCTTGCTATTTACCGAACTTAATAGTg                                                            <  1:2401898/61‑1 (MQ=255)
                 tatTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCaaaaa                                                <  1:1782708/71‑1 (MQ=255)
                  atTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCa                                                        >  1:1982235/1‑62 (MQ=255)
                  atTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCa                                                        >  1:1562264/1‑62 (MQ=255)
                   ttttAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTTAATAGTGCCCAGGCAAAAAt                                               >  1:1066471/1‑70 (MQ=255)
                              cTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACg                                   <  1:555524/71‑1 (MQ=255)
                                      tGAGGTGTTTGGCTTGCTATTTACCGAACTTAATAGTGCCCAGGCAAAAATCATTGCGGAACGTa                                 >  1:2794799/1‑65 (MQ=255)
                                        aGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTaaa                         <  1:2713714/71‑1 (MQ=255)
                                          gtgtTTGGCTTGCTATTTACCGAACTTAATAGTGCCCAGGCAAAAATc                                              <  1:1996680/48‑1 (MQ=255)
                                          gtgtTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAAt                                               <  1:1893308/47‑1 (MQ=255)
                                             tttGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTa                                 >  1:2786313/1‑58 (MQ=255)
                                              ttGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTa                                 >  1:1947028/1‑57 (MQ=255)
                                               tGGCTTGCTATTTACCGAACTTAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTc                   >  1:1144347/1‑70 (MQ=255)
                                                ggCTTGCTATTTACCGAACTTAATAGTGCCCAGGCAAAAATCATTGAGGAACGTATGCGTAAAAATGTCg                  <  1:77419/70‑1 (MQ=255)
                                                ggCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCg                  <  1:2001850/70‑1 (MQ=255)
                                                   ttGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCAt                                            <  1:718309/41‑1 (MQ=255)
                                                     gCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACt              <  1:1453606/69‑1 (MQ=255)
                                                        aTTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAAt                      <  1:2111054/58‑1 (MQ=255)
                                                         tttACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACg                                   <  1:155333/44‑1 (MQ=255)
                                                           tACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGaa                >  1:1728149/1‑61 (MQ=255)
                                                               gAACTTAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCtt  <  1:1563528/71‑1 (MQ=255)
                                                               gAACTTAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCtt  <  1:1560426/71‑1 (MQ=255)
                                                                   tCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCtt  <  1:1797920/67‑1 (MQ=255)
                                                                    |                                                                 
AGCATCCGCCCCGATGATATTTTAGCGCGACTGGAGGGTGAGGTGTTTGGCTTGCTATTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTT  >  W3110S.gb/1873278‑1873411

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