Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 523,681 A→T 41.3% R61S (AGA→AGT moaA → molybdopterin biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE523,6810AT41.3% 21.2 / 41.9 46R61S (AGA→AGTmoaAmolybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (12/15);  new base T (9/10);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCG  >  minE/523627‑523742
                                                      |                                                             
tAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGTCACCGaaaa                                                 <  1:3179033/69‑1 (MQ=255)
tAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGaaaa                                                 <  1:2304756/69‑1 (MQ=255)
tAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGaaaa                                                 <  1:3081576/69‑1 (MQ=255)
tAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGaaaa                                                 <  1:811758/69‑1 (MQ=255)
tAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGaaaa                                                 <  1:3273726/69‑1 (MQ=255)
 aaaGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGaaaa                                                 <  1:2721583/68‑1 (MQ=255)
 aaaGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGaaaa                                                 <  1:2226454/68‑1 (MQ=255)
  aaGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGaaaa                                                 <  1:1302126/67‑1 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:948202/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:2000120/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:1922552/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:398531/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:694937/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:1638441/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:1500313/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:1498949/1‑68 (MQ=255)
    ggCTTTCTGACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGTCTGGGCACCGAAAAAGtc                                             >  1:1408562/1‑68 (MQ=255)
                   gATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAATAGTGCGCCTGACaggag                                >  1:1759962/1‑67 (MQ=255)
                   gATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAATAGTGCGCCTGACaggag                                >  1:1763868/1‑67 (MQ=255)
                   gATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAg                                               >  1:2744425/1‑52 (MQ=255)
                   gATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACaggag                                >  1:746323/1‑67 (MQ=255)
                   gATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACaggag                                >  1:2121168/1‑67 (MQ=255)
                       aaaTTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGc                          >  1:2415635/1‑69 (MQ=255)
                         aTTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGTGGAGAGcc                         <  1:1074566/68‑1 (MQ=255)
                         aTTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGcc                         <  1:2623213/68‑1 (MQ=255)
                         aTTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGcc                         <  1:1990118/68‑1 (MQ=255)
                         aTTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGcc                         <  1:2244711/68‑1 (MQ=255)
                            cGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTc                      >  1:2755459/1‑68 (MQ=255)
                            cGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTc                      >  1:1666902/1‑68 (MQ=255)
                            cGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTc                      >  1:1676520/1‑68 (MQ=255)
                            cGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTc                      >  1:3342189/1‑68 (MQ=255)
                            cGCCGGGTTACGCGAGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTc                      >  1:2177634/1‑68 (MQ=255)
                            cGCCGGATTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTc                      >  1:7102/1‑68 (MQ=255)
                              ccGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTCtt                    <  1:1697714/68‑1 (MQ=255)
                              ccGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGACGTCtt                    <  1:1748/68‑1 (MQ=255)
                               cGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTCtt                    <  1:2122027/67‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTAcct  <  1:47675/69‑2 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:481750/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:637208/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:686398/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:3387233/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:1798732/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:3198594/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:781625/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:2936322/69‑1 (MQ=255)
                                               tCGCCAGTCTGGGCACCGAAAAAGTCCGTCTGACAGGTGGTGAGCCATCTTTACGCCGCGACTTTACCg  <  1:1047594/69‑1 (MQ=255)
                                                      |                                                             
TAAAGGCTTTCTTACCGTCGATGAAATTCGCCGGGTTACGCGCGCCTTCGCCAGACTGGGCACCGAAAAAGTGCGCCTGACAGGAGGAGAGCCGTCTTTACGCCGCGACTTTACCG  >  minE/523627‑523742

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: