Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,484,871 C→A 46.4% S459S (TCG→TCT yidC ← cytoplasmic insertase into membrane protein, Sec system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,484,8710CA46.4% 33.4 / 100.0 97S459S (TCG→TCTyidCcytoplasmic insertase into membrane protein, Sec system
Reads supporting (aligned to +/- strand):  ref base C (20/32);  new base A (20/25);  total (40/57)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 112 total aligned reads displayed.

GACATCTTCTGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATGTAG  >  minE/2484803‑2484937
                                                                    |                                                                  
gACATCTTCTGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGcc                                                                    <  1:3037365/69‑1 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAgg                                                              >  1:2343144/1‑69 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAgg                                                              >  1:1043051/1‑69 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAgg                                                              >  1:2056907/1‑69 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAgg                                                              >  1:758063/1‑69 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAgg                                                              >  1:530042/1‑69 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAgg                                                              >  1:588374/1‑69 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAgg                                                              >  1:2908874/1‑69 (MQ=255)
      ttctGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAg                                                               >  1:2048259/1‑68 (MQ=255)
            atgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTg                                                         <  1:904992/68‑1 (MQ=255)
             tgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTg                                                         <  1:2199643/67‑1 (MQ=255)
             tgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTg                                                         <  1:1058298/67‑1 (MQ=255)
             tgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTg                                                         <  1:710190/67‑1 (MQ=255)
             tgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGa                                                       <  1:2906798/69‑1 (MQ=255)
             tgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGa                                                       <  1:1324521/69‑1 (MQ=255)
             tgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGa                                                       <  1:2950907/69‑1 (MQ=255)
             tgaaGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGa                                                       <  1:648443/69‑1 (MQ=255)
                    cATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCcac                                                  <  1:1068274/67‑1 (MQ=255)
                    cATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCcac                                                  <  1:3271105/67‑1 (MQ=255)
                    cATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCcac                                                  <  1:3306336/67‑1 (MQ=255)
                     aTCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCcac                                                  <  1:1313184/66‑1 (MQ=255)
                           aCGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAAc                                         >  1:1468191/1‑69 (MQ=255)
                           aCGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAAc                                         >  1:199179/1‑69 (MQ=255)
                           aCGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAAc                                         <  1:2292373/69‑1 (MQ=255)
                           aCGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAAc                                         >  1:364074/1‑69 (MQ=255)
                           aCGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAAc                                         <  1:1020553/69‑1 (MQ=255)
                           aCGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAAc                                         >  1:734401/1‑69 (MQ=255)
                           aCGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAAc                                         >  1:807228/1‑69 (MQ=255)
                                cATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGc                                    >  1:2363354/1‑69 (MQ=255)
                                cATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGc                                    >  1:1656009/1‑69 (MQ=255)
                                cATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGc                                    >  1:2301208/1‑69 (MQ=255)
                                cATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGc                                    >  1:2493985/1‑69 (MQ=255)
                                cATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGc                                    >  1:1291004/1‑69 (MQ=255)
                                cATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGc                                    >  1:1220624/1‑69 (MQ=255)
                                cATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGc                                    >  1:2531089/1‑69 (MQ=255)
                                  tCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCcac                                                  <  1:1028559/53‑1 (MQ=255)
                                  tCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCt                                  >  1:1793997/1‑69 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:2693118/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:904504/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:822820/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:2574834/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1003589/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:624578/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:295587/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:493586/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:3047354/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:40060/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:3391429/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:235175/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:2235911/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1062755/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1436066/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1304909/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1220978/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1735860/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1241298/69‑1 (MQ=255)
                                        tCGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1240092/69‑1 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:3424567/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:3344714/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:331383/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:166474/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:3194362/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:1800319/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:2955605/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:2845617/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:105476/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:649475/1‑39 (MQ=255)
                                         cGGCAGGATGTAGTACGGGTCCTGTGCAGACAGGTCGTg                                                         >  1:1996416/1‑39 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACTCAGTTc                        <  1:2352907/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:1294373/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:1691061/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:1179541/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:1971791/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:1278840/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:1986803/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:3118172/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:1766235/69‑1 (MQ=255)
                                            cAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:278702/69‑1 (MQ=255)
                                             aGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:878024/68‑1 (MQ=255)
                                                 tGTAGTACGGGTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:3301278/60‑1 (MQ=255)
                                                          ggTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATc           >  1:1253739/1‑68 (MQ=255)
                                                          ggTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATc           >  1:1000049/1‑68 (MQ=255)
                                                          ggTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATc           >  1:1862961/1‑68 (MQ=255)
                                                          ggTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATc           >  1:2738099/1‑68 (MQ=255)
                                                           gTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTc                        <  1:3252837/54‑1 (MQ=255)
                                                           gTCCTGTGCAGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCa                            <  1:1998421/50‑1 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCAt      >  1:12199/1‑64 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:2409461/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:1616671/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  <  1:2414387/68‑1 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:2426526/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:2637574/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:1222404/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  <  1:1834190/68‑1 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:3174852/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:3244695/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  <  1:3249769/68‑1 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  >  1:3269323/1‑68 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATgtag  <  1:3388542/68‑1 (MQ=255)
                                                                   ccGACAGGTCGTGGATCCACAGTGCAAACGATGCCTGACGCAGTTCAACGGAACCCATCAGCATgta   >  1:2716203/1‑67 (MQ=255)
                                                                    |                                                                  
GACATCTTCTGAATGAAGAACATCGTTACGCCCATCAGGATCGGCAGGATGTAGTACGGGTCCTGTGCCGACAGGTCGTGGATCCACAGTGCAAACGGTGCCTGACGCAGTTCAACGGAACCCATCAGCATGTAG  >  minE/2484803‑2484937

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