Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 796,568 T→C 35.2% Q58Q (CAA→CAG ycgR ← protein involved in flagellar function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE796,5680TC35.2% 57.4 / 51.4 88Q58Q (CAA→CAGycgRprotein involved in flagellar function
Reads supporting (aligned to +/- strand):  ref base T (33/24);  new base C (13/18);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTG  >  minE/796503‑796630
                                                                 |                                                              
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCtt                                                               <  1:2955567/67‑1 (MQ=255)
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCtt                                                               <  1:2659726/67‑1 (MQ=255)
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCtt                                                               <  1:1852368/67‑1 (MQ=255)
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCtt                                                               <  1:1508015/67‑1 (MQ=255)
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCtt                                                               <  1:55212/67‑1 (MQ=255)
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTg                                                              <  1:392570/68‑1 (MQ=255)
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTg                                                              <  1:1549868/68‑1 (MQ=255)
gCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTg                                                              <  1:3207116/68‑1 (MQ=255)
  aCCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACt                                                           >  1:3414654/1‑69 (MQ=255)
  aCCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACt                                                           >  1:265118/1‑69 (MQ=255)
           ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:737545/69‑1 (MQ=255)
           ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:2331606/69‑1 (MQ=255)
           ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:2973364/69‑1 (MQ=255)
           ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:1591505/69‑1 (MQ=255)
           ttCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:41788/69‑1 (MQ=255)
            tCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATc                                                  <  1:2559930/68‑1 (MQ=255)
                    aatggtaatgTGCTTTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:611662/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)
                    aatggtaatgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:2329259/69‑1 (MQ=255)
                    aatggtaatgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAg                                         <  1:892774/69‑1 (MQ=255)
                           atgTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCtt                                                               <  1:2782111/40‑1 (MQ=255)
                             gtgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCgg                                <  1:2762299/69‑1 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2698652/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2433638/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2441536/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:2607288/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:1002754/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:460908/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:539980/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:1454315/1‑69 (MQ=255)
                              tgctgtgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCggg                               >  1:1243786/1‑69 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATc                                   >  1:2358138/1‑65 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:3400566/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:371401/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:44054/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:670775/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:778469/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:782360/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:919085/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:2556711/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:1102682/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:2279073/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:1280398/1‑68 (MQ=255)
                              tgctgcgcCTTTAGCACGGCGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCgga                               >  1:1250888/1‑68 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:2204236/1‑69 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:1731973/1‑69 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:2576353/1‑69 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:3348212/1‑69 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:3238713/1‑69 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:787873/1‑69 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:1007158/1‑69 (MQ=255)
                                   tgcCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTAt                          >  1:911180/1‑69 (MQ=255)
                                     ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCgggg                              >  1:2281006/1‑63 (MQ=255)
                                     ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:1514593/1‑65 (MQ=255)
                                     ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:2085166/1‑65 (MQ=255)
                                     ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:3368666/1‑65 (MQ=255)
                                     ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:2605214/1‑65 (MQ=255)
                                     ccTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGtt                            >  1:2761259/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:1346802/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2006174/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:631069/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:1479420/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:15145/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:1664763/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:38060/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2598855/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2260856/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:3211947/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:271670/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2220600/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2974818/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:3426847/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTTGc            <  1:2818713/69‑1 (MQ=255)
                                                 cGATGTTGTCTTCGGCCTGGCTGCCGAAATCCAGCACCAGTTTATCCGGAGTTATTGCCAGTATTTGGc            <  1:1464929/69‑1 (MQ=255)
                                                   aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:329187/69‑1 (MQ=255)
                                                   aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:1183170/69‑1 (MQ=255)
                                                   aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:766765/69‑1 (MQ=255)
                                                   aTGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          <  1:2099421/69‑1 (MQ=255)
                                                     gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCt           >  1:644168/1‑66 (MQ=255)
                                                     gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:701318/1‑67 (MQ=255)
                                                     gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:3053846/1‑67 (MQ=255)
                                                     gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:2226759/1‑67 (MQ=255)
                                                     gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:3122782/1‑67 (MQ=255)
                                                     gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:1536530/1‑67 (MQ=255)
                                                     gttgtCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTg          >  1:1975355/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)
                                                                 |                                                              
GCACCCTGAGTTTCGGCGGTAATGGTAATGTGCTGTGCCTTTAGCACGGCGATGTTGTCTTCGGCTTGACTGCCGAAATCCAGCACCAGTTTATCCGGGGTTATTGCCAGTAATTTGCTGATCAGCTG  >  minE/796503‑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: