Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,070,061 A→G 38.6% intergenic (‑272/‑119) 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,0610AG38.6% 57.6 / 81.4 119intergenic (‑272/‑119)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 (42/31);  new base G (21/25);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 120 total aligned reads displayed.

AACGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGA  >  minE/2070001‑2070123
                                                            |                                                              
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:930250/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:594081/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:486511/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:3429686/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:198349/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:3339883/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:2154987/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:1359453/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:3054298/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:293407/69‑1 (MQ=255)
aaCGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGtt                                                        <  1:26238/69‑1 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:1853038/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:1046651/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:1887965/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:2439189/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:2459646/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:2626333/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:3273477/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:328131/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:334951/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:3379202/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:933862/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:1634550/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:1261015/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:1304507/1‑67 (MQ=255)
      cATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACa                                                    >  1:157563/1‑67 (MQ=255)
          tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                              >  1:3414828/1‑69 (MQ=255)
          tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                              >  1:1803614/1‑69 (MQ=255)
          tGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGc                                              >  1:2112052/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:3209121/1‑68 (MQ=255)
                              cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                           >  1:3197843/1‑68 (MQ=255)
                              cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                           >  1:2741026/1‑68 (MQ=255)
                              cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                           >  1:1769487/1‑68 (MQ=255)
                              cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                           >  1:2588151/1‑68 (MQ=255)
                              cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                           >  1:1018689/1‑68 (MQ=255)
                              cATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATg                           >  1:1118800/1‑68 (MQ=255)
                                tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                         <  1:1132026/68‑1 (MQ=255)
                                tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                         <  1:3115004/68‑1 (MQ=255)
                                tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                         <  1:393717/68‑1 (MQ=255)
                                tACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGAt                         <  1:59250/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:416602/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                       >  1:2967628/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                       >  1:2944721/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                       >  1:814663/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                       >  1:2221089/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                       >  1:1173887/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                       >  1:2551048/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATaa                       >  1:2224997/1‑68 (MQ=255)
                                  cccGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATa                        >  1:737031/1‑67 (MQ=255)
                                        aaaaaCCGTGCAAAATAACAACAAATGTTAACTGATAGCATTaaa                                        <  1:2991832/45‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:781319/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:3265028/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:826156/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:3402816/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:47514/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1086833/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1682191/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1519008/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1693450/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1632603/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1628737/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1888012/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:2141354/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1562140/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:2316557/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:2356503/66‑1 (MQ=255)
                                                  cAAAATAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGttt         <  1:1481463/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:2977403/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:328332/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:3152735/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:835165/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:2867573/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:195919/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:2002694/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:698645/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:11175/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:2687268/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:2308871/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:1338413/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:2429916/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:1385413/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:253330/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:2554452/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:1215365/1‑69 (MQ=255)
                                                      aTAACAGCAAATGTTAGCAAATAGCATTAAATATTGCGCAAATGATAACCACATTTGTGTTTATCCCGa  >  1:3329117/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)
                                                            |                                                              
AACGTGCATTTGCGCAGCAATCATCAAATCCATACCCGACAAAAACCGTGCAAAATAACAACAAATGTTAACAGATAGCATTAAATATTGCACAAATGATAACCGAATTTGTGTTTATCCCGA  >  minE/2070001‑2070123

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: