Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,740,484 | T→C | 24.7% | 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 | 24.7% | 50.9 / 15.4 | 33 | intergenic (+49/+22) | hutF/PP_5037 | formiminoglutamate deiminase/lipocalin family lipoprotein |
Reads supporting (aligned to +/- strand): ref base T (9/15); new base C (5/3); total (15/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.52e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01 |
TGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGT > NC_002947/5740338‑5740625 | tGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTg < 1:188850/149‑1 (MQ=255) tGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa < 2:240160/149‑1 (MQ=255) gCAATGGGTGGTGCGCCAGCGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATACCCTGGAAGCGCGCCCCGGGGTGGAACGCCAcc > 4:117549/1‑149 (MQ=255) aaTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc < 6:219598/138‑1 (MQ=255) aaTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgc > 5:219598/1‑138 (MQ=255) gtggtgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCCTTTCAGCCCATAGCGTGCAAGCCCGCGAAGGGTTGGAACGCCACCCCaaaa > 8:261403/1‑148 (MQ=255) tgCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTCCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTc < 1:29867/149‑1 (MQ=255) ccAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGAAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTg < 2:22292/149‑1 (MQ=255) aatgcgggcggGCGTTTGCCATGGTGTTGCGCCAGGTTTTGGGCTTAGCTTACGGGGGTTGCGTTCCAGCCGTTCGCTTCCAAGCCCGCTAAGGGTTGAAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGAc < 3:323067/145‑1 (MQ=255) aaaGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGAcg < 7:261403/149‑1 (MQ=255) aGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg < 5:69301/96‑1 (MQ=21) aGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg > 6:69301/1‑96 (MQ=21) cAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGAGCCCATAGCGGGCCAGCCCGCGCACGGTTTGAAACCAACCCCCAAAA‑ACAAAATGGAATTTTTATATTTTGTGATAGGACGCCCCCCAGAAAAAGCGGcg > 8:63876/1‑148 (MQ=255) aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg > 5:268436/1‑149 (MQ=255) gtgtGGCGCCAGGTTTTGGGCTGAGCGTACGGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgtg < 2:213682/149‑1 (MQ=255) ttGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc < 3:63922/148‑1 (MQ=255) cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCCGCCCTTCGCGTGCAAGCCCGCTATGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc < 6:313398/149‑1 (MQ=255) cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc < 1:101499/149‑1 (MQ=255) ggggTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGAcagca < 7:63876/149‑1 (MQ=255) ttgcAACCCTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCAc > 2:201889/1‑149 (MQ=38) aacccTTCGCGGGCTTGCACGCGAAGGGTTGCAACGCAACCCCAAAAAGACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCt > 1:266753/3‑149 (MQ=38) cTTCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg > 3:274696/1‑149 (MQ=39) ttCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGATGGTGTCATACACCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg > 5:266162/1‑149 (MQ=255) ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt < 5:56024/123‑1 (MQ=38) ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt > 6:56024/1‑123 (MQ=38) ttCGCGGGCTTGCACGCGAAGGGCTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg > 8:276170/1‑149 (MQ=39) gCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt < 7:162090/92‑1 (MQ=255) gCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGgt > 8:162090/1‑92 (MQ=255) cgctAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgt < 3:117549/149‑1 (MQ=255) gcgAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtg > 3:184867/1‑149 (MQ=255) cgAAGGGCTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGgtgt < 4:274696/149‑1 (MQ=255) gTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAg > 2:108294/1‑69 (MQ=255) gTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAg < 1:108294/69‑1 (MQ=255) | TGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACGCAACCCCAAAAA‑ACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGCCGGCAGTGTCGGTGT > NC_002947/5740338‑5740625 |
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: ‑ |