| New junction evidence | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | reads (cov) | reads (cov) | score | skew | freq | annotation | gene | product | ||
| * | ? | NC_000913 | = 3588435 | 16 (0.720) | 3 (0.170) | 3/108 | 1.7 | 18.2% | coding (746/1071 nt) | ugpC | sn‑glycerol‑3‑phosphate ABC transporter ATPase |
| ? | NC_000913 | = 3589209 | 14 (0.780) | coding (819/846 nt) | ugpE | sn‑glycerol‑3‑phosphate ABC transporter permease | |||||
| Rejected: Frequency below/above cutoff threshold. | |||||||||||
TGACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGCCGATAAAACTCGCTACAAACAGGCT > NC_000913/3588372‑3588499 | tGACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTAtt > 2:129471/1‑68 (MQ=255) gccGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTc > 1:778044/3‑68 (MQ=255) gACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTgcgcg > 1:727928‑M1/1‑63 (MQ=255) gACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTc > 1:307279/1‑68 (MQ=255) gACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTc > 2:359395/1‑68 (MQ=255) gACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTc > 1:679636/1‑68 (MQ=255) gACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCAGTCCAGTTCGAACTGCGTGCCTTCGTTATTc > 2:149586/1‑68 (MQ=255) cGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTCACGCGGCCTGTc < 2:390165/68‑1 (MQ=255) tttAGCGGCAGCTCAATACCGCCGTCCAGTTCGATATGCGTGCCTTCGTTATTCACGCGGCCTGTctg < 2:514221/68‑3 (MQ=255) ttAGCGGCAGCTCAATGCCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTCACGCGGCCTGTcagc > 1:544704/1‑68 (MQ=255) cGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTCACGCGGCCTGTCAGCAGGt < 1:588630/68‑1 (MQ=255) aCCGCCGTCCAGTTCGAAATGCGTGCCTTCGTgcgcggcctggtcgatagtgagaaataagatggcag > 2:589884‑M1/1‑32 (MQ=255) gccgTCCAGTTCGAAATGCGTGCCTTCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGc > 2:286824/1‑68 (MQ=255) gccgTCCAGTTCGAAATGCGTGCCTTCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGc > 2:498658/1‑68 (MQ=255) cgTCCAGTTCGAAATGCGTGCCTTCGTgcgcggcctggtcgatagtgagaaataagatggcaggactg < 1:589884‑M1/68‑42 (MQ=255) aaaTGCGTGCCTTCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGCCGATAAAACTCGc < 2:709626/68‑1 (MQ=255) gTGCCTTCGTgcgcggcctggtcgatagtgagaaataagatggcaggactgaaattacaggcagtaac < 2:727928‑M1/68‑59 (MQ=255) ccTTCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGCCGATAAAACTCGCTACAAACAg > 1:89654/1‑68 (MQ=255) cTTCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGCCGATAAAACTCGCTACAAACAgg < 1:286825/68‑1 (MQ=255) ttCGTgcgcggcctggtcgatagtgagaaataagatggcaggactgaaattacaggcagtaaccaaaa > 1:531491‑M1/1‑5 (MQ=255) tCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGCCGATAAAACTCGCTACAAACAGGCt < 2:679636/68‑1 (MQ=255) | TGACGGTAGCCACCGTTTAGCGGCAGCTCAATACCGCCGTCCAGTTCGAAATGCGTGCCTTCGTTATTCACGCGGCCTGTCAGCAGGTTCATCGCCGGACTGCCGATAAAACTCGCTACAAACAGGCT > NC_000913/3588372‑3588499 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ 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. |