Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,659,376 A→G 18.4% G35G (GGA→GGG glpE → thiosulfate:cyanide sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,659,3760AG18.4% 106.2 / 13.7 70G35G (GGA→GGGglpEthiosulfate:cyanide sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base G (7/6);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AAAGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACTTT  >  minE/2659326‑2659442
                                                  |                                                                  
aaaGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:936934/69‑1 (MQ=255)
aaaGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:330506/69‑1 (MQ=255)
aaaGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:1860271/69‑1 (MQ=255)
aaaGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:2868697/69‑1 (MQ=255)
 aaGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:1496149/68‑1 (MQ=255)
 aaGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:2535425/68‑1 (MQ=255)
 aaGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:2139971/68‑1 (MQ=255)
   gagGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTc                                                  <  1:1216822/66‑1 (MQ=255)
     ggCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCAttt                                             <  1:58649/69‑1 (MQ=255)
      gCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCAttt                                             <  1:99294/68‑1 (MQ=255)
         gTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAAcc                                         <  1:3235994/69‑1 (MQ=255)
               gTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAAcc                                         >  1:2357366/1‑63 (MQ=255)
               gTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGac                                   >  1:1903726/1‑69 (MQ=255)
               gTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGac                                   >  1:997161/1‑69 (MQ=255)
               gTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGac                                   >  1:214220/1‑69 (MQ=255)
               gTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGac                                   >  1:302850/1‑69 (MQ=255)
               gTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGac                                   >  1:1090152/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCt                              >  1:1220219/1‑68 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:2067607/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:1227589/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:936053/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:726342/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:618030/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:1358864/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCTGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:1913820/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCt                              >  1:2498537/1‑68 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:2340517/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:2068146/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:524992/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:49012/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:41235/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:2897986/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:2300597/1‑69 (MQ=255)
                     atTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTg                             >  1:1638336/1‑69 (MQ=255)
                        cgcgATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTggg                           >  1:644892/1‑68 (MQ=255)
                        cgcgATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGgc                          >  1:2785548/1‑69 (MQ=255)
                        cgcgATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGgc                          >  1:987494/1‑69 (MQ=255)
                        cgcgATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGgc                          >  1:728753/1‑69 (MQ=255)
                        cgcgATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGgc                          >  1:546759/1‑69 (MQ=255)
                        cgcgATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGgc                          >  1:458232/1‑69 (MQ=255)
                        cgcgATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATATAACCAACGACACGCTGGgc                          >  1:502457/1‑69 (MQ=255)
                             tCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCttt                     <  1:2743583/69‑1 (MQ=255)
                             tCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCttt                     <  1:2229586/69‑1 (MQ=255)
                             tCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCttt                     <  1:145777/69‑1 (MQ=255)
                               cacaGAGTTTCGCAATGGGACTTGCGGTGCAGGCTTTCCAttt                                             <  1:948728/43‑1 (MQ=255)
                                 caGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTAt                  >  1:1659580/1‑68 (MQ=255)
                                 caGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGAGCTTTTATg                 >  1:1889184/1‑69 (MQ=255)
                                     gTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTg             <  1:3305477/69‑1 (MQ=255)
                                     gTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTg             <  1:1390352/69‑1 (MQ=255)
                                     gTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTg             <  1:1804686/69‑1 (MQ=255)
                                     gTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTg             <  1:1471235/69‑1 (MQ=255)
                                         cGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCa                                        <  1:2595721/38‑1 (MQ=255)
                                         cGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCa                                        <  1:97085/38‑1 (MQ=255)
                                         cGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCa                                        <  1:870151/38‑1 (MQ=255)
                                         cGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCa                                        <  1:2592454/38‑1 (MQ=255)
                                         cGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCa                                        <  1:1773169/38‑1 (MQ=255)
                                         cGCAATGGGGCATGCGGTGCAGGCTTTCCATTTAACCa                                        <  1:3249597/38‑1 (MQ=255)
                                              tGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACt    <  1:690239/69‑1 (MQ=255)
                                              tGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACt    <  1:2785645/69‑1 (MQ=255)
                                                ggACATGCGTTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:345489/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:3210858/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:2432752/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:764716/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:820108/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:2133329/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:1311322/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:1585205/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:443683/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:2688351/69‑1 (MQ=255)
                                                ggACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACttt  <  1:366963/69‑1 (MQ=255)
                                                  |                                                                  
AAAGAGGCGGTGCTGGTCGATATTCGCGATCCACAGAGTTTCGCAATGGGACATGCGGTGCAGGCTTTCCATTTAACCAACGACACGCTGGGCGCTTTTATGCGTGATAACGACTTT  >  minE/2659326‑2659442

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: