Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,860,168 | A→G | 13.3% | I309I (ATT→ATC) | ligA ← | DNA ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,860,168 | 0 | A | G | 13.3% | 66.5 / 3.4 | 30 | I309I (ATT→ATC) | ligA | DNA ligase |
Reads supporting (aligned to +/- strand): ref base A (12/14); new base G (2/2); total (14/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01 |
GTTGTGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGC > NC_002947/4860021‑4860316 | gttgtGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAAt > 1:250751/1‑150 (MQ=255) ttgtGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATc > 1:186672/1‑150 (MQ=255) gtgCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGc > 1:123136/1‑150 (MQ=255) tgCAGGGTGGCGTTGGAGACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACCGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGcc < 2:525967/150‑1 (MQ=255) tgCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGcc < 1:353177/150‑1 (MQ=255) gTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTc > 2:43825/1‑150 (MQ=255) cacCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa > 2:452595/1‑135 (MQ=255) cacCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa < 1:452595/138‑4 (MQ=255) cacaCCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACTGGCGCGGaa < 2:586942/150‑1 (MQ=255) caaggCAACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGTCACTCCCCCAGCGCGGCTCACGGGCGCGGAAccc > 2:100295/4‑150 (MQ=255) acCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGcc > 1:74195/1‑150 (MQ=255) ggCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATCACAG‑ACCGCGGGCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAg > 1:536962/1‑150 (MQ=255) ggCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGACACCAACGAACTCGACATC‑AACCACCGCGTTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAg > 1:537185/1‑150 (MQ=255) ccTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTcgc > 1:342316/1‑150 (MQ=255) ttCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCGCGAGCTGggcttcc < 2:139191/150‑5 (MQ=255) cgcCACCGGCGTTACTGCACCCGTACGGCCCACCTGGACCTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCAAGCTCGCGCTGGGCGCCCAGGCTGtt < 2:130244/150‑1 (MQ=255) cGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGaa < 1:441997/150‑1 (MQ=255) aCCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTAcccc < 2:186672/150‑1 (MQ=255) aCGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGAt > 2:494206/1‑150 (MQ=255) gCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGAtctc > 2:524968/1‑150 (MQ=255) tGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGc > 1:37406/1‑150 (MQ=255) gAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAg < 2:537185/150‑1 (MQ=255) gAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAg < 2:536962/150‑1 (MQ=255) aaCTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAgg < 1:524968/150‑1 (MQ=255) aCTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa < 2:385391/64‑4 (MQ=255) aCTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCGATTGCCCacaa > 1:385391/1‑61 (MQ=255) cAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTc < 1:115215/150‑1 (MQ=255) gggAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAgg < 2:250751/150‑1 (MQ=255) ggCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGAGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGAt < 1:62296/150‑1 (MQ=255) cAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGc < 2:37406/150‑1 (MQ=255) | GTTGTGCAGGGTGGCGTTGGACACCGTCACACCGGCCACCTTGACCGGCTTCAGCCGCGCCACAGGCGTTACTGCACCCGTACGGCCAACCTGGAACTCGACATC‑CAGCACCTCGGTCAGCTCTTCCATGGCCGGGAACTTGTGGGCAATCGCCCAGCGCGGCTCACGGGCGCGGAAGCCCAGCTCGCGCTGGGCGGCCAGGCTGTTGACCTTGAACACTACCCCGTCGATCTCATACGGCAGGCTGTTACGCCGCTCGCCAATGTCGCGATAGTAGGCCAGGCACTCTTCGATGC > NC_002947/4860021‑4860316 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |