Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 227,540 A→G 18.7% T112A (ACG→GCG)  crl → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE227,5400AG18.7% 75.1 / 9.9 53T112A (ACG→GCG) crlDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (38/5);  new base G (3/7);  total (41/12)
Fisher's exact test for biased strand distribution p-value = 4.54e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

CACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGC  >  minE/227473‑227604
                                                                   |                                                                
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:2705043/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:2958733/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:3035679/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:2059207/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:3420277/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:455622/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:1619916/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:116082/1‑69 (MQ=255)
cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                                 >  1:1334050/1‑69 (MQ=255)
    gAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGacgacg                                                             >  1:665101/1‑69 (MQ=255)
    gAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGacgacg                                                             >  1:2585842/1‑69 (MQ=255)
         ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa                                                        >  1:2990706/1‑69 (MQ=255)
         ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa                                                        >  1:2375904/1‑69 (MQ=255)
         ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGa                                                         >  1:3228813/1‑68 (MQ=255)
         ggTTGAACGACTGGAGCACACCATGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa                                                        >  1:584141/1‑69 (MQ=255)
          gTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAAt                                                       >  1:1690027/1‑69 (MQ=255)
                cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTTCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                                 <  1:83764/69‑1 (MQ=255)
                cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                                 <  1:2715066/69‑1 (MQ=255)
                cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                                 <  1:924377/69‑1 (MQ=255)
                cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                                 <  1:329984/69‑1 (MQ=255)
                cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                                 <  1:596310/69‑1 (MQ=255)
                cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                                 <  1:2940824/69‑1 (MQ=255)
                 gattGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                                 <  1:3125062/65‑1 (MQ=255)
                           cacCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                       >  1:3325121/1‑68 (MQ=255)
                           cacCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                       >  1:210194/1‑68 (MQ=255)
                           cacCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                       >  1:1002125/1‑68 (MQ=255)
                           cacCCTACGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                       >  1:5030/1‑68 (MQ=255)
                                 gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGATTCTGAAGCTGGAACCGGCGgatga                                >  1:1552064/1‑69 (MQ=255)
                                 gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                                >  1:1104842/1‑69 (MQ=255)
                                 gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                                >  1:2500118/1‑69 (MQ=255)
                                 gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                                >  1:2785851/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGaa                                           >  1:3387593/1‑51 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:659103/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:716738/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:2581773/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:2465730/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:1760950/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:2402411/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:1607577/1‑69 (MQ=255)
                                        tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                         >  1:1587135/1‑69 (MQ=255)
                                           cACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAGCt                                               >  1:1873687/1‑44 (MQ=38)
                                           cACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAGCTGGAACCGGCAGATGATTTTCGTGa                       >  1:1628862/1‑68 (MQ=255)
                                           cACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAGCTGGAACCGGCAGATGATTTTCGTGa                       >  1:3347633/1‑68 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:2056010/1‑67 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:164187/1‑67 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:512026/1‑67 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:2166785/1‑67 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:3123686/1‑67 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:1467696/1‑67 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:141730/1‑67 (MQ=255)
                                                    cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg               >  1:2646620/1‑67 (MQ=255)
                                                      gCGTGAACTGCTGACGACGCTGAATGTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGaa           >  1:2398965/1‑69 (MQ=255)
                                                      gCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGaa           >  1:1857766/1‑69 (MQ=255)
                                                         tGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt        <  1:1470458/69‑1 (MQ=255)
                                                         tGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt        <  1:401956/69‑1 (MQ=255)
                                                         tGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt        <  1:2828128/69‑1 (MQ=255)
                                                         tGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGACGATTTTCGTGACGAGCCGGTGAAGtt        <  1:1868453/69‑1 (MQ=255)
                                                          gAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt        <  1:2580942/68‑1 (MQ=255)
                                                             ctgctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACg    >  1:2427149/1‑69 (MQ=255)
                                                                ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc  >  1:1691162/1‑68 (MQ=255)
                                                                ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc  >  1:347347/1‑68 (MQ=255)
                                                                ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc  >  1:2331780/1‑68 (MQ=255)
                                                                ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc  >  1:2919848/1‑68 (MQ=255)
                                                                ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc  >  1:985672/1‑68 (MQ=255)
                                                                   acgacgCTGAATCTGCAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt        <  1:2023072/59‑1 (MQ=255)
                                                                   |                                                                
CACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGC  >  minE/227473‑227604

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