Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,491,045:1 | +C | 21.1% | coding (701/1902 nt) | cysNC → | sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,491,045 | 1 | . | C | 21.1% | 39.7 / 13.5 | 19 | coding (701/1902 nt) | cysNC | sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase |
Reads supporting (aligned to +/- strand): ref base . (8/7); new base C (2/2); total (10/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01 |
CCTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATT > NC_002947/1490896‑1491172 | ccTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGa > 1:328863/1‑150 (MQ=255) cTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTAAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑g < 2:190033/150‑1 (MQ=255) aTCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAAc > 2:327228/1‑150 (MQ=255) cTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCg > 2:253688/1‑150 (MQ=255) cTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCAA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCa < 2:39008/150‑1 (MQ=255) ttCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCAc < 1:140999/128‑1 (MQ=255) ttCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCAc > 2:140999/1‑128 (MQ=255) ggCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTcc > 1:164435/1‑150 (MQ=255) cAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCt < 2:141862/150‑1 (MQ=255) gggTGACAACGTGGTCCACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGaa < 2:131423/150‑1 (MQ=255) aaCCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAAc > 2:258266/1‑129 (MQ=255) aaCCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAAc < 1:258266/129‑1 (MQ=255) aCCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCgg < 2:164435/150‑1 (MQ=255) aTCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTAAACCAGTGCGGCTTCGCGGGCACCTTTGccgt > 1:260051/1‑148 (MQ=255) aTGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGAtt > 1:44463/1‑150 (MQ=255) agccGAAACC‑TCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggt > 2:220405/4‑123 (MQ=255) agccGAAACC‑TCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggt < 1:220405/120‑1 (MQ=255) cTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTc < 1:311489/51‑1 (MQ=255) cTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTc > 2:311489/1‑51 (MQ=255) ccGAAACC‑TCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtg > 1:483970/2‑119 (MQ=255) ccGAAACC‑TCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtg < 2:483970/118‑1 (MQ=255) | CCTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATT > NC_002947/1490896‑1491172 |
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: ‑ |