Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,145,600 A→G 36.4% L198L (TTG→CTG)  chbF ← cryptic phospho‑beta‑glucosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,145,6000AG36.4% 72.9 / 74.5 118L198L (TTG→CTG) chbFcryptic phospho‑beta‑glucosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (43/32);  new base G (16/27);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 5.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 118 total aligned reads displayed.

GACTTGCCATTTATCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCTT  >  minE/1145545‑1145657
                                                       |                                                         
gACTTGCCATTTATCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATaa                                               <  1:3038652/68‑1 (MQ=255)
     gCCATTTATCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAatcatc                                         >  1:3112928/1‑69 (MQ=255)
     gCCATTTATCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAatcatc                                         >  1:1450473/1‑69 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:1005482/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:2128905/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:2134308/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:2323590/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:1476714/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:1418674/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:447004/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:631175/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:1223246/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:931717/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:1141731/1‑64 (MQ=255)
        aTTTACCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATAGATAGATAAatca                                           >  1:2405779/1‑64 (MQ=255)
        aTTTACCAGCAAATCCTTAATGAACACCATATGGTTGAGGCCGAACAGATCGATAGATAAatca                                           >  1:2980166/1‑64 (MQ=255)
           tATCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTc                                   >  1:2124396/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       <  1:791823/69‑1 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       <  1:981392/69‑1 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:2151287/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:1850601/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       <  1:1090864/69‑1 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:2871965/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:3288046/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:1218547/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:915947/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:511740/1‑69 (MQ=255)
                       cTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAg                       >  1:692998/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCTATAGATAAATCATCGCTGTCTTTCAGCATCAGAACAt                  >  1:2604486/1‑68 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:3001162/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:2969998/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:959334/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:25807/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:2487383/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:2300842/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:226193/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:330201/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:1025238/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:1597387/1‑69 (MQ=255)
                             gAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATc                 >  1:1120925/1‑69 (MQ=255)
                                cacCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATcgcg              >  1:2193954/1‑69 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:1038684/69‑1 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:157238/69‑1 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2923797/69‑1 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:3186662/69‑1 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:400858/69‑1 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:1366728/69‑1 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2047759/69‑1 (MQ=255)
                                    atatGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCAGTCTTTCAGCATCAGAACATCGCGAATa          <  1:1363351/69‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:1463985/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:146945/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:866956/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2704993/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2059867/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2575300/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2546943/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2072398/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:221344/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2298552/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:1712008/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATa          <  1:2002931/66‑1 (MQ=255)
                                       tGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAAc       >  1:2595445/1‑69 (MQ=255)
                                       tGGTTGACGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAAc       >  1:2410565/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:2587804/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:2274187/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:2728787/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:2731297/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:517703/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:2312305/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:1566496/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:850138/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACAt     >  1:126214/1‑69 (MQ=255)
                                         gTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATAGCGAATAAACAt     >  1:1694267/1‑69 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1930496/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:723570/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:703093/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1315332/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1345733/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:557160/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:972065/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1239124/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:3323427/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:442594/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1603815/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:267303/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:2506508/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:14441/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1458668/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:3236187/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:2366423/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1761112/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1496800/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1512563/69‑1 (MQ=255)
                                            gAGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAACATCAGAACATCGCGAATAAACATCtt  <  1:1215725/69‑1 (MQ=255)
                                            gAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:772306/69‑1 (MQ=255)
                                            gAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  >  1:2292366/1‑69 (MQ=255)
                                            gAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  >  1:518977/1‑69 (MQ=255)
                                             aGGCCGAACAGATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:3325237/68‑1 (MQ=255)
                                             aGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1663070/68‑1 (MQ=255)
                                               gCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCtt  <  1:1796147/66‑1 (MQ=255)
                                                       |                                                         
GACTTGCCATTTATCAGCACATCCTTAATGAACACCATATGGTTGAGGCCGAACAAATCGATAGATAAATCATCGCTGTCTTTCAGCATCAGAACATCGCGAATAAACATCTT  >  minE/1145545‑1145657

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: