Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 796,571 A→G 36.4% S57S (AGT→AGC ycgR ← protein involved in flagellar function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE796,5710AG36.4% 49.1 / 51.6 85S57S (AGT→AGCycgRprotein involved in flagellar function
Reads supporting (aligned to +/- strand):  ref base A (33/21);  new base G (13/18);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

ACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTG  >  minE/796505‑796630
                                                                  |                                                           
aCCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACt                                                           >  1:265118/1‑69 (MQ=255)
aCCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACt                                                           >  1:3414654/1‑69 (MQ=255)
         ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:41788/69‑1 (MQ=255)
         ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:1591505/69‑1 (MQ=255)
         ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:737545/69‑1 (MQ=255)
         ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:2973364/69‑1 (MQ=255)
         ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:2331606/69‑1 (MQ=255)
          tCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:2559930/68‑1 (MQ=255)
                  aatggtaatgTGCTTTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:611662/69‑1 (MQ=255)
                  aatggtaatgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:892774/69‑1 (MQ=255)
                  aatggtaatgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:2329259/69‑1 (MQ=255)
                  aatggtaatgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:2705993/69‑1 (MQ=255)
                  aatggtaatgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:2587772/69‑1 (MQ=255)
                  aatggtaatgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:1693466/69‑1 (MQ=255)
                           gtgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCgg                                <  1:2762299/69‑1 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2441536/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2433638/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2698652/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2607288/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:1243786/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:539980/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:460908/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:1002754/1‑69 (MQ=255)
                            tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:1454315/1‑69 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATc                                   >  1:2358138/1‑65 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:778469/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:1250888/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:2556711/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:1280398/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:919085/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:1102682/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:782360/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:2279073/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:3400566/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:670775/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:371401/1‑68 (MQ=255)
                            tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:44054/1‑68 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:1007158/1‑69 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:911180/1‑69 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:787873/1‑69 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:3238713/1‑69 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:3348212/1‑69 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:2204236/1‑69 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:2576353/1‑69 (MQ=255)
                                 tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:1731973/1‑69 (MQ=255)
                                   ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCgggg                              >  1:2281006/1‑63 (MQ=255)
                                   ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:2085166/1‑65 (MQ=255)
                                   ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:3368666/1‑65 (MQ=255)
                                   ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:1514593/1‑65 (MQ=255)
                                   ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:2761259/1‑65 (MQ=255)
                                   ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:2605214/1‑65 (MQ=255)
                                   ccTTTAGCAAGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGtt                                       >  1:831047/1‑54 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGACATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:3015201/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTCTTTTGc            <  1:438993/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2006174/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:631069/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:38060/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:3426847/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:1664763/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2220600/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:15145/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:1479420/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2260856/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:3211947/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2974818/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2818713/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:1346802/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:271670/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2598855/69‑1 (MQ=255)
                                               cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTGGc            <  1:1464929/69‑1 (MQ=255)
                                                 aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:766765/69‑1 (MQ=255)
                                                 aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:1183170/69‑1 (MQ=255)
                                                 aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:2099421/69‑1 (MQ=255)
                                                 aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:329187/69‑1 (MQ=255)
                                                   gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCt           >  1:644168/1‑66 (MQ=255)
                                                   gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:701318/1‑67 (MQ=255)
                                                   gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:1536530/1‑67 (MQ=255)
                                                   gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:1975355/1‑67 (MQ=255)
                                                   gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:2226759/1‑67 (MQ=255)
                                                   gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:3122782/1‑67 (MQ=255)
                                                   gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:3053846/1‑67 (MQ=255)
                                                         ttCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTg  <  1:171198/69‑1 (MQ=255)
                                                         ttCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTg  <  1:2513226/69‑1 (MQ=255)
                                                         ttCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTg  <  1:2625729/69‑1 (MQ=255)
                                                         ttCGGCTTGACTGCCGAAATCCAGCACCAGTCTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTg  <  1:541710/69‑1 (MQ=255)
                                                                  |                                                           
ACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTG  >  minE/796505‑796630

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: