Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 227,569 G→A 24.4% A121A (GCG→GCA crl → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE227,5690GA24.4% 64.2 / 15.5 53A121A (GCG→GCAcrlDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (32/8);  new base A (2/11);  total (34/19)
Fisher's exact test for biased strand distribution p-value = 5.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGCGTGAGTGAAATGTGCCGGA  >  minE/227506‑227623
                                                               |                                                      
gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGATTCTGAAGCTGGAACCGGCGgatga                                                   >  1:1552064/1‑69 (MQ=255)
gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                                                   >  1:1104842/1‑69 (MQ=255)
gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                                                   >  1:2785851/1‑69 (MQ=255)
gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                                                   >  1:2500118/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:2465730/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:716738/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:659103/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:2581773/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:2402411/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:1760950/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:1607577/1‑69 (MQ=255)
       tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                                            >  1:1587135/1‑69 (MQ=255)
          cACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAGCTGGAACCGGCAGATGATTTTCGTGa                                          >  1:1628862/1‑68 (MQ=255)
          cACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAGCTGGAACCGGCAGATGATTTTCGTGa                                          >  1:3347633/1‑68 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:2166785/1‑67 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:141730/1‑67 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:2646620/1‑67 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:1467696/1‑67 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:3123686/1‑67 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:2056010/1‑67 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:512026/1‑67 (MQ=255)
                   cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg                                  >  1:164187/1‑67 (MQ=255)
                     gCGTGAACTGCTGACGACGCTGAATGTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGaa                              >  1:2398965/1‑69 (MQ=255)
                     gCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGaa                              >  1:1857766/1‑69 (MQ=255)
                        tGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt                           <  1:401956/69‑1 (MQ=255)
                        tGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt                           <  1:1470458/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:347347/1‑68 (MQ=255)
                               ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc                     >  1:985672/1‑68 (MQ=255)
                               ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc                     >  1:2919848/1‑68 (MQ=255)
                               ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc                     >  1:2331780/1‑68 (MQ=255)
                               ctgACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGc                     >  1:1691162/1‑68 (MQ=255)
                                  acgacgCTGAATCTGCAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGtt                           <  1:2023072/59‑1 (MQ=255)
                                             tCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTACCGGCGTGAGTGAAATGTGc      >  1:1056970/1‑69 (MQ=255)
                                             tCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGCg                    >  1:2054556/1‑55 (MQ=255)
                                             tCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGCGTGAGTGAAATGTGc      >  1:2690388/1‑69 (MQ=255)
                                             tCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGCGTGAGTGAAATGTGc      >  1:634754/1‑69 (MQ=255)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:364027/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:858964/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:81105/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:1258404/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:1963546/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:2751547/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:3114617/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:2999020/69‑1 (MQ=38)
                                             tCTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:2675119/69‑1 (MQ=38)
                                              cTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:542759/68‑1 (MQ=38)
                                              cTGAAGCTGGAACCGGCAGATGATTTTCGTGATGAACCGGTGAAGTTAACGGCGTGAATGAAATGTGc      <  1:1173121/68‑1 (MQ=38)
                                                   gCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGCGTGAGTGAAATGTGCCGGa  <  1:3158986/67‑1 (MQ=255)
                                                   gCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGCGTGAGTGAAATGTGCCGGa  <  1:1656925/67‑1 (MQ=255)
                                                               |                                                      
GCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTTAACGGCGTGAGTGAAATGTGCCGGA  >  minE/227506‑227623

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: