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

TTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCT  >  minE/2608027‑2608164
                                                                     |                                                                     
ttCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTc                                                                      >  1:2190464/1‑71 (MQ=255)
   aGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGt                                                                   <  1:1280541/71‑1 (MQ=255)
   aGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGt                                                                   <  1:200145/71‑1 (MQ=255)
         aTACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttttt                                                            <  1:720811/71‑1 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttttt                                                            >  1:3564136/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttttt                                                            >  1:4035757/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttttt                                                            >  1:3896372/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtgttttt                                                            >  1:1852855/1‑70 (MQ=255)
          tACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGttttttct                                                            <  1:1104533/71‑3 (MQ=255)
           aCTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGCCCGttttttct                                                            >  1:2290163/1‑68 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTaa                                                          <  1:2013794/70‑1 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAag                                                         <  1:927690/71‑1 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAag                                                         <  1:3887649/71‑1 (MQ=255)
            cTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAag                                                         <  1:2343505/71‑1 (MQ=255)
              tCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGtttttt                                                              >  1:3904238/1‑65 (MQ=255)
                cTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAagag                                                       >  1:1310187/1‑69 (MQ=255)
                cTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAagag                                                       >  1:792563/1‑69 (MQ=255)
                   ttAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAgcg                                                  >  1:174008/1‑71 (MQ=255)
                       tAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTaa                                                          <  1:4091003/59‑1 (MQ=255)
                         atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttttt                                                            <  1:3884657/55‑1 (MQ=255)
                         atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGttttt                                                               >  1:2265395/1‑53 (MQ=255)
                         atatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTc                                            >  1:204053/1‑71 (MQ=255)
                          tatTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAgcgc                                                 <  1:152597/65‑1 (MQ=255)
                           atTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCg                                           >  1:1975708/1‑70 (MQ=255)
                            ttGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                         >  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)
                               aTGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAg                                      >  1:365979/1‑71 (MQ=255)
                                tGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAgg                                     >  1:2293183/1‑71 (MQ=255)
                                tGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAgg                                     <  1:1394670/71‑1 (MQ=255)
                                  cTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGt                                   <  1:3978835/71‑1 (MQ=255)
                                  cTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGt                                   <  1:508340/71‑1 (MQ=255)
                                    ggTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGa                                         >  1:1659849/1‑63 (MQ=255)
                                    ggTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTAt                                 >  1:3453483/1‑71 (MQ=255)
                                       gCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAgagga                                                     >  1:2896288/1‑48 (MQ=255)
                                        cGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTa                              <  1:848125/70‑1 (MQ=255)
                                          aCCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTaa                                                          >  1:3623905/1‑40 (MQ=255)
                                           ccACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTa                                  >  1:3061273/1‑63 (MQ=255)
                                            cACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTagagag                         >  1:1620635/1‑71 (MQ=255)
                                            cACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTagagag                         <  1:3776699/71‑1 (MQ=255)
                                             aCTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAgcgc                                                 <  1:1602685/46‑1 (MQ=255)
                                               tGAGGAATCTTTACAATTCACGTCCGTTTTTTTAag                                                         <  1:4105910/36‑1 (MQ=255)
                                               tGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGaa                      >  1:2911873/1‑71 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTTTAAgaggag                                                    >  1:780990/1‑39 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTTTAAgagg                                                      >  1:3320090/1‑37 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAggg                                    >  1:3285699/1‑55 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAAt                     >  1:3516879/1‑70 (MQ=255)
                                                 aGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATc                    <  1:3351166/71‑1 (MQ=255)
                                                  ggAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATc                    <  1:3968652/70‑1 (MQ=255)
                                                   gAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAg                  >  1:862746/1‑71 (MQ=255)
                                                    aaTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAgg                 >  1:2564826/1‑71 (MQ=255)
                                                     aTCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgt                >  1:3333432/1‑71 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAg                  >  1:647437/1‑68 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgtg               >  1:2884662/1‑71 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgtg               >  1:546391/1‑71 (MQ=255)
                                                      tCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGgt                >  1:2083605/1‑70 (MQ=255)
                                                        tttACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCg                                           <  1:88098/41‑1 (MQ=255)
                                                        tttACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGAt                                        >  1:1375928/1‑44 (MQ=255)
                                                        tttACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGtgtg             >  1:3676509/1‑71 (MQ=255)
                                                          tACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGtgtg             <  1:1857284/69‑1 (MQ=255)
                                                          tACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGc            >  1:3250003/1‑70 (MQ=255)
                                                            cAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCt           <  1:3235412/69‑1 (MQ=255)
                                                               ttCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTa                              <  1:96967/47‑1 (MQ=255)
                                                               ttCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTgg      >  1:373527/1‑71 (MQ=255)
                                                               ttCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTgg      >  1:2647913/1‑71 (MQ=255)
                                                                  aCGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGccc   >  1:3366330/1‑71 (MQ=255)
                                                                   cGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTaga                            >  1:2718107/1‑45 (MQ=255)
                                                                   cGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCt  >  1:1765890/1‑71 (MQ=255)
                                                                     tCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCt           <  1:2468835/60‑1 (MQ=255)
                                                                     |                                                                     
TTCAGGCAAATACTTCCTTTTAGTAATATTGATGCTGGTGCGACCACTGAGGAATCTTTACAATTCACGTCCGTTTTTTTAAGAGGAGCGCAACTCGGATAGGGTATTTAGAGAGGAATCAGGTGTGCTAGTGGCCCT  >  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: