Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 643,493 G→C 11.4% A63A (GCC→GCG 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,4930GC11.4% 225.4 / 14.9 79A63A (GCC→GCGompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base G (42/28);  new base C (3/6);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

ATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGAC  >  minE/643427‑643560
                                                                  |                                                                   
atatTCCCACTGACCATAACCGGTCAGATAGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCa                                                                  >  1:2092347/1‑70 (MQ=255)
   ttCCCACTGACCATAACCAGTGAATTAATCAGCGATCTGAGTTTCACCTTTGAAACCGATACGGGCATatg                                                              <  1:1624665/71‑3 (MQ=2)
     cccACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGt                                                             <  1:2949791/70‑1 (MQ=255)
      ccACTGACCATAACCGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATa                                                                <  1:1896485/65‑1 (MQ=255)
        aCTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg                                                         >  1:4718141/1‑71 (MQ=255)
        aCTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg                                                         >  1:4164276/1‑71 (MQ=255)
         cTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCAt                                                                 <  1:888056/62‑1 (MQ=255)
              cATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCa                                                                  >  1:3595559/1‑56 (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:5361736/1‑71 (MQ=255)
                 aaCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTg                                                >  1:236655/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)
                             cGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCa                                                           >  1:114233/1‑48 (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:1169998/1‑71 (MQ=255)
                                 aTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                <  1:4679742/71‑1 (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:335279/39‑1 (MQ=9)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                                                      <  1:979285/39‑1 (MQ=9)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                <  1:670260/45‑3 (MQ=2)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                <  1:4658526/45‑3 (MQ=2)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                <  1:4665175/45‑3 (MQ=2)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                <  1:4325777/45‑3 (MQ=2)
                                            tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                >  1:4732126/1‑42 (MQ=2)
                                            tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                >  1:1392188/1‑42 (MQ=2)
                                            tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                >  1:5123028/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:5041755/1‑50 (MQ=255)
                                                    ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                  >  1:4361800/1‑50 (MQ=255)
                                                    ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                  >  1:4355507/1‑50 (MQ=255)
                                                    ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                                  >  1:3584994/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:2371831/1‑43 (MQ=255)
                                                      aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                       >  1:2859589/1‑43 (MQ=255)
                                                      aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                       >  1:1926999/1‑43 (MQ=255)
                                                       aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAaat           <  1:2786438/70‑1 (MQ=255)
                                                          ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg       >  1:1483044/1‑71 (MQ=255)
                                                          ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg       <  1:3853084/71‑1 (MQ=255)
                                                          ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg       <  1:1073922/71‑1 (MQ=255)
                                                          ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg       <  1:4367416/71‑1 (MQ=255)
                                                            aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaaaa            >  1:1557125/1‑64 (MQ=255)
                                                            aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACACTTGGaaaaa            >  1:1801515/1‑64 (MQ=255)
                                                             aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                       >  1:5147641/1‑36 (MQ=255)
                                                             aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                       >  1:4179006/1‑36 (MQ=255)
                                                             aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                       >  1:3435306/1‑36 (MQ=255)
                                                             aGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCa     <  1:3336347/70‑1 (MQ=255)
                                                              gACGGGCATAGGTCATGTCACCATTGCCACCGTAACTGTTTTCACCGt                          <  1:3296958/48‑1 (MQ=255)
                                                               aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                >  1:1408658/1‑41 (MQ=255)
                                                               aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                >  1:5432919/1‑41 (MQ=255)
                                                               aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                >  1:5410093/1‑41 (MQ=255)
                                                               aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt                                >  1:4928621/1‑41 (MQ=255)
                                                                cGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGAc  >  1:2221897/1‑70 (MQ=255)
                                                                  |                                                                   
ATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGAC  >  minE/643427‑643560

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: