Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB143 | 5,300,676 | 0 | A | C | 41.9% | 14.1 / 10.1 | 31 | intergenic (‑106/+3467) | recC/PP_RS08750 | RecBCD enzyme subunit RecC/tRNA‑Val |
Reads supporting (aligned to +/- strand): ref base A (15/3); new base C (12/1); total (27/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.21e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCG > ACB143/5300545‑5300808 | 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) gCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGTAATGAATCTTGACGATGACCTACTCTCACATGTGGAAACCCCACACTAcc > 8:145218/1‑134 (MQ=255) cGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc > 6:117282/1‑135 (MQ=255) cGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc > 2:1956/1‑135 (MQ=255) cGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGc > 6:350409/1‑135 (MQ=255) cGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACATACTCTCACATGGGGAAACCCCACACTACCATCgg > 2:245964/1‑134 (MQ=255) gaaGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCg < 4:332251/133‑1 (MQ=255) | GGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGG‑CCGTGCCCGCTCCCGCAGCG‑ATTCGCGCAGGTCGA‑GTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCG > ACB143/5300545‑5300808 |
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: ‑ |