Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 643,503 A→T 10.2% M60K (ATG→AAG) ‡ ompF ← outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported643,5030AT10.2% 240.0 / 16.6 88M60K (ATG→AAG) ‡ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (51/28);  new base T (3/6);  total (54/34)
Fisher's exact test for biased strand distribution p-value = 8.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

ACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACC  >  Exported/643435‑643573
                                                                    |                                                                      
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)
                     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: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:4658526/45‑3 (MQ=2)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                             <  1:4665175/45‑3 (MQ=2)
                                   gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct                                                             <  1:4325777/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:2371831/1‑43 (MQ=255)
                                              aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                    >  1:1926999/1‑43 (MQ=255)
                                              aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                    >  1:2859589/1‑43 (MQ=255)
                                               aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAaat                        <  1:2786438/70‑1 (MQ=255)
                                                  ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                    <  1:1073922/71‑1 (MQ=255)
                                                  ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                    >  1:1483044/1‑71 (MQ=255)
                                                  ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                    <  1:3853084/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:3435306/1‑36 (MQ=255)
                                                     aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                    >  1:5147641/1‑36 (MQ=255)
                                                     aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa                                                    >  1:4179006/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:1205247/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt      >  1:3284373/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt      >  1:5667624/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt      <  1:5544781/71‑1 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt      >  1:5531783/1‑71 (MQ=255)
                                                                ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt      >  1:1396457/1‑71 (MQ=255)
                                                                 gTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                    <  1:3724021/56‑1 (MQ=255)
                                                                 gTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg                    <  1:2587379/56‑1 (MQ=255)
                                                                   cATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTAc   >  1:4122511/1‑71 (MQ=255)
                                                                    aTGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTAcc  >  1:648151/1‑71 (MQ=255)
                                                                    |                                                                      
ACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTTTACC  >  Exported/643435‑643573

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: