Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 189,568 G→A 24.8% A355A (GCG→GCA ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE189,5680GA24.8% 134.9 / 37.4 108A355A (GCG→GCAldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base A (15/12);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 8.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

Only 100 of 109 total aligned reads displayed.

CACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTG  >  minE/189511‑189634
                                                         |                                                                  
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCgg                                                                   >  1:2208962/1‑59 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:3321566/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:2464134/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:2468456/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:2573134/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:2631389/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:2939980/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:3065936/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:3241961/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:3263124/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:364129/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:54795/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:624837/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:1332358/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:869181/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:1394297/1‑69 (MQ=255)
cacCCATTTTCATCCGATCTAACAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:1086996/1‑69 (MQ=255)
cacCCATTTTCATCCGATATACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGAt                                                         >  1:834374/1‑69 (MQ=255)
     catttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCg                                                    <  1:2063078/67‑1 (MQ=255)
     aTTTTCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCg                                                    <  1:1658421/69‑1 (MQ=255)
     aTTTTCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCg                                                    <  1:2494850/69‑1 (MQ=255)
     aTTTTCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCg                                                    <  1:3198025/69‑1 (MQ=255)
     aTTTTCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCg                                                    <  1:2784503/69‑1 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:776662/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:940438/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:253841/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:882276/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:260693/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:868209/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:2996070/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:3147853/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:72894/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:1081765/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:336837/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:639226/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:62321/1‑69 (MQ=255)
      ttttCATCCAATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCAGGAAAAGTGATCTTCGa                                                   >  1:1004670/1‑69 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTg                                                           <  1:83950/59‑1 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTg                                                           <  1:1581250/59‑1 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTg                                                           <  1:422737/59‑1 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTg                                                           <  1:417747/59‑1 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTg                                                           <  1:3281932/59‑1 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTg                                                           <  1:2678313/59‑1 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGaaa                                                 >  1:23135/1‑69 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGaaa                                                 >  1:3016047/1‑69 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTCCGGGAAAAGTg                                                           <  1:118681/59‑1 (MQ=255)
        ttCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGAGAGCGTGTTGCGGGAAAAGTg                                                           <  1:3278756/59‑1 (MQ=255)
                 tCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGaaa                                                               >  1:3144112/1‑46 (MQ=255)
                 tCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGaaa                                                               >  1:553503/1‑46 (MQ=255)
                 tCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGaaa                                                               >  1:1785030/1‑46 (MQ=255)
                 tCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGaaa                                                               >  1:232895/1‑46 (MQ=255)
                       aGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCa                                             <  1:242336/58‑1 (MQ=255)
                       aGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCa                                             <  1:2602170/58‑1 (MQ=255)
                       aGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCa                                             <  1:1089521/58‑1 (MQ=255)
                       aGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCa                                             <  1:2069847/58‑1 (MQ=255)
                       aGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCcac                                   >  1:2950196/1‑68 (MQ=255)
                        gggTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaa                                 <  1:1206077/69‑1 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:1662716/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:712502/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:2315135/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:2556654/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:1796988/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:1970603/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:585194/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:1249848/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaaa                               >  1:2711653/1‑69 (MQ=255)
                          gTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACaaa                                >  1:2103265/1‑68 (MQ=255)
                                        aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                      <  1:2115927/48‑1 (MQ=255)
                                        aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                      <  1:3193605/48‑1 (MQ=255)
                                        aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                      <  1:2933852/48‑1 (MQ=255)
                                        aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                      <  1:834198/48‑1 (MQ=255)
                                        aGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGAcc                                      <  1:2366415/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:1338841/69‑1 (MQ=255)
                                           ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:109310/69‑1 (MQ=255)
                                           ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2295949/69‑1 (MQ=255)
                                           ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2221391/69‑1 (MQ=255)
                                           ggcgAGCGTGTTGCAGGAAAAGTGATCTTCGAAACGCAGTCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2190002/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:2319434/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:758558/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2335422/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:1767990/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:3306359/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2263803/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2860337/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:214393/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2201670/57‑1 (MQ=255)
                                                       gCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCg              <  1:2660847/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)
                                                         |                                                                  
CACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCGTGTTGCGGGAAAAGTGATCTTCGAAACGCAATCGACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTG  >  minE/189511‑189634

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: