Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,462,655 IS1 (+) +9 bp 100% intergenic (‑85/‑100) yjgM ← / → yjgN predicted acetyltransferase/conserved inner membrane protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 4462655 =0 (0.000)40 (1.500) 20/278 NT 100% intergenic (‑85/‑108) yjgM/yjgN predicted acetyltransferase/conserved inner membrane protein
?REL606 4521571 = NA (NA)noncoding (768/768 nt) IS1 repeat region
* ? REL606 = 44626630 (0.000)40 (1.510) 25/276 NT 100% intergenic (‑93/‑100) yjgM/yjgN predicted acetyltransferase/conserved inner membrane protein
?REL606 = 4506109 NA (NA)noncoding (1/768 nt) IS1 repeat region

GATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGAACTTACTTATCACTTTATTGAAT  >  REL606/4505964‑4506257
                                                                                                                                                 |                                                                                                                                                    
gaTAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCg                                                                                                                                                  >  2:399324/1‑150 (MQ=31)
     cgcTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGACGCCACTGGATCACCTCATAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACTGGCGtgttg                                                                                                                                             >  1:104871/1‑150 (MQ=255)
            gTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAAc                                                                                                                                      <  2:204893/150‑1 (MQ=40)
                   ggTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCAATGGAGCACCTCAAAAGCACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAg                                                                                                                               <  1:132947/150‑1 (MQ=255)
                     ggagtttGCCGTTACGCACTACCTCGTCAGTGGCTGAAAAGAAGTGACACTTGATCGAAACAGAGGCCACAGGAGTACCTCTAAAACACCATCATACACGAAATCAGTAAGTTGGCAGCATAACCGGCGTGTTGTTGCAACATCGTtgtg                                                                                                                             <  1:335037/146‑1 (MQ=255)
                           tGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGa                                                                                                                       <  2:83183/150‑1 (MQ=255)
                                    gCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTa                                                                                                              >  2:336635/1‑150 (MQ=255)
                                             cGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGcttctt                                                                                                     <  2:164943/150‑1 (MQ=255)
                                                     gCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTAttcgtt                                                                                             <  2:506589/150‑1 (MQ=255)
                                                             ggaggGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGCAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCa                                                                                     >  2:594098/1‑150 (MQ=255)
                                                                  gACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGcc                                                                                >  2:289051/1‑150 (MQ=255)
                                                                         gATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATAttt                                                                         >  1:346533/1‑150 (MQ=255)
                                                                                aCAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTg                                                                  <  1:398172/150‑1 (MQ=255)
                                                                                                  ccTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCGCGTGTTTTTGCAACATCGCAGAGTGTGGATTGCGCCTATATTTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTACTGGGTGCAATTGATTAtttt                                                >  1:171503/1‑150 (MQ=37)
                                                                                                               atcatACACTAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCct                                  >  2:515660/1‑150 (MQ=255)
                                                                                                                     cacTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGa                             >  1:290492/1‑150 (MQ=255)
                                                                                                                         aaaTCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCCTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGa                         >  1:525074/1‑150 (MQ=255)
                                                                                                                                   gTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGAACTTACTTAt               >  1:335179/1‑150 (MQ=255)
                                                                                                                                      ggCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGAACTTACTTATCAc            <  2:291319/150‑1 (MQ=255)
                                                                                                                                             tCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGAACTTACTTATCACTTTATTg     >  2:435632/1‑150 (MQ=255)
                                                                                                                                                ccGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGAACTTACTTATCACTTTATTGAAt  <  2:105176/150‑1 (MQ=255)
                                                                                                                                                ccGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGAACTTACTTATCACTTTATTGAAt  <  2:127158/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCGGCGTGTTGTTGCAACATCGTAGTGTGTGGATTGCGCCTAGAGCTTCTTTATTCGTTCGAGCGCATCGCCAATATTTTCTATTGGGTGCCATTGATTATTTTTGCTTGAAATCCCTCTCGAATGAACTTACTTATCACTTTATTGAAT  >  REL606/4505964‑4506257

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.