Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,819 A→T 16.4% V161V (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,8190AT16.4% 294.7 / 50.0 132V161V (GTT→GTAompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (45/55);  new base T (12/9);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 144 total aligned reads displayed.

CTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCA  >  minE/1439751‑1439888
                                                                    |                                                                     
cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                                     >  1:901969/1‑71 (MQ=255)
cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                                     >  1:3794781/1‑71 (MQ=255)
            tCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGCCCGTCAACCAGACCgaag                                                           >  1:735838/1‑69 (MQ=255)
             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                                      >  1:2246951/1‑57 (MQ=255)
             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                                      >  1:2544921/1‑57 (MQ=255)
             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                                      >  1:2568228/1‑57 (MQ=255)
             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                                        >  1:4258669/1‑71 (MQ=255)
             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                                        >  1:4820904/1‑71 (MQ=255)
             cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                                        >  1:4628469/1‑71 (MQ=255)
              aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                       <  1:4362760/71‑1 (MQ=255)
              aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                       <  1:5050667/71‑1 (MQ=255)
              aCCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                       <  1:447477/71‑1 (MQ=255)
               ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGt                                                       >  1:4219703/1‑70 (MQ=255)
                  gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                                     >  1:1593116/1‑53 (MQ=255)
                  gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                                           <  1:4076419/63‑1 (MQ=255)
                  gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                                   >  1:829005/1‑71 (MQ=255)
                  gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                                   >  1:346052/1‑71 (MQ=255)
                  gATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCCTCAACCAGACCgaag                                                           <  1:3613102/63‑1 (MQ=255)
                    tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                                                 <  1:594962/71‑1 (MQ=255)
                      ggTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                                               <  1:2939076/71‑1 (MQ=255)
                        ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                                     <  1:5095300/63‑1 (MQ=255)
                        ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                                     <  1:2691678/63‑1 (MQ=255)
                          ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAgt                                                   <  1:1222257/61‑1 (MQ=2)
                          ttcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACga                                            >  1:245171/3‑69 (MQ=255)
                           tcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                                           >  1:5514461/2‑68 (MQ=2)
                           ccGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                                            <  1:3448231/53‑1 (MQ=255)
                            cGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                                           >  1:2673864/1‑67 (MQ=9)
                             gttgttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATaa                                         <  1:913249/70‑2 (MQ=255)
                             gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                                        >  1:4051572/1‑71 (MQ=255)
                              tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCaa                                                                      >  1:4609455/1‑40 (MQ=255)
                              tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                       <  1:5418776/71‑1 (MQ=255)
                              tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                                       >  1:3365252/1‑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:2293748/1‑71 (MQ=255)
                                tttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcg                                     >  1:5234665/1‑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:3486189/71‑1 (MQ=9)
                                 ttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgc                                    <  1:5324948/71‑1 (MQ=9)
                                  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)
                                   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)
                                        ggTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg                             >  1:1261773/1‑71 (MQ=255)
                                         gTACTGAACAGCAAAGTTCAGGCCTTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                            <  1:2621710/71‑1 (MQ=255)
                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc                              <  1:1149380/69‑1 (MQ=255)
                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCg                             <  1:4788670/70‑1 (MQ=255)
                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                            >  1:3983863/1‑71 (MQ=255)
                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                            >  1:2312503/1‑71 (MQ=255)
                                         gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                            <  1:798438/71‑1 (MQ=255)
                                          tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCGGTGTTACGGTAGGTCGCGAAGc                               <  1:5030844/67‑1 (MQ=255)
                                          tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt                           <  1:2037482/71‑1 (MQ=255)
                                           aCTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGaa                                 >  1:5538641/1‑64 (MQ=255)
                                            cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                            <  1:2206403/52‑1 (MQ=255)
                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                                            <  1:5617523/35‑1 (MQ=255)
                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaa                                                            <  1:5518907/35‑1 (MQ=255)
                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTc                                                      <  1:3644829/41‑1 (MQ=255)
                                             tGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAg                                                    <  1:1015214/43‑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:1589268/1‑60 (MQ=255)
                                               aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc                      >  1:4472473/1‑71 (MQ=255)
                                                aCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgt                                                   >  1:316611/1‑41 (MQ=255)
                                                 cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc                    <  1:1538613/71‑1 (MQ=255)
                                                 cAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgc                    <  1:5579924/71‑1 (MQ=255)
                                                 aaGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAgg                                        >  1:1010531/2‑51 (MQ=255)
                                                 aaGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAgg                                        >  1:72759/2‑51 (MQ=255)
                                                   gCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt                            >  1:4543191/1‑61 (MQ=255)
                                                     aaaGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACgctgct                >  1:3365846/1‑71 (MQ=255)
                                                       aGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                            <  1:1058974/41‑1 (MQ=255)
                                                       aGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                                            <  1:1747386/41‑1 (MQ=255)
                                                       aGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGtt                                      >  1:4334623/1‑46 (MQ=255)
                                                           cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATg           >  1:3495420/1‑70 (MQ=255)
                                                           cAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATg           >  1:5108051/1‑70 (MQ=255)
                                                             ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg        <  1:3414457/71‑1 (MQ=255)
                                                             ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg        <  1:792965/71‑1 (MQ=255)
                                                             ggCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAg        <  1:5264623/71‑1 (MQ=255)
                                                               ccGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGt       >  1:3921721/1‑70 (MQ=255)
                                                                cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg     <  1:5125337/71‑1 (MQ=255)
                                                                cGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTg     <  1:1414085/71‑1 (MQ=255)
                                                                  tCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGaa         <  1:5215595/65‑1 (MQ=255)
                                                                  tCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt    >  1:1187373/1‑70 (MQ=255)
                                                                   cAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt    >  1:4838349/1‑69 (MQ=255)
                                                                    aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt    >  1:68521/1‑68 (MQ=255)
                                                                    aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt    >  1:1870619/1‑68 (MQ=255)
                                                                    aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGt    >  1:3217950/1‑68 (MQ=255)
                                                                    aaCCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCa  <  1:2939990/70‑1 (MQ=255)
                                                                    |                                                                     
CTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCACGCTGCTGCATGAAGTTGTCA  >  minE/1439751‑1439888

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: