Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,357,755 | 0 | A | T | 53.6% | 12.6 / 41.4 | 28 | P530P (CCA→CCT) | yhcD | fimbrial biogenesis outer membrane usher protein |
Reads supporting (aligned to +/- strand): ref base A (10/3); new base T (0/15); total (10/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.18e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGA > NZ_CP009273/3357667‑3357844 | cGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCAt < 2:768451/90‑1 (MQ=255) aCTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATc > 1:308000/1‑88 (MQ=255) aCTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATc < 2:308000/88‑1 (MQ=255) tGGCAGCATCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAac < 1:92429/90‑1 (MQ=255) tGGCACCATCAACTTAGTGCTTAAGAAAAAAATGAAGAGGTTGGTGATTTTTTCTCCGTAAACATTAGTATTCCTTTTGGAAGATAAAac < 1:569132/90‑1 (MQ=255) tGGCACCACCAACTTAGTGCTTAAAAAAATAATAAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCTTTTTGGCAGAAAAAac < 1:734979/90‑1 (MQ=255) ctCATCTTAGAGCTTCAAAAAAAAAAGAAGACGATGGTGATCTTTTCTCCGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCca < 1:411287/89‑1 (MQ=255) aTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCca < 1:103677/90‑1 (MQ=255) tCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCcac < 1:758988/90‑1 (MQ=255) tCAACTTAGTGCTTCAAAAAATAATAAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCcac < 1:454419/90‑1 (MQ=255) aaCTTAGTGCTTAAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCATTATTCTTTTTGGCAGATCAAACACCGCCacac < 1:376182/90‑1 (MQ=255) cTTAGTGCTTAAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTg < 1:98686/90‑1 (MQ=255) gTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGtata > 2:367390/1‑90 (MQ=255) tCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGc > 1:650580/1‑90 (MQ=255) tCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGc > 2:420665/1‑90 (MQ=255) aaaatGAAGACGAGGGTGATCTTTTCTACGTCAACATCAGTTTTCCTTTTGGAAAACAAAACCCCGCCACACTGTATACTCGCCATGACg < 1:757585/88‑1 (MQ=255) taatGAAGACGATGGTGATCTGTTCTACGTAAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACGa < 2:470664/90‑1 (MQ=255) atGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATc > 2:374104/1‑90 (MQ=255) gACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATaaa < 2:48884/90‑1 (MQ=255) aCGATGGTGTTCTTTTCTACTCAAACATAATTTTTCATTTTGGCAAAAAAAACCCCGCCACACTGTATACTCGCCATGACGATCATaaaa < 1:374586/90‑1 (MQ=255) ttGGTGCTCTGTTCTCCGTCAACCTCAGTTTCCCTTTTGGAACAACAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAAccc < 1:804377/89‑1 (MQ=255) aTGGTGATTTTTTCTACTTAAACATCAGTATTCCTTTTGGAAAATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAccc < 1:348313/90‑1 (MQ=255) gttCTGTTCTACTTAAACATAAGTATTCTTTTTGGAAAATCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTAt < 1:382727/88‑1 (MQ=255) gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc > 2:180536/1‑90 (MQ=255) gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc > 2:50861/1‑90 (MQ=255) tCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTg > 2:22300/1‑90 (MQ=255) cAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgt > 2:605199/1‑90 (MQ=255) gTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgtcg > 2:776246/1‑90 (MQ=255) aTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCaga > 1:363615/1‑90 (MQ=255) | CGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGA > NZ_CP009273/3357667‑3357844 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |