Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,695,634 G→A 11.1% Q126* (CAG→TAG)  yfhQ ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,695,6340GA11.1% 153.8 / 8.9 54Q126* (CAG→TAG) yfhQpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (24/24);  new base A (2/4);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

AACCGCCATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGC  >  minE/1695563‑1695697
                                                                       |                                                               
aaCCGCCATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAAT‑‑CGACATGATAATGGCATTTCTGc                                                                <  1:354437/71‑1 (MQ=255)
         cGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCtt                                                         >  1:2609688/1‑71 (MQ=255)
            cAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCAt                                                      >  1:910072/1‑71 (MQ=255)
               gTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCtt                                                         <  1:208151/65‑1 (MQ=255)
               gTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCtt                                                         <  1:2050098/65‑1 (MQ=255)
               gTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCtt                                                         <  1:3883248/65‑1 (MQ=255)
                 tCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGt                                                  >  1:1400010/1‑70 (MQ=255)
                 tCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGt                                                  >  1:3786611/1‑70 (MQ=255)
                    gCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTc                                                        <  1:2017061/61‑1 (MQ=255)
                    gCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCaa                                               <  1:809149/70‑1 (MQ=255)
                     cGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCtt                                                         <  1:1969227/59‑1 (MQ=255)
                     cGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCtt                                                         <  1:845223/59‑1 (MQ=255)
                      gAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGcc                                            <  1:203985/71‑1 (MQ=255)
                      gAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGcc                                            <  1:2260840/71‑1 (MQ=255)
                        gCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGcc                                            <  1:2245419/69‑1 (MQ=255)
                          tGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAActc                                                           <  1:2157796/52‑1 (MQ=255)
                          tGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCtt                                                         <  1:59455/54‑1 (MQ=255)
                          tGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCtt                                                         <  1:393293/54‑1 (MQ=255)
                             aTTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACgcgc                                     >  1:415618/1‑71 (MQ=255)
                              ttCCGGGTTAGCCGCAATCGCGACATGATATTGGCATTTCTGCAACTCTTCATTGGTCAAGc                                             >  1:1596530/1‑62 (MQ=255)
                              ttCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAAc                                                             <  1:2854336/46‑1 (MQ=255)
                              ttCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAAc                                                             <  1:2293134/46‑1 (MQ=255)
                              ttCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCt                                    >  1:2507616/1‑71 (MQ=255)
                              ttCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCt                                    >  1:522146/1‑71 (MQ=255)
                                ccGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGAc                                         >  1:186003/1‑64 (MQ=255)
                                  gggTTAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTcgcg                                <  1:1483664/71‑1 (MQ=255)
                                     ttAGCCGCAATCGCGACATGATAATGGCATTTCTGCAAc                                                             <  1:1678977/39‑1 (MQ=255)
                                     ttAGCCGCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGAc                              >  1:3296318/1‑70 (MQ=255)
                                      tAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGAcc                             >  1:2952496/1‑70 (MQ=255)
                                       aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGAcc                             >  1:479282/1‑69 (MQ=255)
                                       aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                            <  1:357242/70‑1 (MQ=255)
                                       aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                            >  1:3684619/1‑70 (MQ=255)
                                       aGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                            <  1:3208332/70‑1 (MQ=255)
                                        gCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaa                           <  1:1168564/70‑1 (MQ=255)
                                          cGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTg                                                    >  1:173738/1‑43 (MQ=255)
                                          cGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTc                                                 >  1:3348117/1‑46 (MQ=255)
                                           gCAATCGCGACATGATAATGGCATTTCTACAACTCTTCATTGGTCAAGcc                                            <  1:1154782/50‑1 (MQ=255)
                                             aaTCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                          >  1:3545380/1‑66 (MQ=255)
                                                cgcgACAAGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGa                                          >  1:3191931/1‑47 (MQ=255)
                                                 gcgACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                          >  1:607866/1‑62 (MQ=255)
                                                  cgACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGcc                <  1:1138253/71‑1 (MQ=255)
                                                   gACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGa                                          >  1:11292/1‑44 (MQ=255)
                                                       tGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCACTcgc                                 <  1:3091255/49‑1 (MQ=255)
                                                        gATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCggg          >  1:3590394/1‑71 (MQ=255)
                                                             tGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTAt      <  1:1511481/70‑1 (MQ=255)
                                                              ggCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAgc                   >  1:2251990/1‑56 (MQ=255)
                                                              ggCATTTCTACAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCa               >  1:2741442/1‑60 (MQ=255)
                                                               gCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGTATTTGc  <  1:3131534/71‑1 (MQ=255)
                                                                cATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGc  >  1:3679179/1‑71 (MQ=255)
                                                                 aTTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCa                            <  1:1713659/44‑1 (MQ=255)
                                                                 aTTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAacac                       <  1:2592078/49‑1 (MQ=255)
                                                                    tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCaaa                          >  1:3264870/1‑43 (MQ=255)
                                                                    tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGc  >  1:771226/1‑67 (MQ=255)
                                                                    tCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGc  >  1:2045653/1‑67 (MQ=255)
                                                                       |                                                               
AACCGCCATCGCCAGGTTCAGCGAGCTGTATTCCGGGTTAGCCGCAATCGCGACATGATAATGGCATTTCTGCAACTCTTCATTGGTCAAGCCGACGCGCTCGCGACCAAACACCAGCGCCACCGGGGTATTTGC  >  minE/1695563‑1695697

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: