Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,510,395 | 0 | C | A | 48.8% | 25.5 / 33.0 | 41 | S131R (AGC→AGA) | ydcR | putative DNA‑binding transcriptional regulator and putative aminotransferase |
Reads supporting (aligned to +/- strand): ref base C (16/5); new base A (0/20); total (16/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.44e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTAA > NC_000913/1510249‑1510543 | cGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCaaaa < 1:770617/151‑1 (MQ=255) aTGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGc < 1:342570/151‑1 (MQ=255) cGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGACACAGCGAt > 2:90466/1‑150 (MQ=255) atCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCg > 1:199623/1‑151 (MQ=255) aaCACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTg > 1:886834/1‑150 (MQ=255) cttttATTTTTTTG‑ATATTTTGCAGCCCGGCCGCGACCCGTCGGTC‑TTCCCTTTGCCCCCGGCCTTC‑CCCGCCCGGCTATTTCCCCCCCCCC‑AAAAACAAACCCCCGCGCGGGCCGAGGAAAAAAAAACCGCCACAGCGATGAGc < 2:197792/142‑1 (MQ=255) taattttttGTATTTGTGGCAGGCCCGCCTCGTCCCGTCGGTCTCTCCTTTCC‑CGCCGCCCTTC‑CCCCCCCCGCGTCTTTCCCCCCCCC‑AACAACAAACCCGCTCGCGGGCGCGGGTAACAAAAACCGCCACAGCGATGAGCGTGAt < 2:691030/143‑1 (MQ=255) atatTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGaaa > 1:855880/1‑150 (MQ=255) tttttttt‑ATTTGTTGCAGCCCTGCCGCGTTCGGTCGTTCTTCCCTTTTC‑CCCCGCCCTTC‑CCCGCCCCGCGATTTTCCCCCCCCC‑AAAAAATAACCCGCTCGCTGGCGCAGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTGaaa < 1:816696/147‑1 (MQ=255) tatTTTTG‑ATATGTTGCGGGCCAGTCGCGCTCCGTCGGTTTTTCCGTTTG‑CCTCGGCCTTC‑CCCGCCCCGCGTCTTTCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGCAGGAAAGCAAAACCGCCACAGCGATGAGCGTGATTg < 1:461461/146‑1 (MQ=255) tttttA‑ATTTGTTGCAGGCCAGCCGCGACCCGTCGGTCTTCCCGTTTC‑CCTCGGCCTTC‑CCCGACCCGCGACTTTCCCCCCCCC‑AACAACTAACCCGCTCGCTGGCGCAGATAAACAAAACCGCCACAGCGATGAGCGTGATTGaaa < 1:554986/147‑1 (MQ=255) ttgttttgtggcggcccAGCCGCGACCCGTCGTTCTTCCCTTTTC‑CCCCGCCCTTC‑CCCGCCCCGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCTCGCTGGCGCAGGTAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTAc < 1:232549/135‑1 (MQ=255) atatGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGCTTGAAAACTTACCGcc > 2:628953/1‑151 (MQ=255) tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg > 2:121917/1‑151 (MQ=255) ggCCGGTCGCGACCCGTCGTTCGTCCCTTTTC‑CCCCGGCCTTC‑CCCCACCCGCGATTTTCCCCCCCCA‑AAAAACAAACCCGCCCGCTGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACg < 1:622037/149‑1 (MQ=255) gcgcTCCGTGGTTCTTCCCTTTCG‑CCTGGCCCTTC‑CCCGCCCCGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCTCGGGGGCGCGGAAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGc < 1:21690/150‑1 (MQ=255) cgggcgttttttCCTTTTT‑CCGCGGCTTTT‑CCCCCCCCGCCATTTTCCCCCCCCC‑AAAAACAAAACCGCTCCTGGGCGGGGAAAAGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc < 1:152191/138‑1 (MQ=255) gTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGACGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCt > 2:501305/1‑150 (MQ=255) tttcccccgCCTTTCCCCGCCCCGCGTTTTTCCCCCCCCA‑AACAACAAACCCCCCCGGGGGCGCAGAAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt < 1:925827/132‑1 (MQ=255) tG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTAt > 2:14902/1‑151 (MQ=255) ttt‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAgg > 1:677957/1‑151 (MQ=255) tccccccccccGCGGCTTTCCCCCCCCA‑AACATCTAACCCGCCCGCTGGCGAGGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACa < 1:258119/141‑1 (MQ=255) ccgccccgggcTTTTCCCCCCCCAAAAAACAAACCCCCCCGGGGGCGGAGTAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAg < 1:362037/136‑1 (MQ=255) ccGCGTTTTTCCCCCCCCC‑AACAACAAACCCGCCCGCGGGCGAAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCa < 1:154940/151‑1 (MQ=255) cgcttttcCCCCCCCC‑AACAACAAACCCGCTCGCGGGCGAAGGAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAtt < 1:217906/146‑1 (MQ=255) ccTTTTCCCCCCCCAAACAACAAACCCCCTCGCGGCCGCAGGAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtc < 2:277863/150‑1 (MQ=255) cccccTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat > 2:49922/1‑151 (MQ=255) c‑aacaaCTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg > 1:297717/1‑151 (MQ=255) aacaacTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTg > 2:276127/1‑150 (MQ=255) aacaaaAAACCCGCTCGCTGGCGAAGAAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAAt < 1:821320/148‑1 (MQ=255) aataaacccccccgctGGCGCAGTAAAGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCa < 2:476225/138‑1 (MQ=255) acTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCAtt > 1:736923/1‑151 (MQ=255) accccGCTCTGTGGCGCAGTTAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACt < 2:679352/148‑1 (MQ=255) aacccGCGCGCGGGCGCAGGAAAAAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTa < 1:326225/146‑1 (MQ=255) cgcccgcgGGGGCAGGAAACAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCg < 1:142496/151‑1 (MQ=255) gctcgctGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCgg > 2:514705/1‑151 (MQ=255) ctGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGt > 1:479297/1‑151 (MQ=255) ggCGCAGTAAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGc < 2:522546/147‑1 (MQ=255) cgagggaaaGAAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt < 2:534526/141‑1 (MQ=255) cAGGAAAACAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTg < 2:361978/136‑1 (MQ=255) gCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTaa > 1:543243/1‑150 (MQ=255) | CGAGATGAAGCAGTCGATATCAACACTTATATTTTTG‑ATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTG‑CCTCGGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCC‑AACAACTAAACCGCTCGCTGGCGCAGGTAAGCAAAACCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTAA > NC_000913/1510249‑1510543 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |