Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,740,484 | T→C | 15.1% | intergenic (+49/+22) | hutF → / ← PP_5037 | formiminoglutamate deiminase/lipocalin family lipoprotein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,740,484 | 0 | T | C | 15.1% | 70.2 / 7.1 | 34 | intergenic (+49/+22) | hutF/PP_5037 | formiminoglutamate deiminase/lipocalin family lipoprotein |
Reads supporting (aligned to +/- strand): ref base T (8/19); new base C (3/2); total (13/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.10e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01 |
GGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGA > NC_002947/5740339‑5740631 | ggTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GGTGGAGCCC‑ATAGCGTGGAAGCCCGCGCAGGGttt > 7:193959/1‑147 (MQ=255) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGGGAGCCCGCGCAGGGGTggaa‑c > 7:241067/1‑149 (MQ=255) gtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACt < 3:188151/149‑1 (MQ=255) gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGGTTGC‑GCCC‑TTAGCGGGGGAGCCCGCGCAGGGGTTGAA‑CGCCACCCCaaaaaaaa > 3:281431/1‑145 (MQ=255) gcgcCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTTGAGCCC‑TTAGCGTGGCAGCCCGCGAAGGGTTGGAA‑CGCAACCCCCAAAAAca > 1:91901/1‑146 (MQ=255) ctggggggGAAGGCCCTGTAAGTTCAAGGGGGGGGGGGGTTTGCCTTGGTGTGGCGCCTGTTTTTGGGGTGTGCGTAGCGGGGTGGG‑GTTGCAGCCT‑TTCGGGTCCAAGCCCGCTATGGTTTTAAA‑CGCAACCCCAAAAAACTCAGTGc < 8:116320/147‑1 (MQ=255) ggggCGGCATGCCCAGTATGTTTAACGCGGGGGGGGTTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGCGGTTACGGGGGTTGG‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCTATGGGTTGAAA‑CGCAACCCCAAAAAACTCAGTGCAt < 4:281431/149‑1 (MQ=255) gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTgcaa‑c > 8:72293/1‑119 (MQ=255) gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTgcaa‑c < 7:72293/119‑1 (MQ=255) cAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCTAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc < 3:70629/149‑1 (MQ=255) aGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTCCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAAGCGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGAt < 1:220370/149‑1 (MQ=255) gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGG‑GGTGGAGCCCATACGCG‑GGAAGGCCGCGCAGGGTTTGAA‑CGCCAACCCAAAAAAAGCAGTGGAGTTTTAAATTTTGATGATggggc > 3:26918/1‑147 (MQ=255) gcgggcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa < 5:162650/149‑1 (MQ=255) gcgggcgTTGGCCAAGGTGTTCCGCCAGCTTTTGGGCTGCGCGTACCGGGGTTGC‑GTTGCAGCCC‑TGCGCGTGCAAGCCCGCTAAGGTTTGAAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCa < 4:113214/149‑1 (MQ=255) gcgTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGGGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACt < 5:150461/107‑1 (MQ=255) gTTTGCCAAGGTGTGGCGCCAGCTTTTGGGCTTAGCGTACCGGGGTGGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCaa < 3:281810/149‑1 (MQ=255) tGCCAAGGGTGTGCGCCCGCTTTTGGGCTGAGCGTACCGGGGTTGC‑TTTGCACCCC‑TTCGCGTGCAAGCCCGCTATGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc < 5:176933/148‑1 (MQ=255) ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc < 1:6861/148‑1 (MQ=255) ccAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGCGCAACCCCGCGACGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc < 4:26918/149‑1 (MQ=255) aGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGc < 3:113007/118‑1 (MQ=255) gggTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag < 2:91901/149‑1 (MQ=255) gggTTGC‑GTTCCAGCCC‑TTCGCGTGCAACCCCGCTATGGGTTGAAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagcag < 3:145409/149‑1 (MQ=255) gttgcAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCa > 5:157943/1‑149 (MQ=255) aaccc‑TTCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCa < 8:61840/92‑1 (MQ=21) aaccc‑TTCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCa > 7:61840/3‑94 (MQ=21) ttCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg > 1:84715/1‑149 (MQ=255) ttCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg > 5:186865/1‑133 (MQ=38) ttCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAgg < 6:186865/133‑1 (MQ=38) ttCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg > 7:136618/1‑149 (MQ=39) aaGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAgtg > 2:152333/1‑149 (MQ=255) gCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTc > 4:203417/1‑149 (MQ=255) cccGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCg > 5:156431/1‑149 (MQ=255) cgAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt < 1:152333/149‑1 (MQ=255) gggTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTgcgc < 3:203417/149‑1 (MQ=255) gTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa < 6:157943/149‑1 (MQ=255) ttGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGa > 7:28938/1‑148 (MQ=255) | GGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGC‑GTTGCAGCCC‑TTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGTGCGCGA > NC_002947/5740339‑5740631 |
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: ‑ |