Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 805,245 T→C 9.7% intergenic (+46/+54) dhaR → / ← ycgV predicted DNA‑binding transcriptional regulator, dihydroxyacetone/predicted adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE805,2450TC 9.7% 109.5 / 3.7 41intergenic (+46/+54)dhaR/ycgVpredicted DNA‑binding transcriptional regulator, dihydroxyacetone/predicted adhesin
Reads supporting (aligned to +/- strand):  ref base T (16/21);  new base C (2/2);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCG  >  minE/805177‑805313
                                                                    |                                                                    
ggCAGTTTAAGCGCCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTcc                                                                    <  1:1547581/71‑1 (MQ=255)
ggCAGTTTAAGCGCCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTc                                                                     <  1:2847842/70‑1 (MQ=255)
   aGTTTAAGCGCCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTc                                                                     >  1:311036/1‑67 (MQ=255)
   aGTTTAAGCGCCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGgc                                                                 >  1:2478593/1‑71 (MQ=255)
         aGCGCCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCgg                                                                  >  1:2079058/1‑64 (MQ=255)
          gcgcCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCt                                                          <  1:2718482/71‑1 (MQ=255)
             ccGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCAtt                                                                      >  1:3414283/1‑56 (MQ=255)
             ccGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTaa                                                        >  1:1186053/1‑70 (MQ=255)
              cGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCg                                                                   <  1:920204/58‑1 (MQ=255)
                ggTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAAt                                                    >  1:2592746/1‑71 (MQ=255)
                  tATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAAttt                                                  <  1:1216201/71‑1 (MQ=255)
                    tGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAAtttt                                                 <  1:3679231/70‑1 (MQ=255)
                        aGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGc                                                           >  1:785004/1‑56 (MQ=255)
                        aGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCa                                            >  1:2723250/1‑71 (MQ=255)
                         gACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCAt                                           <  1:1624523/71‑1 (MQ=255)
                           cAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCg                                         >  1:3553938/1‑71 (MQ=255)
                           cAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATCCCGGCGCACGCTAACAAt                                                    >  1:2249574/1‑60 (MQ=255)
                            aGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTc                                                <  1:2685844/63‑1 (MQ=255)
                             gAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATCCCGGCGCACGCTAACAATTTTCAGCATCGtt                                       <  1:519317/71‑1 (MQ=255)
                              aaaCGATTTCTGATACATCAGAGTGATCTGTATTTCATCCCGGCGCACGCTAACaa                                                     >  1:1848209/1‑56 (MQ=255)
                              aaaCGAATTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAg                                              >  1:513614/1‑63 (MQ=255)
                               aaCGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTa                                     <  1:3860467/71‑1 (MQ=255)
                                 cGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTaa                                    <  1:1544203/70‑1 (MQ=255)
                                 cGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAg                                   <  1:2182893/71‑1 (MQ=255)
                                 cGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAg                                   <  1:3423165/71‑1 (MQ=255)
                                 cGATTTCTGATACATCAGAGTGATCTGTATTTCATCCCGGCGCACGCTAACAATTTTCAGCATCGttt                                      <  1:2438540/68‑1 (MQ=255)
                                     ttCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCt                               >  1:1106355/1‑71 (MQ=255)
                                         gATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTc                           >  1:1218459/1‑71 (MQ=255)
                                           tACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACaa                                                     <  1:400324/43‑1 (MQ=255)
                                             cATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGc                                             >  1:1700533/1‑49 (MQ=255)
                                              aTCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATcc                      <  1:1815287/71‑1 (MQ=255)
                                               tCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGc                                             <  1:3248879/47‑1 (MQ=255)
                                                  gagTGATCTGTATTTCATTCCGGCGCACGCTAACaa                                                     >  1:2744192/1‑36 (MQ=255)
                                                   agTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCAc                 <  1:3821612/71‑1 (MQ=255)
                                                         cTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTaa                                    >  1:1561847/1‑46 (MQ=255)
                                                         cTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTaa                                    >  1:292027/1‑46 (MQ=255)
                                                            tATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTaa             <  1:1732625/66‑1 (MQ=255)
                                                            tATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCATTTAAGGGCTTGTCTATCCCGCACTTAAAAg          <  1:2636538/69‑1 (MQ=255)
                                                              tttCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTg       >  1:2210476/1‑70 (MQ=255)
                                                              tttCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGa      <  1:2538620/71‑1 (MQ=255)
                                                                 cATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGc                   <  1:1988599/55‑1 (MQ=255)
                                                                 cATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAg    <  1:3448391/70‑1 (MQ=255)
                                                                   ttCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCt        <  1:2434040/64‑1 (MQ=255)
                                                                   ttCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAg    >  1:250409/1‑68 (MQ=255)
                                                                   ttCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCg  <  1:2057108/70‑1 (MQ=255)
                                                                    |                                                                    
GGCAGTTTAAGCGCCGGGTATGAAAGACAGAAACGATTTCTGATACATCAGAGTGATCTGTATTTCATTCCGGCGCACGCTAACAATTTTCAGCATCGTTTAAGGGCTTGTCTATCCCGCACTTAAAAGCTGAAGCG  >  minE/805177‑805313

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