Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 643,235 G→A 22.6% D149D (GAC→GAT 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,2350GA22.6% 158.8 / 34.6 114D149D (GAC→GATompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base G (24/64);  new base A (18/8);  total (42/72)
Fisher's exact test for biased strand distribution p-value = 1.66e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 124 total aligned reads displayed.

ATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAC  >  minE/643169‑643303
                                                                  |                                                                    
aTCAACCAGACCAAAGATGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:1515352/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:957303/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:1767242/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:3139302/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:3140775/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:1620329/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:1907168/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:3162554/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:1458551/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:3225458/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:3229467/69‑1 (MQ=255)
aTCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTc                                                                    <  1:2240601/69‑1 (MQ=255)
  ccaCCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCAt                                                                  <  1:3372394/67‑1 (MQ=255)
  cAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCAt                                                                  <  1:1765231/69‑1 (MQ=255)
  cAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCAt                                                                  <  1:1264841/69‑1 (MQ=255)
     ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAATCATCGc                                                               <  1:1505310/69‑1 (MQ=255)
     ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAATCATCGc                                                               <  1:658382/69‑1 (MQ=255)
     ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAATCATCGc                                                               <  1:1703154/69‑1 (MQ=255)
     ccAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAATCATCGc                                                               <  1:1340104/69‑1 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGc                                                               >  1:2180665/1‑68 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGc                                                               >  1:1796134/1‑68 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAATCATCGc                                                               <  1:2058103/68‑1 (MQ=255)
      cAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAATCATCGc                                                               <  1:3413845/68‑1 (MQ=255)
       aGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTg                                                             <  1:3302409/69‑1 (MQ=255)
              agaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:3196370/69‑1 (MQ=255)
              agaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:1719640/69‑1 (MQ=255)
              agaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:384384/69‑1 (MQ=255)
              agaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:805092/69‑1 (MQ=255)
              agaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:1583784/69‑1 (MQ=255)
              agaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:1467254/69‑1 (MQ=255)
               gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:282003/68‑1 (MQ=255)
               gaagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg                                                      <  1:353210/68‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:44512/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:3379449/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:18991/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:1591268/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:2795650/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:799334/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:1041511/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:3111775/56‑1 (MQ=255)
                aagTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATc                                                                 <  1:2322069/56‑1 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCAt                                                                  >  1:872173/1‑54 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:2250611/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:908709/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:2696451/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:117640/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:2563540/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:1401063/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:1404393/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:1428926/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:607968/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:55121/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:532162/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:461820/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:1948494/1‑69 (MQ=255)
                 agTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTAt                                                   >  1:2960976/1‑69 (MQ=255)
                                  gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                  >  1:2948228/1‑69 (MQ=255)
                                  gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgc                                  >  1:1377944/1‑69 (MQ=255)
                                  gTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCAACACCAAATTCTGgc                                  >  1:1681201/1‑69 (MQ=255)
                                       aaCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                              <  1:2625664/68‑1 (MQ=255)
                                       aaCGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGgcagca                              <  1:1180211/68‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:10462/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:3350152/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:1584962/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:2349572/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:1574455/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:263293/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:754084/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:792235/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:3039462/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          <  1:1235895/69‑1 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTAt                                                   >  1:688714/1‑44 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCatat                           >  1:3062957/1‑68 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:1366753/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:1537152/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:1277152/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:356647/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:184460/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:2039228/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:3253989/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:3205382/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:269146/1‑69 (MQ=255)
                                          gccgccAACACGACCAACGAAGAAATCATCGCTATATGCAGTATCACCACCAAATTCTGGCAGCATATc                          >  1:1754285/1‑69 (MQ=255)
                                                        cAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg            <  1:2965326/69‑1 (MQ=255)
                                                        cAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTg            <  1:3042991/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:2588326/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:2908484/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:933465/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:1074880/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:2845461/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:1636460/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:1675607/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:2936251/69‑1 (MQ=255)
                                                          aCGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGca          <  1:188649/69‑1 (MQ=255)
                                                                  gTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc  <  1:2662793/69‑1 (MQ=255)
                                                                  gTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc  <  1:1989608/69‑1 (MQ=255)
                                                                  gTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc  <  1:780268/69‑1 (MQ=255)
                                                                  gTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc  <  1:2396980/69‑1 (MQ=255)
                                                                  gTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAc  <  1:3013060/69‑1 (MQ=255)
                                                                  gTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCACAAAc  <  1:1441523/69‑1 (MQ=255)
                                                                  |                                                                    
ATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCACCAAATTCTGGCAGCATATCGGTGTAACCCAGTGCATCATAAAC  >  minE/643169‑643303

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