Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,151,710 | 0 | C | A | 42.9% | 0.7 / 10.9 | 14 | intergenic (+64/+1) | yegJ/yegK | uncharacterized protein/ser/thr phosphatase‑related protein |
Reads supporting (aligned to +/- strand): ref base C (5/3); new base A (6/0); total (11/3) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.09e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTCTTCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCG‑CCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTG‑AAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACC > NC_000913/2151581‑2151914 | gctCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTTCTTTCTTCCA‑CCCCCCTGAGAAAATTACTCTGTCCACAgccc > 1:1069749/3‑146 (MQ=255) gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTCTTCAA‑CCGCCCTGAGAAACTCACTCTGCCCACCGCG‑CCAGTGCCAACGGTTTat > 1:895157/1‑165 (MQ=255) ggTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGA‑TTGTTTTTTTCAA‑ACCGCAAAAGAAAATCACTCTGCCCCCCACGCCCAGGCCCAAC‑TTTTATTATCCTCCGCACGCCCGGTTACCCCCCGGCCTGGCAAAAACTCCGTTAACAAATACGGTAATAAACTGAGTTGTTCCTGCGtt > 2:778232/1‑237 (MQ=255) cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGA‑TTGTTGTGTTCAA‑CCCCCCTGTGAAACTCACCCTGCCCACAGAG‑GCCGTGTCCACGTTTTCTCTTATTACCTACGCTCGCTTTCCCCCCGGCGGGACAAAAAATGCTTTTTCAATTCATGTCATAAATCGAGATGTTTCTGCGGTGCTgccc > 2:478053/1‑237 (MQ=255) aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAA‑TTGTTTTTTTTCA‑CCCCCATGAGAAAATCACTCTGCCCCCAACG‑CCACTGCCCACCTTTTTTTTTCATCCCTAAGCCCCGTTTCCGCCCGGCGGGGAAAAAACTCCGTTTAAAATTTAAGGAATAAAAAGAGTGTTTTCTGCGTtgctgtccccc > 2:681623/1‑233 (MQ=255) acgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑ATGTTGTTTTCAA‑ACGCCATGAGAAAATCACTCCTCCCACAGCG‑CCAGGGCCAACGGTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCCGGACAAAAAATGCTTTAACAATTCCGGTAATAAATCGAGGTGTTCCTGCGGTGGTGCCGCCAGTCCATTGCAAAAACggg > 2:792753/1‑239 (MQ=255) aaaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGA‑ATTTTTTATTCAA‑CCCCCCTGAGAAACTCACTCTCCCCACAACG‑CCAGTGCCCACGTTTTATCATCATCCGTACGCCCGTTTACCGCCCGGCTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCCTTG‑AAAAACGGGGgaaagaaaagta > 2:926455/1‑234 (MQ=255) attgttttccTAAACCGCCTTTTTTTTTTTTTTTTAACCGGGGGGTTTAATTTTTTTTCTCCAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCG‑CCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTG‑AAAAACGGGGTAAAGAAAGGCACAt < 1:681623/230‑1 (MQ=255) acccgCCGTTATTTTTTTTTTTAAATCGGGGGGTTGAT‑TTTTTGGCTTAAC‑CCGCCATGAGAAACTCACTCTGCCCACAGCG‑CCAGTGCCa < 1:1396431/87‑1 (MQ=255) ccGTTATATTTATTGTTTAAACTGGCGGTTTAA‑TTGTTGTCTTCAA‑CCGCCCTGAGAAACTCACTCTGCCCACAGCG‑CCCGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATaaa > 1:1054483/1‑160 (MQ=255) atatTTATTGTTTAAACTGGCGGTTTTA‑ATGTTTTTTTCCA‑CCCCCACGAGAAACTCACTCTTCCCACACCG‑CCAGTGCCAACGttttttt > 2:1396431/1‑88 (MQ=255) ttAAACTGGGGGGTGAA‑TTGGTTTTTTAAACCCGCCATGAGAAACTCACTCTGCCCACAGCG‑CCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGc < 1:796292/110‑1 (MQ=255) tAAACTGGCGGTTGAA‑TTGTTGTCTTCAA‑CCCCCCTGAGAAAATCCCTCTGCCCACAGCG‑CCAGTGCCAAAGTTTTTTCATCATCCGTACGCCCGTTTACCGCCGGGCCGGACAAAAACTGCTTTaa > 1:728152/1‑127 (MQ=255) cGGTTGAA‑TTGTTGTCTTCAA‑CCGCCATGAGAAAAACACCCTGCCCACAACG‑CCAGGGCCCAAGTTTTATTATCATCCGGACGACCGTTTACCGCCGGGGCGGGCAAAAAATGCCTTTACAATTCAAGGTATAAAACGAGGTGTTTCAGCGTTGCTGGCGCCAGTCCCTTT‑GAAAACAGGGGAAAAAAAGGCACATAAGGAGAATTATCCACCACATTTAAAGCCAGTCGCCGGATaac > 2:884262/1‑237 (MQ=255) actcactcTGCCCACAGCG‑CCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTG‑AAAAACgggg < 1:16092/148‑1 (MQ=255) actcactcTGCCCACAGCG‑CCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTG‑AAAAACgggg > 2:16092/1‑148 (MQ=255) | GATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTCTTCAA‑CCGCCATGAGAAACTCACTCTGCCCACAGCG‑CCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTG‑AAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACC > NC_000913/2151581‑2151914 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |