Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,574,079 | 0 | C | A | 42.2% | 75.4 / 20.9 | 45 | G375V (GGT→GTT) | pqqL | putative periplasmic M16 family zinc metalloendopeptidase |
Reads supporting (aligned to +/- strand): ref base C (6/20); new base A (19/0); total (25/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.46e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGCGCAGTGCCGCCAGAAATAGTCTTCAAC > NC_000913/1573937‑1574228 | cAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTTGGTTTTCAATTATCGGGTTCTTCTACAGAAAAAACa > 2:232968/1‑145 (MQ=255) aGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAAATTATCGGGTTTTTTAGCAGAAAAAACa > 1:451694/1‑143 (MQ=255) aGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTAAATCATCGAGTTCTTTTACAGAAAAAACAAGCTg > 1:671367/1‑148 (MQ=255) aCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AAt < 2:697427/151‑1 (MQ=255) aCGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGACATCATCGCGGTTTTTCGCCGAAAAACCAAGCTGAt > 2:486298/1‑147 (MQ=255) ggCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATgg < 1:289495/151‑1 (MQ=255) gCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGt < 2:159835/151‑1 (MQ=255) gTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGGGATTTTGCCTCCACCAGTTCTTTACCACAAAAAACACGCTGAAC‑AAttggt > 2:186588/1‑146 (MQ=255) tattTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCa < 2:58179/150‑1 (MQ=255) tGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACCACATCGAGGTCTTTATCAGAAAAAACAATCTGATT‑AAAGGGTGCCAACACTGcc > 2:494035/1‑151 (MQ=255) tGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTGCATCATCGGGGTTTTCAGCAGAAAAAACACGCTGATC‑AATTGTTGCCCACACTGcc > 2:191065/1‑151 (MQ=255) tGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTGACATCATCGAGTTTTTTAGCAGAAAAAACAAACTGATCAAATGGTGGCC‑ACTCTTCCA‑tt > 1:930013/1‑149 (MQ=255) ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTaa < 1:578155/150‑1 (MQ=255) ggATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAAt < 2:245872/151‑1 (MQ=255) tGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGTCATCAACCGGATTTTTTGCAGAAAAAACACGCCGATC‑AATTGTTTGCCACACTGCCC‑TTTATgt > 1:597081/1‑149 (MQ=255) gCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCAt < 2:316295/151‑1 (MQ=255) ccAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTGTAAACATCGAGGTTTTCCGCAGAAAAAACACGCTGATC‑AATGGGTGCCAACTCTGCCC‑TTAATgtctt > 1:761059/1‑147 (MQ=255) cAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTATATTATCGAGGTATTCAGCCGAAAAAACAAGCCGATT‑AAAGGTTTCCAACTCTGCCC‑TTAATGAAtt > 2:346086/1‑150 (MQ=255) cAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATTGAGGTCTTTAGCAGAAAAAACAAACTGGTT‑AAAGGGTGCCAAATCTTCCC‑TTTAttgatt > 2:583265/1‑145 (MQ=255) gACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACAATATCGCGTTTTTTATCAGAAAAACCACGCTGATC‑AATGGTTGCCCACCCTGCCA‑TTTATGccatctcacc > 2:59140/1‑142 (MQ=255) aCTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTTACATCATCGCGTTCTTCAGCCGAAAAAACATGCTGATC‑AATGGTTGGCAACTCTGCCA‑TTTATGCATTTgcagc > 2:288463/1‑150 (MQ=255) aCTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCgcagca < 1:128961/151‑1 (MQ=255) cGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCAtattat < 1:327202/151‑1 (MQ=255) cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTg < 2:527958/151‑1 (MQ=255) cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTg < 1:73647/151‑1 (MQ=255) aGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCAttt < 1:968279/151‑1 (MQ=255) cTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGATATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCACCTCTGCCA‑TTAATGCATTCGCAGCTTCCTGCATATTATCGTCTCGTGCAttt < 1:224889/149‑1 (MQ=255) cTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTACATCATCGAGGTCTTTAGCAGAAAAAACATGCTGATC‑AATTGTTGCCCACCCTTCCC‑TTTATTCATTTCCAGCATCCTGCATATTTTCGTCTCGTGCCTTTTCACgg > 2:139755/1‑151 (MQ=255) gATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGa < 2:658954/150‑1 (MQ=255) aTCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTCTTCAGCAGAAAAAACAAGCTGATC‑AATGGGTGCCAACTCTTCCC‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGaa > 2:500568/1‑150 (MQ=255) cAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGaaaaa < 2:714105/151‑1 (MQ=255) aaCCGCATTTTTCAGCCAGGTGAGGCGGGGAGAATTTACATCATCCGGTTTTTTTGCAGAAAAAAAAAGCTGGTC‑AATGGGTGGCAAACCTGCCC‑TTTATTCAATCGCCGCCTCCTGCATAATTTTGTTTCCTGCATTTAAACGGgaaaa > 1:128215/1‑150 (MQ=255) aCCGCATTTTTCAGCCAGGTGAGGCGGGTTGATTTTACATCCACCAGATTTTTTGCACAAAAACCAAGCTGATT‑AATGGGTGCCAAATCTTTCATTTATTGCTTTC‑CAGCACCCCTCATATTTTCGTTTCCTGCCTTTTAACGGgaaaa > 1:260749/1‑149 (MQ=255) gCATTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCCTCGAGGTTTTTCTCAGAAAAAACATGCTGATC‑AATGGGTTCCAACTCTGCCC‑TTTATGCATTCCCAGCATCCTGCATATTTTCGTCTCGTGCATTTTCACGGgaaaa > 2:674586/1‑146 (MQ=255) cATTTTTCAGCCAGGTGAGGCGGGTAGATTTTTCATCATCGAGGTTTTCAGCAGAAAAAACATGCTGATC‑AATGGTTGCCAACCCTTCCC‑TTAATGCATTCGCCGCATCCTGCATATTATTGTCTCGTGCATTTTCACGGaaaaaaa > 2:505422/1‑145 (MQ=255) aTTTTTCAGCCAGGTGAGGCGGGGAGATTTTTCATCCTCGAGTTTTTTCACCGAAAAAACAAGCTGATC‑AAAGGTTGCCAACTCTGCCA‑TTTATGCATTTCCAACATCCTGCCTATTTTCGTCTCCTGCATTTTCACGGaaaaa > 1:125596/1‑144 (MQ=255) aTTTTTCAGCCAGGTGAGGCGGGGAGATTTTACATCATCGAGTTCTTTTGCAGAAAAAACAAACTGATC‑AATGGTTGTCCACTCTGCCA‑TTAATGCATTCTCAGCATCCTGCATATTTTCGTCTCGTGCATTTTCACGGGAAAACAGCGa > 2:102292/1‑150 (MQ=255) ccAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCg < 1:805548/151‑1 (MQ=255) ttGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgc > 1:479494/1‑151 (MQ=255) aCATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTgcgcg > 1:151055/1‑151 (MQ=255) tcAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGCGCAGTGCCGCCAg < 2:440093/148‑1 (MQ=255) agcagAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGCGCAGTGCCGCCAGAAATa < 1:148316/151‑1 (MQ=255) cagAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGCGCAGTGCCGCCAGAAATAGt < 2:671853/151‑1 (MQ=255) aaaaaCCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGCGCAGTGCCGCCAGAAATAGTCt < 2:105689/150‑1 (MQ=255) aCCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGCGCAGTGCCGCCAGAAATAGTCTTCAAc > 1:787336/1‑151 (MQ=255) | CAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGCGCGCAGTGCCGCCAGAAATAGTCTTCAAC > NC_000913/1573937‑1574228 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |