Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,948 A→T 14.3% V118V (GTT→GTA ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,9480AT14.3% 257.8 / 31.0 98V118V (GTT→GTAompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (44/40);  new base T (6/8);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 127 total aligned reads displayed.

AAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACA  >  minE/1439880‑1440017
                                                                    |                                                                     
aaGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaaa                                                                      >  1:2516699/1‑70 (MQ=255)
aaGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:3136282/71‑1 (MQ=255)
  gTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaaca                                                                    >  1:4046471/1‑70 (MQ=255)
     gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACg                                                                 <  1:2857174/70‑1 (MQ=255)
     gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                >  1:5192443/1‑71 (MQ=255)
     gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                >  1:5586634/1‑71 (MQ=255)
     gTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGc                                                                >  1:5414727/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:1385996/71‑1 (MQ=255)
               tAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACg                                                      >  1:485216/1‑71 (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)
                         caccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTc                                             >  1:3122079/1‑70 (MQ=255)
                         caccAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTc                                             >  1:3933879/1‑70 (MQ=2)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaa                                                                       >  1:594004/1‑43 (MQ=255)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaa                                                                       >  1:1125920/1‑43 (MQ=255)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaa                                                                       <  1:4747159/43‑1 (MQ=255)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaa                                                                       <  1:2054869/43‑1 (MQ=255)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATaa                                                                       >  1:2982612/1‑43 (MQ=255)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGa                                                     >  1:2560936/1‑61 (MQ=255)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                           <  1:3548694/71‑1 (MQ=255)
                          accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                           <  1:1899483/71‑1 (MQ=255)
                          accGAATTATGGCAGTACGTCGGTACAGGAAGTTACGTCATAaacaa                                                                   >  1:3301902/1‑47 (MQ=255)
                          accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCg                                            >  1:2077768/1‑70 (MQ=2)
                          accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                           <  1:3858876/71‑1 (MQ=2)
                          accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                           <  1:3865498/71‑1 (MQ=2)
                           ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:1343961/44‑1 (MQ=255)
                           ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:5397797/44‑1 (MQ=255)
                           ccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:4828836/44‑1 (MQ=255)
                            caaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                           >  1:2805084/3‑69 (MQ=2)
                            cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:1031628/43‑1 (MQ=255)
                            cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:2420399/43‑1 (MQ=255)
                            cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:2006932/43‑1 (MQ=255)
                            cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:231051/43‑1 (MQ=255)
                            cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:3275574/43‑1 (MQ=255)
                            cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:3477313/43‑1 (MQ=255)
                             gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                        >  1:1664598/1‑55 (MQ=255)
                             gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                        >  1:3116911/1‑55 (MQ=255)
                             gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                        >  1:4875436/1‑55 (MQ=255)
                             gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                        >  1:2565344/1‑55 (MQ=255)
                             gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                        >  1:4681827/1‑55 (MQ=255)
                             gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                        <  1:4759651/71‑1 (MQ=255)
                             aaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                           <  1:1639409/67‑1 (MQ=9)
                             aaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGAAAg                                        >  1:4231383/2‑71 (MQ=2)
                              aaTTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaaca                                                                    >  1:5266472/1‑42 (MQ=255)
                                ttCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:3525080/39‑1 (MQ=255)
                                ttCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:5393229/39‑1 (MQ=255)
                                ttCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATa                                                <  1:1568840/60‑1 (MQ=2)
                                ttCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATa                                                <  1:2916126/60‑1 (MQ=2)
                                  cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAaac                                                                     <  1:3709835/37‑1 (MQ=255)
                                  cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTa                                                <  1:1653329/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)
                                        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)
                                         tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                          <  1:379437/57‑1 (MQ=255)
                                         tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                          <  1:748351/57‑1 (MQ=255)
                                          aCGTCGGTCCAGGAAGTTACGTCATAACCAACGCCGTAGTTa                                                        >  1:2268116/1‑42 (MQ=255)
                                          aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTa                                                        >  1:5438572/1‑42 (MQ=255)
                                          aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAAc                                     >  1:1853459/1‑61 (MQ=255)
                                           cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTgg                          <  1:4442729/71‑1 (MQ=255)
                                           cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTg                           <  1:1887784/70‑1 (MQ=255)
                                           cGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTg                           <  1:4490468/70‑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:3756184/70‑1 (MQ=255)
                                               ggTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAt                       <  1:3516121/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:2205639/1‑68 (MQ=255)
                                                   cAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAttt                     >  1:1061861/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:1089521/1‑38 (MQ=255)
                                                        aGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGaa                        >  1:4401152/1‑60 (MQ=255)
                                                          ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGc            >  1:78193/1‑70 (MQ=255)
                                                          ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCg           <  1:1157656/71‑1 (MQ=255)
                                                          ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCg           <  1:584078/71‑1 (MQ=255)
                                                          ctACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                           <  1:5472907/38‑1 (MQ=17)
                                                           tACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCa                   >  1:3218599/1‑62 (MQ=255)
                                                            aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                        <  1:2195251/40‑1 (MQ=255)
                                                            aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGaa                                      >  1:1402119/1‑42 (MQ=255)
                                                            aCGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGaa                                      >  1:806747/1‑42 (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)
                                                                    |                                                                     
AAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACA  >  minE/1439880‑1440017

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: