Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,848,549 T→C 24.9% P828P (CCA→CCG yjeP ← predicted mechanosensitive channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,848,5490TC24.9% 74.7 / 20.7 64P828P (CCA→CCGyjePpredicted mechanosensitive channel
Reads supporting (aligned to +/- strand):  ref base T (20/28);  new base C (4/12);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGA  >  minE/2848485‑2848613
                                                                |                                                                
tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                                               >  1:360464/1‑68 (MQ=255)
tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                                               >  1:1186792/1‑68 (MQ=255)
tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                                               >  1:3409575/1‑68 (MQ=255)
tGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGTACCGCACCGAGGGTAATCGGc                                                               >  1:1002674/1‑68 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:2492933/69‑1 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:3196451/69‑1 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:807179/69‑1 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:2509663/69‑1 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:3209838/69‑1 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:651918/69‑1 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:1798601/69‑1 (MQ=255)
     aCCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGa                                                         <  1:2755/69‑1 (MQ=255)
       cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                                       <  1:105877/69‑1 (MQ=255)
       cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                                       <  1:2614914/69‑1 (MQ=255)
       cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                                       <  1:567886/69‑1 (MQ=255)
       cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                                       <  1:647648/69‑1 (MQ=255)
       cAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                                       <  1:109800/69‑1 (MQ=255)
        aGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACt                                                       <  1:1411105/68‑1 (MQ=255)
                   gatgatAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTg                                           >  1:196470/1‑69 (MQ=255)
                       atAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTAcc                                       >  1:3224917/1‑69 (MQ=255)
                       atAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTAcc                                       >  1:789629/1‑69 (MQ=255)
                         aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                     <  1:2839088/69‑1 (MQ=255)
                         aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                     <  1:2145824/69‑1 (MQ=255)
                         aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                     <  1:2250571/69‑1 (MQ=255)
                         aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                     <  1:812738/69‑1 (MQ=255)
                         aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                     <  1:448923/69‑1 (MQ=255)
                         aaaCACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGt                                     <  1:225535/69‑1 (MQ=255)
                                 gAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGa                             <  1:3175436/69‑1 (MQ=255)
                                    tGGCTATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:483355/43‑1 (MQ=25)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:1609684/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:3086932/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:1347142/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:62344/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:3355190/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:50455/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:2584699/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:1939040/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:391969/43‑1 (MQ=37)
                                    tGGCAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:2287889/43‑1 (MQ=37)
                                       cAATCAGTACCGCACCGAGGGTAATCGGCTCCAGGCTTTc                                                    <  1:1271763/40‑1 (MQ=37)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:3222526/1‑68 (MQ=255)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:2641474/1‑68 (MQ=255)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:1290977/1‑68 (MQ=255)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:1335482/1‑68 (MQ=255)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:1400146/1‑68 (MQ=255)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:1160019/1‑68 (MQ=255)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:1404143/1‑68 (MQ=255)
                                           cAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCcaca                    >  1:1878512/1‑68 (MQ=255)
                                            aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTTACa                           <  1:1840017/60‑1 (MQ=255)
                                            aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa                           <  1:117581/60‑1 (MQ=255)
                                            aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa                           <  1:1527349/60‑1 (MQ=255)
                                            aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa                           <  1:2338328/60‑1 (MQ=255)
                                            aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa                           <  1:2640119/60‑1 (MQ=255)
                                            aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACa                           <  1:897231/60‑1 (MQ=255)
                                            aGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACACCCTGTACCGTGGAGGTGACa                           <  1:2784814/60‑1 (MQ=255)
                                                  gcaccgAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAata            >  1:854512/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGATGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:2943028/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:2993764/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:965277/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:1214973/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:61820/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:2414044/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:2905367/1‑69 (MQ=255)
                                                            tAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGa  >  1:2968095/1‑69 (MQ=255)
                                                                |                                                                
TGCGCACCAGCTGCGTGGTGATGATAAACACCAGAATGGCAATCAGCACCGCACCGAGGGTAATTGGCTCCAGACTTTCTACGCCCTGTACCGTGGAGGTGACATCCCACAGCGAAATATTTTCGAGGA  >  minE/2848485‑2848613

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: