Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,491,041 | Δ1 bp | 18.2% | coding (697/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,041 | 0 | G | . | 18.2% | 54.0 / 10.0 | 22 | coding (697/1902 nt) | cysNC | sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase |
Reads supporting (aligned to +/- strand): ref base G (9/9); new base . (2/2); total (11/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01 |
CTACCTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATT > NC_002947/1490893‑1491172 | cTACCGGAATTTCGGTGGAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGt < 2:314755/150‑1 (MQ=255) 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) | CTACCTGAAGTTCGCCGAAGCCATCAACCTGACCCCGAGCAGCCTGCACTTCGTGCCGATGTCGGCGCTCAAGGGTGACAACGTGGTCAACCACAGCGAGCAATCGCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGA‑GGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATT > NC_002947/1490893‑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: ‑ |