Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 1,776,434 IS150 (–) +4 bp 56.9% intergenic (‑225/+56) pheS ← / ← pheM phenylalanyl‑tRNA synthetase alpha subunit/phenylalanyl‑tRNA synthetase operon leader peptide

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1776434 =33 (0.390)48 (0.470) 28/48 NT 54.8% intergenic (‑225/+59) pheS/pheM phenylalanyl‑tRNA synthetase alpha subunit/phenylalanyl‑tRNA synthetase operon leader peptide
?REL606 = 2775877 NA (NA)noncoding (1/1443 nt) IS150 repeat region
* ? REL606 = 177643733 (0.390)39 (0.380) 21/48 NT 49.6% intergenic (‑228/+56) pheS/pheM phenylalanyl‑tRNA synthetase alpha subunit/phenylalanyl‑tRNA synthetase operon leader peptide
?REL606 2774435 = NA (NA)noncoding (1443/1443 nt) IS150 repeat region

TTTCAAAACGCGCATACAAAAAAAGCCTCCACTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/1776404‑1776437
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGTACTGACCCCAAAAAGTTGGACAGTT  >  REL606/2774435‑2774462
                                                              
TTTCAAAACGCGCATACAAAAAAAGCCTNNACTGN                             >  1:10900953/1‑35
 TTCAAAACGCGCATACAAAAAAAGCCTCCACTGTG                            >  1:17360307/1‑35
 TTCAAAACGCGCATACAAAAAAAGCCTCCACTGTG                            >  1:4959412/1‑35
    AAAACGCGCATACAAAAAAAGCCTCCACTGTGTAC                         >  1:16680911/1‑35
     AAACGCGCATACAAAAAAAGCCTCCACTGTGTACT                        >  1:8988467/1‑35
     AAACGCGCATACAAAAAAAGCCTCCACTGTGTACT                        >  1:12805942/1‑35
       ACGCGCATACAAAAAAAGCCTCCACTGTGTACTGA                      >  1:20390632/1‑35
       ACGCGCATACAAAAAAAGCCTCCACTGTGTACTGA                      >  1:3738545/1‑35
        CGCGCATACAAAAAAAGCCTCCACTGTGTACTGAC                     >  1:103449/1‑35
         GCGCATACAAAAAAAGCCTCCACTGTGTACTGACC                    >  1:2866294/1‑35
         GCGCATACAAAAAAAGCCTCCACTGTGTACTGACC                    >  1:14164902/1‑35
          CGCATACAAAAAAAGCCTCCACTGTGTACTGACCC                   >  1:20306764/1‑35
          CGCATACAAAAAAAGCCTCCACTGTGTACTGACCC                   >  1:19073783/1‑35
           GCATACAAAAAAAGCCTCCACTGTGTACTGACCCC                  >  1:16001335/1‑35
            CATACAAAAAAAGCCTCCACTGTGTACTGACCCCA                 >  1:779258/1‑35
            CATACAAAAAAAGCCTCCACTGTGTACTGACCCCA                 >  1:196451/1‑35
              TACAAAAAAAGCCTCCACTGTGTACTGACTCCAAA               >  1:12410630/1‑35
              TACAAAAAAAGCCTCCACTGTGTACTGACCCCAAA               >  1:11713723/1‑35
               aaAAAAAAAGCCTCCACTGTGTACTGACCCCAAAA              <  1:10568634/33‑1
                aAAAAAAACCCTCCACTGTGTACTGACCCCAAAAA             <  1:17696431/34‑1
                CAAAAAAAGCCTCCACTGTGTACTGACCCCAAAAA             >  1:6082655/1‑35
                CAAAAAAAGCCTCCACTGTGTACTGACCCCAAAAA             >  1:5825773/1‑35
                 AAAAAAAGCCTCCACTGTGTACTGACCCCAAAAAG            >  1:8259866/1‑35
                  AAAAAAGCCTCCACTGTGTACTGACCCCAAAAAGT           <  1:17336871/35‑1
                  AAAAAAGCCTCCACTGTGTACTGACCCCAAAAAGT           >  1:14045052/1‑35
                  AAAAAAGCCTCCACTGTGTACTGACCCCAAAAAGT           <  1:3235996/35‑1
                    AAAAGCCTCCACTGTGTACTGACCCCAAAAAGTTG         >  1:14929060/1‑35
                    AAAAGCCTCCACTGTGTACTGACCCCAAAAAGTTG         >  1:6394271/1‑35
                     AAAGCCTCCACTGTGTACTGACCCCAAAAAGTTGG        <  1:8019388/35‑1
                      aacCCTCCACTGTGTACTGACCCCAAAAAGTTGGA       <  1:5513946/32‑1
                      AAGCCTCCACTGTGTACTGACCCCAAAAAGTTGGA       >  1:6397845/1‑35
                       AGCCTCAACTGTGTACTGACCCCAAAAAGTTGGAC      <  1:5899527/35‑1
                        nnCTCNNCTGTGTACTGACCCCAAAAAGTTGGACA     <  1:17347595/33‑1
                        cCCTCCACTGTGTACTGACCCCAAAAAGTTGGACA     <  1:17666143/34‑1
                        GCCTCCACTGTGTACTGACCCCAAAAAGTTGGACA     <  1:4777959/35‑1
                         aCTCCACTGTGTACTGACCCCAAAAAGTTGGACAG    <  1:15208444/34‑1
                         CCTCCACTGTGTACTGACCCCAAAAAGTTGGACAG    >  1:7522300/1‑35
                         CCTCCACTGTGTACTGACCCCAAAAAGTTGGACAG    <  1:13429188/35‑1
                          ccCCACTGTGTACTGACCCCAAAAAGTTGGACAGT   >  1:9347595/3‑35
                          CTCCACTGTGTACTGACCCCAAAAAGTTGGACAGT   <  1:13448566/35‑1
                          CTCCACTGTGTACTGACCCCAAAAAGTTGGACAGT   <  1:13475131/35‑1
                          CTCCACTGTGTACTGACCCCAAAAAGTTGGACAGT   <  1:14723222/35‑1
                           TCCACTGTGTACTGACCCCAAAAAGTTGGACAGTT  <  1:10346386/35‑1
                                                              
TTTCAAAACGCGCATACAAAAAAAGCCTCCACTG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/1776404‑1776437
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TGTACTGACCCCAAAAAGTTGGACAGTT  >  REL606/2774435‑2774462

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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.