Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,789 C→T 16.8% Q171Q (CAG→CAA ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,7890CT16.8% 249.1 / 33.6 102Q171Q (CAG→CAAompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (48/36);  new base T (8/9);  total (57/45)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 107 total aligned reads displayed.

CAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAG  >  minE/1439720‑1439858
                                                                     |                                                                     
cAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCt                                                                      <  1:1356938/71‑1 (MQ=255)
cAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCt                                                                      <  1:5696637/71‑1 (MQ=255)
    gCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                                                                  >  1:4927382/1‑71 (MQ=255)
    gCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTa                                                                  >  1:3931361/1‑71 (MQ=255)
      gTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACt                                                                <  1:2349787/71‑1 (MQ=255)
      gTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACt                                                                >  1:4243109/1‑71 (MQ=255)
      gTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACt                                                                <  1:1559443/71‑1 (MQ=255)
      gTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACt                                                                <  1:2968641/71‑1 (MQ=255)
       tCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTgg                                                                    <  1:3752955/66‑1 (MQ=255)
           gACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                             >  1:5640637/1‑69 (MQ=255)
           gACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                             >  1:4105488/1‑69 (MQ=255)
               gttgttAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                        <  1:2955158/70‑1 (MQ=255)
                ttgttAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                        >  1:2907751/1‑69 (MQ=255)
                  gttagttaCGCCACTATTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAg                                                     <  1:4691544/70‑1 (MQ=255)
                  gttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAg                                                     <  1:867747/70‑1 (MQ=255)
                         aCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                             >  1:2776364/1‑71 (MQ=255)
                           gCCACTAGTAAAGCCTTCACAAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGt                                           <  1:2417803/71‑1 (MQ=255)
                               cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                       >  1:4028265/1‑71 (MQ=255)
                               cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                       >  1:3794781/1‑71 (MQ=255)
                               cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                       >  1:901969/1‑71 (MQ=255)
                                tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACa                                                           >  1:5563721/1‑50 (MQ=255)
                                 aGTAAAGCCTTCACAAGATGGGTTGCCGTTTTTACCATgg                                                                    >  1:2312460/1‑40 (MQ=38)
                                  gTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                     <  1:2827934/70‑1 (MQ=255)
                                   tAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAg                                                     <  1:946638/53‑1 (MQ=255)
                                     aaGCGTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGc                                              >  1:4789897/1‑58 (MQ=255)
                                     aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTCAACAGCaaa                                                      >  1:4968357/1‑50 (MQ=255)
                                        ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                        >  1:2616335/1‑61 (MQ=255)
                                          ttCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                             >  1:3829570/1‑70 (MQ=255)
                                           tCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGCCCGTCAACCAGACCgaag                             >  1:735838/1‑69 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCa                                                 >  1:3144134/1‑48 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                        >  1:4472787/1‑57 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                        >  1:2568228/1‑57 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                        >  1:5603419/1‑57 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                        >  1:2246951/1‑57 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                          >  1:4820904/1‑71 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                          >  1:1775763/1‑71 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                          >  1:4258669/1‑71 (MQ=255)
                                            cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                          >  1:4628469/1‑71 (MQ=255)
                                             aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                          <  1:3112795/70‑1 (MQ=255)
                                             aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                         <  1:5050667/71‑1 (MQ=255)
                                             aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                         <  1:4362760/71‑1 (MQ=255)
                                             aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                         <  1:447477/71‑1 (MQ=255)
                                              ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                         >  1:4219703/1‑70 (MQ=255)
                                                aGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                        <  1:3519884/37‑1 (MQ=255)
                                                 gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                       >  1:1593116/1‑53 (MQ=255)
                                                 gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                             <  1:4076419/63‑1 (MQ=255)
                                                 gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                             <  1:2264004/63‑1 (MQ=255)
                                                 gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                     >  1:346052/1‑71 (MQ=255)
                                                 gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                     >  1:829005/1‑71 (MQ=255)
                                                 gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                     >  1:1500305/1‑71 (MQ=255)
                                                 gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCCTCAACCAGACCgaag                             <  1:3613102/63‑1 (MQ=255)
                                                  aTGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTc                                          <  1:3463405/49‑1 (MQ=255)
                                                   tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                   <  1:594962/71‑1 (MQ=255)
                                                     ggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                 <  1:2939076/71‑1 (MQ=255)
                                                       ttGCCGTTTTTACCCTTGTACTGAACAGCAAAGtt                                                   <  1:3747349/35‑1 (MQ=255)
                                                       ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                   <  1:3148972/35‑1 (MQ=255)
                                                       ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                   <  1:2788405/35‑1 (MQ=255)
                                                       ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                   <  1:4288928/35‑1 (MQ=255)
                                                       ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                       <  1:2691678/63‑1 (MQ=255)
                                                       ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                       <  1:5095300/63‑1 (MQ=255)
                                                         ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAgt                     <  1:1222257/61‑1 (MQ=2)
                                                         ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACga              >  1:1921524/3‑69 (MQ=255)
                                                         ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACga              >  1:245171/3‑69 (MQ=255)
                                                         gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                >  1:2116530/1‑36 (MQ=255)
                                                         gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg              <  1:2365396/70‑1 (MQ=255)
                                                          ccGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                              <  1:3448231/53‑1 (MQ=255)
                                                           cGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat             >  1:2673864/1‑67 (MQ=9)
                                                            gttgttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATaa           <  1:913249/70‑2 (MQ=255)
                                                            gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg          >  1:3481802/1‑71 (MQ=255)
                                                             tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                        >  1:4609455/1‑40 (MQ=255)
                                                             tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCg                                <  1:1208302/48‑1 (MQ=255)
                                                             tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt         >  1:3365252/1‑71 (MQ=255)
                                                             tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt         <  1:5418776/71‑1 (MQ=255)
                                                             ttgttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGt         <  1:4163262/68‑1 (MQ=9)
                                                             ttattACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGt         >  1:944137/4‑71 (MQ=255)
                                                             ttattACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGt         >  1:5490431/4‑71 (MQ=255)
                                                              tattACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGtt        >  1:2472725/3‑70 (MQ=255)
                                                               tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                     >  1:2124127/1‑41 (MQ=255)
                                                               tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                     >  1:5206816/1‑41 (MQ=255)
                                                               tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg       >  1:5234665/1‑71 (MQ=255)
                                                               tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg       >  1:2293748/1‑71 (MQ=255)
                                                               gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcg       >  1:2122834/2‑71 (MQ=255)
                                                               gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcg       >  1:1540319/2‑71 (MQ=255)
                                                               gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc        <  1:1903393/69‑1 (MQ=255)
                                                               gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc        <  1:4359271/69‑1 (MQ=255)
                                                                ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc      <  1:5324948/71‑1 (MQ=9)
                                                                ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc      <  1:399180/71‑1 (MQ=9)
                                                                ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc      <  1:3486189/71‑1 (MQ=9)
                                                                ttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                              <  1:4326007/47‑1 (MQ=255)
                                                                ttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                              <  1:755752/47‑1 (MQ=255)
                                                                 tACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc      >  1:1925733/1‑70 (MQ=9)
                                                                 tACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc      >  1:1432098/1‑70 (MQ=9)
                                                                  aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgcg     <  1:1728528/70‑1 (MQ=255)
                                                                  aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg     >  1:2255134/1‑70 (MQ=255)
                                                                  aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg     >  1:77650/1‑70 (MQ=255)
                                                                   cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGCAGAAGTCAGTGTTACGGTAGGTCGCGaa   <  1:4488359/71‑1 (MQ=255)
                                                                   cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGtt                  <  1:1886572/56‑1 (MQ=255)
                                                                   cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgc      <  1:2288497/68‑1 (MQ=255)
                                                                   cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa    <  1:446008/70‑3 (MQ=255)
                                                                    ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg  <  1:1428588/71‑1 (MQ=255)
                                                                     |                                                                     
CAGTGCGTCACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAG  >  minE/1439720‑1439858

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: