Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,241,786 | 0 | C | T | 20.0% | 61.0 / 10.1 | 30 | intergenic (‑42/+41) | ygjV/uxaA | Imp‑YgjV family inner membrane protein/altronate hydrolase |
Reads supporting (aligned to +/- strand): ref base C (11/13); new base T (3/3); total (14/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATG > NC_000913/3241656‑3241917 | aaGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGtttt < 4:103955/149‑1 (MQ=255) gCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGtttgct > 6:152633/1‑146 (MQ=255) cGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGtttgcttttc < 5:152633/149‑8 (MQ=38) ttCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTaac < 8:2774/149‑1 (MQ=255) ttGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTc > 1:76210/1‑149 (MQ=255) aaaaaTGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCTGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGa < 8:132489/148‑1 (MQ=255) ttgtGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTtata < 2:76210/148‑1 (MQ=255) aTACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGt < 2:23379/113‑1 (MQ=39) aTACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGt > 1:23379/1‑113 (MQ=39) cAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTAcgcc < 7:95912/149‑1 (MQ=255) tCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCt > 6:158162/1‑149 (MQ=255) aGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCtt < 5:158162/148‑1 (MQ=255) cccACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCtt > 8:148649/1‑133 (MQ=255) ccaATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGAt < 3:131734/124‑1 (MQ=255) ccaATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGAt > 4:131734/1‑124 (MQ=255) aaTACGCGCTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTa < 1:57619/59‑1 (MQ=255) aaTACGCGCTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTa > 2:57619/1‑59 (MQ=255) tACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGt > 6:72961/1‑149 (MQ=255) aGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACttt > 8:136717/1‑73 (MQ=255) aGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACttt < 7:136717/73‑1 (MQ=255) aTGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTc < 4:89375/118‑1 (MQ=255) aaaaGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAAAAAGACATTACGACTTTAAAGCGTTACGCCGCTTTTGAAAAACGCCAGTAGGCGGAAGTCGTTACGATCGTTAAAGGTTTGCTTACCGTTGGCAAAAGCAACGa > 6:162912/1‑149 (MQ=255) cGGCGTTTCATAATGTGAAGCGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTg > 6:133923/1‑73 (MQ=25) cGGCGTTTCATAATGTGAAGCGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTg < 5:133923/73‑1 (MQ=25) ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTa < 5:31055/120‑1 (MQ=255) ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTa > 6:31055/1‑120 (MQ=255) ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATg > 4:32466/1‑141 (MQ=255) ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATg < 3:32466/141‑1 (MQ=255) tCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGt < 7:90649/107‑1 (MQ=255) tCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGt > 8:90649/1‑107 (MQ=255) | AAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATG > NC_000913/3241656‑3241917 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |