Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,070,074 G→A 51.6% intergenic (‑285/‑106) aer ← / → ygjG fused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component/putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,070,0740GA51.6% ‑2.3 / 87.4 95intergenic (‑285/‑106)aer/ygjGfused signal transducer for aerotaxis sensory component and methyl accepting chemotaxis component/putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (22/27);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

TGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGG  >  minE/2070011‑2070142
                                                               |                                                                    
tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                                                 >  1:2112052/1‑69 (MQ=255)
tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                                                 >  1:1803614/1‑69 (MQ=255)
tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                                                 >  1:3414828/1‑69 (MQ=255)
    cagcaATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTa                                                             <  1:1528062/69‑1 (MQ=255)
      gcaATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTaaa                                                           <  1:3131411/69‑1 (MQ=255)
      gcaATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTaaa                                                           <  1:2998807/69‑1 (MQ=255)
                  tCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAAt                                               <  1:806804/69‑1 (MQ=255)
                  tCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAAt                                               <  1:1392963/69‑1 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:1018689/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:3209121/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:3197843/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:1951680/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:1769487/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:1600007/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:1118800/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:2588151/1‑68 (MQ=255)
                    cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                              >  1:2741026/1‑68 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:59250/68‑1 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:3055744/68‑1 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:393717/68‑1 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:3115004/68‑1 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:1132026/68‑1 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                           <  1:54872/69‑1 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                           <  1:836775/69‑1 (MQ=255)
                      tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                           <  1:3267817/69‑1 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:2551048/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:2944721/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:2948642/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:2967628/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:2221089/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:2224997/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:3211469/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:416602/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:814663/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:1173887/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                          >  1:819269/1‑68 (MQ=255)
                        cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                           >  1:737031/1‑67 (MQ=255)
                              aaaaaCCGTGCAAAATAACAACAAATGTTAACTGATAGCATTaaa                                                           <  1:2991832/45‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:3265028/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:3298674/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:3388965/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:3402816/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:418438/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:47514/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:781319/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:826156/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:987615/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1888012/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1086833/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1128694/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1419299/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1481463/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1519008/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1562140/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1628737/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1632603/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1682191/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1693450/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:2141354/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:2316557/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:245544/66‑1 (MQ=255)
                                        cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:2356503/66‑1 (MQ=255)
                                         aaaaTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:3318837/65‑1 (MQ=255)
                                            aaaacagcAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2623487/3‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATccc                       >  1:1389271/1‑67 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCg                      >  1:2483077/1‑68 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:698645/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2687268/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:1385413/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2867573/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2429916/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:1338413/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2977403/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:1215365/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:3329117/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2554452/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:328332/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:195919/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2002694/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:835165/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:11175/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:3152735/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:253330/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2113267/1‑69 (MQ=255)
                                            aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                     >  1:2308871/1‑69 (MQ=255)
                                            aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:2188078/46‑1 (MQ=255)
                                            aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:2861802/46‑1 (MQ=255)
                                            aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                            <  1:2081250/46‑1 (MQ=255)
                                                    aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg             <  1:137587/69‑1 (MQ=255)
                                                    aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg             <  1:340478/69‑1 (MQ=255)
                                                    aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg             <  1:3354655/69‑1 (MQ=255)
                                                    aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg             <  1:1157035/69‑1 (MQ=255)
                                                    aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg             <  1:2047333/69‑1 (MQ=255)
                                                      aTGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:1959240/52‑1 (MQ=255)
                                                        gTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                            <  1:819363/50‑1 (MQ=255)
                                                               gATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCgg  <  1:314541/69‑1 (MQ=255)
                                                               |                                                                    
TGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGG  >  minE/2070011‑2070142

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: