Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,982,523 | 0 | A | G | 50.0% | 26.1 / 11.4 | 32 | intergenic (‑136/+31) | otsB/araH | trehalose‑6‑phosphate phosphatase, biosynthetic/L‑arabinose ABC transporter permease |
Reads supporting (aligned to +/- strand): ref base A (3/13); new base G (16/0); total (19/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.22e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGG > NC_000913/1982402‑1982666 | tctcTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTCTGGG‑GGGGAAAAAAACTAAAATGTTGCAGaaaaaaaa > 2:126240/1‑149 (MQ=255) ctTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGCATGGG‑GGGGAAAAAAGCTAAAATGTTGCAGAAAAAAACAt > 2:201258/1‑150 (MQ=255) ttCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGGATGGTGGGGAAAAAAACCCAAATT‑TTGCCGaaaaaaagataaga > 2:421871/1‑143 (MQ=255) tCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGAc < 1:162722/151‑1 (MQ=255) ttATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGG‑GGGGAAAAAAACTAAATTGTTTTAGAAAAAAACATCAAACAg > 2:97077/1‑151 (MQ=255) tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGGAAAAAAGCTAAATTGTTGCAGAAAAAAACATCAGACAGTGCGTTTTGCTTTTTTCTTGTAACGGGCGAAGa > 2:172292/1‑151 (MQ=255) tGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGG‑GGGGAAAAAAGCTAAATTGTTGCAGAAAAAAACATCAGACAGTGCGTTTTGCTTTTTGCTTTTAACGGGCGAAGa > 2:78079/1‑151 (MQ=255) tAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTTGGTATGGG‑GGGGAAAAAAGCAAAAATATTTCAGAAAAAAAGATCAGACAGTGCGGTTTGCTTTTTTTTTTTAAAGGGCGaaaaat > 2:41961/1‑147 (MQ=255) aaGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTGTGGT‑GGGGAAAAAAACTAAAATGTTGCAGAAAAAAAAAACAGACAGTGCGGTTTGCTTTTTTTTTGTAACGGGCGAAGATCa > 1:52339/1‑151 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGGTATGGG‑GGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACATCAGACAGTGGGGTTTGCTTTTTTTTTGTAAAGGGCGAAGAACa > 2:353439/1‑150 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTTGCTTTGGG‑GGGGAAAAAAGCTAAAATTTTGCAGAAAAAAACAACAGAAAGTGCGGTTTGCTTTTTTCTTTTAAAGGGGGAAGAACa > 2:85986/1‑150 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGT‑GGGGAAAAACGCTAAATTTTTGCAGAAAAAAACATCAGACAGTGCGTTTCGCTTTTTTCTTTTAACGGTCGAAGATCa > 2:317320/1‑150 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGT‑GGGGAAAAAAACCAAAATTTTTCAGAAAAAAACAACAGAAAATGGGGTTTGCTTTTTTTTTTTAAAGGGGGAAGaaaa > 1:142079/1‑147 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGGAAAAAAGCTAAATTGTTGCAGAAAAAAAAATCAGACAGTGCGGTTTGCTTTTTTCTTGTAAAGGGCGAAGATCa > 2:38691/1‑150 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATGGT‑GGGGAAAAAAACTAAAATTTTTCAGAAAAAAACATCAGAAAGTGCGGTTTGCTTTTTGTTTGTAAAGGGCGGAGAtaa > 1:160673/1‑148 (MQ=255) aGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGGTGGGTATGGG‑GGGGAAAAAAACTAAAATGTTGCAGAAAAAAAGATAAGACAGAGCGGTTTGCTTTTTTCTTGTAAAGGGCggaaaaca > 1:350687/1‑143 (MQ=255) gTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGGTATGGG‑GGGGAAAAAAGCTAAATTGTTGCAGAAAAAAACATCAGACAGTGCGGTTTGCTTTTTTCTTGTTACGGGCGAAGATCAc > 2:280510/1‑150 (MQ=255) ataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTGGAAGATCACCGc < 1:323985/151‑1 (MQ=255) ataGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGc < 2:20258/151‑1 (MQ=255) gtgtGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGG‑GGGGAAAAAAGCTAAATTTTTGCAGAAAAAAACAACAGACAGTGCGGTTCGCTTTTTTTTTGTAACGGTCGGAGAAAACCGCCGCcag > 2:261338/1‑150 (MQ=255) aCGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGAt < 1:439895/151‑1 (MQ=255) gTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGc < 2:222291/150‑1 (MQ=255) cAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGt < 1:332050/151‑1 (MQ=255) aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc < 2:24796/151‑1 (MQ=255) aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc < 2:348573/151‑1 (MQ=255) aCGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTAc < 2:357341/151‑1 (MQ=255) aagaaCAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTg < 1:214461/150‑1 (MQ=255) gTTATGGCGGTTGGCTATGGT‑GGGAAAAAAAGCTAAATTGTTGCAGAAAAAAACATCAGAAAGTGCGGTTCGCTTTTTGCTTGTAAAGGGCGAAGATAACCCCTGCCCGAAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGaaaa > 2:100393/1‑150 (MQ=255) tATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAatat < 2:469871/151‑1 (MQ=255) cGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGc < 1:359525/151‑1 (MQ=255) tGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGt > 1:77025/1‑151 (MQ=255) tggt‑ggGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATgg < 2:350102/151‑1 (MQ=255) | TCTCTTTTCTTATTGTGTTGCGCAGAACGCGCTTTACAGACATGTAAGTATAGACAGTGTGACGGGGGTCGCCATTTGCAACGAAGAACAGCCAACCGTTATGGCGGTTGGCTATGGT‑GGGAAAAAACGCTAAATTGTTGCAGAAAAAAGCATCAGACAGTGCGTTTCGCTTTTTGCTTGTAACGGTCGAAGATCACCGCTGCCAGCAGGATTAAGCCGCGAACCACGTACTGCGCGAAAGGAGAAATATTAAGCAGGTTCATGG > NC_000913/1982402‑1982666 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |