Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 5,300,655 | 0 | A | G | 38.5% | 13.0 / 11.3 | 30 | intergenic (‑85/+3488) | recC/PP_RS08750 | RecBCD enzyme subunit RecC/tRNA‑Val |
Reads supporting (aligned to +/- strand): ref base A (13/3); new base G (10/0); total (25/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.62e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAA > ACB143/5300523‑5300789 | gCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGATAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGAGGGAGc > 7:71341/1‑135 (MQ=255) gCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGc > 7:18963/1‑135 (MQ=255) aCGATCATGAAGCCCGGGTGTAGAGGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGCGGGCAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGt < 3:226359/134‑1 (MQ=255) tGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGGCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGG‑CCCTGCCCCCTcc > 5:326843/1‑136 (MQ=255) gggTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCTCCCGCAGCg > 8:252591/1‑135 (MQ=255) gTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGG‑GCCTGGCCCCTGCCGCCCCG‑At > 3:205248/1‑135 (MQ=255) ggTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGCGCGG‑GCGGGGCCGCGCCCGCCCCG‑ATTagcggccggc > 6:78862/1‑128 (MQ=255) ggTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCCCTCCCGCAGCG‑GTTCGCGCAggg > 1:285860/1‑132 (MQ=255) ggTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGGGCCCGCAACCCCAGCG‑ATTAGCGCAGGtg > 5:287209/1‑133 (MQ=255) tGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCACCCCCCGCG‑AAAAGCCCAGGTCCA‑Ggtt > 1:348269/1‑135 (MQ=255) aGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGG‑CCGTGCCCCCAACCGCAGCG‑CTTAGCGCAAGTCGA‑Ggtttt > 3:332251/1‑135 (MQ=255) aTCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCGCCCGCCGCG‑GATAGCGCAGGGCGA‑Gtttttttc > 1:317709/1‑133 (MQ=255) cGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCCG‑GCGTGGCCCCGCCCCCCCCG‑CTTAGCGCACGTCGAGGTTTTTTGccc > 6:2894/1‑136 (MQ=255) aaGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCcac > 2:235464/1‑134 (MQ=255) aaGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGAGGGAGCGG‑GCGTGCCCGCGGCCGCACCG‑ATTAGCGCAGGTGGA‑GTTTTTTTCCcaaa > 5:245579/1‑133 (MQ=255) ctctGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGG‑CCCGGCCCGCACCCGCCGCG‑CTTAGCGCACGGCGT‑GTTTTTTGGCCACAAAGAc > 5:72092/1‑135 (MQ=255) aGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGG‑CGGTGGCCCCGCCCCCCACG‑AATGGCGCAGGGGGG‑GTTTTTTTGCCCCAAAAGCaaaaaac > 6:145390/1‑133 (MQ=255) cTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCCG‑CCGTGCCCGCTCCCGCAGCG‑ATTAGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTa > 2:304279/1‑134 (MQ=255) cTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCCG‑CCATGCCCGCAACCGCCACG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAACCCAAAACCCCAAc > 2:304281/1‑135 (MQ=255) cTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCCGGGGCGG‑CCCGGGCCCCGCCCGCCCCGCATACG‑GCACGTCGACGTTTTTTG‑CCACAAAAACAAAAcaccaaa > 6:324636/1‑132 (MQ=255) tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGG‑CCGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGCCAAAACCCCAAcc > 8:279592/1‑135 (MQ=255) tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCCTGCCCGCACCCGCAGCG‑CTTCGCGCAGGTCGA‑GTTTTTTGCCCCCAAAGACAAAACCCCAAcc > 8:183580/1‑135 (MQ=255) gtgggCTGAGCGGGCAACGATACCACGCGTACATGAACCGGCCGAATGCCATCAGGCGGCGGGCGCGG‑CGGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCt < 2:118739/132‑1 (MQ=37) aTGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGGGGCCCGGGCCCCCCCCCCCCCCC‑ATAGGCGCCGGGGGG‑GTTTTTTGGCCCAAAAAAAAAAACCCCAAcaac > 5:26672/1‑133 (MQ=255) ccTGAGAGGGCAACTATATCCCGCGTACCTGAAAGGGCCTAAGCCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCTCCCGAAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCtgcgt < 2:275743/134‑1 (MQ=37) agGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGCGCGG‑CCGTGCCCGCTCCCGCCGCG‑CGTCGCGGCGGGGGG‑GGGTTTTTCCCACAAAGACAAAACCCCTACCggggtgcggcg > 2:342216/1‑132 (MQ=255) aCGGCCGAATGCCATCAGGCTGCGGGGGCGG‑CCCGGCCCCCGCCCCCCGCG‑CTTCGCGCCGGTGGA‑GTTTTTTTCCCACCAAGACAAAAACCCAAACAGGGGGCGCAGCTAGGGGGTTTTGGAAATGAATCttta > 2:230681/1‑133 (MQ=255) ccaaTCCCATCTGGCGCCGGGCGCGG‑CGGGCCCCCCCCCCGCCGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATg < 7:268551/133‑1 (MQ=17) tGCCATCAGGCGGGGGTAGCGG‑CGGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCt < 1:56540/135‑1 (MQ=255) gAGCGG‑CCGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTTCCCACAAAGACAAAACCCCTACCTGCGGGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGCCCTACTCGCACATGGGGaaa > 6:89102/1‑136 (MQ=255) | GCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAA > ACB143/5300523‑5300789 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |