Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,998,178 | 0 | G | A | 42.9% | 46.9 / 23.4 | 42 | A66T (GCG→ACG) | uvrD | DNA‑dependent ATPase I and helicase II |
Reads supporting (aligned to +/- strand): ref base G (6/18); new base A (18/0); total (24/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.01e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA > NC_000913/3998030‑3998328 | aGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAgcg < 2:864295/151‑1 (MQ=255) ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAgcggcggc < 1:474322/150‑1 (MQ=255) ccGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGga < 1:323792/151‑1 (MQ=255) cgcgCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGGTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGGGTAAGGTACCCACAAAAcggcggggggg > 1:164052/1‑149 (MQ=255) cGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGGAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGTCGGTTTCCAACAAAacagcgggcggg > 1:609422/1‑139 (MQ=255) agcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATc < 2:950267/151‑1 (MQ=255) cTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGGGGAGGTTTA‑CAAAAAAAAGGCGGCGGGGGGGCG‑GCCTCGGA‑TCgggg > 2:951025/1‑150 (MQ=255) ggTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCggg < 2:894697/148‑1 (MQ=255) gCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACt < 2:295901/150‑1 (MQ=255) ggcgggcgCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGAAGGTTTCCCACAAAACGGCGGGGGGGGTGCG‑TTCTCGTA‑TCCGGGAAatgatg > 2:146995/1‑151 (MQ=255) gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTCCAACAAAACGGCGGCGGGGATGGG‑GCATCGTA‑TCGGGCAACTGATggg > 2:548709/1‑150 (MQ=255) gggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAACGGCGGCGGGGATGGG‑GCCTCGTA‑TCGGGCAACAGAAggg > 2:624164/1‑150 (MQ=255) gcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGa < 1:945603/151‑1 (MQ=255) gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGc < 2:20862/151‑1 (MQ=255) tAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGAAGGTTTCCCACAAAACGGCGGGGGGGATGGG‑GCATCGTA‑TCGGGGAACAGATGGGGACGAGCCCgggggg > 2:293253/1‑150 (MQ=255) aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGt < 2:949435/150‑1 (MQ=255) aaGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGGTTACCAAAAAAAAAGCGGCGGGGAATCGCGCTTCTAA‑TGGGCCAAC‑GATGGGGACGAGCCAggg > 2:966092/1‑146 (MQ=255) aCGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGGTTTTCAACAAAACGGGGGGGGGGGAGGG‑GCATCGTATTGGGGGAACT‑ATGGGGACGAGGCCggggggggtg > 1:635818/1‑146 (MQ=255) cgTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTc < 2:777745/151‑1 (MQ=255) cATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTcc < 1:416662/151‑1 (MQ=255) aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACggg < 1:868320/150‑1 (MQ=255) ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCACCTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc < 1:50026/151‑1 (MQ=255) ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgc < 2:86030/151‑1 (MQ=255) cTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTCCCACAAAACGGCGGCGGGGGAGCG‑TCCTCGTA‑TCGGGCCACTGATGGGGACGGGCCCGGGGGGGGTGGGGGTCGGCACCTTCCCCCGGCCGgcg > 2:828444/1‑151 (MQ=255) ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTCCAACAAAACGGCGGCGGAGGTGCG‑TCATCGTA‑TCCGGCCACCCATGGGGGCGAGCCAAGGGGGGATGGTGGTCGGGCCCCTCCCCCGGGTGGCGcccc > 2:469126/1‑151 (MQ=255) tgatgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGGTTATCCACAAAAAGGCGGCGGGGGGTGG‑TCCTTTTA‑TCGGGGGACCGAAGGGCACGAGCCCGGGGGGGTTGTGGGGCTGGACCCTcccccgggcg > 2:739716/1‑135 (MQ=255) gaGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGttt < 1:441477/150‑1 (MQ=255) gCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTTACGCTTTCCAACCAAAAGGCGGCGGAGCGGCG‑CCCCCGTA‑ACGGGGGACTGAAGGGGGCGAGGCCGGGGGGGGTGTTGGGCGGCACCCTCCCCGGGCCGGGGCAcc > 1:406234/1‑144 (MQ=255) gTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTTTCAACAAAACGGCGGCGGGGATGCG‑TCATCGTA‑TCGGGGAACCGAAGGGGACGAGGCAGGGGGGGATGTGGGGCGGCACCCTCCACCGGCCGGCGCACCGgtt > 2:500296/1‑146 (MQ=255) gAAAACTGCTCGCCATACTCGATTATGGC‑GGTGCCGGTTATCAACAAAACGGCGGCGGGGGAGGG‑GTATCGTA‑ATGGGGAACTGATGGGGACGAGGCCGGGGGGGGTGGGGGGGGGGGCCTTCCACCGGCCGgcggccccgttggt > 1:666113/1‑135 (MQ=255) tGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAACGGCGGCGGGGATGCG‑GCATCGTA‑TCGGGCCACTGATGGGGACGAGCCCGGGCGGGATGTGGGGCGGCACCCTCCCCGGGGTGGCGCACCGTTTTCTGCGTGCGCAcc > 2:21435/1‑151 (MQ=255) cTCGCCATACTCGATTATGGC‑GGTGACGTTTTTCAAAAAAACCGCCGCGGGGGAGCG‑TCATCGTA‑TCTGGCGACTGAATGGGGCGGGCCAGGGGGGGGTGTGGGGTGGGGCCCTCCCCCGGCCGGGGCCCCCGTTGGTGCGGGCGcccc > 1:81483/1‑149 (MQ=255) cGCCATACTCGATTATGGC‑GGTGGCGTTTTACAACAAAACGGGGGCGGGGGAGGG‑TCATCGTA‑TCGGGGAAAAGATGGGGACGAGCCCGGGGGGGATGTGGGGGGGGACCCTCCCCCGGCCGGCGCCCCGgtt > 1:493495/1‑133 (MQ=255) cGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACAata < 2:9001/150‑1 (MQ=255) aTACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAc < 2:772396/151‑1 (MQ=255) aTACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAc < 2:275455/151‑1 (MQ=255) aTACTCGATTATGGC‑GGTGAAGTTTTTCCAAAAAAAGGCGGGGGGGAAGCG‑TCCTCGTA‑TCGGGCAACCGATGGGGACGAGCCCGGGGGGGATGTGGGGCGGGGCCCTCCACCGGCCGGGGCACCCCTTTCtgcgggggcccc > 1:47312/1‑135 (MQ=255) aCTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGc < 1:641625/151‑1 (MQ=255) aTGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTg > 2:250347/1‑150 (MQ=255) tGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGCCTCGCGCACCGTTTTCTGCGTGCGCGCCATATGGACGCCAATCTGc > 2:24931/1‑150 (MQ=255) c‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGc > 2:534165/1‑150 (MQ=255) gTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTc < 2:847344/151‑1 (MQ=255) ccAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTcc > 2:931969/1‑147 (MQ=255) caacaaAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTcc > 2:662790/1‑146 (MQ=255) aacaaAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc > 1:777708/1‑151 (MQ=255) gcggcggcgGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACa < 1:406395/151‑1 (MQ=255) | AGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGC‑GGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCG‑TCATCGTA‑TCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACA > NC_000913/3998030‑3998328 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |