Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 556,511 C→T 42.0% E257E (GAG→GAA moeA ← molybdopterin biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE556,5110CT42.0% 18.3 / 39.6 57E257E (GAG→GAAmoeAmolybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (10/23);  new base T (8/16);  total (18/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATA  >  minE/556449‑556571
                                                              |                                                            
aTCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCAcc                                                         <  1:2690015/68‑1 (MQ=255)
aTCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCAcc                                                         <  1:3302994/68‑1 (MQ=255)
  cGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCa                                                       <  1:698655/68‑1 (MQ=255)
  cGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCa                                                       <  1:2504372/68‑1 (MQ=255)
  cGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCa                                                       <  1:2129134/68‑1 (MQ=255)
  cGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCa                                                       <  1:1649141/68‑1 (MQ=255)
  cGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCa                                                       <  1:1856744/68‑1 (MQ=255)
  cGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCa                                                       <  1:419889/68‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:1898531/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:3017527/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:2990743/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:3185615/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:274678/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:2104373/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:3400862/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:192328/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:1735727/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:731475/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:1414877/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:570353/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:1481796/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:1461643/69‑1 (MQ=255)
     cAGTTTCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:621683/69‑1 (MQ=255)
     cAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGa                                                   <  1:3038713/69‑1 (MQ=255)
     cAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGa                                                   <  1:1116514/69‑1 (MQ=255)
     cAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGa                                                   <  1:1213267/69‑1 (MQ=255)
      agtttCCAGAAGGCAATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCTTCACCCACTGa                                                   <  1:1717554/65‑1 (MQ=255)
             aGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACAcc                                             <  1:422246/67‑1 (MQ=255)
                     aTCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTg                                   <  1:1141994/69‑1 (MQ=255)
                     aTCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTg                                   <  1:2562687/69‑1 (MQ=255)
                     aTCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTg                                   <  1:1870375/69‑1 (MQ=255)
                        tcCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCgg                                        <  1:2137207/61‑1 (MQ=255)
                        tcCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCgg                                        <  1:384394/61‑1 (MQ=255)
                        tcCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCgg                                        <  1:20752/61‑1 (MQ=255)
                        tcCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCgg                                        <  1:1967463/61‑1 (MQ=255)
                        tcCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCgg                                        <  1:385313/61‑1 (MQ=255)
                            ccAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              <  1:2845622/67‑1 (MQ=255)
                                     cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:3206900/1‑58 (MQ=255)
                                     cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:918196/1‑58 (MQ=255)
                                     cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:275273/1‑58 (MQ=255)
                                     cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:261161/1‑58 (MQ=255)
                                     cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:3006625/1‑58 (MQ=255)
                                     cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:1801977/1‑58 (MQ=255)
                                     cAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:1318788/1‑58 (MQ=255)
                                     cAAGAATCGTTTTGGCGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcac                              >  1:134649/1‑58 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGctt                   >  1:1452156/1‑66 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCTTg                  >  1:2465795/1‑69 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCTTg                  >  1:2092400/1‑69 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCTTg                  >  1:75684/1‑69 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCTTg                  >  1:1748377/1‑69 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCTTg                  >  1:1608289/1‑69 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCTTg                  >  1:591483/1‑69 (MQ=255)
                                      aaGAATCGTTTTGGTGTAATCCGCTTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCTTg                  >  1:612180/1‑69 (MQ=255)
                                        gAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcaccac                           >  1:1165985/1‑58 (MQ=255)
                                        gAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGc                >  1:2082795/1‑69 (MQ=255)
                                                      tAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATa  <  1:447892/69‑1 (MQ=255)
                                                      tAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATa  <  1:1850942/69‑1 (MQ=255)
                                                              |                                                            
ATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGGCTGTCGGCTTCAATA  >  minE/556449‑556571

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: