Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,969 G→A 14.6% Y111Y (TAC→TAT ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,9690GA14.6% 264.1 / 25.4 96Y111Y (TAC→TATompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base A (6/8);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGT  >  minE/1439900‑1440039
                                                                     |                                                                      
gtCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCg                                                                        >  1:1314500/1‑70 (MQ=255)
   gccaccaAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGt                                                                    >  1:1029332/2‑71 (MQ=2)
   accaccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAg                                                                     <  1:234209/70‑1 (MQ=255)
     caccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTc                                                                   >  1:3122079/1‑70 (MQ=255)
     caccAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTc                                                                   >  1:3933879/1‑70 (MQ=2)
      accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                                 <  1:3548694/71‑1 (MQ=255)
      accGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                                 <  1:1899483/71‑1 (MQ=255)
      accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCg                                                                  >  1:2077768/1‑70 (MQ=2)
      accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                                 <  1:2183233/71‑1 (MQ=2)
      accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                                 <  1:3865498/71‑1 (MQ=2)
      accAAATTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                                 <  1:3858876/71‑1 (MQ=2)
        caaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                                 >  1:2805084/3‑69 (MQ=2)
         gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGa                                                                 >  1:1224023/1‑68 (MQ=255)
         gAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAg                                                              <  1:4759651/71‑1 (MQ=255)
         aaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                                 <  1:1639409/67‑1 (MQ=9)
         aaaTTCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGAAAg                                                              >  1:4231383/2‑71 (MQ=2)
            ttCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATa                                                                      <  1:2916126/60‑1 (MQ=2)
            ttCTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATa                                                                      <  1:1568840/60‑1 (MQ=2)
              cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTa                                                                      <  1:2263981/58‑1 (MQ=255)
              cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTa                                                                      <  1:1653329/58‑1 (MQ=255)
              cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCCTa                                                                      <  1:4663147/58‑1 (MQ=255)
              cTGGCAGCATATCGGTCCATGCTTCTACGTCATATACGACACCGTAGTTACGACCATAGTCGaaa                                                               >  1:3960016/1‑65 (MQ=2)
                    gTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGa                                                             <  1:4038902/61‑1 (MQ=255)
                    gTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATcc                                                   <  1:556987/71‑1 (MQ=255)
                     tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTATTTACGACCGTAGTCGAAAGAACCCACATcc                                                   <  1:5483669/70‑1 (MQ=255)
                     tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCg                                                                  <  1:379425/55‑1 (MQ=255)
                     tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCg                                                                  <  1:5182644/55‑1 (MQ=255)
                     tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                                                <  1:4545436/57‑1 (MQ=255)
                     tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                                                <  1:379437/57‑1 (MQ=255)
                     tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGaa                                                                <  1:748351/57‑1 (MQ=255)
                      aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAAc                                                           >  1:3098874/1‑61 (MQ=255)
                      aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAAc                                                           >  1:1853459/1‑61 (MQ=255)
                      aCGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTg                                                 >  1:855759/1‑71 (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)
                                      ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCg                                 <  1:1933110/71‑1 (MQ=255)
                                      ctACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                                 <  1:5472907/38‑1 (MQ=17)
                                      ctACGTCATATACGACACCGTAGTTACGACCATAGTCGa                                                                 <  1:107261/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)
                                        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:3226849/71‑1 (MQ=255)
                                                  caacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACgg                     <  1:2057451/71‑1 (MQ=255)
                                                      gCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTcc                 <  1:72899/71‑1 (MQ=255)
                                                       ccGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAAttt                                           <  1:2193276/44‑1 (MQ=255)
                                                        cGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAg               <  1:1995693/71‑1 (MQ=255)
                                                        cGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGGCCAg               <  1:2002255/71‑1 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAAGCCACATCCTGGAATTTc                                          >  1:1462831/1‑40 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTTGAATTTCAGACCTGCGAATGCCacac                       >  1:4641129/1‑59 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGaa                                              >  1:985576/1‑36 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGaa                                              >  1:4183890/1‑36 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGaa                                              >  1:4154746/1‑36 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTc                                          >  1:1967222/1‑40 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCt                                    <  1:3095008/46‑1 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCt                                    <  1:2243771/46‑1 (MQ=255)
                                                            gTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgt           >  1:3962824/1‑71 (MQ=255)
                                                             ttACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGa                                       <  1:3851435/42‑1 (MQ=255)
                                                               aCGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAc                                      >  1:865864/1‑41 (MQ=255)
                                                               aCGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGAc                                      >  1:1098079/1‑41 (MQ=255)
                                                               aCGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgttgt        >  1:3231190/1‑71 (MQ=255)
                                                               aCGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgttgt        >  1:2263199/1‑71 (MQ=255)
                                                               aCGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgttgt        >  1:1771934/1‑71 (MQ=255)
                                                               aCGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgtt          >  1:810727/1‑69 (MQ=255)
                                                                cGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGc                                  <  1:2952794/44‑1 (MQ=255)
                                                                cGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAgttgt        <  1:1134603/70‑1 (MQ=255)
                                                                 gACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca                          >  1:4494731/1‑51 (MQ=255)
                                                                 gACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCca                          >  1:5175015/1‑51 (MQ=255)
                                                                   ccGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACgg                     >  1:1495117/1‑54 (MQ=255)
                                                                     gTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttttcgt  <  1:5010771/71‑1 (MQ=255)
                                                                     |                                                                      
GTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGT  >  minE/1439900‑1440039

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: