Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,462,881 | 0 | C | . | 22.2% | 62.7 / 17.3 | 27 | intergenic (+235/‑131) | fadL/yfdF | long‑chain fatty acid outer membrane transporter/uncharacterized protein |
Reads supporting (aligned to +/- strand): ref base C (9/12); new base . (3/3); total (12/15) | |||||||||||
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.80e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCAT > NC_000913/2462746‑2463021 | cAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACGGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCt < 2:51057/149‑1 (MQ=255) aCAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTAtt < 5:155318/149‑1 (MQ=255) gCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCAt < 1:169538/148‑1 (MQ=255) ccTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTc > 1:209084/1‑113 (MQ=255) ccTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTc < 2:209084/113‑1 (MQ=255) cTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATATTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAAtt > 6:151854/1‑148 (MQ=255) aaaaCAATAAGCGGTATACACAGACAAAACTGAAGGGAAAGGCATATATTCAGTCGTTCTGCACGACCTTCTCAAAGAGTTTACTTTTCTGCATATCCAGGATACTCCACCCCCCAGCGTATTGTACGCTCAAACAATCAAATTAAAg > 6:151400/1‑148 (MQ=255) acaGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGc > 7:94754/1‑149 (MQ=255) aaaCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑CCCCCCCCTGGCTATTGTGCGCTc < 1:213146/105‑1 (MQ=255) aaaCCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑CCCCCCCCTGGCTATTGTGCGCTc > 2:213146/1‑105 (MQ=255) cTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGAt < 8:94754/149‑1 (MQ=255) tGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATtg > 2:48299/1‑94 (MQ=255) tGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATtg < 1:48299/94‑1 (MQ=255) aaaGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCGGGCTATTGTGCGCTCATCCACTCAAATTAAAGATAGGTTCTGAATAAATGCGCGTTTTTTGATAGTCTAt > 1:59850/1‑149 (MQ=255) cAGGCCTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATAATCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTaa < 3:98495/148‑1 (MQ=255) cTCGTCCTTCTCAAAGAGTTTACTTTTCTTCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATATGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGt < 5:151400/148‑1 (MQ=255) ctcAAAGAGTTTACTTTTCTGCATTTCCAGGATACT‑CCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATCGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAaatcaa < 5:151854/148‑1 (MQ=255) agagTTTTCTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGGTATTGTGCGCGCATACACTCAAATTAAAGATAGGTTCTAAATAAATGCGCGTTTTTTGATAGGCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCaaa < 2:59850/149‑1 (MQ=255) agagTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaa > 6:16988/1‑71 (MQ=255) agagTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaa < 5:16988/71‑1 (MQ=255) tttACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATCCACTAAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTGTATTAAATTAGGTAATATAAATTTGTAAAAAATCAATCAAAATg > 4:117141/1‑148 (MQ=255) cTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCGCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAAt > 1:127313/1‑148 (MQ=255) tGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAATATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATc < 3:117141/148‑1 (MQ=255) cAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaaa > 5:85211/1‑50 (MQ=255) cAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTaaa < 6:85211/50‑1 (MQ=255) cAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTAcc < 2:127313/148‑1 (MQ=255) aGGATACTCCCCCCCCCTGGCGATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATAATGCTACCAt > 1:86453/1‑149 (MQ=255) gATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCAtt > 1:46814/1‑92 (MQ=255) gATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCAtt < 2:46814/92‑1 (MQ=255) | CAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTTTCAGGCGTTCTGCTCGTCCTTCTCAAAGAGTTTACTTTTCTGCATTTCCAGGATACTCCCCCCCCCTGGCTATTGTGCGCTCATACACTCAAATTAAAGATAGGTTCTAAATAAATGAGCGTTTTTTGATAGTCTATTTCATTAGGTAATATATATTTGTAACAAATCAATCAAAATGGAATAAAATCATGCTACCAT > NC_000913/2462746‑2463021 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |