Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,608,096 T→C 35.6% intergenic (‑206/+104) 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,0960TC35.6% 35.1 / 38.8 45intergenic (‑206/+104)pitA/nikRphosphate transporter, low‑affinity/DNA‑binding transcriptional regulator, Ni‑binding
Reads supporting (aligned to +/- strand):  ref base T (15/14);  new base C (11/5);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

TCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCTT  >  minE/2608028‑2608165
                                                                    |                                                                      
tCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGccc                                                                      >  1:3930310/1‑71 (MQ=255)
 cAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCg                                                                     >  1:1374382/1‑71 (MQ=255)
 cAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCg                                                                     >  1:857431/1‑71 (MQ=255)
 cAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCg                                                                     >  1:1915611/1‑71 (MQ=255)
   ggCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGtt                                                                   >  1:3893897/1‑71 (MQ=255)
      aaaTACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGttttt                                                                >  1:644422/1‑71 (MQ=255)
        aTACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt                                                               >  1:4923204/1‑70 (MQ=255)
        aTACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt                                                               >  1:2118391/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:4228536/1‑70 (MQ=255)
                        atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAAcg                                              >  1:631526/1‑70 (MQ=255)
                         tatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCg                                            <  1:2676509/71‑1 (MQ=255)
                           ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                          >  1:5236651/1‑71 (MQ=255)
                           ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                          >  1:1057640/1‑71 (MQ=255)
                           ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAAcgtg                                            <  1:2997344/70‑4 (MQ=255)
                             gATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                          >  1:2577980/1‑69 (MQ=255)
                                    gTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAaga                                                         <  1:4810235/47‑1 (MQ=255)
                                    gTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAaga                                                         <  1:874692/47‑1 (MQ=255)
                                    gTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAaga                                                         <  1:3617268/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:5499600/54‑1 (MQ=255)
                                              tGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATa                                        <  1:64158/54‑1 (MQ=255)
                                              tGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATa                                        <  1:1322613/54‑1 (MQ=255)
                                               gAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAACGTGGAAAggtt                                    >  1:1562825/1‑56 (MQ=255)
                                                    aTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                          >  1:3811757/1‑46 (MQ=255)
                                                    aTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                          >  1:5631197/1‑46 (MQ=255)
                                                    aTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                          >  1:827686/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)
                                                                    |                                                                      
TCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCTT  >  minE/2608028‑2608165

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: