Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,501,865 | 0 | A | C | 57.1% | 37.0 / 56.3 | 42 | L123V (TTG→GTG) | fecD | Fe(3+) dicitrate ABC transporter permease subunit FecD |
Reads supporting (aligned to +/- strand): ref base A (15/3); new base C (0/24); total (15/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.27e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-09 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAG > NZ_CP009273/4501782‑4501945 | ggCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGc > 2:325421/1‑90 (MQ=255) gCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATgcc > 2:38704/1‑90 (MQ=255) gCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATgcc > 2:627685/1‑90 (MQ=255) gCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCACCTTCAGCAATATCAACCCCCCCATGccgccc > 1:772360/1‑90 (MQ=255) aaaGCCCCCCGCGGGGGAGCGCCAGCTTCACCGGCTGGGGGGTCTTCGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCaa < 1:910085/90‑1 (MQ=255) aGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAg > 2:461237/1‑90 (MQ=255) aGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAg > 2:375913/1‑90 (MQ=255) cACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCcag > 1:497919/1‑90 (MQ=255) cgcgGGTGAGCCCCAGCTTCATCGGCGGGGGGGTTTTTGCCAGCATCTTCAGAAATATCACCCCCGCCATGCCGCCCGCAAAGGCcagca < 1:202170/90‑1 (MQ=255) gTGAGCCCCAGCTTCCCGGGCTGGGGGGTCTTTGCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgc < 1:471457/90‑1 (MQ=255) gAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCCCCATGCCGCCCGCAAAGGCCAGCAGCGgcag > 1:1210624/1‑90 (MQ=255) aGCGCCAGCTTCATCGGCGGGGGGGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgcagc < 2:595851/90‑1 (MQ=255) gcgcCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGgcagca > 1:121095/1‑90 (MQ=255) cgcCTGTTTCCTCGGGTGGGGGTTTTTTCCAAGCTTCTCAAGCAATACAACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCAc < 1:1378408/90‑1 (MQ=255) cTTCATCGGCTGGGGGGTCTTTGCCAGCATCTCCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCAc < 1:490280/90‑1 (MQ=255) ccggcggGTGGGTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCa < 1:929119/89‑1 (MQ=255) tggtggGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc > 2:928114/1‑90 (MQ=255) tggtggGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc > 2:737886/1‑90 (MQ=255) tggtggGGCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc > 1:1006450/1‑90 (MQ=255) tgggggGTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGAc < 1:807729/90‑1 (MQ=255) gggggTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACgg < 2:724971/88‑1 (MQ=255) gggTTTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCa < 2:759755/90‑1 (MQ=255) gTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 2:1247989/90‑1 (MQ=255) gTCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 1:386752/90‑1 (MQ=255) gTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATa < 1:38704/90‑1 (MQ=255) tCTTTGCCAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaa < 2:182724/90‑1 (MQ=255) cTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaag > 1:570318/1‑90 (MQ=255) ttGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATaagaa > 1:1224924/1‑90 (MQ=255) tGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTg > 1:615971/1‑77 (MQ=255) tGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTg < 2:615971/77‑1 (MQ=255) gCCAGCATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGt < 1:927193/90‑1 (MQ=255) ccAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTa < 1:223119/90‑1 (MQ=255) ccAGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTa < 1:220620/90‑1 (MQ=255) ggCATCTGAACCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 1:508550/89‑1 (MQ=255) aGCATTTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 2:1115035/90‑1 (MQ=255) aGCATCTTCAGCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTaga < 1:12099/90‑1 (MQ=255) cATCTTCAGCAATACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGc < 1:506526/90‑1 (MQ=255) gaaaTACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 2:972107/88‑1 (MQ=255) gaaaTACCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 1:1175958/88‑1 (MQ=255) gCAATATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 2:389777/90‑1 (MQ=255) gCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 1:737886/90‑1 (MQ=255) aaaaTATCACCCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACag < 1:279355/88‑1 (MQ=255) | GGCCCAGCATGCAGAAAGCGCCACGCCGGTGAGCGCCAGCTTCATCGGCTGGTGGGTCTTTGCCAGCATCTTCAGCAATATCAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAG > NZ_CP009273/4501782‑4501945 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |