Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,592 A→G 40.8% Q363Q (CAA→CAG ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,5920AG40.8% 27.5 / 53.6 76Q363Q (CAA→CAGldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base A (15/30);  new base G (15/16);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAA  >  minE/189534‑189645
                                                          |                                                     
aGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCcac                                              >  1:2278314/1‑68 (MQ=255)
aGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCcac                                              >  1:2950196/1‑68 (MQ=255)
 gggTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaa                                            <  1:1206077/69‑1 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:1970603/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:2556654/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:2711653/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:1796988/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:1662716/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:2533226/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:585194/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:2315135/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:1249848/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                                          >  1:712502/1‑69 (MQ=255)
   gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaa                                           >  1:2103265/1‑68 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:2115927/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:2366415/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:2041397/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:834198/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:2933852/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:3193605/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:1464669/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                                 <  1:701176/48‑1 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTa                            >  1:574182/1‑69 (MQ=255)
                 aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTa                            >  1:1347004/1‑69 (MQ=255)
                    ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2295949/69‑1 (MQ=255)
                    ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:109310/69‑1 (MQ=255)
                    ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2221391/69‑1 (MQ=255)
                    ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2190002/69‑1 (MQ=255)
                    ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:1338841/69‑1 (MQ=255)
                     gcgAGCGTGTTGCAGGAAAAGTGATCTTCGTAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATct                         <  1:1281597/68‑2 (MQ=255)
                     gcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2710542/68‑1 (MQ=255)
                                gCGGGAAATGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:3239273/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2335422/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2860337/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2660847/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:1767990/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:214393/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2319434/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:3306359/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2263803/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:758558/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:2201670/57‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg             <  1:951077/69‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg             <  1:1707793/69‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg             <  1:3082972/69‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg             <  1:2872832/69‑1 (MQ=255)
                                gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTg             <  1:2247072/69‑1 (MQ=255)
                                 cGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:1951951/56‑1 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATc                          >  1:346900/2‑54 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:2230900/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:2192034/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:63278/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:536430/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:357579/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:3309554/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:2366258/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:2369701/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:3182008/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:3084641/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:2527156/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:1711881/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:1823945/2‑68 (MQ=255)
                                  aggAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGa            >  1:2907357/2‑68 (MQ=255)
                                     aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCca        <  1:1480924/69‑1 (MQ=255)
                                     aaaaGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCca        <  1:214569/69‑1 (MQ=255)
                                         gTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGAGGCTTCGCTGATCCACATTa   <  1:2105721/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCTCTGATCCACATTaa  <  1:1236701/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:2843994/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:858484/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:1076258/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:2914234/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:550792/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:35972/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:2319149/69‑1 (MQ=255)
                                           gATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCGTCGCTGATCCACATTaa  <  1:3085861/69‑1 (MQ=255)
                                              cTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg                         <  1:65019/43‑1 (MQ=255)
                                                          |                                                     
AGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAA  >  minE/189534‑189645

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