Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,070,092 A→G 47.1% intergenic (‑303/‑88) 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,0920AG47.1% 5.6 / 92.5 104intergenic (‑303/‑88)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 (29/26);  new base G (22/27);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGC  >  minE/2070029‑2070159
                                                               |                                                                   
tCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAAt                                                                <  1:1392963/69‑1 (MQ=255)
tCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAAt                                                                <  1:806804/69‑1 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:1018689/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:3209121/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:1600007/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:3197843/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:1769487/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:1951680/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:2741026/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:2588151/1‑68 (MQ=255)
  cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                                                               >  1:1118800/1‑68 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:59250/68‑1 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:3055744/68‑1 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:1132026/68‑1 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:393717/68‑1 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:3115004/68‑1 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                                            <  1:54872/69‑1 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                                            <  1:3267817/69‑1 (MQ=255)
    tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                                            <  1:836775/69‑1 (MQ=255)
      cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                                           >  1:2948642/1‑68 (MQ=255)
      cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                                           >  1:2967628/1‑68 (MQ=255)
      cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                                           >  1:2551048/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:1173887/1‑68 (MQ=255)
      cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                                           >  1:814663/1‑68 (MQ=255)
      cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                                                           >  1:819269/1‑68 (MQ=255)
      cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                                                            >  1:737031/1‑67 (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:781319/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:245544/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:2316557/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:826156/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1086833/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1128694/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1693450/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:2141354/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1682191/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1632603/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1628737/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1419299/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1481463/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1562140/66‑1 (MQ=255)
                      cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:1519008/66‑1 (MQ=255)
                       aaaaTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt                                             <  1:3318837/65‑1 (MQ=255)
                          aaaacagcAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2623487/3‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCg                                       >  1:2483077/1‑68 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:3152735/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:11175/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:3329117/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:1385413/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:1338413/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:698645/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:1215365/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:835165/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:195919/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2429916/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:253330/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2977403/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2113267/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2554452/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2002694/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2867573/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2687268/1‑69 (MQ=255)
                          aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa                                      >  1:2308871/1‑69 (MQ=255)
                          aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:2188078/46‑1 (MQ=255)
                          aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:2081250/46‑1 (MQ=255)
                          aTAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                                                             <  1:2861802/46‑1 (MQ=255)
                                  aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                              <  1:3354655/69‑1 (MQ=255)
                                  aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                              <  1:2047333/69‑1 (MQ=255)
                                  aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                              <  1:1157035/69‑1 (MQ=255)
                                  aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                              <  1:137587/69‑1 (MQ=255)
                                  aaaTGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTcgcg                              <  1:340478/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:1401037/69‑1 (MQ=255)
                                                         tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc       <  1:1130996/69‑1 (MQ=255)
                                                         tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc       <  1:3369027/69‑1 (MQ=255)
                                                         tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc       <  1:3136503/69‑1 (MQ=255)
                                                         tatTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATc       <  1:1216766/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:703928/68‑1 (MQ=255)
                                                               acaAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGc  <  1:352043/68‑1 (MQ=255)
                                                               acaAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGc  <  1:341707/68‑1 (MQ=255)
                                                               |                                                                   
TCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGATTTTCGCGATCGCAGCCGGAGTGGCGCAATCCCTGC  >  minE/2070029‑2070159

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: