Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,103,141 A→G 32.8% A204A (GCT→GCC agaR ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,103,1410AG32.8% 121.1 / 83.0 149A204A (GCT→GCCagaRDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (27/73);  new base G (20/29);  total (47/102)
Fisher's exact test for biased strand distribution p-value = 9.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 150 total aligned reads displayed.

TTGAACTTACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCA  >  minE/2103073‑2103208
                                                                    |                                                                   
ttGAACTTACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGa                                                                     >  1:1726472/1‑69 (MQ=255)
 tGAACTTACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAg                                                                    <  1:874337/69‑1 (MQ=255)
     cTTACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGcttc                                                                <  1:2939556/69‑1 (MQ=255)
     cTTACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGcttc                                                                <  1:1102047/69‑1 (MQ=255)
     cTTACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGcttc                                                                <  1:2321653/69‑1 (MQ=255)
      ttACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGcttc                                                                <  1:2471624/68‑1 (MQ=255)
      ttACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGcttc                                                                <  1:3054437/68‑1 (MQ=255)
                 ggTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    >  1:1728034/1‑69 (MQ=255)
                 ggTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    <  1:25679/69‑1 (MQ=255)
                 ggTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    >  1:2219260/1‑69 (MQ=255)
                 ggTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    <  1:338724/69‑1 (MQ=255)
                 ggTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    <  1:1898237/69‑1 (MQ=255)
                 ggTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    >  1:3408199/1‑69 (MQ=255)
                 ggTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    >  1:1010492/1‑69 (MQ=255)
                  gTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTAtg                                                    <  1:2250055/68‑1 (MQ=255)
                       cACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCt                                              <  1:2252355/69‑1 (MQ=255)
                       cACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCt                                              <  1:3200379/69‑1 (MQ=255)
                       cACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCt                                              <  1:1636383/69‑1 (MQ=255)
                       cACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCt                                              <  1:1824745/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:1265679/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:2522510/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:2931917/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:862318/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:2228283/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:2192808/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:2080674/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:626210/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:1522510/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:469188/69‑1 (MQ=255)
                         cGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGa                                            <  1:1357133/69‑1 (MQ=255)
                                             tCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAgatcgat                         >  1:1855277/1‑68 (MQ=255)
                                             tCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAgatcgat                         >  1:361340/1‑68 (MQ=255)
                                             tCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAgatcgat                         >  1:1600727/1‑68 (MQ=255)
                                               atccaTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGc                      <  1:727527/66‑1 (MQ=255)
                                               atccaTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGc                      <  1:1933694/66‑1 (MQ=255)
                                               atccaTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGc                      <  1:596629/66‑1 (MQ=255)
                                               atccaTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGc                      <  1:373298/66‑1 (MQ=255)
                                               atccaTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGc                      <  1:177892/66‑1 (MQ=255)
                                               atccaTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGc                      <  1:1747303/66‑1 (MQ=255)
                                                cacaTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGatcgatcg                       <  1:1572798/67‑1 (MQ=255)
                                                  caTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGAcgc                                         >  1:25481/1‑47 (MQ=255)
                                                  caTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGAcgc                                         >  1:1384921/1‑47 (MQ=255)
                                                  caTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGAcgc                                         >  1:1026456/1‑47 (MQ=255)
                                                  caTCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGAcgc                                         >  1:1246093/1‑47 (MQ=255)
                                                  caTCCGAGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGAcgc                                         >  1:2406546/1‑46 (MQ=39)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCa                                >  1:1072151/1‑52 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:3197591/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:496859/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:3089891/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:3259067/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:3377850/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:891736/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:695217/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:720468/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:2460784/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:2572461/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:2078523/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:1375257/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:1372482/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAg                               >  1:2346141/1‑53 (MQ=37)
                                                      cGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTac                >  1:3334103/1‑68 (MQ=255)
                                                      cGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTac                >  1:641552/1‑68 (MQ=255)
                                                      cGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTac                >  1:681932/1‑68 (MQ=255)
                                                      cGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTac                >  1:23350/1‑68 (MQ=255)
                                                          ggTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGc                      <  1:1459222/58‑1 (MQ=255)
                                                            tttAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGa         <  1:1445303/69‑1 (MQ=255)
                                                            tttAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGa         <  1:2682094/69‑1 (MQ=255)
                                                            tttAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGa         <  1:3311860/69‑1 (MQ=255)
                                                            tttAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGa         <  1:1052667/69‑1 (MQ=255)
                                                            tttAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGa         <  1:3341766/69‑1 (MQ=255)
                                                            tttAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGa         <  1:3356325/69‑1 (MQ=255)
                                                             ttAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGa         <  1:2661737/68‑1 (MQ=255)
                                                              tAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGaaa       <  1:2310351/69‑1 (MQ=255)
                                                              tAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGaaa       <  1:1574048/69‑1 (MQ=255)
                                                              tAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGaaa       <  1:554498/69‑1 (MQ=255)
                                                              tAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGaaa       <  1:1007121/69‑1 (MQ=255)
                                                              tAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGaaa       <  1:842724/69‑1 (MQ=255)
                                                              tAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGaaa       <  1:324672/69‑1 (MQ=255)
                                                              tAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGaaa       <  1:135974/69‑1 (MQ=255)
                                                                 aCGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAGATCGATCGCATCTACACCAAGAAACa     <  1:1142688/68‑1 (MQ=255)
                                                                 aCGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAGATCGATCGCATCTACACCAAGAAACa     <  1:595449/68‑1 (MQ=255)
                                                                 aCGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAGATCGATCGCATCTACACCAAGAAACa     <  1:1333328/68‑1 (MQ=255)
                                                                 aCGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAGATCGATCGCATCTACACCAAGAAACa     <  1:1665351/68‑1 (MQ=255)
                                                                 aCGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAGATCGATCGCATCTACACCAAGAAACa     <  1:1662008/68‑1 (MQ=255)
                                                                 aCGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAGATCGATCGCATCTACACCAAGAAACa     <  1:3291636/68‑1 (MQ=255)
                                                                 aCGGGCTTCATCTTCATTATGCGTGCTGATGCCGCACTCCAGATCGATCGCATCTACACCAAGAAACa     <  1:3099546/68‑1 (MQ=255)
                                                                  cGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAg    <  1:2807827/68‑1 (MQ=255)
                                                                  cGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGc   <  1:456448/69‑1 (MQ=255)
                                                                  cGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGc   <  1:2157430/69‑1 (MQ=255)
                                                                  cGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGc   <  1:2205522/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:975683/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:1024027/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:2817396/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:2449542/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:773543/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:3049416/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:2183564/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:672676/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:659086/69‑1 (MQ=255)
                                                                   gAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCa  <  1:1634288/69‑1 (MQ=255)
                                                                    |                                                                   
TTGAACTTACTGGAATCGGTGACCACGATGATCCGTTCCGCAACTTCGCACATCCGACGGTTTAAACGAGCTTCATCTTCATTATGTGTGCTGACGCCGCGCTCCAGATCGATCGCATCTACACCAAGAAACAGCA  >  minE/2103073‑2103208

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: