Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,801 A→C 15.8% A167A (GCT→GCG ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,8010AC15.8% 297.5 / 51.2 133A167A (GCT→GCGompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (52/60);  new base C (12/9);  total (64/69)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 139 total aligned reads displayed.

GTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGC  >  minE/1439735‑1439870
                                                                  |                                                                     
gttgttAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                                    <  1:2955158/70‑1 (MQ=255)
 ttgttAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                                    >  1:2907751/1‑69 (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)
                 tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACa                                                                       >  1:5563721/1‑50 (MQ=255)
                      aaGCGTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGc                                                          >  1:4789897/1‑58 (MQ=255)
                         ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                    >  1:2616335/1‑61 (MQ=255)
                           ttCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                         >  1:3829570/1‑70 (MQ=255)
                             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCa                                                             >  1:3144134/1‑48 (MQ=255)
                             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                    >  1:2246951/1‑57 (MQ=255)
                             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                    >  1:2544921/1‑57 (MQ=255)
                             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                      >  1:1775763/1‑71 (MQ=255)
                             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                      >  1:4628469/1‑71 (MQ=255)
                             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                      >  1:4820904/1‑71 (MQ=255)
                              aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                     <  1:447477/71‑1 (MQ=255)
                              aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                     <  1:5050667/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)
                                  gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                 >  1:829005/1‑71 (MQ=255)
                                  gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                 >  1:346052/1‑71 (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:194325/35‑1 (MQ=255)
                                        ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                               <  1:4288928/35‑1 (MQ=255)
                                        ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                               <  1:3148972/35‑1 (MQ=255)
                                        ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                               <  1:2788405/35‑1 (MQ=255)
                                        ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                   <  1:5095300/63‑1 (MQ=255)
                                          ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAgt                                 <  1:1222257/61‑1 (MQ=2)
                                          ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACga                          >  1:245171/3‑69 (MQ=255)
                                          gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                            >  1:2116530/1‑36 (MQ=255)
                                          gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                          <  1:2365396/70‑1 (MQ=255)
                                           tcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                         >  1:5514461/2‑68 (MQ=2)
                                            cGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                         >  1:2673864/1‑67 (MQ=9)
                                             gttgttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATaa                       <  1:913249/70‑2 (MQ=255)
                                             gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                      >  1:4051572/1‑71 (MQ=255)
                                             gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                      >  1:3481802/1‑71 (MQ=255)
                                              tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCg                                            <  1:1208302/48‑1 (MQ=255)
                                              tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                     <  1:5418776/71‑1 (MQ=255)
                                              tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                     >  1:3365252/1‑71 (MQ=255)
                                              ttgttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGt                     <  1:4163262/68‑1 (MQ=9)
                                              ttattACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGt                     >  1:944137/4‑71 (MQ=255)
                                                tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                                 >  1:5206816/1‑41 (MQ=255)
                                                tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCa                                                 >  1:2124127/1‑41 (MQ=255)
                                                tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                   >  1:5234665/1‑71 (MQ=255)
                                                tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                   >  1:2293748/1‑71 (MQ=255)
                                                gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcg                   >  1:1540319/2‑71 (MQ=255)
                                                gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcg                   >  1:2122834/2‑71 (MQ=255)
                                                gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc                    <  1:1903393/69‑1 (MQ=255)
                                                gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc                    <  1:4359271/69‑1 (MQ=255)
                                                 ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                  <  1:399180/71‑1 (MQ=9)
                                                 ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                  <  1:5324948/71‑1 (MQ=9)
                                                 ttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                          <  1:755752/47‑1 (MQ=255)
                                                 ttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                          <  1:4326007/47‑1 (MQ=255)
                                                  tACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                  >  1:1925733/1‑70 (MQ=9)
                                                  tACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                  >  1:1432098/1‑70 (MQ=9)
                                                   aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                 >  1:2255134/1‑70 (MQ=255)
                                                   aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg                 >  1:77650/1‑70 (MQ=255)
                                                   aCCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa                >  1:4607544/1‑69 (MQ=9)
                                                    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)
                                                       tgatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGc             >  1:1718842/4‑70 (MQ=9)
                                                        ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                             >  1:2497059/1‑53 (MQ=255)
                                                        ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg           >  1:1261773/1‑71 (MQ=255)
                                                        gatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAGTCAGTATTACGATAAGTCGCGACGcca           >  1:3367110/3‑70 (MQ=9)
                                                         gTACTGAACAGCAAAGTTCAGGCCTTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt          <  1:2621710/71‑1 (MQ=255)
                                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc            <  1:888923/69‑1 (MQ=255)
                                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc            <  1:4409602/69‑1 (MQ=255)
                                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg           <  1:4788670/70‑1 (MQ=255)
                                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt          <  1:4839570/71‑1 (MQ=255)
                                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt          >  1:1097808/1‑71 (MQ=255)
                                                          tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCGGTGTTACGGTAGGTCGCGAAGc             <  1:5030844/67‑1 (MQ=255)
                                                          tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt         <  1:5506495/71‑1 (MQ=255)
                                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                          <  1:5617523/35‑1 (MQ=255)
                                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                          <  1:5518907/35‑1 (MQ=255)
                                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                    <  1:2209405/41‑1 (MQ=255)
                                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                    <  1:3644829/41‑1 (MQ=255)
                                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                  <  1:1015214/43‑1 (MQ=255)
                                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGc             <  1:4635692/64‑1 (MQ=255)
                                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt         <  1:1523523/68‑1 (MQ=255)
                                                              gAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                         <  1:3798414/35‑1 (MQ=255)
                                                              gAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                   <  1:4521908/41‑1 (MQ=255)
                                                               aaCAGCAAAGTTCAGGCCGTGAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa               >  1:4705293/1‑60 (MQ=255)
                                                               aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa               >  1:4595357/1‑60 (MQ=255)
                                                               aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc    >  1:4472473/1‑71 (MQ=255)
                                                                aCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                 >  1:316611/1‑41 (MQ=255)
                                                                 cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc  <  1:5579924/71‑1 (MQ=255)
                                                                 cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc  <  1:1538613/71‑1 (MQ=255)
                                                                 aaGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAgg                      >  1:5216332/2‑51 (MQ=255)
                                                                 aaGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAgg                      >  1:72759/2‑51 (MQ=255)
                                                                 aaGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAgg                      >  1:3733531/2‑51 (MQ=255)
                                                                 aaGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAgg                      >  1:3497126/2‑51 (MQ=255)
                                                                 aaGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAgg                      >  1:1010531/2‑51 (MQ=255)
                                                                  |                                                                     
GTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGC  >  minE/1439735‑1439870

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: