Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 938,456 C→T 28.6% G162G (GGG→GGA adhP ← alcohol dehydrogenase, 1‑propanol preferring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE938,4560CT28.6% 99.9 / 39.7 95G162G (GGG→GGAadhPalcohol dehydrogenase, 1‑propanol preferring
Reads supporting (aligned to +/- strand):  ref base C (22/45);  new base T (8/19);  total (30/65)
Fisher's exact test for biased strand distribution p-value = 8.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

CATTCTTCGCGTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGTAATGCTG  >  minE/938392‑938519
                                                                |                                                               
cATTCTTCGCGTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTg                                                             <  1:146260/69‑1 (MQ=255)
cATTCTTCGCGTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTg                                                             <  1:1990399/69‑1 (MQ=255)
cATTCTTCGCGTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTg                                                             >  1:2583898/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGaa                                                       >  1:1083469/1‑68 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGaa                                                       >  1:350034/1‑68 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:2321036/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:2820856/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:1407791/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:137788/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:2111067/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:1172527/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:1134827/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:2528674/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:974042/1‑69 (MQ=255)
       cgcgTATTGCAGGGCGAGGTTACCCAGACCGACAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAt                                                      >  1:1973448/1‑69 (MQ=255)
         cgTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttt                                                    >  1:1468169/1‑69 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             >  1:535319/1‑59 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTGGACGAAttt                                                    >  1:1208521/1‑68 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTGGACGAAttt                                                    >  1:2572913/1‑68 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTGGACGAAttt                                                    >  1:200366/1‑68 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTGGACGAAttt                                                    >  1:779863/1‑68 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTGGACGAAttt                                                    >  1:1945055/1‑68 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTGGACGAAttt                                                    >  1:1642973/1‑68 (MQ=255)
          gTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGAATCCACTGTCCTGGACGAAttttt                                                  >  1:983111/1‑69 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:1566337/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:2750883/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:2691349/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:909546/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:2137052/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:429037/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:262085/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:476809/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:611363/69‑1 (MQ=255)
           tATTGCAGGGCGAGGTTACCCAGACCGCCAAGACAGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:3176692/69‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCTAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:2660392/55‑1 (MQ=39)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:1006662/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:2538546/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:494181/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:358373/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:329443/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:264940/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:312950/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:1713392/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:2225743/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:1529872/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:1419878/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:139004/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:1185268/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:2220517/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:1896085/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:1945707/55‑1 (MQ=255)
              tGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:137608/55‑1 (MQ=255)
               gCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGTCCTg                                                             <  1:814285/54‑1 (MQ=255)
                          ttACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAAttttt                                                  <  1:1010669/54‑1 (MQ=255)
                           tACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCtt                                  >  1:2224801/1‑69 (MQ=255)
                            aCCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCtt                                  >  1:2064359/1‑68 (MQ=255)
                            aCCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCtt                                  >  1:1076852/1‑68 (MQ=255)
                            aCCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCtt                                  >  1:97199/1‑68 (MQ=255)
                            aCCCAGACCGCCAAGACCGTAGATAGCAATCAACTGCCCTGGACGAATTTTTGACAGCTTAACGGCtt                                  >  1:315111/1‑68 (MQ=255)
                              ccAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTg                                <  1:1289002/68‑1 (MQ=255)
                              ccAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTg                                <  1:2307561/68‑1 (MQ=255)
                              ccAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTg                                <  1:21350/68‑1 (MQ=255)
                              ccAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTg                                <  1:2380034/68‑1 (MQ=255)
                              ccAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTg                                <  1:2026224/68‑1 (MQ=255)
                               cAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:2947647/69‑1 (MQ=255)
                               cAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:3007560/69‑1 (MQ=255)
                               cAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:2680059/69‑1 (MQ=255)
                               cAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:3317301/69‑1 (MQ=255)
                               cAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:3190572/69‑1 (MQ=255)
                                  aCCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:1200136/66‑1 (MQ=255)
                                  aCCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:3091288/66‑1 (MQ=255)
                                  aCCGCCAAGACCGTAGATAGCAATCCACCGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTa                              <  1:1006155/66‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:564047/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:2538330/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:2457158/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:907251/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:980460/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:1273352/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:2710520/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:2848733/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:1864484/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:1297826/69‑1 (MQ=255)
                                     gCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAAAGGCTTTGTAggtggt                        <  1:884614/69‑1 (MQ=255)
                                      ccAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAggtggt                        <  1:1522803/68‑1 (MQ=255)
                                          gACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTc                   <  1:2876861/69‑1 (MQ=255)
                                               tataTAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGca              <  1:1141745/66‑1 (MQ=255)
                                                   tAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAgg          <  1:613388/69‑1 (MQ=255)
                                                   tAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAgg          <  1:2227608/69‑1 (MQ=255)
                                                    aGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGt         <  1:575230/69‑1 (MQ=255)
                                                    aGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGt         <  1:790533/69‑1 (MQ=255)
                                                    aGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGt         <  1:3086601/69‑1 (MQ=255)
                                                           cACTGGCCTGGACGTATTTTTGACAGCTTATCGGGTTTGTa                              <  1:972925/41‑1 (MQ=25)
                                                           cACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGTAAtgctg  >  1:2780884/1‑69 (MQ=255)
                                                           cACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGTAAtgctg  >  1:1834224/1‑69 (MQ=255)
                                                           cACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGTAAtgctg  >  1:2421145/1‑69 (MQ=255)
                                                                |                                                               
CATTCTTCGCGTATTGCAGGGCGAGGTTACCCAGACCGCCAAGACCGTAGATAGCAATCCACTGCCCTGGACGAATTTTTGACAGCTTAACGGCTTTGTAGGTGGTGACTCCCGCACAGGTAATGCTG  >  minE/938392‑938519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: