Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 942,108 C→T 33.6% V127V (GTC→GTT osmC → osmotically inducible, stress‑inducible membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE942,1080CT33.6% 64.8 / 40.8 80V127V (GTC→GTTosmCosmotically inducible, stress‑inducible membrane protein
Reads supporting (aligned to +/- strand):  ref base C (43/10);  new base T (14/13);  total (57/23)
Fisher's exact test for biased strand distribution p-value = 9.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGCATTGCCGGATGACG  >  minE/942045‑942176
                                                               |                                                                    
gCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCa                                                                 >  1:3247765/1‑69 (MQ=255)
gCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCa                                                                 >  1:1554039/1‑69 (MQ=255)
          gATGCCTGTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                                                       >  1:2107983/1‑69 (MQ=255)
          gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                                                       >  1:2065727/1‑69 (MQ=255)
          gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                                                       >  1:2613825/1‑69 (MQ=255)
          gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                                                       >  1:2299671/1‑69 (MQ=255)
          gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                                                       >  1:3402224/1‑69 (MQ=255)
          gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                                                       >  1:2154319/1‑69 (MQ=255)
          gATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGaaa                                                       >  1:2970330/1‑69 (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:1628393/1‑68 (MQ=255)
             gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                                                     >  1:3226815/1‑68 (MQ=255)
             gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                                                     >  1:1043949/1‑68 (MQ=255)
             gCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGc                                                     >  1:2891870/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:409635/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:2544812/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:2583818/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:1997681/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:1979962/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:2484610/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:958563/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:1126371/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:777412/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:1738066/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:1263710/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:2955928/1‑48 (MQ=255)
                        tGACGTCATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGt                                                              >  1:294075/1‑48 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCGGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:2636444/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:2878484/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:1141863/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:1459944/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:1616533/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:1644402/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:1648286/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:1831163/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:529730/67‑1 (MQ=255)
                            gtcATTATCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:2726488/67‑1 (MQ=255)
                               attaTCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:3186729/63‑1 (MQ=255)
                                ttaTCCAGAAAGCAAAAGCAGGATGCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:1468020/63‑1 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa                                                  >  1:528217/1‑50 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa                                                  >  1:2099402/1‑50 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa                                                  >  1:1255292/1‑50 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa                                                  >  1:1457877/1‑50 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa                                                  >  1:3298472/1‑50 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa                                                  >  1:1021899/1‑50 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGa                                                  >  1:2927043/1‑50 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCa                                >  1:3014303/1‑68 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:657816/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:1070177/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:8947/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:2930617/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:1703167/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:176656/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:2826058/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:2804578/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:2654865/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:2053298/1‑69 (MQ=255)
                                  aTCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAg                               >  1:2151546/1‑69 (MQ=255)
                                          agcaaaagcaTGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTg                                        >  1:2777378/1‑52 (MQ=255)
                                          agcaaaagcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTg                                        >  1:589608/1‑52 (MQ=255)
                                          agcaaaagcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTg                                        >  1:2269081/1‑52 (MQ=255)
                                          agcaaaagcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTg                                        >  1:1565959/1‑52 (MQ=255)
                                          agcaaaagcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTg                                        >  1:1444445/1‑52 (MQ=255)
                                                agcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGAtt                                    <  1:1883774/50‑1 (MQ=255)
                                                 gcaGGATGCCCGGTCTCTCAGGTTCTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGc                <  1:1838105/69‑1 (MQ=255)
                                                 gcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGc                <  1:3072711/69‑1 (MQ=255)
                                                 gcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGc                <  1:2915870/69‑1 (MQ=255)
                                                 gcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGc                <  1:1944080/69‑1 (MQ=255)
                                                 gcaGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGc                <  1:1795435/69‑1 (MQ=255)
                                                     gATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGCATTg            >  1:2808290/1‑69 (MQ=255)
                                                     gATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGCATTg            >  1:2212418/1‑69 (MQ=255)
                                                        gCCCGGTTTCTCAGGTACTGAAAGCGGAAATTACGCTGGAt                                     <  1:2534695/41‑1 (MQ=255)
                                                             gTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGc                <  1:2227487/57‑1 (MQ=255)
                                                               ctctcAGGTACTGGAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGCATTGCCGGATGAcg  >  1:1845513/1‑69 (MQ=255)
                                                               ctctcAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGCATTGCCGGATGAcg  >  1:1032638/1‑69 (MQ=255)
                                                               |                                                                    
GCCGGGTATTGATGCCTCTACCTTTGACGGCATAATCCAGAAAGCAAAAGCAGGATGCCCGGTCTCTCAGGTACTGAAAGCGGAAATTACGCTGGATTACCAGTTGAAATCGTAAAGCATTGCCGGATGACG  >  minE/942045‑942176

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: