Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 643,487 A→C 13.7% L65L (CTT→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,4870AC13.7% 162.0 / 19.7 68L65L (CTT→CTGompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (33/23);  new base C (3/6);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATG  >  minE/643427‑643555
                                                            |                                                                    
atatTCCCACTGACCATAACCGGTCAGATAGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCa                                                             >  1:2092347/1‑70 (MQ=255)
  atTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACg                                                                 >  1:3054175/1‑64 (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)
         cTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCaa                                                                     <  1:680037/53‑1 (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)
                  aCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCaa                                                                     >  1:2854609/1‑44 (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: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:39490/70‑1 (MQ=255)
                                          aGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa                   <  1:5263873/70‑1 (MQ=255)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                                                 <  1:979285/39‑1 (MQ=9)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg                                                 <  1:335279/39‑1 (MQ=9)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                           <  1:670260/45‑3 (MQ=2)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                           <  1:4325777/45‑3 (MQ=2)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                           <  1:4665175/45‑3 (MQ=2)
                                           gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                           <  1:4658526/45‑3 (MQ=2)
                                            tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                           >  1:5123028/1‑42 (MQ=2)
                                            tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                           >  1:1392188/1‑42 (MQ=2)
                                            tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                           >  1:4732126/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:1441848/54‑1 (MQ=255)
                                                   tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTc                          <  1:754256/54‑1 (MQ=255)
                                                    ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                             >  1:4361800/1‑50 (MQ=255)
                                                    ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                             >  1:4355507/1‑50 (MQ=255)
                                                    ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt                             >  1:5041755/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:2859589/1‑43 (MQ=255)
                                                      aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                  >  1:2371831/1‑43 (MQ=255)
                                                      aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                  >  1:1926999/1‑43 (MQ=255)
                                                       aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAaat      <  1:2786438/70‑1 (MQ=255)
                                                          ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg  <  1:3853084/71‑1 (MQ=255)
                                                          ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg  >  1:1483044/1‑71 (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)
                                                            |                                                                    
ATATTCCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATG  >  minE/643427‑643555

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: