Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 227,533 A→G 13.5% E109E (GAA→GAG crl → DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE227,5330AG13.5% 127.5 / 8.8 73E109E (GAA→GAGcrlDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (58/5);  new base G (3/7);  total (61/12)
Fisher's exact test for biased strand distribution p-value = 4.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

AGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTT  >  minE/227470‑227598
                                                               |                                                                 
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                                                              >  1:1052905/1‑69 (MQ=255)
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                                                              >  1:31800/1‑69 (MQ=255)
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                                                              >  1:3434018/1‑69 (MQ=255)
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                                                              >  1:1688300/1‑69 (MQ=255)
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                                                              >  1:3172577/1‑69 (MQ=255)
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                                                              >  1:864117/1‑69 (MQ=255)
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgct                                                              >  1:2736199/1‑69 (MQ=255)
aGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgc                                                               >  1:2987621/1‑68 (MQ=255)
 gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                                                             >  1:3028918/1‑69 (MQ=255)
 gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                                                             >  1:303026/1‑69 (MQ=255)
 gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                                                             >  1:481535/1‑69 (MQ=255)
 gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                                                             >  1:2099838/1‑69 (MQ=255)
 gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                                                             >  1:436515/1‑69 (MQ=255)
 gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                                                             >  1:1665575/1‑69 (MQ=255)
 gACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAActgctg                                                             >  1:920498/1‑69 (MQ=255)
   cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                           >  1:3035679/1‑69 (MQ=255)
   cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                           >  1:2958733/1‑69 (MQ=255)
   cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                           >  1:2059207/1‑69 (MQ=255)
   cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                           >  1:2705043/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:1334050/1‑69 (MQ=255)
   cacTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTgac                                                           >  1:116082/1‑69 (MQ=255)
       gAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGacgacg                                                       >  1:665101/1‑69 (MQ=255)
       gAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGacgacg                                                       >  1:2585842/1‑69 (MQ=255)
            ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa                                                  >  1:2375904/1‑69 (MQ=255)
            ggTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGaa                                                  >  1:2990706/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:924377/69‑1 (MQ=255)
                   cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                           <  1:329984/69‑1 (MQ=255)
                   cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                           <  1:2940824/69‑1 (MQ=255)
                   cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                           <  1:596310/69‑1 (MQ=255)
                   cGATTGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                           <  1:2715066/69‑1 (MQ=255)
                    gattGGAGCACACCCTGCGTGAATTTCACGAGAAGCTGCGTGAGCTGCTGGCGACGCTGAATCTGAAg                                           <  1:3125062/65‑1 (MQ=255)
                              cacCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                 >  1:1002125/1‑68 (MQ=255)
                              cacCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                 >  1:3325121/1‑68 (MQ=255)
                              cacCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                 >  1:210194/1‑68 (MQ=255)
                              cacCCTACGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACcgg                                 >  1:5030/1‑68 (MQ=255)
                                    gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGATTCTGAAGCTGGAACCGGCGgatga                          >  1:1552064/1‑69 (MQ=255)
                                    gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                          >  1:2785851/1‑69 (MQ=255)
                                    gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                          >  1:1104842/1‑69 (MQ=255)
                                    gCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGgatga                          >  1:2500118/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGaa                                     >  1:3387593/1‑51 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:1760950/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:1607577/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:1587135/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:659103/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:2402411/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:2465730/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:716738/1‑69 (MQ=255)
                                           tttCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGt                   >  1:2581773/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:164187/1‑67 (MQ=255)
                                                       cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg         >  1:2056010/1‑67 (MQ=255)
                                                       cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg         >  1:512026/1‑67 (MQ=255)
                                                       cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg         >  1:3123686/1‑67 (MQ=255)
                                                       cTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCgg         >  1:2166785/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: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)
                                                               |                                                                 
AGACACTGAAGTGGTTGAACGACTGGAGCACACCCTGCGTGAGTTTCACGAGAAGCTGCGTGAACTGCTGACGACGCTGAATCTGAAGCTGGAACCGGCGGATGATTTTCGTGACGAGCCGGTGAAGTT  >  minE/227470‑227598

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: