Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,873,630 C→A 38.0% intergenic (+774/+599) ygcE → / ← ygcF predicted kinase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,873,6300CA38.0% 22.1 / 13.7 21intergenic (+774/+599)ygcE/ygcFpredicted kinase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base A (5/3);  total (12/9)
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.43e-01

GCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCG  >  minE/1873568‑1873698
                                                              |                                                                    
gCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCgg                                                              >  1:3425924/1‑71 (MQ=255)
     cgcgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAac                                                         <  1:573945/71‑1 (MQ=255)
       cgGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAacac                                                       <  1:2474164/71‑1 (MQ=255)
         gggAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGGTTATCCCCGCTGGCGCGGGGAACACgg                                                     <  1:645242/71‑1 (MQ=255)
            aaCACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGgcg                                                  >  1:414778/1‑71 (MQ=255)
              cacCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCa                                                >  1:237225/1‑71 (MQ=255)
               acCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCAc                                               >  1:3322285/1‑71 (MQ=255)
                 cTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACg                                              >  1:187127/1‑70 (MQ=255)
                            aaTCTATGGACTATTGCTATTCGGTTTATCCCCGCTg                                                                    >  1:2545313/1‑37 (MQ=255)
                            aaTCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGcc                                  >  1:515874/1‑71 (MQ=255)
                             aTCTATGGACTATTGCTATTCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCg                                 >  1:3605298/1‑71 (MQ=255)
                                    gACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCt                          <  1:956894/71‑1 (MQ=255)
                                     aCTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTg                         <  1:27050/71‑1 (MQ=255)
                                         ttGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTc                      <  1:285464/70‑1 (MQ=255)
                                             tATTCGGTTTATCCCCGATGGCGCGGGGAACACGGgcg                                                  >  1:846845/1‑38 (MQ=255)
                                               ttCGGTTTATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTAt                            >  1:1601538/1‑58 (MQ=255)
                                                   gTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATccc            >  1:2197886/1‑70 (MQ=255)
                                                    tttATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTg                         <  1:2471406/56‑1 (MQ=255)
                                                      tATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGc         <  1:2925278/70‑1 (MQ=255)
                                                      tATCCCCGATGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGc         <  1:3488144/70‑1 (MQ=255)
                                                            cGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGcgcg  >  1:3590826/1‑71 (MQ=255)
                                                              |                                                                    
GCTGGCGCGGGGAACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCG  >  minE/1873568‑1873698

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: