Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,290,935 A→C 20.0% G549G (GGA→GGC ppc → phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,290,9350AC20.0% 138.1 / 33.6 85G549G (GGA→GGCppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (30/38);  new base C (2/15);  total (32/53)
Fisher's exact test for biased strand distribution p-value = 2.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCGGG  >  minE/2290869‑2291003
                                                                   |                                                                    
tCGTGGCCTGATTCAGGGCAATCAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:630629/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:1058671/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:2113569/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:2452272/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:2961514/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:693089/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:1234691/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:2859504/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:3255397/69‑1 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    >  1:1933733/1‑69 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    >  1:3339698/1‑69 (MQ=255)
tCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGt                                                                    <  1:3413625/69‑1 (MQ=255)
  gTGGCCTTATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTg                                                                   <  1:488932/68‑1 (MQ=255)
  gTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTg                                                                   <  1:2914875/68‑1 (MQ=255)
  gTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTg                                                                   <  1:1590965/68‑1 (MQ=255)
  gTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTg                                                                   <  1:1573422/68‑1 (MQ=255)
  gTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGa                                                                  >  1:1491142/1‑69 (MQ=255)
  gTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGa                                                                  >  1:1616263/1‑69 (MQ=255)
  gTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGa                                                                  >  1:1201861/1‑69 (MQ=255)
  gTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGa                                                                  >  1:1710608/1‑69 (MQ=255)
  gTGGCCTGATTAAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGa                                                                  >  1:2455581/1‑69 (MQ=255)
            tCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTGAGCAAAAGATGCGGGCGTGATGGCAGCTTc                                                        >  1:2893279/1‑69 (MQ=255)
            tCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGCGTGATGGCAGCTTc                                                        >  1:1357342/1‑69 (MQ=255)
                    aaCAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGgcgc                                                >  1:2776545/1‑69 (MQ=255)
                      cAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCaa                                              >  1:2850923/1‑69 (MQ=255)
                      cAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCaa                                              >  1:1300005/1‑69 (MQ=255)
                      cAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCa                                               >  1:1735414/1‑68 (MQ=255)
                      cAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCa                                               >  1:2973762/1‑68 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTc                                                        >  1:3152908/1‑47 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGcaca                                   >  1:2267085/1‑68 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  >  1:1578803/1‑69 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  <  1:1759328/69‑1 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  <  1:1799534/69‑1 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  <  1:44522/69‑1 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  >  1:541765/1‑69 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  <  1:2665561/69‑1 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  <  1:607732/69‑1 (MQ=255)
                                  atTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  <  1:891391/69‑1 (MQ=255)
                                   ttGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAg                                  <  1:759187/68‑1 (MQ=255)
                                    tGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGTCAGCTTCCTGGGCGCAATATCAGGCACAGGa                                >  1:1955730/1‑69 (MQ=255)
                                    tGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGca                                     <  1:2641209/64‑1 (MQ=255)
                                    tGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAgg                                 >  1:399671/1‑68 (MQ=255)
                                    tGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGa                                >  1:1089568/1‑69 (MQ=255)
                                    tGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGa                                >  1:1468879/1‑69 (MQ=255)
                                    tGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGa                                >  1:1802725/1‑69 (MQ=255)
                                       cTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATg                              <  1:902941/68‑1 (MQ=255)
                                       cTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATg                              <  1:508499/68‑1 (MQ=255)
                                       cTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATg                              <  1:2054437/68‑1 (MQ=255)
                                       cTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATg                              <  1:3433769/68‑1 (MQ=255)
                                       cTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATg                              <  1:3282842/68‑1 (MQ=255)
                                       cTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATg                              <  1:3217457/68‑1 (MQ=255)
                                       cTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATg                              <  1:1696148/68‑1 (MQ=255)
                                         aatCCGACTCAGCAAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCa                            <  1:456184/67‑1 (MQ=255)
                                         aTTCCGACTCAGCAAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCa                            <  1:2278713/69‑1 (MQ=255)
                                         aTTCCGACTCAGCAAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCa                            <  1:506038/69‑1 (MQ=255)
                                         aTTCCGACTCAGCAAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCa                            <  1:378089/69‑1 (MQ=255)
                                         aTTCCGACTCAGCAAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCa                            <  1:1723156/69‑1 (MQ=255)
                                         aTTCCGACTCAGCAAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCa                            <  1:2038451/69‑1 (MQ=255)
                                         aTTCCGACTCAGCAAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCa                            <  1:2709150/69‑1 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:2897210/1‑68 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:1717884/1‑68 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:2164818/1‑68 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:1202429/1‑68 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:1291402/1‑68 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:1366470/1‑68 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:1473509/1‑68 (MQ=255)
                                            ccGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTa                         >  1:3167680/1‑68 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCTTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:3140394/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:91838/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:2099077/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:2653157/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:1405220/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:1206136/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:1439194/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:3278132/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:3297945/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:3057532/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:1577548/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:2976431/68‑1 (MQ=255)
                                               aCTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:2349855/68‑1 (MQ=255)
                                                cTCAGCAAAAGATGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:2570749/67‑1 (MQ=255)
                                                 tCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCaaa                   >  1:3268747/1‑69 (MQ=255)
                                                             tGCGGGCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:584791/55‑1 (MQ=255)
                                                                gggCGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATc                      <  1:1381952/52‑1 (MQ=255)
                                                                gggAGTGATGGCAGCTTCCTGGGCGCAATATCAGGca                                     <  1:1125711/37‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:1797601/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:2964054/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:901193/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:1117849/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:1217788/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:2983189/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:514487/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:1315135/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:2065532/69‑1 (MQ=255)
                                                                   aGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCggg  <  1:2021812/69‑1 (MQ=255)
                                                                   |                                                                    
TCGTGGCCTGATTCAGGGCAAACAGATGGTGATGATTGGCTATTCCGACTCAGCAAAAGATGCGGGAGTGATGGCAGCTTCCTGGGCGCAATATCAGGCACAGGATGCATTAATCAAAACCTGCGAAAAAGCGGG  >  minE/2290869‑2291003

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: