Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,874,182 A→G 20.0% Y143C (TAC→TGC)  yeaJ → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,874,1820AG20.0% 117.8 / 24.8 55Y143C (TAC→TGC) yeaJpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (9/2);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAT  >  W3110S.gb/1874114‑1874252
                                                                    |                                                                      
tgtgGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTGCg                                                                      >  1:805491/1‑71 (MQ=255)
 gtgGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACgg                                                                     >  1:424133/1‑71 (MQ=255)
     aaGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGtt                                                                   <  1:3099386/69‑1 (MQ=255)
      aGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAg                                                                <  1:4235728/71‑1 (MQ=255)
      aGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAg                                                                <  1:1585004/71‑1 (MQ=255)
             gcaCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTa                                                                        >  1:2877172/1‑56 (MQ=255)
             gcaCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTgg                                                         >  1:3921494/1‑71 (MQ=255)
             gcaCTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTgg                                                         <  1:3145163/71‑1 (MQ=255)
                cTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTGCGGTTCAGGTCTTGGGa                                                       >  1:548318/1‑70 (MQ=255)
                cTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTGCGGTTCAGGTCTTGGGa                                                       >  1:3045783/1‑70 (MQ=255)
                cTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTGCGGTTCAGGTCTTGGGa                                                       >  1:3682977/1‑70 (MQ=255)
                  aTTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATgg                                                    >  1:2759292/1‑71 (MQ=255)
                     aCCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGcc                                                  >  1:4197906/1‑70 (MQ=255)
                     aCCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGcc                                                  >  1:132282/1‑70 (MQ=255)
                      cccGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCg                                                 >  1:1670796/1‑70 (MQ=255)
                      cccGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCg                                                 >  1:1611585/1‑70 (MQ=255)
                      cccGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCg                                                 >  1:3887574/1‑70 (MQ=255)
                       ccGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGaa                                               <  1:3340990/71‑1 (MQ=255)
                       ccGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGaa                                               <  1:1231452/71‑1 (MQ=255)
                       ccGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGaa                                               <  1:1992890/71‑1 (MQ=255)
                        cGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTGCGGTTCAGGTCTTGGGATGGCCg                                                 >  1:1787164/1‑68 (MQ=255)
                        cGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGaa                                               <  1:3965565/70‑1 (MQ=255)
                         gCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGa                                                       >  1:3616785/1‑61 (MQ=255)
                         gCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGaa                                               >  1:4394330/1‑69 (MQ=255)
                          cTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAagag                                            <  1:4475463/71‑1 (MQ=255)
                             ttGATCAGTCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAgg                                                               >  1:156843/1‑49 (MQ=255)
                             ttGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTa                                                                        >  1:2139802/1‑40 (MQ=255)
                                       aGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACtat                                >  1:3316534/1‑70 (MQ=255)
                                          ggAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTAtc                            >  1:1583428/1‑71 (MQ=255)
                                           gAAGGCAACCGCCACCAGGATGATTGCGGTTCAGGTCTTgg                                                         >  1:3783933/1‑41 (MQ=255)
                                           gAAGGCAACCGCCACCAGGATGATTGCGGTTCAGGTCTTgg                                                         >  1:3692006/1‑41 (MQ=255)
                                           gAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGaa                                               >  1:3028148/1‑51 (MQ=255)
                                                  aCCGCCACCAGGATGATTGCGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtct                     >  1:2804919/1‑70 (MQ=255)
                                                  aCCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAtctc                    <  1:450478/71‑1 (MQ=255)
                                                   ccGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCa                   <  1:2656255/71‑1 (MQ=255)
                                                    cGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAAtt                                       >  1:78905/1‑50 (MQ=255)
                                                    cGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCaa                  >  1:3378141/1‑71 (MQ=255)
                                                     gCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAAt                 <  1:3183324/71‑1 (MQ=255)
                                                        accaGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCaa                  >  1:1071814/1‑67 (MQ=255)
                                                        accaGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCaa                  >  1:1325215/1‑67 (MQ=255)
                                                        accaGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGAc              <  1:374875/71‑1 (MQ=255)
                                                        accaGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGAc              <  1:3138371/71‑1 (MQ=255)
                                                          caGGATGATTACGGTTCAGGTCTTGGGATGGCCGa                                                <  1:701696/35‑1 (MQ=255)
                                                          caGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAt                       >  1:2660167/1‑60 (MQ=255)
                                                          caGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAt                       >  1:4920366/1‑60 (MQ=255)
                                                          caGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATATCACTATTATCTCGAt                       >  1:2505926/1‑60 (MQ=255)
                                                           aGGATGATTACGGTTCAGGTCTTTGGATGGCCGAAGAGAAATTTCACta                                 >  1:1174594/1‑49 (MQ=255)
                                                           aGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGaaa                                         <  1:4271932/41‑1 (MQ=255)
                                                           aGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGaaa                                         <  1:5081759/41‑1 (MQ=255)
                                                           aGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACta                                 >  1:4692868/1‑49 (MQ=255)
                                                            ggATGATTGCGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGc           >  1:3958004/1‑70 (MQ=255)
                                                             gatgatTACGGTTCAGGTCTTGGGATGGCGGAAGAGAAATTTCAc                                   <  1:761343/45‑1 (MQ=255)
                                                             gatgatTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCAc                                   <  1:4942196/45‑1 (MQ=255)
                                                               tgatTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGAt                       >  1:1166765/1‑55 (MQ=255)
                                                                gatTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACtattat                             >  1:4423350/1‑48 (MQ=255)
                                                                 atTGCGGTTCAGGTCTTGGGATGGCCTAAGAGAAATTTCACTATTAtc                            <  1:5140272/48‑1 (MQ=255)
                                                                 atTGCGGGTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTAtc                            <  1:1590775/48‑1 (MQ=255)
                                                                    aCGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAt  >  1:2464542/1‑71 (MQ=255)
                                                                    |                                                                      
TGTGGAAGCTGCAGCACTATTACCCGCTTTTGATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTAT  >  W3110S.gb/1874114‑1874252

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: