Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,720,956 G→A 17.8% H178H (CAC→CAT yfhA ← predicted DNA‑binding response regulator in two‑component system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,720,9560GA17.8% 58.2 / 7.2 28H178H (CAC→CATyfhApredicted DNA‑binding response regulator in two‑component system
Reads supporting (aligned to +/- strand):  ref base G (8/15);  new base A (2/3);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

AGCAATTGCTCGGGCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCGTTAATCAAAACGCTGACGTCTGA  >  minE/1720888‑1721024
                                                                    |                                                                    
aGCAATTGCTCGGGCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTtgtg                                                                    >  1:1671051/1‑71 (MQ=255)
 gCAATTGCTCTGGCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTATgg                                                                   <  1:3455550/71‑1 (MQ=255)
  cAATTGCTCGGGCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGa                                                                  <  1:2087760/71‑1 (MQ=255)
       gCTCGGGCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGc                                                              <  1:2212707/70‑1 (MQ=255)
        cTCGGGCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGcc                                                             <  1:2766054/70‑1 (MQ=255)
           gggCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTg                                                           <  1:3326045/69‑1 (MQ=255)
           gggCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTg                                                           <  1:2494214/69‑1 (MQ=255)
           gggCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTg                                                           <  1:3930426/69‑1 (MQ=255)
            ggCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGc                                                        >  1:2427582/1‑71 (MQ=255)
              cAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCg                                                       >  1:852124/1‑70 (MQ=255)
                 tGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGaaaa                                                   >  1:2747360/1‑71 (MQ=255)
                 tGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGaaa                                                    >  1:3056527/1‑70 (MQ=255)
                   cgcCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTtgt                                                                     >  1:1122353/1‑51 (MQ=255)
                   cgcCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTtgt                                                                     >  1:2345403/1‑51 (MQ=255)
                    gcCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAAtct                                                <  1:550442/71‑1 (MQ=255)
                     ccACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGa                                                      <  1:2782887/64‑1 (MQ=255)
                       acaGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCtt                                             >  1:1170876/1‑71 (MQ=255)
                       acaGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTATGGATAGCCTGGGCGAAAATCTCtt                                             >  1:1053973/1‑71 (MQ=255)
                        caGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCttt                                            <  1:2302133/71‑1 (MQ=255)
                           ttAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTg                                                           <  1:1683393/53‑1 (MQ=255)
                           ttAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTATGGATAGCCTGGGCGAAAATCTCTTTccg                                         <  1:213230/71‑2 (MQ=255)
                             aaTAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTc                                           >  1:879842/1‑67 (MQ=255)
                             aaTAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTATGGATAGCCTGGGCGAAAATCTCTTTCCCg                                        <  1:633212/70‑1 (MQ=255)
                                    taaaaTGGTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGaaaa                                                   <  1:1354685/49‑1 (MQ=255)
                                       aaTGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCt                                                            >  1:3707849/1‑40 (MQ=255)
                                              tGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGctct                              >  1:2126804/1‑63 (MQ=255)
                                                cTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCGTTaa                     <  1:2584206/70‑1 (MQ=255)
                                                   ttGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCg                         <  1:3826467/63‑1 (MQ=255)
                                                   ttGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCg                         <  1:3830962/63‑1 (MQ=255)
                                                        cgGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCGTTAATCAAAACGc            <  1:3452738/71‑1 (MQ=255)
                                                          ggCTGGCGTTATGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCGTTAATCAAAACGCt           >  1:2160308/1‑70 (MQ=255)
                                                                   tgtgGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCGTTAATCAAAACGCTGACGTCTGa  >  1:2700444/1‑70 (MQ=255)
                                                                    gtgGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCGTTAATCaaaa               <  1:2366512/56‑1 (MQ=255)
                                                                    |                                                                    
AGCAATTGCTCGGGCAATGCGCCACAGTTAATAGCAATAAATGGTTTGCTGTTGCGCGGGCTGGCGTTGTGGATAGCCTGGGCGAAAATCTCTTTCCCGGTGCCGCTCTGACCGTTAATCAAAACGCTGACGTCTGA  >  minE/1720888‑1721024

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: