Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,070,071 A→G 43.7% intergenic (‑282/‑109) 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,0710AG43.7% 25.9 / 89.7 112intergenic (‑282/‑109)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 A (41/22);  new base G (22/27);  total (63/49)
Fisher's exact test for biased strand distribution p-value = 3.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

CATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCG  >  minE/2070007‑2070131
                                                                |                                                            
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTaa                                                              >  1:2738272/1‑65 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:1046651/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:933862/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:3379202/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:334951/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:328131/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:3273477/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:2626333/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:2459646/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:1853038/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:1727033/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:1634550/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:157563/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:1261015/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:1304507/1‑67 (MQ=255)
cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                            >  1:1323085/1‑67 (MQ=255)
    tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                                      >  1:2112052/1‑69 (MQ=255)
    tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                                      >  1:3414828/1‑69 (MQ=255)
    tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                                      >  1:1803614/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)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:1018689/1‑68 (MQ=255)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:1951680/1‑68 (MQ=255)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:2588151/1‑68 (MQ=255)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:1769487/1‑68 (MQ=255)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:2741026/1‑68 (MQ=255)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:1600007/1‑68 (MQ=255)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:3197843/1‑68 (MQ=255)
                        cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                   >  1:1118800/1‑68 (MQ=255)
                          tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                 <  1:3055744/68‑1 (MQ=255)
                          tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                 <  1:3115004/68‑1 (MQ=255)
                          tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                 <  1:59250/68‑1 (MQ=255)
                          tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                 <  1:1132026/68‑1 (MQ=255)
                          tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                <  1:836775/69‑1 (MQ=255)
                          tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                <  1:3267817/69‑1 (MQ=255)
                          tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                <  1:54872/69‑1 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:2944721/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:416602/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:819269/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:814663/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:2221089/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:2224997/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:2948642/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:3211469/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                               >  1:2551048/1‑68 (MQ=255)
                            cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                >  1:737031/1‑67 (MQ=255)
                                  aaaaaCCGTGCAAAATAACAACAAATGTTAACTGATAGCATTaaa                                                <  1:2991832/45‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1562140/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1086833/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1628737/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1682191/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:3298674/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1519008/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:3388965/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:3265028/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1128694/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1481463/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:3402816/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1693450/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:418438/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:47514/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1888012/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:1419299/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:245544/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:2356503/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:2316557/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:781319/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:2141354/66‑1 (MQ=255)
                                            cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                 <  1:826156/66‑1 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATccc            >  1:1389271/1‑67 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCg           >  1:2483077/1‑68 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:1385413/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:3329117/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:835165/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:11175/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:698645/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:3152735/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:2002694/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:195919/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:2308871/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:2429916/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:1215365/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:253330/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:2554452/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:2687268/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:2867573/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:2977403/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:1338413/1‑69 (MQ=255)
                                                aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa          >  1:328332/1‑69 (MQ=255)
                                                aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                 <  1:2861802/46‑1 (MQ=255)
                                                aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                 <  1:2188078/46‑1 (MQ=255)
                                                aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                 <  1:2081250/46‑1 (MQ=255)
                                                        aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg  <  1:137587/69‑1 (MQ=255)
                                                        aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg  <  1:1157035/69‑1 (MQ=255)
                                                        aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg  <  1:340478/69‑1 (MQ=255)
                                                        aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg  <  1:3354655/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)
                                                                |                                                            
CATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCG  >  minE/2070007‑2070131

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: