Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,177,748 | 0 | C | A | 36.6% | 33.1 / 15.0 | 41 | E55* (GAG→TAG) | nudF | ADP‑ribose pyrophosphatase |
Reads supporting (aligned to +/- strand): ref base C (15/11); new base A (0/15); total (15/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.41e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGG > NC_000913/3177599‑3177888 | acTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTc > 1:307935/1‑150 (MQ=255) cACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAAGCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGTTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCaaaaa < 1:80820/150‑1 (MQ=255) cACCCTCTTCAATCAGCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCaaaa < 1:142285/149‑1 (MQ=255) cccTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAAttt < 2:428174/151‑1 (MQ=255) atcCCGACAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCTGCAATCCGAATCTGCTCAATCAGCACACCTTCATCACACACTGGTTCAAA‑GGGTAACAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCAc < 1:69558/149‑1 (MQ=255) cGGCAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATg > 1:597740/1‑151 (MQ=255) ggCAACCACCTCCAGTAGCAAAGGGGTTTCT‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCCCTGGGTAAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGa < 2:95340/151‑1 (MQ=255) ctcCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCAtt < 2:254365/148‑1 (MQ=255) ccAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTCCGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCCCGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGCCC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGcc < 2:155505/150‑1 (MQ=255) ccAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTAAACCAGCACAACTTCATCCCGCACTGGGTAAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTg < 2:468633/148‑1 (MQ=255) cAAGGGATTTCGCCGGGTGTCGTCCGCGGCAATCCAAATCTCCACAACCAGAACATCGTCAACACGCATTGGGGCAAG‑GGGTACCAAGCCTGCGGGGTGACG‑GCCCAAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGaa < 1:414965/150‑1 (MQ=255) ggggTTTCG‑CTGGTGTCGTACGCGGCAATCCGCATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAg > 1:357837/1‑149 (MQ=255) ggggTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGAc > 2:547222/1‑151 (MQ=255) tttCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCACCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGAt < 2:45039/150‑1 (MQ=255) cGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat > 2:368111/1‑150 (MQ=255) cGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTAAAA‑GGGTAGCAAGACTGCGGCGTGCCC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtat < 2:421304/150‑1 (MQ=255) gTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCtata < 2:506102/150‑1 (MQ=255) aCGCGGCAATCCAAATCTGCCCAACCAGCACAATTTACCCCCGCCATGGTTCAAA‑GGGCAGCAACACTCGGTCTTGTCC‑CCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACTATGACGAAATCTATAAAg < 1:227893/151‑1 (MQ=255) cgGCAATCCGAATCTGCTCAATCAGCACAACTTCACCACCCACTGGGTCAAA‑GGGTAGAAAGCCTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATc < 2:574887/151‑1 (MQ=255) cgGCAATCCGAATCTGCCCAATCAGCACAACTCCACCACTCACGGGGTCAAA‑GGGTAGCAAGACGGCCGCGTTACC‑GCGCTAAAAAGTTTCCCGCCGCACCTCATGACTCATGTGCCCGTTGAAGAGACGATGACGAAATCTATaa < 1:498505/146‑1 (MQ=255) aaTCCAAATTTGCGAAATCAGCAAAACTTCACCCCGCACTGGCAAAAA‑GGCTAGCAAACCGGCGGCGGAACC‑GCTCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAAt < 2:416314/151‑1 (MQ=255) aTCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACA‑GCGCTCAAAAATTTCCCGCCGCCCCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCt > 1:138705/1‑147 (MQ=255) tCCGAATCTGCTCAATCCGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATg > 1:325913/1‑150 (MQ=255) ccGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGa > 2:103710/1‑150 (MQ=255) tCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa > 2:269469/1‑150 (MQ=255) gCTCAATCAGCACAACTTCATCCCGCACGGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGaaaaaa < 2:359735/147‑1 (MQ=255) gCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCAGTTGAATAGACGATGACGACATCTATAAAGATCTAATGAAAAAAAGcc > 1:520820/1‑151 (MQ=255) ccAGCACAACTTCATCACGCACT‑GGTCAAAGGGGTAGCAAGACTGCGG‑GTGGCCGGCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGAt < 1:626843/149‑1 (MQ=255) acaACTTCATCACGCACTGGGTAAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgt < 2:434033/150‑1 (MQ=255) cagcTTCAGCCCGCACGGGGTCAAG‑GGGTAGCAAGACTGCGGGGTGCCC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAgtgt < 2:395771/148‑1 (MQ=255) atcaCGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGt > 2:552321/1‑151 (MQ=255) cGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaa < 2:649546/151‑1 (MQ=255) gCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaa > 1:463896/1‑150 (MQ=255) cACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaata < 2:126981/151‑1 (MQ=255) gggTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAAttt > 2:203460/1‑151 (MQ=255) gggTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCaataa < 2:73122/148‑1 (MQ=255) gCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGCTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCACTAATTTCTACATCGtttt > 1:518941/1‑150 (MQ=255) cggcgTGACC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaaa < 2:405028/148‑1 (MQ=255) ggcgTGCCC‑GCGCTAAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCaaa < 2:449662/147‑1 (MQ=255) ggcgTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGt > 1:37782/1‑150 (MQ=255) gTGACC‑GCCCTAAAAAATTCCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGt < 1:445577/147‑1 (MQ=255) tGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACgg > 1:589446/1‑151 (MQ=255) | ACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCG‑CTGGTGTCGTACGCGGCAATCCGAATCTGCTCAATCAGCACAACTTCATCACGCACTGGGTCAAA‑GGGTAGCAAGACTGCGGCGTGACC‑GCGCTCAAAAATTTCCCGCCGCACCTCATGACTCATTTGCCCGTTGAATAGACGATGACGAAATCTATAAAGATCTAATGAAAAAAAGCCGCGATAAAGTGTTTCTCGTGCAATAATTTCTACATCGTTTTTGCCAAATGTAACGG > NC_000913/3177599‑3177888 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |