Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,440,194 A→G 14.3% G36G (GGT→GGC ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,440,1940AG14.3% 248.6 / 21.9 100G36G (GGT→GGCompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (46/36);  new base G (10/2);  total (61/39)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 101 total aligned reads displayed.

AAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAG  >  minE/1440125‑1440262
                                                                     |                                                                    
aaGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTa                                                                      >  1:3048947/1‑70 (MQ=255)
aaGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                                     >  1:5072382/1‑71 (MQ=255)
  gACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   >  1:4437036/1‑71 (MQ=255)
   aCGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   >  1:4789128/1‑70 (MQ=255)
    cGCATGTATGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                                 <  1:4228482/71‑1 (MQ=255)
      cATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                               <  1:5384074/71‑1 (MQ=255)
        tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   <  1:4823278/65‑1 (MQ=255)
        tGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGa                                                             <  1:1611298/71‑1 (MQ=255)
          tAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATc                                                           >  1:757451/1‑71 (MQ=255)
          tAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATc                                                           >  1:371044/1‑71 (MQ=255)
            ggTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                                     >  1:243293/1‑59 (MQ=255)
            ggTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                                     >  1:595585/1‑59 (MQ=255)
              tCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                                       <  1:106870/71‑1 (MQ=255)
               cTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGTCCGTCTACTTTACCGTACAGATCTAAtt                                                      >  1:46530/1‑71 (MQ=255)
                  gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCGGTCTACTTTACCGTACAGATCTAAttt                                                     <  1:2837429/69‑1 (MQ=255)
                  gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAttt                                                     <  1:5245665/69‑1 (MQ=255)
                  gTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   <  1:4054538/71‑1 (MQ=255)
                   tCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgtt                                                  <  1:3972393/71‑1 (MQ=255)
                   tCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgtt                                                  <  1:4415298/71‑1 (MQ=255)
                      ccATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGcc                                               <  1:624218/71‑1 (MQ=255)
                      ccATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGcc                                               >  1:3906280/1‑71 (MQ=255)
                       cATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCg                                              >  1:689984/1‑71 (MQ=255)
                         tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                              <  1:5167924/53‑1 (MQ=255)
                         tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                              <  1:5394859/53‑1 (MQ=255)
                         tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAg                                                              <  1:5521815/53‑1 (MQ=255)
                         tCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCAGTCTACTTTACCGTACAGATCTAAtt                                                      <  1:3704217/61‑1 (MQ=255)
                          cTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                           >  1:2346539/1‑71 (MQ=255)
                          cTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                           >  1:3663035/1‑71 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                                 >  1:4845125/1‑48 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTa                                                                 >  1:78174/1‑48 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:2906427/1‑71 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:4393193/1‑71 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:431957/1‑71 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:2625383/1‑71 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:1934871/1‑71 (MQ=255)
                           tACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          >  1:3347306/1‑71 (MQ=255)
                             cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCTTCTttg                                        <  1:180053/71‑1 (MQ=255)
                             cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCt                                           <  1:2752969/68‑1 (MQ=255)
                             cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttg                                        <  1:3816773/71‑1 (MQ=255)
                             cATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttg                                        <  1:696502/71‑1 (MQ=255)
                              aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   <  1:914361/43‑1 (MQ=255)
                              aTCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCg                                                                   <  1:4857748/43‑1 (MQ=255)
                               tCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                                       <  1:780484/54‑1 (MQ=255)
                               tCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   >  1:4338540/1‑58 (MQ=255)
                                cTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTAc                                                                     >  1:5208711/1‑39 (MQ=255)
                                   tgttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTa                                   <  1:4519547/70‑1 (MQ=255)
                                   tgttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTa                                   >  1:2978238/1‑70 (MQ=255)
                                    gttgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTaa                                  <  1:3901941/70‑1 (MQ=255)
                                     ttatCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATaa                                  >  1:4584766/4‑69 (MQ=2)
                                     ttatCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:5075339/4‑71 (MQ=9)
                                      tgtCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTttgtt                                      >  1:5530487/1‑64 (MQ=255)
                                       atCAGAAAAATAATGCAGGCCATCAACTTTGCCATACAGATCTAATttgt                                                   >  1:1072543/2‑50 (MQ=2)
                                        tCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACtt                              >  1:323894/1‑70 (MQ=255)
                                        tCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACtt                              <  1:2968300/70‑1 (MQ=255)
                                          agagAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                               >  1:3978740/1‑35 (MQ=255)
                                          agagAAATAGTGCAGGCCGTCTACTTTACCGTACa                                                               >  1:3062996/1‑35 (MQ=255)
                                          agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAAt                                                       <  1:4981965/43‑1 (MQ=255)
                                          agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   >  1:63663/1‑47 (MQ=255)
                                          agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATttgt                                                   >  1:1637962/1‑47 (MQ=255)
                                          agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                           >  1:2106970/1‑71 (MQ=255)
                                          agagAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcag                           >  1:4428153/1‑71 (MQ=255)
                                            aaaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAAttt                                                     <  1:3358368/41‑1 (MQ=2)
                                             gaAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTcagcag                        >  1:1661661/1‑71 (MQ=255)
                                             aaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGt                                             >  1:5335882/2‑50 (MQ=14)
                                             aaaaTAATGCAGGCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTcggctg                        >  1:1430497/2‑69 (MQ=9)
                                              aaaTAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTc                             >  1:3737092/1‑65 (MQ=255)
                                                  aGTGCAGGCCGTCTACTTTTCCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                   >  1:1521748/1‑71 (MQ=255)
                                                  aGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                            >  1:4284700/1‑62 (MQ=255)
                                                   gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTaaa                                 >  1:2672177/1‑56 (MQ=255)
                                                   gTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtg                  <  1:3563717/71‑1 (MQ=255)
                                                   gTGCAGGCCGTCTACTTTAACGTACAGATCTAATTTGTTGCCg                                              >  1:1951569/1‑43 (MQ=255)
                                                      cAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                            <  1:4235278/58‑1 (MQ=255)
                                                       aGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTtgctgc              <  1:482148/71‑1 (MQ=255)
                                                        ggCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcg             <  1:767531/71‑1 (MQ=255)
                                                        ggCCATCAACTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTCGGCTGCATGtgcagca             <  1:2463377/71‑2 (MQ=12)
                                                            gTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCtt                                          <  1:1041773/38‑1 (MQ=255)
                                                             tCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTg         <  1:3581500/70‑1 (MQ=255)
                                                               tACTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGCCAGc        <  1:3016962/69‑1 (MQ=2)
                                                               tACTTTACCGTACAGATCTAATTTGTTGCCGTCTttgtt                                      >  1:605422/1‑39 (MQ=255)
                                                               tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTca                            >  1:5107266/1‑49 (MQ=255)
                                                               tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTc                             >  1:2609717/1‑48 (MQ=255)
                                                               tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGttt                   <  1:3279653/58‑1 (MQ=255)
                                                               tACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa      >  1:1528568/1‑71 (MQ=255)
                                                               aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:4446776/2‑45 (MQ=17)
                                                               aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:5362326/2‑45 (MQ=17)
                                                               aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:1064332/2‑45 (MQ=17)
                                                               aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAAc                                >  1:1861550/2‑45 (MQ=17)
                                                               aaCTTTGCCATACAGATCTAATTTGTTGCCGTCTTTATTATAAACTTcg                            >  1:2065895/2‑48 (MQ=17)
                                                                aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:165270/1‑65 (MQ=2)
                                                                aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:470653/1‑65 (MQ=2)
                                                                aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:2294653/1‑65 (MQ=2)
                                                                aCTTTTCCGTACAAATCCAGTTTATTACCATCTTTATTATAGACTTCAGCCGCATTTGCTGCGcc           >  1:1285980/1‑65 (MQ=2)
                                                                 cTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTa      <  1:4598928/69‑1 (MQ=255)
                                                                  tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTgcgc            <  1:4836368/62‑1 (MQ=255)
                                                                  tttACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTAc     >  1:4481721/1‑69 (MQ=255)
                                                                    tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:1937160/1‑70 (MQ=255)
                                                                    tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:4389517/1‑70 (MQ=255)
                                                                    tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:1621446/1‑70 (MQ=255)
                                                                    tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:520484/1‑70 (MQ=255)
                                                                    tACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACcag  >  1:496206/1‑70 (MQ=255)
                                                                     |                                                                    
AAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACAGATCTAATTTGTTGCCGTCTTTGTTGTAAACTTCAGCAGCGTTTGCTGCGCCTGCTACCAG  >  minE/1440125‑1440262

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: