Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,608,096 T→C 7.2% 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,0960TC 7.2% 197.9 / 6.3 69intergenic (‑206/+104)pitA/nikRphosphate transporter, low‑affinity/DNA‑binding transcriptional regulator, Ni‑binding
Reads supporting (aligned to +/- strand):  ref base T (41/23);  new base C (3/2);  total (44/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

TTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCT  >  minE/2608027‑2608164
                                                                     |                                                                     
ttCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTc                                                                      >  1:2190464/1‑71 (MQ=255)
   aGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGt                                                                   <  1:1280541/71‑1 (MQ=255)
   aGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGt                                                                   <  1:200145/71‑1 (MQ=255)
         aTACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt‑t                                                            <  1:720811/71‑1 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt‑t                                                            >  1:3564136/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt‑t                                                            >  1:4035757/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt‑t                                                            >  1:3896372/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtgtttt‑t                                                            >  1:1852855/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGttttttct                                                            <  1:1104533/71‑3 (MQ=255)
           aCTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGttttttct                                                            >  1:2290163/1‑68 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑Taa                                                          <  1:2013794/70‑1 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAag                                                         <  1:927690/71‑1 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAag                                                         <  1:3887649/71‑1 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAag                                                         <  1:2343505/71‑1 (MQ=255)
              tCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt                                                              >  1:3904238/1‑65 (MQ=255)
                cTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAagag                                                       >  1:1310187/1‑69 (MQ=255)
                cTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAagag                                                       >  1:792563/1‑69 (MQ=255)
                   ttAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAgcg                                                  >  1:174008/1‑71 (MQ=255)
                       tAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑Taa                                                          <  1:4091003/59‑1 (MQ=255)
                         atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt‑t                                                            <  1:3884657/55‑1 (MQ=255)
                         atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttt                                                               >  1:2265395/1‑53 (MQ=255)
                         atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTc                                            >  1:204053/1‑71 (MQ=255)
                          tatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAgcgc                                                 <  1:152597/65‑1 (MQ=255)
                           atTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCg                                           >  1:1975708/1‑70 (MQ=255)
                            ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGa                                         >  1:1209400/1‑71 (MQ=255)
                            ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAGCGCAAcgtg                                           <  1:2123149/70‑4 (MQ=255)
                              gATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAgcgc                                                 >  1:1202427/1‑62 (MQ=255)
                              gATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGTTTTTTCTAAGAGGAgcgc                                                 >  1:253055/1‑62 (MQ=255)
                               aTGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAg                                      >  1:365979/1‑71 (MQ=255)
                                tGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAgg                                     >  1:2293183/1‑71 (MQ=255)
                                tGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAgg                                     <  1:1394670/71‑1 (MQ=255)
                                  cTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGt                                   <  1:3978835/71‑1 (MQ=255)
                                  cTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGt                                   <  1:508340/71‑1 (MQ=255)
                                    ggTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGa                                         >  1:1659849/1‑63 (MQ=255)
                                    ggTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTAt                                 >  1:3453483/1‑71 (MQ=255)
                                       gCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAgagga                                                     >  1:2896288/1‑48 (MQ=255)
                                        cGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTa                              <  1:848125/70‑1 (MQ=255)
                                          aCCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑Taa                                                          >  1:3623905/1‑40 (MQ=255)
                                           ccACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTa                                  >  1:3061273/1‑63 (MQ=255)
                                            cACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTagagag                         >  1:1620635/1‑71 (MQ=255)
                                            cACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTagagag                         <  1:3776699/71‑1 (MQ=255)
                                             aCTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAgcgc                                                 <  1:1602685/46‑1 (MQ=255)
                                               tGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAag                                                         <  1:4105910/36‑1 (MQ=255)
                                               tGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGaa                      >  1:2911873/1‑71 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAgaggag                                                    >  1:780990/1‑39 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAgagg                                                      >  1:3320090/1‑37 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAggg                                    >  1:3285699/1‑55 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAAt                     >  1:3516879/1‑70 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATc                    <  1:3351166/71‑1 (MQ=255)
                                                  ggAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATc                    <  1:3968652/70‑1 (MQ=255)
                                                   gAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAg                  >  1:862746/1‑71 (MQ=255)
                                                    aaTCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAgg                 >  1:2564826/1‑71 (MQ=255)
                                                     aTCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgt                >  1:3333432/1‑71 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAg                  >  1:647437/1‑68 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgtg               >  1:2884662/1‑71 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgtg               >  1:546391/1‑71 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgt                >  1:2083605/1‑70 (MQ=255)
                                                        tttACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCg                                           <  1:88098/41‑1 (MQ=255)
                                                        tttACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGAt                                        >  1:1375928/1‑44 (MQ=255)
                                                        tttACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGtgtg             >  1:3676509/1‑71 (MQ=255)
                                                          tACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGtgtg             <  1:1857284/69‑1 (MQ=255)
                                                          tACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGc            >  1:3250003/1‑70 (MQ=255)
                                                            cAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCt           <  1:3235412/69‑1 (MQ=255)
                                                               ttCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTa                              <  1:96967/47‑1 (MQ=255)
                                                               ttCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTgg      >  1:373527/1‑71 (MQ=255)
                                                               ttCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTgg      >  1:2647913/1‑71 (MQ=255)
                                                                  aCGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGccc   >  1:3366330/1‑71 (MQ=255)
                                                                   cGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTaga                            >  1:2718107/1‑45 (MQ=255)
                                                                   cGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCt  >  1:1765890/1‑71 (MQ=255)
                                                                     tCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCt           <  1:2468835/60‑1 (MQ=255)
                                                                     |                                                                     
TTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTT‑TAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCT  >  minE/2608027‑2608164

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: ‑