Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,740,467 | 0 | T | G | 26.8% | 39.0 / 24.5 | 41 | intergenic (+32/+39) | hutF/PP_5037 | formiminoglutamate deiminase/lipocalin family lipoprotein |
Reads supporting (aligned to +/- strand): ref base T (13/17); new base G (7/4); total (20/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.61e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGC > NC_002947/5740320‑5740610 | gCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGt < 5:79058/148‑1 (MQ=255) aGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCaa < 3:164539/149‑1 (MQ=255) tgATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAgg < 7:219188/148‑1 (MQ=255) gATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCATAGCGTGCAAGCCCGCGTACGGt > 1:213523/1‑149 (MQ=255) aTGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCACCCCATCGCGTGCAAGCCCGCGAAGGGtt > 7:122582/1‑149 (MQ=255) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa‑c > 2:223257/1‑149 (MQ=255) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGcttgc > 6:5453/1‑128 (MQ=255) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGcttgc < 5:5453/132‑5 (MQ=255) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg > 6:266742/1‑145 (MQ=38) gAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTg < 5:266742/145‑1 (MQ=38) aaTGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGACAGCGGGCGGGCGTTTGCCAAGGGGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTgcaa‑cg < 7:361673/149‑1 (MQ=255) aaTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaa‑c < 1:59364/141‑1 (MQ=38) aaTGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAACCCTTCGCGGGCTTGCACGCGAAGGGCTgcaa‑c > 2:59364/1‑141 (MQ=38) gcATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAACCGCAACCCCAAAAAACTCAGTGCATTTTTa < 1:223257/149‑1 (MQ=255) cAGGAAGATGACCGCGGGCGGGCTTTTGCCAAGGGGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGc < 2:213523/149‑1 (MQ=255) aaGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGAt < 5:110433/149‑1 (MQ=255) gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt < 6:289084/98‑1 (MQ=255) gAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGtt > 5:289084/1‑98 (MQ=255) cgggcgggCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTTGCGCCCATAGCGGGGAAGCCCCCGCCGGGTTGGCC‑CGCAACCCCCAAAAACGCAGGGGGTTTTTTAGATTTTCTGAGGGGGcgcccg > 2:37127/1‑149 (MQ=255) ggcgggctTTTGACAAGGTGTTGCGCAAGATTTTGGGCTGAGCGTACCGGGTTTGCTTTGAAGCCCTTCGCGTTCAATCCCGCTAAGGGTTGAAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGAcgcccgcc < 5:205070/149‑1 (MQ=255) aGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTGGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTgg < 1:37127/149‑1 (MQ=255) gtTGCGCCAGCTTTGGGGCGGAGCGTACCGGGGTTCCGTTGCACCCCTTCGCGTGCAAGCCCGCTAAGGGTTGAAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc < 7:153156/149‑1 (MQ=255) cgcCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAc < 5:356833/149‑1 (MQ=255) ttttGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAccc < 7:184191/144‑1 (MQ=255) gggCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCTAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGcc < 5:121478/148‑1 (MQ=255) gCTGAGCGTACCGGGGTTGCGTTTCAGCCCATAGCGTTCAAGCCCGCCAACGGTTTCAA‑CGCCACCCCAAAAAAAAAAGTGCAATTTTAATATTTTCTGATCGGACCCCCGCCACATAAAACCGCTGGTGTAAAAACCCAGCTGCcccg > 2:55647/1‑149 (MQ=255) ggtttcAGACCTTCGCGGGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtct < 6:377684/147‑1 (MQ=255) cgttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACACCTGCTGCCGCGCCCTGGACAGCAGGTtct > 1:83767/1‑148 (MQ=255) gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAcccc < 4:149154/119‑1 (MQ=255) gttgcAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATAcccc > 3:149154/1‑119 (MQ=255) cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt < 4:134323/129‑1 (MQ=255) cAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCt > 3:134323/1‑129 (MQ=255) aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTc > 2:314774/1‑107 (MQ=255) aGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTc > 5:52206/1‑149 (MQ=255) gCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTcc < 6:52206/149‑1 (MQ=255) cTTCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg > 5:58087/1‑149 (MQ=39) cTTCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAg > 6:183469/1‑149 (MQ=39) ttCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg > 4:311551/1‑149 (MQ=255) ttCGCGGGCTTGCACGCGAAGGGCTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAgg > 5:203087/1‑149 (MQ=39) tCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGc > 7:42311/1‑149 (MQ=255) gcgTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc < 5:22626/133‑1 (MQ=255) gcgTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCAAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTtc > 6:22626/1‑133 (MQ=255) | GCCAGGTGCGCGATGTGATGGTGAATGGGCAATGGGTGGTGCGCCAGGGGCGGCATGCCCAGGAAGATGAAAGCGGGCGGGCGTTTGCCAAGGTGTTGCGCCAGCTTTTGGGCTGAGCGTACCGGGGTTGCGTTGCAGCCCTTCGCGTGCAAGCCCGCGAAGGGTTGCAA‑CGCAACCCCAAAAAACTCAGTGCATTTTTACTATTTGCTGATCGGACGCCCGCCAGATCAAGCGGCTGGTGTCATACCCCTGCTGCCGCGCCCTGGACAGCAGGTTCTCACGCTCCCAGGC > NC_002947/5740320‑5740610 |
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: ‑ |