Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,954 G→A 13.7% G116G (GGC→GGT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,9540GA13.7% 285.7 / 24.3 102G116G (GGC→GGTompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (47/41);  new base A (6/8);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 105 total aligned reads displayed.

GTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCC  >  minE/1439885‑1440024
                                                                     |                                                                      
gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACg                                                                        <  1:2857174/70‑1 (MQ=255)
gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                       >  1:5414727/1‑71 (MQ=255)
gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                       >  1:5161832/1‑71 (MQ=255)
gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                       >  1:77581/1‑71 (MQ=255)
gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                       >  1:5586634/1‑71 (MQ=255)
gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                       >  1:5192443/1‑71 (MQ=255)
   aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGt                                                                    <  1:3815406/71‑1 (MQ=255)
   aGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGt                                                                    <  1:4179201/71‑1 (MQ=255)
       ccGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGt                                                                 >  1:3274189/1‑70 (MQ=255)
         gTAGGTGTCGCCACCGAATTCTGGCAGAACGTCGGTCCAGGAAGTTACGTCGTAAACAACGCCGTAGTTAc                                                              >  1:5286585/1‑71 (MQ=255)
          tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGt                                                                    <  1:2183552/64‑1 (MQ=255)
          tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTAc                                                              <  1:5600500/70‑1 (MQ=255)
          tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACg                                                             >  1:485216/1‑71 (MQ=255)
          tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACg                                                             <  1:1385996/71‑1 (MQ=255)
               gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGtt                                                                >  1:2382272/1‑63 (MQ=255)
               gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCg                                                         >  1:1314500/1‑70 (MQ=255)
                  gccaccaAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGt                                                     >  1:1029332/2‑71 (MQ=2)
                  accaccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAg                                                      <  1:234209/70‑1 (MQ=255)
                    caccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTAc                                                              >  1:2317450/1‑60 (MQ=255)
                    caccAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTc                                                    >  1:3933879/1‑70 (MQ=2)
                     accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGa                                                            >  1:2560936/1‑61 (MQ=255)
                     accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGAcc                                                          >  1:5221540/1‑63 (MQ=255)
                     accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGAcc                                                          >  1:5682205/1‑63 (MQ=255)
                     accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                  <  1:1899483/71‑1 (MQ=255)
                     accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                  <  1:3548694/71‑1 (MQ=255)
                     accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCg                                                   >  1:2077768/1‑70 (MQ=2)
                     accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                  <  1:3865498/71‑1 (MQ=2)
                     accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                  <  1:2183233/71‑1 (MQ=2)
                     accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                  <  1:3858876/71‑1 (MQ=2)
                       caaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                  >  1:2805084/3‑69 (MQ=2)
                        gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                               >  1:1664598/1‑55 (MQ=255)
                        gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                               >  1:4681827/1‑55 (MQ=255)
                        gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                               >  1:2565344/1‑55 (MQ=255)
                        gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                               >  1:4875436/1‑55 (MQ=255)
                        gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                               >  1:3116911/1‑55 (MQ=255)
                        gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                  >  1:1224023/1‑68 (MQ=255)
                        gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                               <  1:4759651/71‑1 (MQ=255)
                        aaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                  <  1:1639409/67‑1 (MQ=9)
                           ttCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTAc                                                              <  1:259777/53‑1 (MQ=255)
                           ttCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATa                                                       <  1:2916126/60‑1 (MQ=2)
                           ttCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATa                                                       <  1:1568840/60‑1 (MQ=2)
                             cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTa                                                                   >  1:357278/1‑46 (MQ=255)
                             cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTa                                                       <  1:1653329/58‑1 (MQ=255)
                             cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTa                                                       <  1:2263981/58‑1 (MQ=255)
                             cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCCTa                                                       <  1:4663147/58‑1 (MQ=255)
                             cTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGaaa                                                >  1:3960016/1‑65 (MQ=2)
                                   gTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGt                                                                 <  1:5207718/42‑1 (MQ=255)
                                   gTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTAc                                                              <  1:3922148/45‑1 (MQ=255)
                                   gTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGa                                              <  1:4038902/61‑1 (MQ=255)
                                   gTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATcc                                    <  1:556987/71‑1 (MQ=255)
                                    tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTATTTACGACCGTAGTCGAAAGAACCCACATcc                                    <  1:5483669/70‑1 (MQ=255)
                                    tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCg                                                   <  1:5182644/55‑1 (MQ=255)
                                    tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCg                                                   <  1:379425/55‑1 (MQ=255)
                                    tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                                 <  1:379437/57‑1 (MQ=255)
                                    tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                                 <  1:748351/57‑1 (MQ=255)
                                    tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                                 <  1:4545436/57‑1 (MQ=255)
                                     aCGTCGGTCCAGGAAGTTACGTCATAACCAACGCCGTAGTTa                                                               >  1:2268116/1‑42 (MQ=255)
                                     aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGt                                                                 <  1:3183425/40‑1 (MQ=255)
                                     aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                               >  1:5438572/1‑42 (MQ=255)
                                     aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAAc                                            >  1:1853459/1‑61 (MQ=255)
                                     aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAAc                                            >  1:3098874/1‑61 (MQ=255)
                                      cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTgg                                 <  1:4442729/71‑1 (MQ=255)
                                      cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTg                                  <  1:4490468/70‑1 (MQ=255)
                                      cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTg                                  <  1:1887784/70‑1 (MQ=255)
                                        tCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                               <  1:1407141/55‑1 (MQ=255)
                                        tCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                               <  1:1064637/55‑1 (MQ=255)
                                        tCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTg                                  <  1:5367911/68‑1 (MQ=255)
                                          ggTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAtt                             <  1:4948393/71‑1 (MQ=255)
                                          ggTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAt                              <  1:3516121/70‑1 (MQ=255)
                                          ggTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAt                              <  1:3756184/70‑1 (MQ=255)
                                           gTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATc                                     <  1:3583432/62‑1 (MQ=255)
                                           gTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAttt                            >  1:803002/1‑71 (MQ=255)
                                              cAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                               >  1:4526654/1‑49 (MQ=255)
                                              cAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAttt                            >  1:1061861/1‑68 (MQ=255)
                                              cAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAttt                            >  1:2205639/1‑68 (MQ=255)
                                                 gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGa                        >  1:1133195/1‑69 (MQ=255)
                                                 gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAcc                      >  1:5313212/1‑71 (MQ=255)
                                                 gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAcc                      >  1:2011459/1‑71 (MQ=255)
                                                   aGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGt                                                     >  1:3406862/1‑38 (MQ=255)
                                                   aGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGt                                                     >  1:1089521/1‑38 (MQ=255)
                                                   aGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGaa                               >  1:4401152/1‑60 (MQ=255)
                                                     ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGc                   >  1:78193/1‑70 (MQ=255)
                                                     ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCg                  <  1:584078/71‑1 (MQ=255)
                                                     ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCg                  <  1:1157656/71‑1 (MQ=255)
                                                     ctACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                  <  1:107261/38‑1 (MQ=17)
                                                     ctACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                  <  1:5472907/38‑1 (MQ=17)
                                                      tACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCa                          >  1:3218599/1‑62 (MQ=255)
                                                       aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                               <  1:2195251/40‑1 (MQ=255)
                                                       aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGaa                                             >  1:806747/1‑42 (MQ=255)
                                                       aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGaa                                             >  1:1402119/1‑42 (MQ=255)
                                                       aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGa                 <  1:5567178/70‑1 (MQ=255)
                                                             tAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGc                   >  1:5628671/1‑62 (MQ=255)
                                                             tAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca           >  1:5311221/1‑70 (MQ=255)
                                                              aaaCAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCcaca         >  1:773917/1‑71 (MQ=255)
                                                               aacaacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAAt               >  1:4470062/1‑64 (MQ=255)
                                                                acaacGACGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGa                        >  1:735530/1‑54 (MQ=255)
                                                                 caacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACgg      >  1:1844581/1‑71 (MQ=255)
                                                                 caacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACgg      <  1:2057451/71‑1 (MQ=255)
                                                                 caacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACgg      <  1:3226849/71‑1 (MQ=255)
                                                                     gCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTcc  <  1:72899/71‑1 (MQ=255)
                                                                     |                                                                      
GTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCC  >  minE/1439885‑1440024

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