Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 555,561 G→C 33.9% T162T (ACC→ACG moeB ← molybdopterin synthase sulfurylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE555,5610GC33.9% 55.7 / 39.0 56T162T (ACC→ACGmoeBmolybdopterin synthase sulfurylase
Reads supporting (aligned to +/- strand):  ref base G (29/8);  new base C (5/14);  total (34/22)
Fisher's exact test for biased strand distribution p-value = 3.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTT  >  minE/555496‑555611
                                                                 |                                                  
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:571485/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:1709407/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:637982/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:2660852/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:2501468/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:1342046/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:1578065/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:197783/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTg                                                 >  1:2177682/1‑69 (MQ=255)
gCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTACTGATAAGTAAAGACGGTg                                                 >  1:2951479/1‑69 (MQ=255)
 cATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGa                                                >  1:2282100/1‑69 (MQ=255)
 cATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGa                                                >  1:1366356/1‑69 (MQ=255)
 cATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGa                                                >  1:1268109/1‑69 (MQ=255)
 cATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGa                                                >  1:1238294/1‑69 (MQ=255)
        aCCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGAc                                          <  1:293081/68‑1 (MQ=255)
        aCCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGAc                                          <  1:174719/68‑1 (MQ=255)
              caaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCAt                                   >  1:2024660/1‑69 (MQ=255)
              caaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCAt                                   >  1:1143387/1‑69 (MQ=255)
              caaacGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCAt                                   >  1:1724050/1‑69 (MQ=255)
                  cGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGa                               <  1:1477269/69‑1 (MQ=255)
                  cGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGa                               <  1:2935658/69‑1 (MQ=255)
                     cTAAGGCAGTGATAGCACGGTTCGCCGTCCTGATAAGTGAAGACCGTGAt                                               <  1:1574206/50‑1 (MQ=37)
                              tgATAGCACGGTTCGCCGTCCTGATAAGTGAAGACCGTGAt                                               <  1:1862351/40‑1 (MQ=25)
                                        gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:2971881/1‑68 (MQ=255)
                                        gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:1764776/1‑68 (MQ=255)
                                        gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:126520/1‑68 (MQ=255)
                                        gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:2772746/1‑68 (MQ=255)
                                        gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:3085591/1‑68 (MQ=255)
                                        gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:3173403/1‑68 (MQ=255)
                                        gTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGc          >  1:3436078/1‑68 (MQ=255)
                                          tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:1342001/69‑1 (MQ=255)
                                          tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:804407/69‑1 (MQ=255)
                                          tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:3396860/69‑1 (MQ=255)
                                          tCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGt       <  1:739646/69‑1 (MQ=255)
                                              cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3286226/36‑1 (MQ=37)
                                              cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3297646/36‑1 (MQ=37)
                                              cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3074155/36‑1 (MQ=37)
                                              cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:3055848/36‑1 (MQ=37)
                                              cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:77179/36‑1 (MQ=37)
                                              cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:837120/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)
                                              cGTCCTGATAAGTGAAGACCGTGATTTGCCCTTCCa                                    <  1:1175198/36‑1 (MQ=37)
                                               gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:3045963/1‑69 (MQ=255)
                                               gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:194711/1‑69 (MQ=255)
                                               gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:696539/1‑69 (MQ=255)
                                               gTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCtt  >  1:1367404/1‑69 (MQ=255)
                                                  cTGATAAGTAAAGACGGTGATTTGACCTTCCATACg                                >  1:1641317/1‑36 (MQ=255)
                                                            aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCa            >  1:1923185/1‑46 (MQ=39)
                                                            aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:2550149/1‑47 (MQ=39)
                                                            aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:2341024/1‑47 (MQ=39)
                                                            aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:1393604/1‑47 (MQ=39)
                                                            aaGACCGTGATTTGCCCTTCCATACGAATTGCCGCGCCGGAAACCAg           >  1:1314686/1‑47 (MQ=39)
                                                                 |                                                  
GCATTTTCACCAAACAAACGGCTAAGGCAGCGATAGCACGGTTCACCGTCCTGATAAGTAAAGACGGTGATTTGACCTTCCATACGAATTGCCGCGCCGGAAACCAGCGGTACCTT  >  minE/555496‑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: