Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,872,347 A→G 7.7% K83K (AAA→AAG yeaI → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,872,3470AG 7.7% 134.9 / 4.4 52K83K (AAA→AAGyeaIpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (29/19);  new base G (2/2);  total (31/21)
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.50e-01

TCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATTTTTCT  >  W3110S.gb/1872280‑1872415
                                                                   |                                                                    
tCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAGCt                                                                    <  1:3483190/70‑1 (MQ=255)
 cAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTc                                                                   <  1:2227938/70‑1 (MQ=255)
      tATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCca                                                             >  1:2711647/1‑71 (MQ=255)
       aTTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaaaa                                                                      >  1:2650195/1‑61 (MQ=255)
          gCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACt                                                                    >  1:2976029/1‑60 (MQ=255)
           cTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGa                                                         <  1:1225227/70‑1 (MQ=255)
            tGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAGCTCATCCCACAGAt                                                        >  1:4021087/1‑70 (MQ=255)
             gACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATg                                                       >  1:3376441/1‑70 (MQ=255)
             gACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATg                                                       >  1:211037/1‑70 (MQ=255)
              aCGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGaa                                                     >  1:1960336/1‑71 (MQ=255)
              aCGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGaa                                                     >  1:1840993/1‑71 (MQ=255)
                gTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAAtt                                                   <  1:2356288/71‑1 (MQ=255)
                gTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAAt                                                    >  1:1173044/1‑70 (MQ=255)
                   aaTTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTaaaa                                                                      >  1:3721355/1‑49 (MQ=255)
                       ttattttGCATTACTACCTTATATAACGAAAACACAAATGTaaaa                                                                      >  1:3182321/1‑45 (MQ=255)
                       ttattttGCATTACTACCTTATATAACGAAAACACAAATGTaaaa                                                                      >  1:1765983/1‑45 (MQ=255)
                           tttGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCAt                                                                 <  1:4037840/46‑1 (MQ=255)
                             tGCATTACTACCTTATATAACGAAAACACAAATGTaaaa                                                                      <  1:1075902/39‑1 (MQ=255)
                             tGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGgttgt                                      >  1:913767/1‑71 (MQ=255)
                             tGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGgttgt                                      >  1:1392705/1‑71 (MQ=255)
                              gCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGgttgt                                      <  1:1307904/70‑1 (MQ=255)
                                aTTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTtgttgt                                   <  1:3553710/71‑1 (MQ=255)
                                 ttACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATg                                                       <  1:4272609/50‑1 (MQ=255)
                                 ttACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGa                                                      <  1:1788818/51‑1 (MQ=255)
                                   actacCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAAt                                                    >  1:4176323/1‑51 (MQ=255)
                                    ctacCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTa                                                  >  1:365861/1‑52 (MQ=255)
                                    ctacCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATCAATTAcc                                                >  1:2630371/1‑54 (MQ=255)
                                      acCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAAtt                                                   <  1:214167/49‑1 (MQ=255)
                                       ccTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCttt                             >  1:3834559/1‑70 (MQ=255)
                                       ccTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCtt                              <  1:1446915/69‑1 (MQ=255)
                                        cTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGt                           <  1:3463320/71‑1 (MQ=255)
                                        cTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGt                           >  1:1847393/1‑71 (MQ=255)
                                         ttATATAACGAAAACACAAATGTAAAACTCATCCcac                                                            <  1:605896/37‑1 (MQ=255)
                                         ttATATAACGAAAACACAAATGTAAAACTCATCCCACAGATg                                                       <  1:3557266/42‑1 (MQ=255)
                                          tatataACGAAAACACAAATGTAAAACTCATCCCACAGa                                                         >  1:215322/1‑39 (MQ=255)
                                           atataACGAAAACACAAATGTAAAACTCATCCcac                                                            >  1:1862949/1‑35 (MQ=255)
                                           atataACGAAAACACAAATGTAAAACTCATCCCACAg                                                          >  1:2631508/1‑37 (MQ=255)
                                           atataACGAAAACACAAATGTAAAACTCATCCCACAg                                                          >  1:2081449/1‑37 (MQ=255)
                                           atataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGtttt                        >  1:3151463/1‑71 (MQ=255)
                                           atataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGtttt                        >  1:3216790/1‑71 (MQ=255)
                                            tataACGAAAACACAAATGTAAAGCTCATCCCACAGa                                                         <  1:1362668/37‑1 (MQ=255)
                                            tataACGAAAACACAAATGTAAAACTCATCCCACAGATGAAt                                                    >  1:627316/1‑42 (MQ=255)
                                            tataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGttttt                       >  1:3444702/1‑71 (MQ=255)
                                             ataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTc                      >  1:3127374/1‑71 (MQ=255)
                                             ataACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTc                      >  1:2336305/1‑71 (MQ=255)
                                              taACGAAAACACAAATGTAAAACTCATCCCACAGATGa                                                      <  1:2680697/38‑1 (MQ=255)
                                               aaCGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCt                    >  1:1326364/1‑71 (MQ=255)
                                                  gAAAACACAAATGTAAAGCTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTaa                 >  1:531785/1‑71 (MQ=255)
                                                  gAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTaa                 >  1:3046632/1‑71 (MQ=255)
                                                   aaaaCACATATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTc                      <  1:171181/65‑1 (MQ=255)
                                                   aaaaCACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGgt                                         <  1:2063157/46‑1 (MQ=255)
                                                    aaaCACAAATGTAAAACTCATCCCACAGATGAATTa                                                  >  1:1086164/1‑36 (MQ=255)
                                                     aaCACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCt                               >  1:4238936/1‑54 (MQ=255)
                                                     aaCACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGc              <  1:3505670/71‑1 (MQ=255)
                                                     aaCACAAATGTAAAACTCATCACACAGATGAATTACCTGATGgttgt                                      <  1:2551265/47‑1 (MQ=255)
                                                      acacAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTtgttgt                                   >  1:2831845/1‑49 (MQ=255)
                                                               tAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCg            <  1:1377697/63‑1 (MQ=255)
                                                                 aaaCTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATttttct  <  1:1536871/71‑1 (MQ=255)
                                                                 aaaCTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATttttct  <  1:3922566/71‑1 (MQ=255)
                                                                 aaaCTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATttttct  <  1:1784619/71‑1 (MQ=255)
                                                                   aCTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGc              >  1:1182157/1‑57 (MQ=255)
                                                                   |                                                                    
TCAGTTTATTGCTGACGTTAATTTTATTTTGCATTACTACCTTATATAACGAAAACACAAATGTAAAACTCATCCCACAGATGAATTACCTGATGGTTGTTGTGGCTTTGTTTTTCCTTAACGCCGTCATTTTTCT  >  W3110S.gb/1872280‑1872415

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