| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_pae_tpiA | 1,446,357 | 1 | . | T | 25.0% | 31.7 / 14.5 | 16 | intergenic (+151/+21) | ydbL/feaR | DUF1318 family protein/transcriptional activator for tynA and feaB |
| Reads supporting (aligned to +/- strand): ref base . (7/5); new base T (2/2); total (9/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCGCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGG > NC_000913_3_pae_tpiA/1446215‑1446505 | ccGCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATg > 1:98168/1‑150 (MQ=255) tCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagtt > 1:102092/1‑126 (MQ=255) tCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagtt < 2:102092/126‑1 (MQ=255) tATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCt > 3:55933/1‑150 (MQ=255) aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt > 3:58039/1‑150 (MQ=255) tCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCTACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt < 4:63227/150‑1 (MQ=255) tttGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt < 2:98168/150‑1 (MQ=255) ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaa > 1:31626/1‑150 (MQ=255) tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGaaaa > 1:142065/1‑150 (MQ=255) cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATg > 3:83456/1‑150 (MQ=255) atatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGAt < 2:31626/150‑1 (MQ=255) aaCTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGtt > 4:47994/1‑150 (MQ=255) ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt < 1:80559/80‑1 (MQ=37) ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt > 2:80559/1‑80 (MQ=37) ttGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGCTACTCGCCTGGCATAACCCCGAACGGTTGCTTTAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTg < 3:8735/150‑1 (MQ=255) cA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCgg < 4:58039/150‑1 (MQ=255) | CCGCAATATCGGGATAACGCTACAGGAGAAGGCGATTGATCTATTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCAGAAAATCCCCAATGAAAGCCGATGCCTGCCAGTTTTTCATCATCTGCGG > NC_000913_3_pae_tpiA/1446215‑1446505 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |