Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,159 | G→A | 15.8% | intergenic (+65/‑82) | purH → / → purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,485,159 | 0 | G | A | 15.8% | 83.6 / 8.2 | 38 | intergenic (+65/‑82) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base G (14/18); new base A (2/4); total (16/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01 |
CTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCA > NC_002947/5485021‑5485302 | cTGAAGTCATCGCTGCTGCCGACGAGGCCTGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATGCGGCGGTGAATCCAccatc < 2:155609/149‑4 (MQ=38) tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCAccgcc > 1:216128/1‑149 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCC‑GCGAATAAGAATGTGGATTAACACCCCCATTagc > 2:73660/1‑149 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTcgc < 2:66335/149‑1 (MQ=255) tgcCGACGAGGCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTcgc > 1:212329/1‑141 (MQ=255) tgcCGACGAGGCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTcgc < 2:212329/141‑1 (MQ=255) ccGACGAGGCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc < 1:30329/147‑1 (MQ=25) aGGCCGGCATCGC‑GTTTGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCAGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTcc < 1:73660/149‑1 (MQ=39) gCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGAAtc < 8:162718/116‑2 (MQ=255) gCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGAAtc > 7:162718/1‑115 (MQ=255) gCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGAAtc < 4:202317/116‑2 (MQ=255) gCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGAAtc > 3:202317/1‑115 (MQ=255) c‑gATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATACCCCGTAACCCGAGGCATGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgc < 3:43037/149‑1 (MQ=38) gTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGTATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTGGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCaa < 1:39764/148‑1 (MQ=255) cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCC‑GCGGATACGATGGGGGATTAACCGCCCCCTTAGCGGGCAAGCCCGCGGCCCCAGTGTTCGGCGCAAGCCCCg > 5:72366/1‑149 (MQ=255) gCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCACGCTCCCACAGTGTTCGGCGCAAGCCACGGGAt < 2:216128/148‑1 (MQ=255) gCATCCGCCACTTCCGCCACTAATTCCGCGTATCCCGAGGCAAGTCCAGATCCTGTGGGAGCGGGCTTGCCC‑GCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGAt < 2:258135/148‑1 (MQ=38) gcgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCggg < 4:123755/105‑1 (MQ=255) ccACTTCCGCCACTAATTACGCGCATCCCGTGGCAAGGGGAGATCGTGTGGGTGCGGGCTTGCCG‑GCTAATACGATGGTGGATTCCCGGCCGCATGGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGaa < 2:237479/149‑1 (MQ=255) cACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt < 5:260702/149‑1 (MQ=255) cctcggCAAGTCGAGATCCTGGGGGCGCGGGTTGGCCC‑GCTAATACGATGGTGGATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt < 7:32952/142‑1 (MQ=255) acctaaGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGat‑cat < 5:88087/145‑1 (MQ=255) tggcggCGGGGTTGCCC‑GCTAATCCGAGGGTGGTTCCCCCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGa < 2:285620/144‑1 (MQ=255) cTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCcctg > 2:200653/1‑149 (MQ=255) ttGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCcctgg > 7:45611/1‑149 (MQ=255) ccc‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCTCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCcctggcct > 1:225069/1‑149 (MQ=255) ccggctaaTACTAGGGGGGATTCCCCGCGGCATTCGCGGGCATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCcctggcct < 4:108989/143‑1 (MQ=255) cc‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGat‑cat < 6:92046/116‑1 (MQ=255) cc‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGat‑cat > 5:92046/1‑116 (MQ=255) cc‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCAc < 1:153212/138‑1 (MQ=255) cc‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCAc > 2:153212/1‑138 (MQ=255) cc‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:167511/1‑149 (MQ=255) cc‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 2:215174/1‑149 (MQ=255) gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg > 2:74217/1‑148 (MQ=255) gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa < 5:167511/149‑1 (MQ=255) gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa > 3:261067/1‑149 (MQ=255) gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa > 5:168831/1‑147 (MQ=255) tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc < 1:200653/149‑1 (MQ=255) gtggTTTCAACGCCGGATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGGGCAATCCAAGGGTTCACTTGTATCGCTGTTTTGACATGAAAGTTTTGATCCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa < 1:74217/148‑1 (MQ=255) | CTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGC‑GATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCC‑GCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGAT‑CATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCA > NC_002947/5485021‑5485302 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |