Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,190,401 G→A 61.0% Q411Q (CAG→CAA yeaJ → predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,190,4010GA61.0% 27.0 / 41.6 64Q411Q (CAG→CAAyeaJpredicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (18/7);  new base A (15/24);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

AAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTC  >  minE/1190335‑1190464
                                                                  |                                                               
aaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTc                                                               <  1:272281/69‑1 (MQ=255)
aaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTc                                                               <  1:118807/69‑1 (MQ=255)
aaatGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTc                                                               <  1:3002779/69‑1 (MQ=255)
aaatGACACCATCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTc                                                               <  1:1449223/69‑1 (MQ=255)
                catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGa                                                >  1:2151070/1‑68 (MQ=255)
                catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGAt                                               >  1:1880800/1‑69 (MQ=255)
                catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGAt                                               >  1:2762757/1‑69 (MQ=255)
                catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGAt                                               <  1:1182617/69‑1 (MQ=255)
                catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGAt                                               <  1:2796204/69‑1 (MQ=255)
                catcaGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGAt                                               >  1:3353291/1‑69 (MQ=255)
                      gaggGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGcc                                         <  1:397427/69‑1 (MQ=255)
                         ggggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCAt                                       >  1:2377079/1‑68 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:3273345/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:252662/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:2544066/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:599388/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:206130/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:2052137/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:2652899/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:1962113/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:1955846/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:536006/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:1668059/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:2721611/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:1145476/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:3403988/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:1194568/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:1327498/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:30348/42‑1 (MQ=255)
                           ggATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAATc                                                               <  1:3169309/42‑1 (MQ=255)
                               ttAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACt                                 >  1:2280137/1‑68 (MQ=255)
                               ttAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTc                                >  1:2253851/1‑69 (MQ=255)
                               ttAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTc                                >  1:1068727/1‑69 (MQ=255)
                               ttAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTc                                >  1:1078944/1‑69 (MQ=255)
                               ttAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTc                                >  1:1121252/1‑69 (MQ=255)
                               ttAGCGATTAAGTTATTAGCTCAGGCGATTAAACAGTCGAtt                                                           >  1:2931254/1‑42 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTc       >  1:96189/1‑68 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:1580808/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:2019275/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:1157716/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:3315473/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:1373139/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:897080/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:683094/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:2423494/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:2928029/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:2806113/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:1552835/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:888032/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:2668677/1‑69 (MQ=255)
                                                         gATTAAACAATCGATTCGTAAAAGTGACTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCt      >  1:154661/1‑69 (MQ=255)
                                                           ttAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    <  1:1560398/69‑1 (MQ=255)
                                                           ttAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    <  1:272538/69‑1 (MQ=255)
                                                           ttAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    <  1:1296675/69‑1 (MQ=255)
                                                           ttAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    <  1:1579912/69‑1 (MQ=255)
                                                           ttAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    <  1:3420787/69‑1 (MQ=255)
                                                           ttAAACAATCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    <  1:2241791/69‑1 (MQ=255)
                                                            tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCg                                    >  1:1070329/1‑36 (MQ=255)
                                                            tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    >  1:3023858/1‑68 (MQ=255)
                                                            tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    >  1:1420589/1‑68 (MQ=255)
                                                            tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    >  1:963980/1‑68 (MQ=255)
                                                            tAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTg    >  1:982916/1‑68 (MQ=255)
                                                             aaaCAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGt   >  1:2757913/1‑68 (MQ=255)
                                                             aaaCAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTc  >  1:1450354/1‑69 (MQ=255)
                                                                  |                                                               
AAATGACACCCTCGGTCATCAGGAGGGGGATTTAGCGATTACGTTATTAGCTCAGGCGATTAAACAGTCGATTCGTAAAAGTGATTATGCCATCCGACTCGGTGGCGATGAATTCTGCATCATTCTTGTC  >  minE/1190335‑1190464

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