Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | NC_000913 | 16,490 | 0 | T→G | 100.0% | 0.9 | 5 | I349S (ATT→AGT) | insL1 | IS186 transposase |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/3): ref base (2/0): total (2/3) | ||||||||||
Fisher's exact test strand distribution p-value = 1.00e-01 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
GAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCA > NC_000913/16435‑16781 | gAACCTGAATGCGCGAAAGCGTGGATATTTGCGAATCTAGTCGCCGCATTTTTAAGTGAGTACATAATCCAGCCATAGCTGGAGTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATc < 2:234649/299‑1 (MQ=2) aaaGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAAGTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGctc < 2:112158/301‑1 (MQ=21) aataaagaCTCGCCTCAGTTTGAAGTGACGAAAAAATCAACGCATTCTGGGATTGCCACCCCGGCAGAGACGGATCCGAAAAGAAGAACGATGTCGTGGTGGAGAATAACAAAAATGCTCATCTGTAGCTCACAGGTGGCCATTCGTGGGACAGTATCCCTTACAGCCTACAAAAAGCAATAGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGTGGCGCCGCATTCGCCAGAAGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCAGTTGCTCTTACTCGGATTCGTaa < 1:639069/294‑1 (MQ=14) aCTCGCCGCATTTTTAATTGCCGACATAATCCAGCCATCGCTGGATTTCCCCCACAGAAGTGCCGGATCAGACAAGAAGACCTCACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGCAAGAATCCCTCACCGCCTACAACACGCCATTGAACACCGCGAGGCATCGACTTAACGAGGACCCGAAGCGTCGCCTTCTTCACATGGTTCAACCCTTAAGTTACCGCTTATGTGACCACCACCCGCCGTTCCTCTGCCCCGGCTTCggcccgttga > 2:639043/1‑291 (MQ=11) cATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGACGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGGGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCTTTGCTCTTCCTCGGCTTCTTAATCCTTGAAAACAGCa > 2:670065/1‑301 (MQ=255) | GAACCTGAACTCGCGAAAGCGTGGATATTTGCTAATCTACTCGCCGCATTTTTAATTGACGACATAATCCAGCCATCGCTGGATTTCCCCCCCAGAAGTGCCGGATCCGAAAAGAAGAACTAACTCGTTGTGGAGAATAACAAAAATGGTCATCTGGAGCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGATCACTCCCCGCCGTTGCTCTTACTCGGATTCGTAAGCCGTGAAAACAGCA > NC_000913/16435‑16781 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
Reads not counted as support for junction |
read_name Not counted due to insufficient overlap past the breakpoint. |
read_name Not counted due to not crossing MOB target site duplication. |