| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 287,172 | 0 | G | C | 44.4% | 61.7 / 18.5 | 45 | P128P (CCG→CCC) | yagH | CP4‑6 prophage; putative xylosidase/arabinosidase |
| Reads supporting (aligned to +/- strand): ref base G (5/20); new base C (20/0); total (25/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.69e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 | |||||||||||
GGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACC > NC_000913/287032‑287322 | ggTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTGACCCCCCCCTGt > 1:564603/1‑148 (MQ=255) ggCTGCTCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCa‑cgacgac < 1:513102/149‑1 (MQ=255) gctgctCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTTACCCCTCCCTGTTcc > 2:554507/1‑140 (MQ=255) ctgctCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTTACCCCTCCCCGTTCCA‑CGACGACGATg > 2:419875/1‑151 (MQ=255) tgctCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTTgccccgcccccttcccccac > 1:695866/1‑133 (MQ=255) tctACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTTGCCCCG‑CCCCGTTccccccccaccaaaggc > 1:294931/1‑135 (MQ=255) tACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTTACCCCCCCCTGTTCCA‑CGACCACGATgggc > 2:526469/1‑146 (MQ=255) tACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGcc < 2:377744/148‑1 (MQ=255) acacCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTTACCCCTCCCTGGTCCA‑CCACGACGATGGcc > 2:104921/1‑147 (MQ=255) cacCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTTACCCCTCCCCGTTCCA‑CGACGACGATgg > 2:126978/1‑144 (MQ=255) ccGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATa < 2:682859/151‑1 (MQ=255) ccGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTGCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAAt < 1:416759/150‑1 (MQ=255) gTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCCTCCCCGTTCCA‑CCACGACGATGGCCGCAAATACta > 2:237121/1‑149 (MQ=255) tGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTGGAACCTCCCCGGTCCATCGA‑GACCATGGCGGCAAAAACgata > 2:623035/1‑150 (MQ=255) tGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTGACCCCTCCCCGTTCCC‑CGACCACGAGGGCCGCCAATACTaaattttcccc > 2:447760/1‑144 (MQ=255) cGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCg < 2:529779/151‑1 (MQ=255) gTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTGACCCGCCCCCGTTCCC‑CCACCACGAAGGCCCCACAAACTATATCTACCCCCCCCggggg > 1:321995/1‑148 (MQ=255) gAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCgcgc < 2:690740/151‑1 (MQ=255) aaaaaCGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTGACCCCTCCCCGTTCCC‑CGCCCACGCAGGCCCCCAATACCATATCTACCCCCCCCGGGGGCCCCGcc > 1:320817/1‑151 (MQ=255) aaCGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGGCCCCTCCCCTTTCCA‑ACACGACCACGGCCGCACATACTCTTTCTAACGCCCGTTGGGGCCGCGccccc > 2:294939/1‑151 (MQ=255) cGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTTACCCCCCCCCGTTCCA‑CGCCCGCGATGGCCCCCAAAACTATATCTACCCCCCCCGGGGGCCGCGccacc > 1:491171/1‑149 (MQ=255) gCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTTACCCCTCCCCGTTCCC‑CGACCACGAAGGGCGCCAAAACTATATCTACCCCCCGCGGGGGGCGCGccacccc > 2:161195/1‑147 (MQ=255) cGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTGACCCCTCCCTGTTCCA‑CGACGACGATGGCCGCCAAAACTATATCTACCGCCCGTGGGGGCCGCGCCACcaca > 2:17913/1‑148 (MQ=255) gCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCaa < 1:343221/151‑1 (MQ=255) ccTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGca < 2:606787/150‑1 (MQ=255) tCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGCGGGGTTTACCCCGCCCCGGTTCC‑‑CGCCCACGACGGCCGCAAAAACCATATCTACCCCCCCCGGGGGGCGCcccccccccaccacccccc > 1:458650/1‑135 (MQ=255) ccGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCCGTTCCA‑CCACCACGATGGCCGCCAAAACTATATTTACCCCCCCTGGGGGCCCGGCCACCACACCAACCCCCCCAAAAcc > 2:554215/1‑149 (MQ=255) gcCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCgtg < 2:683628/151‑1 (MQ=255) cccTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTTGCCCCCCCCTCTTCCC‑CCACGACGCTGGGCGCCAAAACTATATTTACCCCCCCGGGGGGGCGCGCCCCCCCACCAACCCCCCCAAAACCCTCgtg > 1:363128/1‑150 (MQ=255) ccATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGGTTTACCCCGCCCCGTTCCA‑CGACGACGATGGCCCCAAAAACCCTATCTACCCCCCCTGGGGGGCGCGCCCCCACAGCAACCCCCCCCACACCCCCCTGGTACAAGCGTTTTAcc > 1:9120/1‑151 (MQ=255) gagcgagcCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAg < 1:599948/151‑1 (MQ=255) cgagcCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGAc < 1:186148/150‑1 (MQ=255) cgagcCAATCCCGATGGGCAACGGC‑GGGTTTGACCCCTCCCCGTTTCC‑CGACCACGATGGGCGCAAAAACTATATCTACCCCCCCCGGGGGCCCCGCCCCCACAGCAACCCGCACAACACCATCGCGTTACAGGCGTTTGACCCCCa > 1:176474/1‑147 (MQ=255) gagcCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCg < 2:400262/151‑1 (MQ=255) cgcCAATCCCGATGGGCAACGGC‑GGGGTTTACCCCTCCCCGTTCCC‑CGACGAAGGTGGCCGCCAAAACTATATCTACCGCCCCCGGGGGGCGGGCCCCCCCAACAACCCCCCCACCACCCTCCTGTTACAGGCGTTTGACCCGCCGAcc > 1:238633/2‑149 (MQ=255) cAATCCCGATGGGCAACGGC‑GGGTTTTACCCCTCCCCGTTCCA‑CCACGACGACGGCCGCCAAAACTCTATATACCGCCCGTGGGGGCCGCGCCACCCCAGCAAACCCCCCAACACCATCGGGTTACAGGCGGTTTACCCGCAGACCGGCAc > 2:53394/1‑151 (MQ=255) aaTCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACg < 2:634036/151‑1 (MQ=255) cccGATGGGCAACGGC‑GGGTTTGTCCCCTCCCCGTTCCA‑CGACCACGATGGCCGCAAATACTATATCTACCCCCCCTGGGGGGCGCCCCCCCACAGCAACCCCCACAACACCCTCGTGTTACAGGCGTTTTACCCCCCGACCGGCACGccc > 1:689647/1‑149 (MQ=255) gggCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGccc < 2:677828/151‑1 (MQ=255) gCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGa < 1:237322/151‑1 (MQ=255) cAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGa < 2:517813/150‑1 (MQ=255) aaCGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAgc < 1:447723/151‑1 (MQ=255) ggc‑ggGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTCCCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCa < 2:678054/151‑1 (MQ=255) gggTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCaaa > 2:164745/1‑150 (MQ=255) gACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGt < 1:526431/151‑1 (MQ=255) cGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTAc < 2:569583/151‑1 (MQ=255) gTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTAcc < 1:629768/151‑1 (MQ=255) | GGTAAATTCTGGCTGCTCTACACCGACGTGAAGATTGTCGACTCGCCGTGGAAAAACGGCCGCAACTTCCTCGTCACCGCGCCCTCCATCGAGGGGCCATGGAGCGAGCCAATCCCGATGGGCAACGGC‑GGGTTTGACCCGTCCCTGTTCCA‑CGACGACGATGGCCGCAAATACTATATCTACCGCCCGTGGGGGCCGCGCCACCACAGCAACCCGCACAACACCATCGTGTTACAGGCGTTTGACCCGCAGACCGGCACGCTCTCGCCCGAGCGCAAAACGCTGTTTACC > NC_000913/287032‑287322 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |