New junction evidence | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
* | ? | NZ_CP009273 | 2808944 = | 2 (0.090) | 22 (1.020) | 15/164 | 0.3 | 91.7% | coding (855/1557 nt) | gshA | glutamate‑‑cysteine ligase |
? | NZ_CP009273 | 3576788 = | NA (NA) | coding (391/417 nt) | BW25113_RS17825 | YrhA family protein |
GGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAA > NZ_CP009273/2808855‑2809031 | ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGc < 2:194178/90‑1 (MQ=255) gACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCctc < 1:273295/90‑1 (MQ=255) gacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCg < 2:97982‑M1/10‑1 (MQ=255) agctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCGTCt > 2:372097‑M1/78‑90 (MQ=255) agctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCGTCt > 2:98848‑M1/78‑90 (MQ=255) agctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCGTCt > 2:507693‑M1/78‑90 (MQ=255) tAACGCGTTTTGGACGAATTGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactc < 2:645179‑M2/90‑11 (MQ=255) tAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactc > 1:151945‑M2/1‑80 (MQ=255) aCGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactcca < 1:388200‑M2/90‑13 (MQ=255) aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactccaacttattgata > 2:388200‑M2/1‑67 (MQ=255) cGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactccaacttattgatag > 1:93022‑M2/1‑66 (MQ=255) gcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTActc > 1:458287‑M1/52‑90 (MQ=255) tACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactccaacttattgatagtgttttatgttc > 2:560470‑M2/1‑54 (MQ=255) aaaacaccatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATg > 2:98909‑M1/44‑90 (MQ=255) aaaacaccatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATg > 2:634670‑M1/44‑90 (MQ=255) gTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactccaacttattgatagtgttttatgttcagat < 1:225452‑M2/90‑41 (MQ=255) accatcatacactaaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGtt < 1:8083‑M1/52‑1 (MQ=255) tCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatg < 1:440939‑M2/48‑7 (MQ=255) tCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatg > 2:440939‑M2/1‑42 (MQ=255) acactaaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGc > 2:237839‑M1/31‑90 (MQ=255) cactaaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc > 2:145951‑M1/30‑90 (MQ=255) gCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGgtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactt > 1:615198‑M2/1‑35 (MQ=255) aaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc > 2:165791‑M1/26‑86 (MQ=255) aaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc < 1:165791‑M1/61‑1 (MQ=255) aaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGtt < 2:92754‑M1/65‑1 (MQ=255) aaatcagtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGtt < 1:175351‑M1/65‑1 (MQ=255) gtaagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGc < 2:396458‑M1/71‑1 (MQ=255) taagttggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGcc > 2:453102‑M1/19‑90 (MQ=255) tggcagcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa > 2:175351‑M1/14‑90 (MQ=255) gcatcaccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTAc > 1:396458‑M1/9‑90 (MQ=255) tttGCAGCCTCTTACCGgtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccg < 1:217421‑M2/90‑74 (MQ=255) accCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc < 2:19766‑M1/61‑1 (MQ=255) accCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc > 1:19766‑M1/4‑64 (MQ=255) ccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaa < 1:354527‑M1/88‑1 (MQ=255) ccCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaa > 1:157475‑M1/3‑90 (MQ=255) | GGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAA > NZ_CP009273/2808855‑2809031 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ 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. |