Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,340,089 | 0 | G | C | 20.9% | 71.0 / 13.8 | 44 | intergenic (+74/‑257) | ribAB‑II/PP_3814 | 3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base G (14/20); new base C (2/7); total (16/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.46e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCAC > NC_002947/4339942‑4340234 | aaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGa < 8:107728/149‑1 (MQ=255) aTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGAc < 6:154374/149‑1 (MQ=255) aCGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACg > 2:25030/1‑148 (MQ=255) ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGTCGGTGGTGTCACCCTCGCATTcgc < 5:125053/149‑1 (MQ=25) ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGTCGGTGGTGTCACCCTCGCATTcgc < 5:161494/149‑1 (MQ=25) gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGTCGGTGGTGTCACCCTCGCATTcgcg < 3:344844/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCGGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGTCGGTGGTGTGACCCTCGCATTCGCgg < 8:26315/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGTCGGTGGTGTCACCCTCGCATTCGCgg < 5:363192/149‑1 (MQ=25) cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCgg < 1:25030/149‑1 (MQ=255) ggCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAAc < 4:178571/149‑1 (MQ=255) tgcgCTGGCGGGGGTGCCGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGCACACACACGCAATGGATTTAGCGGGGCCCGGTAGGAGC‑GGGGTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCAGTCGCGGG‑TAAACCCGc < 3:243792/145‑1 (MQ=255) tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCt < 5:368126/149‑1 (MQ=255) gAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTTCCCGCGAAT‑GCGAGGGTGACACCCCCGCCGCAATCGCCGG‑GTAACCCGcgc > 2:254189/1‑147 (MQ=255) cacagCATCCCGTTCGCCCCGGGTTTGGAGGGCGCATAACACACAAGCATTTTATTTCGCGGGGCCGTGTCGCCGG‑GGGTTACCCCGCTAAT‑GCGCGGGTCACACCACCGACGCATTCGCGGT‑TAAACCCGCTCCTACAGGGCCC‑‑TTg < 2:87656/145‑1 (MQ=255) cAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTcgc > 6:12643/1‑120 (MQ=255) cAGAGCATCCCGGTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTcgc < 5:12643/120‑1 (MQ=255) ccGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTc < 7:312/116‑1 (MQ=255) gggATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCGTGTAGCTGC‑GGGTTAACCCGCGATT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTg < 1:254189/149‑1 (MQ=38) gggATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTTCCCGCGAAT‑GCGAGGGTGACACCACCGCCGCCTTACCGGG‑TAAACCCGCGCCGAAACGGCCC‑‑GTGGGGCGTTTTACAGGGCGGAGTg > 5:114277/1‑149 (MQ=37) acacaaAACAAATTGATTCCGCGGGGCCCGGTCGCCGC‑GGGTAACCCCGCTAAG‑GCGCGGGTCCCCCCCCCCCCCCATGCGGGGG‑AAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCAt < 1:246986/149‑1 (MQ=255) acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGGTTTAACCCCGAAT‑GGGGGGGGGAACCAACCCCCCCATAGCCGGG‑TAAACCCCAACCAAAAGGGCCG‑‑TTGGGGGTTTTGAAAGGGGGGGGTGTTGGCCGGGGGCCCCAAAc > 2:117404/1‑148 (MQ=255) aGCAATTGATTTAGCGGGGCCCTTTAGGCGC‑GGGTTTACCCGCGAATCGCTCGCTTCACACCACCGACCCAATCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACctct < 1:117404/149‑1 (MQ=255) ccAATTGATTTAGCGGGCCCCTGTAGGCGC‑GGGTTTCCCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTg < 7:1384/148‑1 (MQ=255) cgCGGGGCCCTGTAGGAGC‑GGGTTAACCCGATATT‑GGGTGGGTCACACCCCCGACGCATTCGCGGT‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCtt < 6:114277/147‑1 (MQ=38) gCGGGGCCCTGTAGGAGC‑GGGTTAACCCGCGAAT‑GCGTGGGTTACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg < 3:295018/148‑1 (MQ=255) cccTGTCGGCGG‑GGGTTTACCCGCTAAT‑GCGAGGGTTACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAcc < 5:315704/149‑1 (MQ=39) c‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata > 5:292782/1‑149 (MQ=255) c‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGata > 7:87268/1‑149 (MQ=255) ccGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 7:73813/1‑149 (MQ=255) ccGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTACTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa > 3:195367/1‑149 (MQ=255) cgcgAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCCGATTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCATGGAATACCATAATATGATATCCCGTAg > 1:104012/1‑149 (MQ=255) gcgAAT‑GCGTCGGTGGTGTCACCCTCGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc > 1:315613/1‑149 (MQ=25) gcgAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGa > 5:139392/1‑136 (MQ=255) gcgAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGa < 6:139392/136‑1 (MQ=255) aT‑GCGCGGGTTACACCCCCGCCGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGTATACCATAATATGATATCCCGTAGACCTTg < 1:217065/148‑1 (MQ=255) cggcgcGGGCGACACCCCCGACCCCTTCGCGGG‑TAACCCCGCTCCTACACGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGGCGCCCCATACCTCTTCCCAAACTCTTGGAATACCATAATATGATATCCCGTAGACCTTGa < 8:73813/147‑1 (MQ=255) aGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCCG‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCt > 1:272955/1‑149 (MQ=255) aGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTACTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTg < 4:195367/149‑1 (MQ=255) ggggGACACCACC‑CCGGATTTGGGGGTTAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGc < 6:221111/149‑1 (MQ=255) gACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGc > 1:254833/1‑149 (MQ=255) acacCACCCACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGcc < 1:9069/149‑1 (MQ=255) accaccGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCgg > 4:124882/1‑149 (MQ=255) caccGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCa > 2:261839/1‑149 (MQ=255) accGACGCTTTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCAc < 6:292782/149‑1 (MQ=255) | AATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGC‑GGGTTTACCCGCGAAT‑GCGAGGGTGACACCACCGACGCATTCGCGGG‑TAAACCCGCTCCTACAGGGCCC‑‑TTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGTCTGCAGGCCGGCAC > NC_002947/4339942‑4340234 |
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: ‑ |