Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 555,570 A→C 37.0% G159G (GGT→GGG moeB ← molybdopterin synthase sulfurylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE555,5700AC37.0% 22.3 / 37.5 46G159G (GGT→GGGmoeBmolybdopterin synthase sulfurylase
Reads supporting (aligned to +/- strand):  ref base A (15/14);  new base C (5/12);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

ACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTT  >  minE/555504‑555611
                                                                  |                                         
aCCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGAc                                          <  1:293081/68‑1 (MQ=255)
aCCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGAc                                          <  1:174719/68‑1 (MQ=255)
      caaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCAt                                   >  1:1143387/1‑69 (MQ=255)
      caaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCAt                                   >  1:2024660/1‑69 (MQ=255)
      caaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCAt                                   >  1:1724050/1‑69 (MQ=255)
          cGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGa                               <  1:2935658/69‑1 (MQ=255)
          cGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGa                               <  1:1477269/69‑1 (MQ=255)
                                gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:2772746/1‑68 (MQ=255)
                                gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:2971881/1‑68 (MQ=255)
                                gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:3085591/1‑68 (MQ=255)
                                gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:3173403/1‑68 (MQ=255)
                                gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:1764776/1‑68 (MQ=255)
                                gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:3436078/1‑68 (MQ=255)
                                gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:126520/1‑68 (MQ=255)
                                  tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:3396860/69‑1 (MQ=255)
                                  tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:739646/69‑1 (MQ=255)
                                  tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:804407/69‑1 (MQ=255)
                                  tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:1342001/69‑1 (MQ=255)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3297646/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3286226/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:77179/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3074155/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3055848/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:837120/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:1175198/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:2664772/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:2655228/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:256270/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:2519359/36‑1 (MQ=37)
                                      cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:2471984/36‑1 (MQ=37)
                                       gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:194711/1‑69 (MQ=255)
                                       gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:1367404/1‑69 (MQ=255)
                                       gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:696539/1‑69 (MQ=255)
                                       gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:3045963/1‑69 (MQ=255)
                                          cTGATAAGTAAAGACGGTGATTTGACCTTCCATACg                                >  1:1641317/1‑36 (MQ=255)
                                                    aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCa            >  1:1923185/1‑46 (MQ=39)
                                                    aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:2341024/1‑47 (MQ=39)
                                                    aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:2550149/1‑47 (MQ=39)
                                                    aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:1393604/1‑47 (MQ=39)
                                                    aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:1314686/1‑47 (MQ=39)
                                                                tGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:2411290/39‑1 (MQ=255)
                                                                tGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:2829814/39‑1 (MQ=255)
                                                                tGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:361841/39‑1 (MQ=255)
                                                                tGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:15043/39‑1 (MQ=255)
                                                                tGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:752060/39‑1 (MQ=255)
                                                                tGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:894173/39‑1 (MQ=255)
                                                                  |                                         
ACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTT  >  minE/555504‑555611

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: