Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,444,406 T→C 34.5% intergenic (+27/+173) rcsB → / ← rcsC DNA‑binding response regulator in two‑component regulatory system with RcsC and YojN/hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,444,4060TC34.5% 62.7 / 45.9 81intergenic (+27/+173)rcsB/rcsCDNA‑binding response regulator in two‑component regulatory system with RcsC and YojN/hybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base T (36/17);  new base C (22/6);  total (58/23)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC  >  minE/1444350‑1444468
                                                        |                                                              
gTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTa                                                       <  1:986533/66‑1 (MQ=255)
gTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTa                                                       <  1:2848618/66‑1 (MQ=255)
gTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTa                                                       <  1:403321/66‑1 (MQ=255)
gTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTa                                                       <  1:454940/66‑1 (MQ=255)
gTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTa                                                       <  1:898595/66‑1 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:2147584/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:664206/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:2247491/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:226069/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:199965/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:2924986/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:1821220/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:382107/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:2789444/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:2636795/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:2586216/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:1501203/1‑53 (MQ=255)
       tAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTCAgt                                                             >  1:14595/1‑53 (MQ=255)
        aaGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCAttt                                             >  1:330064/1‑68 (MQ=255)
        aaGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTg                                            >  1:1526406/1‑69 (MQ=255)
        aaGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTg                                            >  1:545999/1‑69 (MQ=255)
        aaGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTg                                            >  1:621959/1‑69 (MQ=255)
                  gATAAAGAGTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:2597892/1‑69 (MQ=255)
                  gATAAAGAGTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:2601903/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:1046368/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:2933300/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:2944601/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:2027446/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:1901866/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:436451/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:1720620/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:525597/1‑69 (MQ=255)
                  gATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:806417/1‑69 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGc                                  >  1:2441978/1‑67 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                <  1:3165062/69‑1 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                <  1:3255908/69‑1 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                >  1:361808/1‑69 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                >  1:377829/1‑69 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                >  1:461247/1‑69 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                <  1:2732643/69‑1 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                >  1:2670922/1‑69 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                <  1:2612080/69‑1 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                <  1:2273227/69‑1 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                <  1:2226829/69‑1 (MQ=255)
                    tAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCg                                >  1:2076557/1‑69 (MQ=255)
                             aTCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATgcggc                        >  1:1760216/1‑68 (MQ=255)
                               cACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     <  1:1904732/69‑1 (MQ=255)
                               cACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     <  1:109368/69‑1 (MQ=255)
                               cACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     <  1:798335/69‑1 (MQ=255)
                               cACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     <  1:1771922/69‑1 (MQ=255)
                               cACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     <  1:2036891/69‑1 (MQ=255)
                                aCCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     >  1:3227628/1‑68 (MQ=255)
                                aCCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     >  1:299598/1‑68 (MQ=255)
                                aCCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     >  1:1190377/1‑68 (MQ=255)
                                aCCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     >  1:981841/1‑68 (MQ=255)
                                aCCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTa                     >  1:2699582/1‑68 (MQ=255)
                                aCCTGTAGGCCAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTaa                    <  1:333855/69‑1 (MQ=17)
                                aCCTGTAGGCCAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTaa                    <  1:562383/69‑1 (MQ=17)
                                aCCTGTAGGCCAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTaa                    <  1:1820579/69‑1 (MQ=17)
                                aCCTGTAGGCCAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTaa                    <  1:3330116/69‑1 (MQ=17)
                                aCCTGTAGGCCAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTaa                    <  1:3302397/69‑1 (MQ=17)
                                aCCTGTAGGCCAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTaa                    <  1:325906/69‑1 (MQ=17)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:1059327/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:854162/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:963078/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:3174937/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:74186/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:1400349/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:1404360/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:1815596/1‑48 (MQ=255)
                                  cTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCg                                       >  1:1923458/1‑48 (MQ=255)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCa                                          >  1:3272209/1‑37 (MQ=255)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTAAACGCCTTATc          >  1:1499699/1‑69 (MQ=14)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTAAACGCCTTATc          >  1:1764188/1‑69 (MQ=14)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTAAACGCCTTATc          >  1:1246244/1‑69 (MQ=14)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTAAACGCCTTATc          >  1:2994096/1‑69 (MQ=14)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTAAACGCCTTATc          >  1:2451674/1‑69 (MQ=14)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTAAACGCCTTATc          >  1:187617/1‑69 (MQ=14)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCACGGATTGCCGGATGCGGCGTAAACGCCTTATc          >  1:2530442/1‑69 (MQ=14)
                                          cAGATAAGACGCGTCAGTGTCTTATCTGGCATTTGCAGGATTGCCGGATGCGGCGTAAACGCCTTATcc         >  1:2449126/1‑69 (MQ=14)
                                                   cgcgTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc  <  1:3260300/68‑1 (MQ=255)
                                                        |                                                              
GTGACCTTAAGTCCGGCAGATAAAGACTAATCACCTGTAGGCCAGATAAGACGCGTTAGTGTCTTATCTGGCATTTGCACCGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAC  >  minE/1444350‑1444468

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: