Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,462,878 | 2 | . | C | 40.0% | ‑2.5 / 14.0 | 10 | intergenic (+232/‑134) | fadL/yfdF | long‑chain fatty acid outer membrane channel/bacteriophage T2 receptor/protein YfdF |
Reads supporting (aligned to +/- strand): ref base . (2/4); new base C (1/3); total (3/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 = 1.18e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑CTCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATAT > NC_000913/2462749‑2462975 | aCAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑ccc > 2:428329/1‑131 (MQ=255) gTTGGCAGCTAGCCCTTTCAAAACAATAAACGGTATCCCCAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCCTCCTTCTCAAAAAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑cccccccccc < 2:231831/133‑1 (MQ=255) cAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑Cccccccccccctgg < 1:620082/133‑5 (MQ=255) gCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑cccccccc > 2:623322/1‑124 (MQ=255) gCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑cccccccc < 1:623322/124‑1 (MQ=255) gCTTGCCCTTTCAAAACAAAAAGGGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCCCGTCCTTCTCAAAGAGTTTTCTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑Ccccccccccccctggc < 2:437524/133‑7 (MQ=255) ccTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GATT‑‑C‑‑‑CCCCCCCCCCC‑‑TGGCTATtgtg < 1:412138/133‑1 (MQ=255) ttttAAATCAATAAGCGGTTTCCCCAGGAAAACCTTAAGGGAAAGGGATTTTTTCAGGGG‑TTTCGCCCGTCCCTCTCAAAGGGTTTTCTTTTTTTCCTTTCCAGG‑GA‑TCCC‑‑‑CCCCCCCCCCC‑‑TGGCTATtgtg < 2:69958/133‑1 (MQ=255) tATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑ACC‑CCCCCCCCCCC‑‑TGGCTTTGGTGCGCACTTACAGTCAAATTa > 1:253216/1‑133 (MQ=255) tATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑CCCCCCCCCCCCGTGGCTATTGTGGCCTCATACCTTCAAATTa > 1:218733/1‑133 (MQ=255) ccACAGGCAAACCTGAAGGGAAAGGCATATTTTCGGGGG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCG‑G‑GA‑T‑‑G‑‑‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGAt < 2:588249/133‑1 (MQ=255) gACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑ACC‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAgg > 2:431304/1‑133 (MQ=255) cAAACCCGAAGGGAAACGGATATTTTGCAGCG‑TTTTTCTCGTCCTTCCCAAAGAGGTTTCTTTTATGGCTTTTCA‑G‑GG‑T‑‑ACC‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGtt < 2:396488/133‑1 (MQ=255) aaaGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑AC‑‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGc > 2:467078/1‑133 (MQ=255) aaaGGAATTTTTTCATGCGTTTCTG‑TCCATCTGCTCAAAAAGTTTTCTTTTCTGCATTTCCA‑G‑GA‑TTCC‑‑‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAg < 2:253216/133‑1 (MQ=255) ttcTCGTCCTTCTCAAAAAGTTTTTTTTTCTGCATTTCCA‑GAGT‑TC‑C‑‑‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCa < 2:245039/131‑1 (MQ=255) tGCTCGTCCTTCTCAAAGAATTTTTTTTTTTTCATTTCCA‑G‑AA‑T‑‑ACC‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCa < 2:218733/133‑1 (MQ=255) cTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑ACC‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCAtt < 2:355792/133‑1 (MQ=255) tCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑ACC‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTa < 1:428329/133‑1 (MQ=255) gTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑ACC‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg < 1:431304/133‑1 (MQ=255) tCCTTCTCCAAGAGTCTTCTTTCCTGCATTACCA‑G‑GA‑T‑‑ACCCCCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAgg < 1:159894/133‑1 (MQ=255) ccaaaaGATTTTTCTTTTTTTCATTTCCA‑G‑AA‑TCCC‑‑‑CCCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAatat < 2:149211/130‑1 (MQ=255) | ACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCG‑TTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCA‑G‑GA‑T‑‑A‑‑‑CTCCCCCCCCC‑‑TGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATAT > NC_000913/2462749‑2462975 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |