Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 973,326 | C→G | 33.3% | G26R (GGC→CGC) | elyC ← | envelope biogenesis factor; DUF218 superfamily protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 973,326 | 0 | C | G | 33.3% | 74.0 / 55.0 | 54 | G26R (GGC→CGC) | elyC | envelope biogenesis factor; DUF218 superfamily protein |
Reads supporting (aligned to +/- strand): ref base C (19/17); new base G (5/13); total (24/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.45e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTC > NC_000913/973179‑973466 | aggcatttacgggagaaaatatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgccaggcc < 8:302529‑M2/56‑1 (MQ=255) tGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCa < 3:98343/148‑1 (MQ=255) gaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATca < 2:288304‑M2/70‑1 (MQ=255) aGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATca < 4:91916/147‑1 (MQ=255) aGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATcaa > 6:286024/1‑148 (MQ=255) ggTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAAAAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATcaa < 4:35205/147‑1 (MQ=255) atcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAacaaca > 7:190160‑M2/73‑148 (MQ=255) gggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCgg > 5:7953/1‑147 (MQ=255) aaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa < 8:190160‑M2/84‑1 (MQ=255) ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGa < 1:180842/147‑1 (MQ=255) ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGa < 4:62061/147‑1 (MQ=255) ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa < 5:286024/147‑1 (MQ=255) ggACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa < 6:7953/147‑1 (MQ=255) aacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCgg > 7:11862‑M2/64‑146 (MQ=255) gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTgg > 4:145359‑M2/59‑147 (MQ=255) gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTgg > 4:140572‑M2/59‑147 (MQ=255) aaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGgcagc > 1:64944‑M2/55‑147 (MQ=255) atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAAc < 2:64944‑M2/96‑1 (MQ=255) gTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa < 3:253876/133‑1 (MQ=255) gTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa > 4:253876/1‑133 (MQ=255) gTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACata > 5:163633/1‑148 (MQ=255) ggCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTg < 2:244383/148‑1 (MQ=255) ccgagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCacttat > 2:121181‑M2/35‑146 (MQ=255) gagaagaatttccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt > 1:212028‑M2/33‑148 (MQ=255) aaCTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt > 2:9042/1‑148 (MQ=255) aCTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt > 2:40453/1‑147 (MQ=255) aaTAACAAGAGCTCCAGACACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAAGAACTTCTGCGGTAGTGGCAGCAACATATTGCCTATCACTTTTTTCAgt < 8:83806/147‑1 (MQ=255) tAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAgtgt > 7:306366/1‑147 (MQ=255) aCAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaa < 1:17181/147‑1 (MQ=255) ccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCTCAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaaaa < 3:140572‑M2/126‑1 (MQ=255) ccccgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTaaaa < 3:145359‑M2/126‑1 (MQ=255) cgcttattcgcaccttccCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTACACCACAGCAGATCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGACAATCACTTTTTTCAGTGTAAAAAGCa < 5:39181‑M2/130‑1 (MQ=255) gcgcCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCAt < 1:9042/148‑1 (MQ=255) agccaCCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCttt < 5:104762/148‑1 (MQ=255) gccaCCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTc < 6:163633/148‑1 (MQ=255) ccaCCCAATGCTAATAAAAATCCTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCt > 7:44345/1‑148 (MQ=255) ccAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGT‑AAAAGCATCCtt > 8:235919/1‑140 (MQ=255) ccAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGT‑AAAAGCATCCtt < 7:235919/140‑1 (MQ=255) tGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtt > 2:320679/1‑148 (MQ=255) cTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACttttt < 5:263825/115‑1 (MQ=255) cTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACttttt > 6:263825/1‑115 (MQ=255) cTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttt > 8:141549/1‑148 (MQ=255) aataaAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCttttt > 1:18158/1‑147 (MQ=255) ataaAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGAGTAAAAAGCATCCTTTCTGGTTCCTTTTTTa < 1:40453/148‑1 (MQ=255) ataaAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttttt > 2:12115/1‑147 (MQ=255) aaaaaTCTTCCCGGGTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTACTGTTTGAAc > 6:223336/1‑148 (MQ=255) aaaTCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCC‑TTTTTAACCat < 7:141549/148‑1 (MQ=255) cTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCata > 2:153364/1‑146 (MQ=255) ggTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAgg < 8:306366/148‑1 (MQ=255) ttttGAAAACGGCTAAACCACAGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGAt < 7:286277/147‑1 (MQ=255) aaCGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGc > 4:50976/1‑147 (MQ=255) cTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACCTCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTCAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATAtt < 5:223336/148‑1 (MQ=255) cTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATAtt < 8:11862/148‑1 (MQ=255) acaGCAGCGCCAGGCGAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAg < 1:320679/148‑1 (MQ=255) acaGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAg < 1:153364/148‑1 (MQ=255) gccaggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTc < 3:50976/148‑1 (MQ=255) ggccagCGCCAATAATCAACAACATCAGCGGGAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGata < 1:135457/127‑1 (MQ=255) ggccagCGCCAATAATCAACAACATCAGCGGGAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGata > 2:135457/1‑127 (MQ=255) | TATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGATCGCCAGCGATATTACACCAGACGGGTTC > NC_000913/973179‑973466 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |