Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 643,504 T→G 10.2% M60L (ATG→CTG) ‡ ompF ← outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE643,5040TG10.2% 238.5 / 18.5 88M60L (ATG→CTG) ‡ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base T (52/27);  new base G (3/6);  total (55/33)
Fisher's exact test for biased strand distribution p-value = 7.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCG  >  minE/643435‑643574
                                                                     |                                                                      
aCTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg                                                                       >  1:4718141/1‑71 (MQ=255)
aCTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg                                                                       >  1:4164276/1‑71 (MQ=255)
        tAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc                                                                     >  1:3717188/1‑65 (MQ=255)
        tAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc                                                                     >  1:844067/1‑65 (MQ=255)
         aaCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCa                                                                 >  1:2067888/1‑68 (MQ=255)
         aaCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTg                                                              >  1:236655/1‑71 (MQ=255)
         aaCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTg                                                              >  1:5361736/1‑71 (MQ=255)
          aCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGc                                                             <  1:5061059/71‑1 (MQ=255)
            cGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCa                                                           <  1:5320295/71‑1 (MQ=255)
             ggTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCAc                                                          >  1:4015265/1‑71 (MQ=255)
             ggTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCAc                                                          >  1:1943712/1‑71 (MQ=255)
                 aGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGt                                                       <  1:3138920/70‑1 (MQ=255)
                   aTCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCAt                                                                        <  1:4611718/51‑1 (MQ=255)
                    tCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCa                                                           >  1:2894089/1‑63 (MQ=255)
                     cGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc                                                                     >  1:1219940/1‑52 (MQ=255)
                       gAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                                <  1:2414017/71‑1 (MQ=255)
                        aaTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGttt                                               <  1:1447785/71‑1 (MQ=255)
                         aTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                              <  1:4679742/71‑1 (MQ=255)
                         aTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                              >  1:1169998/1‑71 (MQ=255)
                           tGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                     >  1:1289310/1‑62 (MQ=255)
                           tGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTc                                             >  1:4538312/1‑70 (MQ=255)
                           tGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCa                                            >  1:1560731/1‑71 (MQ=255)
                            gATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCAc                                           >  1:966891/1‑71 (MQ=255)
                                tGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGt                                        >  1:2483133/1‑70 (MQ=255)
                                tGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGt                                        <  1:2997788/70‑1 (MQ=255)
                                  aGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa                                      <  1:5263873/70‑1 (MQ=255)
                                  aGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa                                      <  1:39490/70‑1 (MQ=255)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                                                                    <  1:979285/39‑1 (MQ=9)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                                                                    <  1:335279/39‑1 (MQ=9)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                              <  1:4325777/45‑3 (MQ=2)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                              <  1:4658526/45‑3 (MQ=2)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                              <  1:4665175/45‑3 (MQ=2)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                              <  1:670260/45‑3 (MQ=2)
                                    tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                              >  1:4732126/1‑42 (MQ=2)
                                    tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                              >  1:5123028/1‑42 (MQ=2)
                                    tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                              >  1:1392188/1‑42 (MQ=2)
                                    tttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCAt                                                                <  1:742732/42‑1 (MQ=255)
                                     tgccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGt                                                       <  1:4086236/48‑1 (MQ=255)
                                          cTTTGAATCCGATACGGGCATATGACATGTCGCCAt                                                                <  1:797672/36‑1 (MQ=255)
                                          cTTTGAAACCGATACGGGCATATGACATGTCGCCAt                                                                <  1:5139789/36‑1 (MQ=16)
                                          cTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTg                                                  <  1:4465853/50‑1 (MQ=255)
                                          cTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTAcc                                    >  1:3924379/1‑64 (MQ=255)
                                           tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTg                                                  <  1:3289042/49‑1 (MQ=255)
                                           tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTc                                             <  1:754256/54‑1 (MQ=255)
                                           tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTc                                             <  1:1441848/54‑1 (MQ=255)
                                            ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                                >  1:3584994/1‑50 (MQ=255)
                                            ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                                >  1:5041755/1‑50 (MQ=255)
                                            ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                                >  1:4355507/1‑50 (MQ=255)
                                            ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                                >  1:4361800/1‑50 (MQ=255)
                                            ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaa                            >  1:1325758/1‑70 (MQ=255)
                                             tAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaaaa                          >  1:3530018/1‑71 (MQ=255)
                                             tAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaaa                           >  1:3906047/1‑70 (MQ=255)
                                              aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                     >  1:1926999/1‑43 (MQ=255)
                                              aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                     >  1:2371831/1‑43 (MQ=255)
                                              aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                     >  1:2859589/1‑43 (MQ=255)
                                               aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAaat                         <  1:2786438/70‑1 (MQ=255)
                                                  ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                     <  1:4367416/71‑1 (MQ=255)
                                                  ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                     <  1:3853084/71‑1 (MQ=255)
                                                  ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                     <  1:1073922/71‑1 (MQ=255)
                                                  ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                     >  1:1483044/1‑71 (MQ=255)
                                                    aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaaaa                          >  1:1557125/1‑64 (MQ=255)
                                                    aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACACTTGGaaaaa                          >  1:1801515/1‑64 (MQ=255)
                                                     aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                     >  1:4179006/1‑36 (MQ=255)
                                                     aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                     >  1:3435306/1‑36 (MQ=255)
                                                     aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                     >  1:5147641/1‑36 (MQ=255)
                                                     aGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCa                   <  1:3336347/70‑1 (MQ=255)
                                                      gACGGGCATAGGTCATGTCACCATTGCCACCGTAACTGTTTTCACCGt                                        <  1:3296958/48‑1 (MQ=255)
                                                       aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                              >  1:4928621/1‑41 (MQ=255)
                                                       aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                              >  1:1408658/1‑41 (MQ=255)
                                                       aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                              >  1:5432919/1‑41 (MQ=255)
                                                       aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                              >  1:5410093/1‑41 (MQ=255)
                                                        cGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGAc                >  1:2221897/1‑70 (MQ=255)
                                                          gggaTAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCaa             <  1:4682230/68‑1 (MQ=255)
                                                           gCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCaa             >  1:1352152/1‑70 (MQ=255)
                                                            cATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGa                 >  1:4410730/1‑65 (MQ=255)
                                                               aGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCt        >  1:505377/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa                                      >  1:4322509/1‑40 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa                                      >  1:3942761/1‑40 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa                                      >  1:2099510/1‑40 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt       >  1:3284373/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt       >  1:1396457/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt       >  1:5531783/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt       <  1:5544781/71‑1 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt       >  1:5667624/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt       >  1:1205247/1‑71 (MQ=255)
                                                                 gTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                     <  1:2587379/56‑1 (MQ=255)
                                                                 gTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                     <  1:3724021/56‑1 (MQ=255)
                                                                   cATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTAc    >  1:4122511/1‑71 (MQ=255)
                                                                    aTGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTAcc   >  1:648151/1‑71 (MQ=255)
                                                                     tGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCg  >  1:5200102/1‑71 (MQ=255)
                                                                     |                                                                      
ACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACCG  >  minE/643435‑643574

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: