Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,586,057 | 1 | . | C | 32.3% | 29.4 / 42.8 | 31 | intergenic (+140/+75) | PP_4061/PP_4063 | hypothetical protein/long‑chain fatty acid‑‑CoA ligase |
Reads supporting (aligned to +/- strand): ref base . (10/11); new base C (4/6); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGG > NC_002947/4585926‑4586202 | cACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑t < 4:178084/149‑1 (MQ=39) ggCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑Tgcg > 8:114541/1‑149 (MQ=39) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTg < 1:362081/149‑1 (MQ=39) aaCAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTg < 5:1821/149‑1 (MQ=39) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCggg < 2:365210/147‑1 (MQ=38) cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGt < 7:14794/148‑1 (MQ=38) ctCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAcc < 7:114541/148‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGt < 8:253881/132‑1 (MQ=38) tCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGt > 7:253881/1‑132 (MQ=38) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGc > 1:27156/1‑118 (MQ=37) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGc < 2:27156/118‑1 (MQ=37) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc > 6:103520/1‑139 (MQ=38) aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTAc < 5:103520/139‑1 (MQ=38) gTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCgg < 1:26232/129‑1 (MQ=38) gTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCgg > 2:26232/1‑129 (MQ=38) gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAATTGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCcgc < 1:252142/131‑1 (MQ=37) tcgcgggTAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACg > 6:221355/6‑149 (MQ=255) tAAACCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCg > 1:199429/1‑148 (MQ=255) aCCCGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAAc > 6:93196/1‑149 (MQ=255) aCCCGCTCCTACCAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcg < 8:132759‑M2/148‑77 (MQ=255) cccGCTCCTACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGACCt > 3:251065/1‑149 (MQ=255) gCTCCTACCAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcg < 3:55649‑M2/144‑77 (MQ=255) cTACCAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcag > 4:49121‑M2/1‑64 (MQ=255) cTACCAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcag < 3:49121‑M2/148‑85 (MQ=255) tACCAGGGACC‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCAcc > 6:361488/1‑149 (MQ=255) tACCAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagct < 2:311515‑M2/149‑87 (MQ=255) aCCAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatc > 3:375667‑M2/1‑62 (MQ=255) aCCAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatc < 4:375667‑M2/115‑54 (MQ=255) ccAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaatggagcagcgcagatcagct > 7:132759‑M2/1‑61 (MQ=255) ccAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaa < 2:357388‑M2/127‑67 (MQ=255) ccAGGGA‑CCGTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAagaggccggcccgatcaatgcaagatcttcaccttgccaagcacccggtcatcctcaacctttcaa > 1:357388‑M2/1‑61 (MQ=255) gggacc‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGaa < 6:201906/129‑1 (MQ=255) gggacc‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGaa > 5:201906/6‑134 (MQ=255) gacc‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCcgg > 2:100080/4‑149 (MQ=255) | CACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGC‑GGTAAACCCGCTCCTACCAGGGA‑C‑GTGCCGGCCTTGAA‑TGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGG > NC_002947/4585926‑4586202 |
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: ‑ |