Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 17,531 A→G 19.0% E92E (GAA→GAG nhaR → DNA‑binding transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE17,5310AG19.0% 181.0 / 29.9 120E92E (GAA→GAGnhaRDNA‑binding transcriptional activator
Reads supporting (aligned to +/- strand):  ref base A (60/37);  new base G (7/16);  total (67/53)
Fisher's exact test for biased strand distribution p-value = 9.47e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 122 total aligned reads displayed.

TATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCG  >  minE/17463‑17592
                                                                    |                                                             
tATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGaa                                                               >  1:1303360/1‑69 (MQ=255)
    gCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCa                                                           >  1:904486/1‑69 (MQ=255)
    gCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCa                                                           >  1:759494/1‑69 (MQ=255)
    gCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCa                                                           >  1:654680/1‑69 (MQ=255)
    gCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCa                                                           >  1:1321340/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:449/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:2495164/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:2440175/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:2745596/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:282697/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:321514/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:3221010/1‑69 (MQ=255)
        tGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAAttt                                                       >  1:3348601/1‑69 (MQ=255)
               aaaaTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGttt                                                <  1:581186/69‑1 (MQ=255)
               aaaaTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGttt                                                <  1:911872/69‑1 (MQ=255)
               aaaaTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGttt                                                <  1:2421095/69‑1 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:982055/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:296940/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:1245176/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:3369960/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:2731951/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:1791557/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:2691554/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:450359/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:454430/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:10182/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:1592777/1‑69 (MQ=255)
                  aTGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGAc                                             >  1:2523716/1‑69 (MQ=255)
                     ttCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGt                                                  <  1:1163663/61‑1 (MQ=255)
                     ttCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGt                                                  <  1:1380311/61‑1 (MQ=255)
                     ttCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGt                                                  <  1:3151314/61‑1 (MQ=255)
                     ttCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGt                                                  <  1:3213289/61‑1 (MQ=255)
                     ttCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGt                                                  <  1:473756/61‑1 (MQ=255)
                     ttCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGt                                                  <  1:1514369/61‑1 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:752650/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:3102917/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:1344/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:1352163/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        <  1:2812594/69‑1 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:237336/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        <  1:351558/69‑1 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:241261/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:4963/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:534776/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:1358114/1‑69 (MQ=255)
                       cACCTTAAGCCAGGAAATGCAGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTgg                                        >  1:1868314/1‑69 (MQ=255)
                             aaGCCAGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTATTGTTTGATGTTggtg                                      >  1:2619558/1‑63 (MQ=255)
                             aaGCCAGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTATTGTTTGATGTTGgtgt                                     >  1:2420452/1‑66 (MQ=255)
                             aaGCCAGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTATTGTTTGATGTTGgtgt                                     >  1:2369861/1‑66 (MQ=255)
                             aaGCCAGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTATTGTTTGATGTTGgtgt                                     >  1:1649069/1‑66 (MQ=255)
                             aaGCCAGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTATTGTTTGATGTTGgtgt                                     >  1:2783050/1‑66 (MQ=255)
                             aaGCCAGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTATTGTTTGATGTTGgtgt                                     >  1:1600373/1‑66 (MQ=255)
                             aaGCCAGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTATTGTTTGATGTTGgtgt                                     >  1:1442770/1‑66 (MQ=255)
                               gCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCt                                 <  1:106029/68‑1 (MQ=255)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:2238413/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:3417086/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:3389723/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:1846245/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:1764042/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:1928140/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:514980/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:3138493/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:3087142/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:591681/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:662040/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:125542/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:818491/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:876141/44‑1 (MQ=25)
                                  aGGAAATGCTGGATATCGTCAACTATCGCAAAGAGTCCAACTTa                                                      <  1:2437220/44‑1 (MQ=25)
                                   ggAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTTATGc                            <  1:933011/69‑1 (MQ=255)
                                   ggAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGc                            <  1:3120778/69‑1 (MQ=255)
                                   ggAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGc                            <  1:2941013/69‑1 (MQ=255)
                                   ggAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGc                            <  1:2802791/69‑1 (MQ=255)
                                   ggAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGc                            <  1:1796518/69‑1 (MQ=255)
                                   ggAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGc                            <  1:929411/69‑1 (MQ=255)
                                    gAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGc                            <  1:1047464/68‑1 (MQ=255)
                                            ggATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCa                    <  1:2400306/68‑1 (MQ=255)
                                            ggATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCa                    <  1:1924973/68‑1 (MQ=255)
                                            ggATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCa                    <  1:2374750/68‑1 (MQ=255)
                                            ggATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCa                    <  1:2740315/68‑1 (MQ=255)
                                            ggATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCa                    <  1:2244093/68‑1 (MQ=255)
                                            ggATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCa                    <  1:2248130/68‑1 (MQ=255)
                                                   gtgAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCt                                 >  1:1698521/1‑48 (MQ=255)
                                                   gtgAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCt                                 >  1:1404511/1‑48 (MQ=255)
                                                   gtgAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCt                                 >  1:1273806/1‑48 (MQ=255)
                                                   gtgAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCt             >  1:640259/1‑68 (MQ=255)
                                                   gtgAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTg            >  1:10599/1‑69 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGta         >  1:2213903/1‑66 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:2560332/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:2780630/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:770103/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:1068474/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:2139345/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:1472617/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:2001375/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTc         >  1:1903329/1‑67 (MQ=255)
                                                        cTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAg       <  1:2596029/69‑1 (MQ=255)
                                                             gCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCg  <  1:3063286/69‑1 (MQ=255)
                                                             gCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCg  <  1:3426918/69‑1 (MQ=255)
                                                                    aTCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCa                    <  1:3350203/44‑1 (MQ=255)
                                                                    |                                                             
TATCGCTATGCCGATAAAATGTTCACCTTAAGCCAGGAAATGCTGGATATTGTGAACTATCGCAAAGAATCCAATTTATTGTTTGACGTTGGCGTGGCTGATGCACTTTCCAAACGCCTGGTCAGTAGCG  >  minE/17463‑17592

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