Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,682,468 | 0 | C | . | 20.0% | 78.7 / 15.5 | 30 | coding (310/720 nt) | rstA | DNA‑binding transcriptional regulator RstA |
Reads supporting (aligned to +/- strand): ref base C (13/11); new base . (2/4); total (15/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.51e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGA > NC_000913/1682338‑1682598 | cATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGccc > 2:141016/1‑133 (MQ=255) gACCATTTGTCGTGATTTCCGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCct < 1:211237/133‑1 (MQ=255) tGATTTACGCACAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACG‑CCCCTGCTGTTTTGCTAg > 1:608185/1‑133 (MQ=255) gACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACG‑CCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGt < 1:79311/133‑1 (MQ=255) aTTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCGTGCTGTTTTGCTCGAGGGTTTACGTGTGCATGTGCGTCAGa > 2:290410/1‑133 (MQ=255) aTTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGa > 2:4805/1‑133 (MQ=255) aTTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGa > 2:131634/1‑133 (MQ=255) aTTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACG‑CCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGaa > 2:18027/1‑133 (MQ=255) ctAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAg < 1:63781/133‑1 (MQ=255) aCCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCca > 2:534919/1‑133 (MQ=255) ctctctCGCTAGCGATATCAACGACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCcaca < 1:258141/133‑1 (MQ=255) tctcGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAGAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTg > 1:530890/1‑133 (MQ=255) gATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCa < 1:534919/133‑1 (MQ=255) gATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACG‑CCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCaa < 2:608185/133‑1 (MQ=255) aGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAg > 1:145393/1‑133 (MQ=255) aGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAg > 1:425751/1‑133 (MQ=255) cGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGTGTCAGAAAGAGCAAGCGACACTGACCAAAGGt > 2:172371/1‑133 (MQ=255) gATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTCTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTc > 1:589124/1‑133 (MQ=255) aCCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGa < 2:145393/133‑1 (MQ=255) aTCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCGAATGTCTTCTGGAAACGt > 1:230635/1‑133 (MQ=255) tGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTctct < 1:131634/133‑1 (MQ=255) gCACTGGATATGGGTGCCTGCGTCTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGa < 1:202962/133‑1 (MQ=255) cACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACG‑CCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACt < 1:141016/133‑1 (MQ=255) cTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTc < 2:589124/133‑1 (MQ=255) tGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTcc > 2:332296/1‑133 (MQ=255) gAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCt > 1:463369/1‑133 (MQ=255) aaTGGGTGCCTGCGACTATATTCTCAAAACGACG‑CCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACa < 1:18027/133‑1 (MQ=255) tatTCTCAGAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGc < 2:530890/133‑1 (MQ=255) tatTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGc < 1:332296/133‑1 (MQ=255) cgacgCCCACTGCTGCTTCGCTCGCGCGTTGACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCAt < 1:172371/133‑1 (MQ=255) cgCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGa > 1:612071/1‑133 (MQ=255) | CATGACCATTTGTCGTGATTTACGCGCAAAGTGGTCTGGACCGATTGTTCTTCTAACCTCTCTCGATAGCGATATGAACCACATCCTGGCACTGGAAATGGGTGCCTGCGACTATATTCTCAAAACGACGCCCCCTGCTGTTTTGCTAGCGCGTTTACGTTTGCATTTGCGTCAGAATGAGCAAGCCACACTGACCAAAGGTCTTCAGGAAACGTCTCTGACTCCCTACAAAGCCCTGCATTTCGGCACGTTGACCATCGA > NC_000913/1682338‑1682598 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |