Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,921,519 | 0 | A | T | 57.7% | 2.1 / 27.7 | 26 | intergenic (+488/‑69) | ppnN/sdaC | nucleotide 5'‑monophosphate nucleosidase PpnN/HAAAP family serine/threonine permease SdaC |
Reads supporting (aligned to +/- strand): ref base A (5/6); new base T (15/0); total (20/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.01e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAAAGCTGAATTATTTGCATTCCTCCAGGAGAAATAGATGGAAACGACTC > NZ_CP009273/2921436‑2921600 | cATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTTCGGTCa > 1:167671/1‑90 (MQ=255) cATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTTCGGTCa > 1:196648/1‑90 (MQ=255) aCTTGTGTGTGTCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAATGCATCCCCGGAAAACAAAATTCATTTTTTCGGTTCGGCa > 1:386768/1‑90 (MQ=255) tgtgtgTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTTATTTTTACGGTCAGGCACCt > 1:16856/1‑90 (MQ=255) tgtgtgTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTTCGGTCAGGCACCt > 1:274211/1‑90 (MQ=255) tgtgtgTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTACGGTCAGGGACCt > 1:21313/1‑90 (MQ=255) tgtgtgTATCTTTCCGCCCTCAAATTATTACGGCGGGAAAAGATTAAACCAACGCCCGTAGACAGATTTTATTTTTTCGGGCAGGGAccc > 2:321738/1‑89 (MQ=255) tgtgtgTATCTTTCCGCCCTCAAATTATTACGGCGGGAAAAAATTAAACCATCGCCGATAAACAAATTTTATTTTTTCCGGCAGGCACCt > 2:249749/1‑90 (MQ=255) gtgtATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTTCGGGCAGGCACCTTcc > 1:116022/1‑90 (MQ=255) aTCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTTCGGTCAGGCACCTTCCCggg > 1:169767/1‑90 (MQ=255) aTCTTTCCGCCCTCAAATTATTACGGCGGTAAAAGATTAAGCCATCGCCGATAGACAGATTTTTTTTTTTCGGTCAGGCACCTTCCCggg > 2:375942/1‑90 (MQ=255) tCTTTCCGCCCTCAAATTATTACGGCGGAAAAAGGTTAAGCCATTGCCGATCGAAAGACTTCATTGTTTAGGGCACAGACCCTCCCGGGc > 1:55463/1‑90 (MQ=255) ttCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTTATTTTTTCGGTCAGGGACCTTCCCCGGCTGa > 2:125438/1‑90 (MQ=255) ttCCGCCCTCAAATTATTACGGCGGTAAATGATTAAACCATCCCCGATAGACAGATTTTTTTTTTTTGGTCAGGGACCCTCCCcgggcgg > 2:108999/1‑86 (MQ=255) cGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACt < 1:321738/90‑1 (MQ=255) tCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTa < 2:89718/90‑1 (MQ=255) taCGGCGGTAAATGATTAAGCCATCGCCGGTAGACAGATTTTATTTTTTCGGGCAGGCACCTTCCCGGGCTGAACTGGCTAAAAACTGaa > 2:110125/1‑90 (MQ=255) cggcggTAAAAGATTAAACCATCGCCGATAAACAGATTTTATTTTTTTGGTCAGGCACCTTCCCCGGGTGAACTGGCTAAAAACTGAatt > 2:323322/1‑90 (MQ=255) gcggTAAATGATTAAGCCATCGCCGATAGACAGATTTTTTTTTTACGGGCAGGGACCCTCCCGGGCTGAACTGGCTAAAAACTGAattat > 2:43120/1‑90 (MQ=255) ggTAAATGATTAAGCCATCGCCGGTAGACAGATTTTATTTTTTCGGGCAGGCACCTTTCCGGGCTGAACTGGCTAAAAGCTGAATTAttt > 2:381048/1‑90 (MQ=255) aaTGATTAAGCCATCGCCGATAGACAGATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAATGCTGAATTATTTGCAt < 2:386768/90‑1 (MQ=255) tAAGCCATCACCGATAGACAGATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAAAGCTGAATTATTTGCATtcctcc > 2:189564/1‑90 (MQ=255) aaGCCATCGCCGATAGACAGATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAAAGCTGAATTATTTGCATTCCTCCa < 1:108999/90‑1 (MQ=255) gATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAAAGCTGAATTATTTGCATTCCTCCAGGAGAAATAGATGGAAACg > 1:490115/1‑90 (MQ=255) tCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAAAGCTGAATTATTTGCATTCCTCCAGGAGAAATAGATGGAAACGACTc < 2:179479/90‑1 (MQ=255) tCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAAAGCTGAATTATTTGCATTCCTCCAGGAGAAATAGATGGAAACGACTc < 1:125438/90‑1 (MQ=255) | CATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTAAGCCATCGCCGATAGACAGATTTCATTTTTACGGTCAGGCACCTTCCCGGGCTGAACTGGCTAAAAGCTGAATTATTTGCATTCCTCCAGGAGAAATAGATGGAAACGACTC > NZ_CP009273/2921436‑2921600 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |