Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,589,325 A→G 16.0% intergenic (‑1/+798) gadE ← / ← hdeD DNA‑binding transcriptional activator/acid‑resistance membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,589,3250AG16.0% 302.7 / 45.5 119intergenic (‑1/+798)gadE/hdeDDNA‑binding transcriptional activator/acid‑resistance membrane protein
Reads supporting (aligned to +/- strand):  ref base A (60/40);  new base G (16/3);  total (76/43)
Fisher's exact test for biased strand distribution p-value = 6.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 121 total aligned reads displayed.

TCGTGATTATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGAATCGT  >  minE/2589257‑2589395
                                                                    |                                                                      
tCGTGATTATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAAc                                                                      >  1:2241166/1‑71 (MQ=255)
tCGTGATTATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAAc                                                                      >  1:2598225/1‑71 (MQ=255)
 cGTGATTATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATGAc                                                                      <  1:1821156/70‑1 (MQ=255)
 cGTGATTATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACt                                                                     >  1:65058/1‑71 (MQ=255)
 cGTGATTATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACt                                                                     >  1:4002113/1‑71 (MQ=255)
      ttATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTc                                                                >  1:4249509/1‑71 (MQ=255)
      ttATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTc                                                                >  1:3742350/1‑71 (MQ=255)
        aTCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGc                                                                  >  1:1768673/1‑67 (MQ=255)
        aTCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCt                                                              >  1:1099925/1‑71 (MQ=255)
         tCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCt                                                              >  1:4519097/1‑70 (MQ=255)
             tCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCtt                                                             <  1:3315743/67‑1 (MQ=255)
             tCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGcc                                                         >  1:673607/1‑71 (MQ=255)
             tCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGcc                                                         >  1:1394879/1‑71 (MQ=255)
             tCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGcc                                                         >  1:3460820/1‑71 (MQ=255)
               aaCTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGt                                                       <  1:3951234/71‑1 (MQ=255)
               aaCTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGt                                                       <  1:3817865/71‑1 (MQ=255)
                aCTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGt                                                       <  1:841832/70‑1 (MQ=255)
                     cAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGt                                                 >  1:1191877/1‑71 (MQ=255)
                      aaaaGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATGACTTGCTCCTTAGCCGTTATCGtt                                                <  1:4264437/71‑1 (MQ=255)
                         aGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTc                                              >  1:476914/1‑70 (MQ=255)
                         aGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTc                                              >  1:2559597/1‑70 (MQ=255)
                         aGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTc                                              >  1:435992/1‑70 (MQ=255)
                           cccTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGaa                                           <  1:1070963/71‑1 (MQ=255)
                            ccTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAat                                          <  1:4629090/71‑1 (MQ=255)
                            ccTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAat                                          <  1:4071160/71‑1 (MQ=255)
                             cTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGaa                                           <  1:207882/69‑1 (MQ=255)
                               gTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGt                                                       <  1:2813413/55‑1 (MQ=255)
                               gTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGt                                                       <  1:2591457/55‑1 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAATATg                                       <  1:1721525/70‑1 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACt                                                                     >  1:505798/1‑40 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGtt                                                      <  1:3537784/55‑1 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGtt                                                      <  1:2151354/55‑1 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGaa                                           >  1:1832402/1‑66 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGaa                                           >  1:133893/1‑66 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAatat                                        <  1:193688/69‑1 (MQ=255)
                                tAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAatat                                        <  1:2396186/69‑1 (MQ=255)
                                 aaaagaaaagAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCt                                                              >  1:786134/1‑46 (MQ=255)
                                 aaaagaaaagAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCg                                                        >  1:3671317/1‑52 (MQ=255)
                                 aaaagaaaagAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAatat                                        <  1:3266966/68‑1 (MQ=255)
                                  aaagaaaagaATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCg                                                        >  1:569840/1‑51 (MQ=255)
                                  aaagaaaagaATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCg                                                        >  1:4641000/1‑51 (MQ=255)
                                  aaagaaaagaATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCa                                    >  1:3737949/1‑71 (MQ=255)
                                  aaagaaaagaATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCa                                    >  1:1140254/1‑71 (MQ=255)
                                  aaagaaaagaATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCa                                    >  1:3651576/1‑71 (MQ=255)
                                  aaagaaaagaATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCa                                    >  1:299929/1‑71 (MQ=255)
                                  aaagaaaagaATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCa                                    >  1:3743323/1‑71 (MQ=255)
                                   aagaaaagaaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCg                                                  <  1:3882483/56‑1 (MQ=255)
                                    agaaaagaaTCTTTCGTCATGAGAAAAATCATAACTTGCt                                                                 >  1:1719082/1‑40 (MQ=255)
                                    agaaaagaaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:1830804/1‑64 (MQ=255)
                                    agaaaagaaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:1413552/1‑64 (MQ=255)
                                    agaaaagaaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATc                                  <  1:1477385/71‑1 (MQ=255)
                                       aaagaaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCtt                               <  1:1511895/71‑1 (MQ=255)
                                        aagaaTCTTTCGTCATGAGAAAAATCATGACTTGCTCCTTAGGCGTTATCGttt                                               >  1:293767/1‑54 (MQ=255)
                                        aagaaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCttt                              >  1:665273/1‑71 (MQ=255)
                                        aagaaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCttt                              >  1:2174603/1‑71 (MQ=255)
                                           aaTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCg                                             <  1:4466108/53‑1 (MQ=255)
                                            aTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCat                          >  1:1690084/1‑71 (MQ=255)
                                            aTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCat                          <  1:2283152/71‑1 (MQ=255)
                                            aTCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCat                          >  1:3599007/1‑71 (MQ=255)
                                               tttCGTCATGAGAAAAATCATGACTTGCTCCTTAGCCGTTa                                                     >  1:4001482/1‑41 (MQ=255)
                                               tttCGTCATGAGAAAAATCATGACTTGCTCCTTAGCCGTTa                                                     >  1:882981/1‑41 (MQ=255)
                                               tttCGTCATGAGAAAAATCATGACTTGCTCCTTAGCCGTTa                                                     >  1:3268399/1‑41 (MQ=255)
                                               tttCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGaa                                           <  1:80760/51‑1 (MQ=255)
                                               tttCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCtttt                             <  1:2144644/65‑1 (MQ=255)
                                               tttCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCatta                        >  1:4269030/1‑70 (MQ=255)
                                               tttCGTAATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCatta                        >  1:4762768/1‑70 (MQ=255)
                                                ttCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCg                                                        <  1:3556116/37‑1 (MQ=255)
                                                 tCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATc                                  >  1:1212899/1‑58 (MQ=255)
                                                 tCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCatta                        <  1:48057/68‑1 (MQ=255)
                                                 tCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCatta                        <  1:3335921/68‑1 (MQ=255)
                                                 tCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCatta                        <  1:4309165/68‑1 (MQ=255)
                                                  cGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATc                                  >  1:2810527/1‑57 (MQ=255)
                                                      aTGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATc                >  1:1750738/1‑71 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:651938/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:691947/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:4700605/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:1999319/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:1332746/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:2032167/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:4340860/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:3321657/1‑44 (MQ=255)
                                                        gagaAAAATCATGACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:4139348/1‑44 (MQ=255)
                                                        gagaAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCattat                       >  1:1810308/1‑62 (MQ=255)
                                                          gaAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTg            >  1:3241235/1‑71 (MQ=255)
                                                             aaaTCATAACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:4625579/1‑39 (MQ=255)
                                                             aaaTCATAACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:2704462/1‑39 (MQ=255)
                                                             aaaTCATAACTTGCTCCTTAGCCGTTATCGTTTCGAata                                         >  1:139028/1‑39 (MQ=255)
                                                             aaaTCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTcc         >  1:1602844/1‑71 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAatat                                        <  1:339916/37‑1 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAatat                                        <  1:2024410/37‑1 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAatat                                        <  1:3020044/37‑1 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGaa      >  1:1964900/1‑71 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGaa      >  1:3129989/1‑71 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGaa      >  1:1719268/1‑71 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGaa      >  1:3094167/1‑71 (MQ=255)
                                                                tCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGaa      <  1:2092745/71‑1 (MQ=255)
                                                                  aTAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGt                                      >  1:3470094/1‑37 (MQ=255)
                                                                  aTAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGAATc    >  1:1517834/1‑71 (MQ=255)
                                                                    aaCTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGAATCgt  >  1:1885792/1‑71 (MQ=255)
                                                                    aaCTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGAATCgt  >  1:2233255/1‑71 (MQ=255)
                                                                    |                                                                      
TCGTGATTATCTTTCAACTGCCAAAAGCCCTGTAAAAGAAAAGAATCTTTCGTCATGAGAAAAATCATAACTTGCTCCTTAGCCGTTATCGTTTCGAATATGTCATCCTTTTCATTATTTACATCCTTGTCCGAATCGT  >  minE/2589257‑2589395

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: