Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,214,118 | 0 | C | A | 40.5% | 43.0 / 17.2 | 37 | A422A (GCC→GCA) | ygjG | putrescine aminotransferase |
Reads supporting (aligned to +/- strand): ref base C (18/4); new base A (0/15); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.01e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATAC > NZ_CP009273/3213993‑3214253 | gTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCa > 1:207007/1‑139 (MQ=255) gAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AAc > 1:681666/1‑139 (MQ=255) aGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCg > 1:54892/1‑139 (MQ=255) aGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCg > 1:59522/1‑139 (MQ=255) cgcgTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGc > 2:637506/1‑138 (MQ=255) cgcgTGGTAAAGGGATGTTGATGGCGATTGAGTTTGATGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGc > 1:88216/1‑138 (MQ=255) aggggTTGTGGATGGCGTTGGTGTTTGTTAAT‑ACGAAAACCGATTTAACCTTTC‑CCAGCGACTTCTCCCGCGTGCGCGTGGTGTCGGCGGAACCGCTAAAAAACACAAAACCGA‑CCCGCATAG‑AACCGCcactgacact < 2:757026/137‑1 (MQ=255) aggggATGTTGAGGGCGATTGATTTTGTTGAAAGCGAAATCGGATAAAA‑CTTTC‑CGAGAGAAATGTCCCCCCGGGGCGTGTGGGTGGCCGGAACGTTAAATAACGCAAAAGCAA‑TCCGCATTG‑AACCGCcactgacact < 1:366107/137‑1 (MQ=255) agacagtGTTGATGGCGTTTGATTTTTTTAAAAACAAAATCGGTTATAA‑CTTTG‑CCAGCGAGATTTCCCCCCAGCGCGTACTGGTGGCCGGAACGCTAAAAACCCCAAAAACGA‑CCCGCATTG‑AACCGCcactgacact < 1:54146/133‑1 (MQ=255) aaGGGTTGTTGGTGGCGTTTGATTTTGTTAAAAACAAAACCGTTTATAA‑TTTTG‑CCAGGGGGTTTTCCCGCCCGCGCTTACTGGGGGCCGGAACGCTCAAAAACCCAAAAACGA‑TCCGCATTG‑AACCGCcactgacact < 1:30322/139‑1 (MQ=255) aaGGGATGTTGATGGCTATTGATTTTTTTGATAACAAAATCGGTTATAT‑CTTTG‑CCAGCGAGTTTTTCCGCCAGCGCGTCCTGGTGCCCGAAACGCTCAAAAACCCAAAAACGC‑CCCGCTTTG‑AACCGCcactgacact < 1:515707/139‑1 (MQ=255) aaGGGATGTTGATGGCGATTGAGTTTGTTAATAACAAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAACCGA‑TCCGCATTG‑AACCGCcactgacact < 2:442461/139‑1 (MQ=255) ggatttttGAGGGCGTTTGAGTTTTTTGAAAACAAAATCGGTTATAA‑TTTTG‑CAAGGGAGATGTCCCGCCAGCGCGTACTGGTGGCGGAACCGCTCAATAACGCAAAACCGA‑TCCGCATTG‑AACCGCcactgacactga < 1:856746/133‑1 (MQ=255) ttGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGGCACTGACCAtt > 2:771049/1‑138 (MQ=255) ttGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCAtt > 1:260957/1‑138 (MQ=255) ttGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCAtt > 1:757026/1‑138 (MQ=255) ttGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCAtt > 1:485430/1‑138 (MQ=255) ttGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCAtt > 1:694279/1‑138 (MQ=255) gggggggattgattTTGTTGAAAACGAAAGGGGTTAT‑A‑CTTTGCCCAGGGGGTTTTCCCCCCGGCGCGTACGGGTGGCCGAACCGCCCAATACCCCAAAACCGT‑CCCTCATTA‑AACCGCCACTGACACTGACCATTGaa < 2:341910/133‑1 (MQ=255) aagggggATTGTGTTTGTTAAAAACGAAATCGGTTATAA‑TTTTG‑CCAGGGAGTTGTCCCCCCGGCGCGTCCTGGGGCCCGGAACGCAAAATACCGCAAAAACGA‑CCCGCATTG‑AACCGCCACTGACACTGACCATTGa < 2:694279/132‑1 (MQ=255) ggggtgATTGTTTTTGTTAATAACGAAACTGGTTAAAA‑CTTTC‑CCAGAGGGTTTTCCCGCCGGGGCTTCCGGGGGGCCGAAACCCTAAATAACCCAAAAACCC‑TCCGTTTTG‑AACCGCCACTGACACTGACCATTGaa < 2:601909/133‑1 (MQ=255) agggCGATTGATTTTTTTGATAACGAAACCGGCTATAA‑TTTTG‑CCAGCGAGATGTCCCGCCAGCGCGTACTGGTGGCCGGAACCCCAAATACCGCAAAAACGA‑CCCGTATTG‑AACCGCCACTGACACTGACCATTGaa < 1:446159/136‑1 (MQ=255) aTGGCGTTTGTGTTTTTTGATAACGAAATCGGTTATAC‑CTTTC‑CCAGGGAGATGTCCCCCCAGCGCGTACTGGTGCCCGAAACGCTCAATAACCCAAAACCGA‑CCCGCATTA‑AACCGCCACTGACACTGACCATTGaa < 1:611418/138‑1 (MQ=255) aTGGCGATTGATTTTGTTGAAAACGAAATCGGTTATAA‑CTTTG‑CCAGGGAGATGTCCCGCCGGCGCGTACGGGGGGCGGAAACGCTAAATAACGCAAAAACGA‑CCCGCATTA‑AACCGCCACTGACACTGACCATTGaa < 2:107689/138‑1 (MQ=255) aTTGAGTTTTTAGAAAACAAAAGCGATTATAT‑TTTTG‑CCAGGGAGGTGTCCCGCCGGCCCGTACTGGTGGCGGAAACGCAAAAAAACGCAAAAACGA‑CCCG‑ATTGAAACCGCCACTGACACTGACCATTGAACAGtgtg < 2:54892/139‑1 (MQ=255) aTTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAAcctgtct > 2:54146/1‑133 (MQ=255) gattttGTTAAAAACGAAATGGGCTATAT‑TTTCC‑CCAGCGAGATGTCCCGCCAGCGCTTACTGTGGGCGGAAACTCACAATACCCCCAAAACCA‑TCCACATTC‑AACCGCCACTGACACTGACCATTGAACAGTGTGaa < 2:485430/135‑1 (MQ=255) tG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTggcg > 2:338010/1‑139 (MQ=255) aGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCa > 1:657533/1‑139 (MQ=255) cGCCAGCGGTTCCGGGGGGCCGGACCGCTCAATAACGCAAAAACGA‑CCCGCATGG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAgtg < 2:46913/139‑1 (MQ=255) cgggTACTGGGGGCCGAAACGCTAAAAAACCACAAACCGC‑CCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCgaag < 2:327426/136‑1 (MQ=255) ggggCCGGAACGCTAAATAAC‑CCAAAACCACCCCGAATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGt < 1:713161/134‑1 (MQ=255) gCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAatat > 2:402653/1‑138 (MQ=255) gCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAatat > 1:42334/1‑138 (MQ=255) ccgtcaaaTAACGCAAAAACGA‑TCCGCATGA‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCg < 2:260957/133‑1 (MQ=255) aaTAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGc > 1:829346/1‑139 (MQ=255) cGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATac > 1:872982/1‑139 (MQ=255) | GTACAGGAAGCGCGTGGTAAAGGGATGTTGATGGCGATTGAGTTTGTTGATAACGAAATCGGCTATAA‑CTTTG‑CCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGA‑TCCGCATTG‑AACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATAC > NZ_CP009273/3213993‑3214253 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |