Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,134,407 | 0 | A | C | 51.7% | 30.6 / 44.9 | 29 | G229G (GGT→GGG) | asmA | outer membrane assembly protein AsmA |
Reads supporting (aligned to +/- strand): ref base A (12/2); new base C (0/15); total (12/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.75e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-07 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATC > NZ_CP009273/2134323‑2134494 | gTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTg > 1:550627/1‑90 (MQ=255) gTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTg > 2:842476/1‑90 (MQ=255) tGAAAGGCGTTTATGTGCCTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTg < 1:922223/90‑1 (MQ=255) tGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTg < 2:794141/90‑1 (MQ=255) aaGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCa > 2:421694/1‑90 (MQ=255) tcccACTGGGCCTGAAAACCCCCCTGCCCCTGAATCCCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAgcc < 2:902306/88‑1 (MQ=255) tGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAgcc < 2:219894/90‑1 (MQ=255) tGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGtt > 2:817194/1‑90 (MQ=255) gggCCTGAAACCCCCCCTCCCCCTAAATCCTTTTTTTTGGAAAATCGGCCCCCTGCAACTGCCATTTAATTTGTTAAATAGCCGCCGTTa < 1:245541/90‑1 (MQ=255) gAAACCCCCCTGCCCCGAAATGCTTTGTTTGGGAAAAACCCCCCCCTACAACTGAAATTAAATTTGTTCAATGGCCGCCGTTAAATCATg < 1:679550/90‑1 (MQ=255) cTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGa > 2:679550/1‑90 (MQ=255) ctgcccctgAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATc > 1:42774/1‑90 (MQ=255) gcccctgAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAg > 1:664757/1‑90 (MQ=255) caccctaAATCCTTTTTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTAAATAGCCGCCGTTAAATCATGCGGATAATCAg < 2:491980/88‑1 (MQ=255) ccctgAATCCCTTTTTTTGGAAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGaa < 2:459342/90‑1 (MQ=255) ccctaAATCCCTTTTTTTGGTAACCCGGCCCCCGGCAACTGCCAGTTAATTTGTCCAATAGCCGCCGTTAAATCATGCGGATAATCAGaa < 1:522354/90‑1 (MQ=255) cctgAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAg > 1:297692/1‑90 (MQ=255) aTACTTTTTTTTGGTAAATCGCCCCCCTCCAACTGCCAGTTATTTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCt < 1:667241/90‑1 (MQ=255) tttttttGGTAAATCGGCCCCCTGCAACTGCCAGTTATTTTTTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgt < 2:819644/87‑1 (MQ=255) tttttttGGAAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgt < 2:792192/87‑1 (MQ=255) ttGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtacctg > 1:491980/1‑88 (MQ=255) tGTTTTGGTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgta < 2:550627/90‑1 (MQ=255) ttttGGTAATTCGGCCCCCTGAACCCCCCAGTTAATTTTTTAAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCtaccgtacc < 2:284968/90‑1 (MQ=255) ggTAAATCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTa < 2:450560/90‑1 (MQ=255) ggTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTa > 1:697757/1‑90 (MQ=255) aTCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGa > 1:780906/1‑90 (MQ=255) tCGGCCCCCTGCAACTGCCATTTAATTTTTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGat < 1:61440/90‑1 (MQ=255) tCGGCCCCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGat < 2:38056/90‑1 (MQ=255) gCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATc > 2:351723/1‑90 (MQ=255) | GTTAAATGAAAGGCGTTTATGTGACTCCTGCCACTGGGCCTGGAAACTCCCCTGCCCCTGAATACCTTGTTTTGGTAAATCGGCACCCTGCAACTGCCAGTTAATTTGTTCAATAGCCGCCGTTAAATCATGCGGATAATCAGAAGCATCTACCGTACCGTTAAGGGATATC > NZ_CP009273/2134323‑2134494 |
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: ‑ |