| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,510,395 | 0 | C | A | 42.9% | 49.3 / 27.5 | 42 | S131R (AGC→AGA) | ydcR | putative DNA‑binding transcriptional regulator and putative aminotransferase |
| Reads supporting (aligned to +/- strand): ref base C (19/5); new base A (0/18); total (19/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.06e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA > NC_000913/1510250‑1510542 | gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑c < 2:898536/151‑1 (MQ=255)gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑c > 1:652401/1‑151 (MQ=255) agaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑cc > 1:819767/1‑151 (MQ=255) gaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑cc < 1:900816/150‑1 (MQ=255) gaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCg < 2:968781/151‑1 (MQ=255) aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGc > 1:514136/1‑151 (MQ=255) tGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGcc > 2:47054/1‑151 (MQ=255) cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGcgcgcg > 1:938590/1‑149 (MQ=255) cacTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGAtt > 2:869870/1‑151 (MQ=255) acTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTg > 1:410805/1‑151 (MQ=255) ttttttttttGATTTTTTGCAGGCCAGCCGCGATCCGTCGGTCTTTCCTTTTCCCCCGGCCTTCCCCGCCCCGCGACTTTCCCCCCCCC‑AACAACAAAC‑CCGC‑TCGCTGGCGCAGGTAAGAAAAA‑CCGCCACAGCGATGAGCGTg < 2:122159/140‑1 (MQ=255) aatatTTTTGATATGTTGCAGGCCAGTCGCGAGACGTGTGTCGTTCCGTTTGCCACTGCCTTTCCTGACCCGTGACTTTTCCCCCTCC‑AATAACTAAC‑CCGC‑TGGCTGACGAAGGTAAGCAAAA‑CCCCGACAGCGATGAGCTTCATTGaa > 1:467370/2‑150 (MQ=255) atatTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGaa > 1:777923/1‑149 (MQ=255) tttttGATTTTTTGCGGCCCAGCCGCACCCCGGGGTTAGTTCCTTTTGCCCCGCCTTCCCCCGACCCGC‑ATTTTTCCCCCCCCAAACAACAAAC‑CCGC‑TCGCGGGCCGAGGAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGaaaa < 2:216888/148‑1 (MQ=255) agccgcgACCCGTGGGTCGTCCCGTTTCCCCCGGCTTTCCCCGACCCGC‑ACCTTTCCCCCCCC‑AAAAACAAAACCCGC‑TCGCGGGCGAGGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGaa < 2:331432/147‑1 (MQ=255) aGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAAc > 2:501246/1‑151 (MQ=255) gccgcgtccgggcggtcgTCCCGTTGCCCGCGGCCTCCCCCCCCCCCCGCCTTTCCCCCCCCC‑AACAACAAAC‑CCCC‑CCGCGGGCGCAGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACt < 2:514004/139‑1 (MQ=255) gtcggtcgTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt > 2:385559/1‑151 (MQ=255) tgggtcgTCCGTTTTCCCCCGGCTTTCCCCGCCCCGC‑ATTTTTCCCCCCCCAAACAACAAAC‑CCGC‑TCGCGGGCGAAGTAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCt < 2:109629/148‑1 (MQ=255) tgttggcgctggcccTTTCCCGCCCCCCG‑CTTTTCCCCCCCA‑CACAATCAAA‑CCGCAGGGCGGGCGGAGAAAGGAAAAG‑CCGCCACAGCGAGGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAgga < 1:551654/127‑2 (MQ=255) cgttttGCCTCGCCTTTCCCCGCCCCGCGTTTTTCCCCCCCCA‑AAAAACAAAC‑CCGC‑CCGGGGGCGCGGAAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCt < 2:375535/144‑1 (MQ=255) ccttttCCCTCGCCCTTCCCCCCCCCGC‑ATTTTTCCCCCCCC‑AAAAAAAAAACCCCC‑TCGCGGGCGCGGAAAAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGt < 2:802372/132‑1 (MQ=255) tttcccccgCCCTTCCCCCCCCGGGGCTTTTCCCCCCCCAAACAACAAAC‑CCGC‑CCGCGGGCGAAGAAAAGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtcgtca < 1:362628/141‑2 (MQ=255) gCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATg > 1:714201/1‑151 (MQ=255) ttCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCCCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAggg > 2:282505/1‑151 (MQ=255) tcccccgccccccgcTTTTCCCCCCCCAAACAACAAAC‑CCCG‑CCGGTGGCGAGGGAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTAt < 1:422216/127‑1 (MQ=255) ccccgccccgcGTTTTTCCCCCCCCC‑AAAAAAAAAC‑CCGC‑CCGCGGGCGGGGGAAAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAgg < 2:912058/148‑1 (MQ=255) ccccccGCGTCTTTCCCCCCCCA‑AACCAC‑AAA‑CCTC‑CCGGCGGGGGGGGTAAAAAAAATCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAggg < 1:823344/144‑1 (MQ=255) ccccgagtgttccccccccc‑AAAAACTAAA‑CCCC‑TACGTGGCGGAGGGAAACAAAATCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAc < 2:251587/139‑1 (MQ=255) ccGCGACTTTTCCCCCTCC‑ATCAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc > 2:950278/1‑150 (MQ=255) gcgACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAtt > 2:887909/1‑151 (MQ=255) aacaacTAAA‑CCGC‑TCGCTCGTGCATTTAAGCAAAA‑CCGCCACAGCCATTAGCGTGATTGACAACTTCCCGCCAGGAAACGCAGAACTGTGTCAGGCTATTGCACTTCGCTAAGCCTTACAGGGCATCACCTTTTCTCCtgatg > 1:78357/1‑144 (MQ=255) aacaacAAAC‑CCGC‑TCGCGGGCGCAGAAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAAt < 2:356927/148‑1 (MQ=255) aaaaacAAAC‑CCGC‑GCGGGGGCGGAGAAAAGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg < 1:566307/147‑1 (MQ=255) aacccGC‑TCGTGGGCGCAGAAAAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCa < 2:361164/143‑1 (MQ=255) aacccGC‑TCGCGGCCACAGAAAAACAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTAc < 1:693712/147‑1 (MQ=255) cgggccatgaaaaGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCg < 2:430866/139‑1 (MQ=255) tGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGACGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGt > 1:502331/1‑150 (MQ=255) gcgcAGGTAAGCAAAA‑CCGCCACGGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGt > 1:803896/1‑148 (MQ=255) gcAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTag > 1:360288/1‑149 (MQ=255) gaaaaGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTag < 1:605786/143‑1 (MQ=255) aaGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTa < 1:402575/151‑1 (MQ=255) | GAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTTCCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGC‑TCGCTGGCGCAGGTAAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTA > NC_000913/1510250‑1510542 |
| 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: ‑ |