Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,070,105 G→A 72.9% intergenic (‑316/‑75) 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,1050GA72.9% 98.4 / 37.4 92intergenic (‑316/‑75)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 (8/17);  new base A (40/27);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 2.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATACTTAAATCGG  >  minE/2070051‑2070173
                                                      |                                                                    
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:3265028/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:1888012/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:2141354/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:2316557/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:2356503/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:245544/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:1086833/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:1519008/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:1481463/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:1562140/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:1419299/66‑1 (MQ=255)
cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                                           <  1:1128694/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: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:698645/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:2867573/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:2977403/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:3152735/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:1338413/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:328332/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:3329117/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:195919/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:2308871/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:1385413/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:2113267/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:11175/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:2002694/1‑69 (MQ=255)
    aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                                    >  1:835165/1‑69 (MQ=255)
            aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                                            <  1:137587/69‑1 (MQ=255)
            aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                                            <  1:1157035/69‑1 (MQ=255)
            aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                                            <  1:3354655/69‑1 (MQ=255)
            aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                                            <  1:340478/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)
                        aTAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCgg                                 <  1:760456/68‑1 (MQ=255)
                                   tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGGTCGCAGCCGGAGTGGCGCAATc                     <  1:320808/69‑1 (MQ=255)
                                   tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCaa                       >  1:2511302/1‑67 (MQ=255)
                                   tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc                     <  1:1216766/69‑1 (MQ=255)
                                   tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc                     <  1:3369027/69‑1 (MQ=255)
                                   tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc                     <  1:3136503/69‑1 (MQ=255)
                                   tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc                     <  1:1130996/69‑1 (MQ=255)
                                   tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc                     <  1:1401037/69‑1 (MQ=255)
                                     ttGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGgcgc                         >  1:113945/1‑63 (MQ=255)
                                     ttGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATccc                   >  1:2243890/1‑69 (MQ=255)
                                     ttGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATccc                   >  1:2533204/1‑69 (MQ=255)
                                       gCACAAATGATAACCGAATTTGTGTTTATCCCGAtt                                                  >  1:1790987/1‑36 (MQ=255)
                                       gCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCTCGATCGCAGCCg                                  >  1:3361191/1‑52 (MQ=255)
                                       gCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCt                  >  1:1434946/1‑68 (MQ=255)
                                       gCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCt                  >  1:2698016/1‑68 (MQ=255)
                                         acaAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGc                <  1:352043/68‑1 (MQ=255)
                                         acaAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGc                <  1:341707/68‑1 (MQ=255)
                                         acaAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGc                <  1:703928/68‑1 (MQ=255)
                                           aaaTGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGc                <  1:1043684/66‑1 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:569191/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:1829677/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:1865510/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:1878785/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:2054511/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:205769/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:2350546/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:756208/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:27214/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:2784632/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:2910474/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:466310/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:3019093/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:36178/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:3238568/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:341037/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATa            >  1:3332988/1‑68 (MQ=255)
                                             aTGATAACCACATTTGTGTTAATCCCGATTTTCGCGATAGCAGCCGGAGTgg                            >  1:3283915/1‑52 (MQ=37)
                                                      gAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATACTTAAATCgg  <  1:3367560/69‑1 (MQ=255)
                                                      |                                                                    
CAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGCAATACTTAAATCGG  >  minE/2070051‑2070173

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: