Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 769,085 G→C 38.0% P206P (CCC→CCG phoP ← DNA‑binding response regulator in two‑component regulatory system with PhoQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE769,0850GC38.0% 84.5 / 111.2 121P206P (CCC→CCGphoPDNA‑binding response regulator in two‑component regulatory system with PhoQ
Reads supporting (aligned to +/- strand):  ref base G (23/52);  new base C (13/33);  total (36/85)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

Only 100 of 128 total aligned reads displayed.

TAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTT  >  minE/769020‑769145
                                                                 |                                                            
tAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGa                                                           <  1:92014/69‑1 (MQ=255)
tAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGa                                                           <  1:2180462/69‑1 (MQ=255)
             aGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGaa                                              >  1:3216073/1‑69 (MQ=255)
             aGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGaa                                              >  1:841645/1‑69 (MQ=255)
                 cAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAAtttt                                          >  1:450414/1‑69 (MQ=255)
                 cAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAAtttt                                          >  1:2899487/1‑69 (MQ=255)
                 cAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAAtttt                                          >  1:66060/1‑69 (MQ=255)
                 cAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAAtttt                                          >  1:1334838/1‑69 (MQ=255)
                         aCAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCa                                  <  1:3292267/69‑1 (MQ=255)
                         aCAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCa                                  <  1:2550300/69‑1 (MQ=255)
                         aCAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCa                                  <  1:1914890/69‑1 (MQ=255)
                              tAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:2373325/62‑1 (MQ=255)
                              tAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:2245461/62‑1 (MQ=255)
                              tAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:52336/62‑1 (MQ=255)
                              tAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:2561219/62‑1 (MQ=255)
                              tAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:3232205/62‑1 (MQ=255)
                              tAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:1594639/62‑1 (MQ=255)
                               aGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTc                            <  1:733260/69‑1 (MQ=255)
                               aGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTc                            <  1:2129151/69‑1 (MQ=255)
                                gCCCTGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTc                            <  1:81381/68‑1 (MQ=255)
                                gCCCTGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTc                            <  1:1983105/68‑1 (MQ=255)
                                  ccTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTccc                          >  1:1805642/1‑68 (MQ=255)
                                  ccTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTccc                          >  1:71608/1‑68 (MQ=255)
                                  ccTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTccc                          >  1:3059380/1‑68 (MQ=255)
                                  ccTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTccc                          >  1:2838306/1‑68 (MQ=255)
                                  ccTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTccc                          >  1:1385228/1‑68 (MQ=255)
                                  ccTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTccc                          >  1:617186/1‑68 (MQ=255)
                                    tGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATc                       >  1:1284288/1‑69 (MQ=255)
                                    tGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATc                       >  1:1144749/1‑69 (MQ=255)
                                    tGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATc                       >  1:403691/1‑69 (MQ=255)
                                    tGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATc                       >  1:2204395/1‑69 (MQ=255)
                                    tGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATc                       >  1:1255141/1‑69 (MQ=255)
                                    tGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATc                       >  1:1296311/1‑69 (MQ=255)
                                    tGGCCGCGAACGGTGGTAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATc                       >  1:383309/1‑69 (MQ=255)
                                        cgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:2795117/67‑1 (MQ=255)
                                        cgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:3131844/67‑1 (MQ=255)
                                        cgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:963873/67‑1 (MQ=255)
                                        cgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:841226/67‑1 (MQ=255)
                                        cgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:760302/67‑1 (MQ=255)
                                        cgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:738869/67‑1 (MQ=255)
                                        cgcgAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:615136/67‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGGTATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:582280/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:2080718/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:1087502/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:3025915/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:1472607/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:3336168/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:2632800/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:1003385/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:1888430/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:398558/69‑1 (MQ=255)
                                                 tggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:220950/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:2569099/68‑1 (MQ=255)
                                                  ggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:307436/68‑1 (MQ=255)
                                                  ggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGt         <  1:1852661/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGt         <  1:3302993/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGt         <  1:1856299/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGt         <  1:2687543/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGt         <  1:2358283/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGt         <  1:1585121/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGt         <  1:1938018/69‑1 (MQ=255)
                                                  ggtAATCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGt         <  1:2954127/69‑1 (MQ=255)
                                                   gtAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:2288295/56‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:564564/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:2610251/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:879067/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:3340984/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:602270/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:3030217/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:1828664/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:1761211/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:1581640/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:1973277/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:2090318/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:2094254/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:115180/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:2611443/68‑1 (MQ=255)
                                                     aaTCACTTCTTGCGGATATTGTGCCTGTATTTTTTTGCGCAGACGTCCCATCAGTACATCTATGGTAt       <  1:2750028/68‑1 (MQ=255)
                                                      aTCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:1317235/64‑1 (MQ=255)
                                                      aTCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:958041/64‑1 (MQ=255)
                                                      aTCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:2669333/64‑1 (MQ=255)
                                                      aTCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:2753986/64‑1 (MQ=255)
                                                      aTCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATgg          <  1:623681/64‑1 (MQ=255)
                                                        cACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:1454273/36‑1 (MQ=255)
                                                        cACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:2007051/36‑1 (MQ=255)
                                                        cACTTCTTGGGGATATTGTGCCTGAATTTTTTTgcg                                    <  1:1768823/36‑1 (MQ=255)
                                                         aCTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGc    >  1:57169/1‑67 (MQ=255)
                                                         aCTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCtt  >  1:3117542/1‑69 (MQ=255)
                                                         aCTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCtt  >  1:1030563/1‑69 (MQ=255)
                                                         aCTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCtt  >  1:2192798/1‑69 (MQ=255)
                                                           ttcttGGGGATATTGTGCCTGAATTTTTTTGCGCAGCCGt                             >  1:2615277/1‑40 (MQ=255)
                                                           ttcttGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGt                             >  1:1912938/1‑40 (MQ=255)
                                                           ttcttGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGt                             >  1:240928/1‑40 (MQ=255)
                                                           ttcttGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGt                             >  1:1107893/1‑40 (MQ=255)
                                                                ggggATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:1912549/43‑1 (MQ=255)
                                                                ggggATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:622286/43‑1 (MQ=255)
                                                                ggggATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:1931821/43‑1 (MQ=255)
                                                                ggggATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:2129696/43‑1 (MQ=255)
                                                                ggggATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:1429558/43‑1 (MQ=255)
                                                                ggggATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAg                     <  1:2698657/43‑1 (MQ=255)
                                                                 |                                                            
TAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTT  >  minE/769020‑769145

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: