Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,135,605 A→G 28.9% H40R (CAT→CGT)  yniA → predicted phosphotransferase/kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,135,6050AG28.9% 72.5 / 33.1 76H40R (CAT→CGT) yniApredicted phosphotransferase/kinase
Reads supporting (aligned to +/- strand):  ref base A (45/9);  new base G (12/10);  total (57/19)
Fisher's exact test for biased strand distribution p-value = 1.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAAC  >  minE/1135546‑1135670
                                                           |                                                                 
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:917187/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:2271148/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:2356534/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:2564526/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:2642964/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:2900390/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:3036280/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:380447/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:1505013/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActtttt                                                          >  1:620147/1‑69 (MQ=255)
gCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACGTGActttt                                                           >  1:2191622/1‑68 (MQ=255)
gCGAAATGAACTGCCAGCCGGACAAGTCCAAGCCGCATGGCATTTGCCCTATGCAGGACGTGac                                                               >  1:2267227/1‑63 (MQ=255)
   tAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTc                                                         <  1:416106/67‑1 (MQ=255)
   tAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTc                                                         <  1:628818/67‑1 (MQ=255)
     aTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCa                                                     >  1:51233/1‑69 (MQ=255)
     aTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCa                                                     >  1:1661540/1‑69 (MQ=255)
     aTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCa                                                     >  1:3012271/1‑69 (MQ=255)
     aTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCa                                                     >  1:519737/1‑69 (MQ=255)
     aTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCa                                                     >  1:2682628/1‑69 (MQ=255)
     aTGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCa                                                     >  1:11866/1‑69 (MQ=255)
          cTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAAtgt                                                <  1:1902996/69‑1 (MQ=255)
                    gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAggg                                      >  1:763161/1‑69 (MQ=255)
                    gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAggg                                      >  1:1256235/1‑69 (MQ=255)
                    gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAggg                                      >  1:2023253/1‑69 (MQ=255)
                    gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAggg                                      >  1:2724539/1‑69 (MQ=255)
                    gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAggg                                      >  1:3409389/1‑69 (MQ=255)
                    gagaAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAg                                        >  1:1119952/1‑67 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     <  1:698706/69‑1 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     <  1:2371893/69‑1 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     >  1:1661484/1‑69 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     <  1:895528/69‑1 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     >  1:2808530/1‑69 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     <  1:2532767/69‑1 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     <  1:2728667/69‑1 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     >  1:1711427/1‑69 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAAGGAACTGCTTCCCGGtttt                     <  1:195293/69‑1 (MQ=255)
                                     tGGCATTTGCGCTATGCAGGACATGATTTTTTCGTAAAATGTGATGAAAGGGAACTGCTTCCCGGtttt                     >  1:2138835/1‑69 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:2037528/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:2591521/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:1078972/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:2499015/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:1686203/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:2877966/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:614019/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:2935166/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:1341129/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTc                             >  1:3028816/1‑59 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAc                   >  1:1646751/1‑69 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAc                   >  1:135070/1‑69 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAc                   >  1:1871984/1‑69 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAc                   >  1:1967232/1‑69 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAc                   >  1:2406747/1‑69 (MQ=255)
                                       gCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCAGGTTTTAc                   >  1:1531159/1‑69 (MQ=255)
                                            tGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccgc                >  1:1318512/1‑67 (MQ=255)
                                            tGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccg                 >  1:405604/1‑66 (MQ=255)
                                                   tGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTAccg                 >  1:1761094/1‑59 (MQ=255)
                                                   tGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc         >  1:863714/1‑67 (MQ=255)
                                                   tGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc         >  1:2341113/1‑67 (MQ=255)
                                                   tGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc         >  1:973063/1‑67 (MQ=255)
                                                   tGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc         >  1:367947/1‑67 (MQ=255)
                                                   tGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc         >  1:316424/1‑67 (MQ=255)
                                                   tGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGcc         >  1:2988537/1‑67 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:726717/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:2428326/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:876642/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:873648/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:848388/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:1787746/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:30401/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:3227701/43‑1 (MQ=255)
                                                     cAGGACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:334798/43‑1 (MQ=255)
                                                     cAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAcc     >  1:645440/1‑69 (MQ=255)
                                                     cAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAcc     >  1:282799/1‑69 (MQ=255)
                                                     cAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGAcc     >  1:2048665/1‑69 (MQ=255)
                                                        gACGTGACTTTTTCGTCAAATGTGATGAAAGGGAACTGCt                               <  1:490855/40‑1 (MQ=255)
                                                        gACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAAc  >  1:484641/1‑69 (MQ=255)
                                                           |                                                                 
GCGTAATGAACTGCCTGGCGGAGAAGTCCATGCCGCATGGCATTTGCGCTATGCAGGACATGATTTTTTCGTCAAATGTGATGAAAGGGAACTGCTTCCCGGTTTTACCGCCGAAGCCGACCAAC  >  minE/1135546‑1135670

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: