Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,415,710 IS150 (–) +3 bp 100% coding (848‑850/1413 nt) cycA → D‑alanine/D‑serine/glycine transporter

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3893554 =NA (NA)38 (1.420) 27/278 NT 100% noncoding (1/1443 nt) IS150 repeat region
?REL606 4415710 = 0 (0.000)coding (848/1413 nt) cycA D‑alanine/D‑serine/glycine transporter
* ? REL606 = 3894996NA (NA)40 (1.510) 27/276 NT 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 = 4415712 0 (0.000)coding (850/1413 nt) cycA D‑alanine/D‑serine/glycine transporter

TCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGT  >  REL606/3893409‑3893699
                                                                                                                                                 |                                                                                                                                                 
tgaaggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtaCAGGCTGTAc                                                                                                                                               <  1:323353‑M2/10‑1 (MQ=255)
tCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTAc                                                                                                                                               <  2:443188/150‑1 (MQ=255)
tCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTAc                                                                                                                                               <  2:37268/150‑1 (MQ=255)
tCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTACGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTAc                                                                                                                                               <  1:31594/150‑1 (MQ=255)
    ggatgctgttacggaatatattgaatactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtgtttaactgtccaactttttggggtcagtaCagactgtctctta                                                                                                                                           >  2:137840‑M2/137‑144 (MQ=255)
     aaaCGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCAc                                                                                                                                          <  2:10990/150‑1 (MQ=255)
     aaaCGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCAc                                                                                                                                          <  2:467372/150‑1 (MQ=255)
          cccGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCAtt                                                                                                                                     <  1:480272/150‑1 (MQ=255)
                gTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTttgttg                                                                                                                               >  2:633515/1‑150 (MQ=255)
                   tttGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGAc                                                                                                                            <  2:447006/150‑1 (MQ=255)
                       cGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATg                                                                                                                        >  2:117003/1‑150 (MQ=255)
                             tactacaacagcagaagaattagcctgaaattaaaaggtctgactccaattgaatatcggaatcagacctatatgcctcgtttttaactgtccaactttttgtggtcagtaCAGGCTTTACTGCACCCATTTTGTTGGACGATGAAATgg                                                                                                                  >  1:312834‑M2/112‑150 (MQ=255)
                                   aaGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGATGATGAAATGGAATAGc                                                                                                            >  2:169747/1‑150 (MQ=255)
                                        gTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa                                                                                                       >  2:598636/1‑150 (MQ=255)
                                        gTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTAACGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTa                                                                                                       >  1:166143/1‑150 (MQ=255)
                                         tACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTaa                                                                                                      <  1:35100/150‑1 (MQ=255)
                                          aCTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAat                                                                                                     <  2:166617/150‑1 (MQ=255)
                                          aCTTGCTGCTGCAACGTAATGCGCCGCGGGCACCAGATCCATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCAGTTTGTTGGACGATGAAATGGAATAGCCCCTAat                                                                                                     <  1:454615/150‑1 (MQ=255)
                                            ttGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACAGTACTGCACCCATTTTGTTGGACGATGAACTGGAATAGCCCCTAatat                                                                                                   >  1:99036/1‑150 (MQ=255)
                                               ctgctgcAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTc                                                                                                <  1:410819/150‑1 (MQ=255)
                                                           aaTGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATa                                                                                    >  1:229521/1‑150 (MQ=255)
                                                                       gCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGaaaa                                                                        >  2:17012/1‑150 (MQ=255)
                                                                        cACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAg                                                                       <  2:224656/150‑1 (MQ=255)
                                                                                   gtgtataagagacagtatatgcctcgtgtttaactgtccaactttttggggtcagtaCAGGCTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGt                                                            <  2:340858‑M2/93‑1 (MQ=255)
                                                                                                  tatatgcctcgtgtttaactgtccaactttttggggtcagtaCAGGCTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAaagcgccttgaagtctgtctcttatacac                                             >  1:339882‑M2/43‑136 (MQ=255)
                                                                                                        gTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAAc                                       <  2:99328/150‑1 (MQ=255)
                                                                                                        ataagagacaggggTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAAc                                       <  2:676608/139‑1 (MQ=255)
                                                                                                            gtggttaactgtccaactttttggggtcagtaCAGGCTGTACTGCCCCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACtgat                                   >  1:209703‑M2/33‑150 (MQ=255)
                                                                                                                   gggTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAAcctgtctcttat                            >  1:674622/1‑139 (MQ=255)
                                                                                                                   gggTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGtt                            >  2:19758/1‑150 (MQ=255)
                                                                                                                        cTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAgg                       <  1:631659/150‑1 (MQ=255)
                                                                                                                          cAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGat                     <  2:230172/150‑1 (MQ=255)
                                                                                                                          agggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGat                     <  2:666122/148‑1 (MQ=255)
                                                                                                                            ggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGatct                   >  1:664161/1‑149 (MQ=255)
                                                                                                                            ggTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGatca                   >  2:535829/1‑150 (MQ=255)
                                                                                                                              tCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAAATGATGATGGTTACAGGatcatc                 <  2:347263/150‑1 (MQ=255)
                                                                                                                                            CAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTgg   <  1:19695‑M2/150‑1 (MQ=255)
                                                                                                                                             AGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTttgat  >  2:556115‑M2/1‑148 (MQ=255)
                                                                                                                                             AGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGACACGTTTTGGt  >  1:19327‑M2/1‑150 (MQ=255)
                                                                                                                                             AGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGAGGGTTACAGGATCATCTCGGCACGTTTTGGt  >  1:224235‑M2/1‑150 (MQ=255)
                                                                                                                                             AGACTGTACTGCACCCATTTTGTTGGACAATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAATGGATGATGGTTACAGAACCATCTCGGCACGTTTTGGt  >  1:376337‑M2/1‑150 (MQ=255)
                                                                                                                                             AGACTGTACTACACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATAGGTCAAAGCCAAAATACCCTTTTTAACAGCGCCTTGAAGTAGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGcacgtttggtt  >  1:622459‑M2/1‑146 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCGGCAAACGCCCGTGGTATTTGCGTCTGCGCGGCAAGCTGTACTTGCTGCTGCAACGTAATGCGCTGCGTGCACCAGATCAATTTGAAATCCCGCCAAACAGGGTTATCGAACTGGGTACTCAGGTCGAAATCTAACGCCAGACTGTACTGCACCCATTTTGTTGGACGATGAAATGGAATAGCCCCTAATATGTCAAAGCCAAAATACCCTTTTGAAAAGCGCCTTGAAGTCGTGAATCACTACTTCACAACTGATGATGGTTACAGGATCATCTCGGCACGTTTTGGT  >  REL606/3893409‑3893699

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ 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.