Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,608,105:1 +C 29.7% intergenic (‑215/+95) pitA ← / ← nikR phosphate transporter, low‑affinity/DNA‑binding transcriptional regulator, Ni‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,608,1051.C29.7% 39.1 / 54.0 37intergenic (‑215/+95)pitA/nikRphosphate transporter, low‑affinity/DNA‑binding transcriptional regulator, Ni‑binding
Reads supporting (aligned to +/- strand):  ref base . (11/15);  new base C (6/5);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCTTTGTTGATCG  >  minE/2608036‑2608174
                                                                      |                                                                     
aTACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt                                                                        >  1:2118391/1‑70 (MQ=255)
aTACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt                                                                        >  1:4923204/1‑70 (MQ=255)
 tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttttt                                                                      >  1:4453845/1‑70 (MQ=255)
 tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt                                                                        >  1:403954/1‑69 (MQ=255)
 tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTa                                                                     <  1:4770038/71‑1 (MQ=255)
      ccTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAaga                                                                  <  1:3109219/69‑1 (MQ=255)
      ccTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAagag                                                                 >  1:305275/1‑71 (MQ=255)
      ccTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAaga                                                                  <  1:673437/70‑1 (MQ=255)
                atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTc                                                      >  1:3669960/1‑71 (MQ=255)
                atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAaga                                                                  >  1:4937412/1‑60 (MQ=255)
                atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAAcg                                                       >  1:631526/1‑70 (MQ=255)
                atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAAcg                                                       >  1:4228536/1‑70 (MQ=255)
                 tatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCg                                                     <  1:2676509/71‑1 (MQ=255)
                   ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                                   >  1:1057640/1‑71 (MQ=255)
                   ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                                   >  1:5236651/1‑71 (MQ=255)
                   ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAAcgtg                                                     <  1:2997344/70‑4 (MQ=255)
                     gATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                                   >  1:2577980/1‑69 (MQ=255)
                            gTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAaga                                                                  <  1:3617268/47‑1 (MQ=255)
                            gTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAaga                                                                  <  1:4810235/47‑1 (MQ=255)
                            gTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAaga                                                                  <  1:874692/47‑1 (MQ=255)
                            gTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTAtt                                          <  1:2300249/71‑1 (MQ=255)
                            gTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAggttt                                            >  1:2578678/1‑67 (MQ=255)
                            gTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAggttt                                            <  1:4867585/70‑4 (MQ=255)
                                gACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTaga                                      <  1:5421076/71‑1 (MQ=255)
                                gACCACTGAGGAATCTTTACAATTCAAGCCCGTTTTTTCTAAgagga                                                               <  1:4595341/47‑1 (MQ=255)
                                  ccACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTagaga                                    <  1:485510/71‑1 (MQ=255)
                                      tGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATa                                                 <  1:1322613/54‑1 (MQ=255)
                                      tGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATa                                                 <  1:5499600/54‑1 (MQ=255)
                                      tGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATa                                                 <  1:64158/54‑1 (MQ=255)
                                       gAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAggtt                                             >  1:1562825/1‑56 (MQ=255)
                                            aTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                                   >  1:3811757/1‑46 (MQ=255)
                                            aTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                                   >  1:827686/1‑46 (MQ=255)
                                            aTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                                   >  1:5631197/1‑46 (MQ=255)
                                              cTTTACAATTCACGCCCGTTTTTTCTAAGAGGAgcgc                                                           <  1:3927858/37‑1 (MQ=255)
                                               tttACAATTCACGTCCGTTTTTTTAAGAGGAGCGCaa                                                         <  1:3927597/37‑1 (MQ=255)
                                               tttACAATTCACGTCCGTTTTTTTAAGAGGAGCGCaa                                                         <  1:1736691/37‑1 (MQ=255)
                                                ttACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGc                      >  1:3951402/1‑71 (MQ=255)
                                                        cACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGcc              >  1:4621154/1‑71 (MQ=255)
                                                        cACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGcc              >  1:815398/1‑71 (MQ=255)
                                                           gTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCtt           >  1:3232550/1‑71 (MQ=255)
                                                               gTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCTttgtt       <  1:4538357/71‑1 (MQ=255)
                                                                tttttttAAGAGGAGCTCAACTCGGATAGGGTATTTAGAGAgg                                  <  1:2453423/43‑1 (MQ=255)
                                                                tttttttAAGAGGAGCGCAACTCGGATAGGGTATTTa                                        <  1:1442689/37‑1 (MQ=255)
                                                                tttttttAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGTAATCAGGTGTGCTAGTGGCCCttt          <  1:2812645/67‑1 (MQ=255)
                                                                  tttttAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCTTTGTTGAt    >  1:1875525/1‑71 (MQ=255)
                                                                    tttAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCTTTGTTGATCg  >  1:2484219/1‑71 (MQ=255)
                                                                      |                                                                     
ATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCTTTGTTGATCG  >  minE/2608036‑2608174

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: