Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 942,078 A→T 41.0% I117I (ATA→ATT osmC → osmotically inducible, stress‑inducible membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE942,0780AT41.0% 41.7 / 73.0 78I117I (ATA→ATTosmCosmotically inducible, stress‑inducible membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/16);  new base T (14/18);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 6.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGAT  >  minE/942030‑942141
                                                |                                                               
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:1686853/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:2549573/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:1046865/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:1570973/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:417730/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:1126523/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:3244411/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:3233528/68‑1 (MQ=255)
gagTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:2577642/68‑1 (MQ=255)
gagTGAAGTGGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAg                                              <  1:115243/68‑1 (MQ=255)
   tGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATg                                          >  1:288943/1‑69 (MQ=255)
   tGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATg                                          >  1:2166612/1‑69 (MQ=255)
   tGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATg                                          >  1:2652564/1‑69 (MQ=255)
   tGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATg                                          >  1:2846349/1‑69 (MQ=255)
   tGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATg                                          >  1:517631/1‑69 (MQ=255)
   tGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATg                                          >  1:1039784/1‑69 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:282641/64‑1 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:2675746/64‑1 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:1029490/64‑1 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:752162/64‑1 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:2260118/64‑1 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:962412/64‑1 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:1010657/64‑1 (MQ=255)
   cgAAGTTGCGGTGCCGGGTATTGATGCCTCCACCTTTGACGTCATTATCCAGAAAGCAAAAGCAg                                              <  1:1616622/64‑1 (MQ=255)
     aaGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGcc                                        >  1:1111204/1‑69 (MQ=255)
     aaGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGcc                                        >  1:3076243/1‑69 (MQ=255)
     aaGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGcc                                        >  1:144570/1‑69 (MQ=255)
     aaGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGcc                                        >  1:1491820/1‑69 (MQ=255)
     aaGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGcc                                        >  1:2349940/1‑69 (MQ=255)
     aaGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGcc                                        >  1:2011076/1‑69 (MQ=255)
       gTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCg                                      <  1:2746963/69‑1 (MQ=255)
       gTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCg                                      <  1:88688/69‑1 (MQ=255)
        ttGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCg                                      <  1:1997914/68‑1 (MQ=255)
               gCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCa                              >  1:1554039/1‑69 (MQ=255)
               gCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCa                              >  1:3247765/1‑69 (MQ=255)
                         gATGCCTGTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:2107983/1‑69 (MQ=255)
                         gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:2970330/1‑69 (MQ=255)
                         gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:2613825/1‑69 (MQ=255)
                         gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:3402224/1‑69 (MQ=255)
                         gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:2299671/1‑69 (MQ=255)
                         gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:2154319/1‑69 (MQ=255)
                         gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:2065727/1‑69 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGaaagcaaaagca                                               >  1:2273610/1‑39 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGAAAGCATAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                  >  1:14000/1‑68 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                    >  1:192252/1‑66 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                  >  1:3000943/1‑68 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                  >  1:2891870/1‑68 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                  >  1:1043949/1‑68 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                  >  1:1628393/1‑68 (MQ=255)
                            gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                  >  1:3226815/1‑68 (MQ=255)
                               tctACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa               <  1:796414/68‑1 (MQ=255)
                                ctACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGaa              >  1:2392442/1‑68 (MQ=255)
                                   ccTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAAt            <  1:2246348/67‑1 (MQ=255)
                                   ccTTTGACGGCATAATACAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAAt            <  1:74027/67‑1 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAgg                            >  1:1143033/1‑47 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:1738066/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:1126371/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:409635/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:777412/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:2955928/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:294075/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:1263710/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:958563/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:2583818/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:2544812/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:2484610/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:1979962/1‑48 (MQ=255)
                                       tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                           >  1:1997681/1‑48 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCGGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:2636444/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:2726488/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:529730/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:1141863/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:2878484/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:1648286/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:1459944/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:1616533/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:1831163/67‑1 (MQ=255)
                                           gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt  <  1:1644402/67‑1 (MQ=255)
                                                |                                                               
GAGTGAAGTTGCGGTGCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGAT  >  minE/942030‑942141

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: