Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 294,671 | 0 | A | C | 27.5% | 123.5 / 10.5 | 51 | Y659D (TAT→GAT) | paoC | aldehyde oxidoreductase molybdenum‑binding subunit PaoC |
Reads supporting (aligned to +/- strand): ref base A (37/0); new base C (0/14); total (37/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.74e-13 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTC > NZ_CP009273/294582‑294756 | tttGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATa > 2:278978/1‑90 (MQ=255) tCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCa > 1:441725/1‑90 (MQ=255) ggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGAtt > 2:41826/1‑90 (MQ=255) gaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTa > 2:165980/1‑90 (MQ=255) tggggggTGGTTATCACCCAGAAAACCCCCCTCTTTTTTGGGT‑ATATCCGAAGGACCCGCCCCCCAATCCCCCGCCATATCGTGATTAAc < 2:457578/87‑1 (MQ=255) tGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGa > 2:177985/1‑90 (MQ=255) tGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGa > 2:177986/1‑90 (MQ=255) gggTCGGTACCACCCAGAAAACCCCCCTCTTTTTTTGGG‑ATACCCAAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAg < 2:383075/90‑1 (MQ=255) ggTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGt > 2:88526/1‑90 (MQ=255) tCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAg > 1:116414/1‑90 (MQ=255) tCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAg > 1:6520/1‑90 (MQ=255) tCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAg > 1:493671/1‑90 (MQ=255) tCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAg > 1:447511/1‑90 (MQ=255) tCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAg > 1:441885/1‑90 (MQ=255) tCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAg > 2:488880/1‑90 (MQ=255) cGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 1:292479/1‑90 (MQ=255) cGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 1:489414/1‑90 (MQ=255) cGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 2:200637/1‑90 (MQ=255) cGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGc > 2:227975/1‑90 (MQ=255) atcaccCAGAAAAACCCCCCCTTTTTTTGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa < 1:359999/90‑1 (MQ=255) cacccAGGAAAACCCCCTCTTGTTTTGGG‑ATATCCGAATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc < 1:200637/87‑1 (MQ=255) cacccAGAAAAACCCCCTCCTGTTTGGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc < 1:165980/87‑1 (MQ=255) atcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg > 2:158848/1‑90 (MQ=255) atcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACg > 2:465529/1‑90 (MQ=255) cAGAAAAATCCCCTCCTTTTTGGGG‑ATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc < 1:397724/90‑1 (MQ=255) ggAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt > 1:75910/1‑90 (MQ=255) ggAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt > 1:85580/1‑90 (MQ=255) ggAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt > 2:100858/1‑90 (MQ=255) ggAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt > 2:242879/1‑90 (MQ=255) ggAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAt > 2:369835/1‑90 (MQ=255) aaaaCCCCCCCCTTTTTGGGG‑ATACCCGCAGGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc < 1:127107/90‑1 (MQ=255) aTCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAcc > 2:190895/1‑90 (MQ=255) aTCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAcc > 1:286017/1‑90 (MQ=255) cctcttTTTTGGGG‑ATACCCGAATGACCCGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa < 1:488880/90‑1 (MQ=255) ctcctTTTTTGGG‑ATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAg < 1:141323/90‑1 (MQ=255) tcttttttGGGAATACCCGCAAGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGc < 2:181329/86‑1 (MQ=255) tttttgGGG‑ATACCCGCAAGACCCGCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc < 1:494894/88‑1 (MQ=255) ttttGGGAATACCCGCATGACCCGCCCCCCCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcct < 1:430763/90‑1 (MQ=255) tGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCAt > 2:154413/1‑90 (MQ=255) tGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCAt > 2:57215/1‑90 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 2:285347/1‑90 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 2:405833/1‑90 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 1:46075/1‑90 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 1:469636/1‑90 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 2:155367/1‑90 (MQ=255) atCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc > 2:131152/1‑90 (MQ=255) cATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGCCATCCACCGCCAGCTCCTCCATCAGTGCCgcgc > 1:223711/1‑90 (MQ=255) taacccGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc < 1:61400/86‑1 (MQ=255) ggCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct < 1:142726/64‑1 (MQ=255) ggCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct > 2:142726/1‑64 (MQ=255) cTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCAGCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt > 2:170349/1‑90 (MQ=255) tCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc > 2:454806/1‑90 (MQ=255) | TTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTC > NZ_CP009273/294582‑294756 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |