Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,574,051 | 0 | G | T | 42.2% | 62.2 / 16.5 | 45 | V384V (GTC→GTA) | pqqL | putative periplasmic M16 family zinc metalloendopeptidase |
Reads supporting (aligned to +/- strand): ref base G (9/17); new base T (19/0); total (28/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.00e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGC > NC_000913/1573906‑1574199 | gTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatc < 1:11308/151‑1 (MQ=255) tttCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAAAATTGTca > 1:259016/1‑148 (MQ=255) tCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGAATTTAc > 1:105649/1‑145 (MQ=255) gAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGAcatca > 1:658742/1‑147 (MQ=255) ctctTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGGTTTTTcatccaccggtttttttg > 1:439935/1‑132 (MQ=255) cGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGGAGATTTGAAATCATCGAGGTCTTTAGCAGaaaaaac > 2:360573/1‑148 (MQ=255) 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 < 1:73647/151‑1 (MQ=255) cGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTg < 2:527958/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) | GTTTCGAAAGTTGATATGTCTCTTCCGGCGACAAGAACGGCGTATTATTTAATGAGCTGGATGCCAGGCGACTGGTCAGCATACGTAAATCACGTTCAGCTTGCTGATCAACCGCATTTTTCAGCCAGGTGAGGCGGGTAGATTTGACATCATCGAGTTCTTCAGCAGAAAAACCATGCTGATC‑AATGGTTGCCAACTCTGCCA‑TTAATGCATTCGCAGCATCCTGCATATTATCGTCTCGTGCATTTACACGGAAAAACAGCGACTGATAATCGGGTGCAATTTTAACGCTGC > NC_000913/1573906‑1574199 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |