Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 1,125,562 IS186 (+) +6 bp coding (131‑136/2544 nt) mdoH → glucosyltransferase MdoH

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 1125562 =1 (0.030)40 (1.000) 29/146 0.3 97.4% coding (131/2544 nt) mdoH glucosyltransferase MdoH
?REL606 = 1858406 NA (NA)noncoding (1343/1343 nt) IS186 repeat region
* ? REL606 = 11255671 (0.030)45 (1.140) 34/144 0.1 97.6% coding (136/2544 nt) mdoH glucosyltransferase MdoH
?REL606 1857064 = NA (NA)noncoding (1/1343 nt) IS186 repeat region

GAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTATTCCAGAA  >  REL606/1858333‑1858481
                                                                         |                                                                           
gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGt                                                                           <  2:1290217/76‑1 (MQ=32)
     aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACNCTTAAGTTAGCGCTTATGGGGTAAAcc                                                                      <  1:1318587/76‑1 (MQ=255)
         tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCg                                                                  <  1:1090804/76‑1 (MQ=255)
           ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCGNTANGTTAGCGCTTATGGNGTAAACCCCCggg                                                                >  2:1190773/1‑76 (MQ=255)
            tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGgc                                                               <  1:1042198/76‑1 (MQ=255)
            tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGgc                                                               >  2:435099/1‑76 (MQ=255)
               cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGcgcg                                                            >  2:1317439/1‑76 (MQ=255)
                gAGGCACCGAGGCGTCGCATTCNTCAGATGGTTAGACCNTTAAGTGAGCGCTGATGGCGCAAACGCCCGGNCGCGt                                                           >  1:21565/1‑76 (MQ=14)
                 aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGtt                                                          <  1:957921/76‑1 (MQ=255)
                 aGGCACCGAGGCGTCGCATTCTTCAGACGGTTCAACCCTTAAGANAGNGCTTATGGGGTAAANCCCCGGGCGCGtt                                                          >  2:1186106/1‑76 (MQ=18)
                    cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCt                                                       >  1:193425/1‑76 (MQ=255)
                     aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCtt                                                      >  2:551126/1‑76 (MQ=255)
                       cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTa                                                    >  2:65214/1‑76 (MQ=255)
                        gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCttttt                                                   >  1:867827/1‑74 (MQ=255)
                        gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTAt                                                   <  1:964816/76‑1 (MQ=255)
                            cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAg                                               <  1:947604/76‑1 (MQ=255)
                             gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGa                                              >  2:363315/1‑76 (MQ=255)
                              tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGAt                                             <  1:1084934/76‑1 (MQ=255)
                              tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGAt                                             <  2:511868/76‑1 (MQ=255)
                               cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATg                                            <  2:386034/76‑1 (MQ=255)
                                    tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGtt                                       >  2:279493/1‑76 (MQ=255)
                                       tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGa                                    >  1:867809/1‑76 (MQ=255)
                                       tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGa                                    >  1:128161/1‑76 (MQ=255)
                                        cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGaa                                   <  1:515599/76‑1 (MQ=255)
                                               ttCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCNGGGCGCGGTCCTTTAGCAAGATGACGTTCGAATANGCCg                            >  1:157499/1‑76 (MQ=38)
                                                tCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCgg                           >  2:1217385/1‑76 (MQ=255)
                                                 cAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGc                          <  2:313864/76‑1 (MQ=255)
                                                   aCCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAg                        >  2:709356/1‑76 (MQ=255)
                                                       ttAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGgcgc                    <  2:999246/76‑1 (MQ=255)
                                                         aaGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCNAGNTGACGTTCGAATATNCCGGCAGGCGCAg                  >  2:1189085/1‑76 (MQ=255)
                                                          aGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAgt                 >  1:262405/1‑76 (MQ=255)
                                                           gTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAgtt                >  1:1264667/1‑76 (MQ=255)
                                                              aGNGCTAATGGAGTAAACCCCCGGGTGCGTTCCTTTATCAAGATGACGTTCGAATATGGCGGCAGGCGCAgtttgt             <  2:1236881/76‑1 (MQ=255)
                                                                  cTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTAt         <  2:877809/76‑1 (MQ=255)
                                                                     aTGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTATTcc      <  1:723766/76‑1 (MQ=255)
                                                                     aTGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTATTcc      <  1:735951/76‑1 (MQ=255)
                                                                     aTGGGGTAAACCCCCGGGCGCGTTCCTNTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTATTcc      >  1:1307934/1‑76 (MQ=255)
                                                                        gggTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTATTCCAGa   <  2:953111/76‑1 (MQ=255)
                                                                         ggTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTATTCCAGaa  <  2:626112/76‑1 (MQ=255)
                                                                         |                                                                           
GAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGGTAAACCCCCGGGCGCGTTCCTTTATCAAGATGACGTTCGAATATGCCGGCAGGCGCAGTTTGTTTATTCCAGAA  >  REL606/1858333‑1858481

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