Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,294,657 A→G 33.5% S196S (AGT→AGC nac ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,294,6570AG33.5% 55.8 / 39.7 77S196S (AGT→AGCnacDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (34/16);  new base G (13/12);  total (48/29)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

TTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATCAACGCTTT  >  minE/1294590‑1294715
                                                                   |                                                          
tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                                                           <  1:242047/69‑1 (MQ=255)
tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                                                           <  1:2462302/69‑1 (MQ=255)
tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                                                           <  1:2904942/69‑1 (MQ=255)
    aCCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGt                                                        <  1:2770705/68‑1 (MQ=255)
    aCCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGt                                                        <  1:1560295/68‑1 (MQ=255)
       aaTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                                                    <  1:2582921/69‑1 (MQ=255)
       aaTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                                                    <  1:3008200/69‑1 (MQ=255)
        aTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                                                    <  1:3279490/68‑1 (MQ=255)
           aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                                                       >  1:2392970/1‑62 (MQ=255)
           aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                                                       >  1:2288677/1‑62 (MQ=255)
           aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                                                       >  1:1292873/1‑62 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:80109/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1837619/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:315232/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:622801/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1148902/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:3314151/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1239962/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:429097/67‑1 (MQ=255)
             ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:508264/67‑1 (MQ=255)
             ttttGCCGTGAGCAGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTCACAGCGCTGTAATCACTg                                                <  1:193745/67‑1 (MQ=255)
             ttttGCCGAGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:1413046/67‑1 (MQ=255)
             ctttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg                                                <  1:2061566/66‑1 (MQ=255)
                     tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:1540754/1‑67 (MQ=255)
                     tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:509537/1‑67 (MQ=255)
                     tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:2658558/1‑67 (MQ=255)
                     tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:2829861/1‑67 (MQ=255)
                     tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:3129163/1‑67 (MQ=255)
                     tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:1410566/1‑67 (MQ=255)
                     tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa                                      >  1:1667411/1‑67 (MQ=255)
                                   gAAAAAGCCTCATCAACACGAAGTCTAATAGCTCTGTAATCACTGGGGAGAAAGAGGTTCATTTgcgc                         <  1:200157/68‑1 (MQ=255)
                                   gAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTgcgc                         <  1:581660/68‑1 (MQ=255)
                                                caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:326691/1‑68 (MQ=255)
                                                caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:979945/1‑68 (MQ=255)
                                                caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:224224/1‑68 (MQ=255)
                                                caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:135880/1‑68 (MQ=255)
                                                caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:1094098/1‑68 (MQ=255)
                                                caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTcaca            >  1:1141663/1‑68 (MQ=255)
                                                caACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGATCATTTGCGCAATAGCATTcac             >  1:2908018/1‑67 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:3147351/1‑37 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:1741398/1‑37 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:2288255/1‑37 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:1399485/1‑37 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:2781195/1‑37 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGgaga                                          >  1:1179890/1‑37 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAgg                                    >  1:2883717/1‑43 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCCCATc          >  1:1026622/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1175499/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:1815171/69‑1 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:514258/69‑1 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3385205/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1088258/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3308906/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:932475/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:3226331/69‑1 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:960593/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2986692/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:2321251/69‑1 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:2440204/69‑1 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2808529/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1990888/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2877775/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:1938604/69‑1 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:2322847/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:1932477/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3055528/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3074455/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          <  1:3074761/69‑1 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:3085988/1‑69 (MQ=255)
                                                 aaCACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATc          >  1:311959/1‑69 (MQ=255)
                                                         gTCTAACAGCTCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:1194068/1‑69 (MQ=255)
                                                         gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:561007/1‑69 (MQ=255)
                                                         gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:1431482/1‑69 (MQ=255)
                                                         gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:612716/1‑69 (MQ=255)
                                                         gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:3227566/1‑69 (MQ=255)
                                                         gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:2058155/1‑69 (MQ=255)
                                                         gTCTAACAGCGCTGTAATCACTGGGGAGAAAAAGGTTCATTTGCGCAATAGCATTAACATCAACGCttt  >  1:2281485/1‑69 (MQ=255)
                                                                   |                                                          
TTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGAGGTTCATTTGCGCAATAGCATTCACATCAACGCTTT  >  minE/1294590‑1294715

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: