Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,134,407 | 0 | A | C | 55.2% | 24.8 / 37.3 | 29 | G229G (GGT→GGG) | asmA | outer membrane assembly protein AsmA |
Reads supporting (aligned to +/- strand): ref base A (10/3); new base C (0/16); total (10/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.43e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGAT > NZ_CP009273/2134321‑2134491 | tGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCAccc > 1:524587/1‑90 (MQ=255) gTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTGCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCACCCTg > 1:817705/1‑90 (MQ=255) gTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCACCCTg > 2:209379/1‑90 (MQ=255) ggCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTGCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGt < 2:817705/90‑1 (MQ=255) tcccACTGGGCCTGAAACCCCCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAgcc < 2:501714/88‑1 (MQ=255) tcccACGGGGCCTGGAACCCCCCCTGCCCCTAAATCCT‑TTTTTTTGGAAAATCGGCCCCCTGCACCTGCCAGTTAATTTGTTCAATAgcc < 1:858964/88‑1 (MQ=255) tGCCTCGGGGCCTGGAACCTCCCCTGCCCCTGAATTCC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgcc < 1:638/90‑1 (MQ=255) tGCCCCTGGGCTTGGAAACCCCCCTCCCCCTGAATCCC‑TTTTTTTGGTAAATCGGCCCCCTTCAACTGCCAGTTAATTTGTTAAATAgcc < 2:12962/90‑1 (MQ=255) tGCCCCTGGGCCTGGAAACCCCCCTGCCCCTGAATCCC‑TTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAgcc < 2:613442/90‑1 (MQ=255) tGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAgcc < 2:524587/90‑1 (MQ=255) cTGGGCCTGGAAACTCCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGt > 1:667304/1‑90 (MQ=255) gaaaaCCCCCCTCCCCCTAAATCCC‑TTGTTTTGGAAAATCGCCCCCCTGCAACTGCCATTTAATTTGTTCAATAGCCGCCGTTAAATCAt < 2:376320/88‑1 (MQ=255) cccccctccccctaAATCCT‑TTTTTTTGGTAAATCGGCCCCCGGCAACTGCCAGTTATTTTGTTAAATAGCCGCCGTTAAATCATGCGGa < 1:390541/88‑1 (MQ=255) cccccctccccctaAATCCC‑TTGTTTTGGTAAACCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAAAAGCCGCCGTTAAATCATGCGGa < 2:87164/88‑1 (MQ=255) ctgcccctgAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATc > 1:798958/1‑90 (MQ=255) ccccttAATTCCTTTTTTTTGGAAAATCGCCCCCCTGCACCTGCCATTTATTTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAg < 1:633523/90‑1 (MQ=255) cccctgAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGa > 2:372540/1‑90 (MQ=255) cctgAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAg > 2:615837/1‑90 (MQ=255) tgAATTCT‑TTTTTTTGAAAAATCGGCCCCCGGCAACTGCCAGTTATTTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCa < 1:436749/90‑1 (MQ=255) ttttttgGGAAAATCGCCCCCCTGCAACTCCCATTTAATTTGTAAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgt < 1:178255/87‑1 (MQ=255) ttGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgt > 1:865282/1‑90 (MQ=255) ttttttGGAAAACCGGCCCCTGGCAACTGCCATTTATTTTTTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgta < 1:303842/88‑1 (MQ=255) tGTTTTGGAAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgta < 2:798958/90‑1 (MQ=255) ttttGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtacc > 1:146328/1‑90 (MQ=255) gTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTaa > 1:702693/1‑90 (MQ=255) tAAATCGCCCCCCTGCAACGGCCATTTATTTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAg < 1:158732/90‑1 (MQ=255) aaaaCCGGCCCCCTGCAACTGCCAGTTATTTTTTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAg < 2:762654/89‑1 (MQ=255) gcGGCCCCCTGCAACTGCCAGTTAATTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGat < 2:667304/89‑1 (MQ=255) ccGGCCCCCTGCAACTGCCAGTTAATTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGat < 2:762628/89‑1 (MQ=255) | TGGTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACC‑TTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGAT > NZ_CP009273/2134321‑2134491 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |