Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 2,238,948 Δ1 bp :: IS186 (–) +6 bp :: Δ2 bp coding (337‑342/1521 nt) mglA ← D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2238948 =0 (0.000)15 (0.590) 15/482 0.6 100% coding (342/1521 nt) mglA D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit
?NC_000913 2514275 = NA (NA)noncoding (3/1345 nt) IS186 repeat region
* ? NC_000913 = 22389530 (0.000)18 (0.710) 18/480 0.4 100% coding (337/1521 nt) mglA D‑galactose/methyl‑galactoside ABC transporter ATP binding subunit
?NC_000913 = 2515616 NA (NA)noncoding (1344/1345 nt) IS186 repeat region

ATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATAAAATGGTGCCGGAGTCTTTTTGATAAATACCAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATT  >  NC_000913/2238707‑2239193
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
gcagtgacttcacgtaatgacatccccccggggccgtaagccagccccagacgtagcagagttgcagcatcacgaatttcgcggcggcgggttagagccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttcaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGcccc                                                                                                                                                                                                                                              <  1:164555‑M1/10‑1 (MQ=255)
                gAACAGATGATTGACCTCTTTTTCGGTTAACGAATAAGTCTGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTTGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaa                                                                                                                                                                                                                               <  1:222781‑M1/250‑20 (MQ=255)
                   cAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTAAAAAATCGCTTTGGTTTAGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggt                                                                                                                                                                                                                           >  1:92395‑M1/1‑228 (MQ=255)
                                     ttCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGTATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaat                                                                                                                                                                                                          <  1:207060‑M1/250‑41 (MQ=255)
                                                        cagagttgcagcatcacgaatttcgcggcggcggggtagagccccggcattacgtgccgaagtatccagttcttcgggcgtaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTAc                                                                                                                                                                                      >  2:20423‑M1/186‑251 (MQ=255)
                                                               gcagaacaacgaatttcgcggctgcgggttagagccccggcattacgggccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacggtaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGa                                                                                                                                                                                >  2:154117‑M1/179‑250 (MQ=255)
                                                                    atcacgaatttcgcggcggcgggttagagccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATAcc                                                                                                                                                                          >  1:4515‑M1/174‑251 (MQ=255)
                                                                                 cgTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCCGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGGTATCGcccataagcgctaacttaagggttgaaccctctgcagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttca                                                                                                                                                              >  1:288820‑M1/1‑166 (MQ=255)
                                                                                           aaaCGCTTTGGCGATTTCGATCATCTGCATTTGCGACCCGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAACATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacacctcggtgcctcgttaggacgatgcctcgcgttcttcaattgcgttttg                                                                                                                                                   >  2:269579‑M1/1‑156 (MQ=255)
                                                                                                 ccccggcattacgtgccgaagtatccagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGTAATACCATTTTCCAGGGCTTCTTTTGCAGAATgg                                                                                                                                              <  1:20656‑M1/106‑1 (MQ=255)
                                                                                                            cgtgccgaagtatcctgttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATctct                                                                                                                                  <  2:35946‑M1/118‑1 (MQ=255)
                                                                                                            cGATCATCGGCATTTGCGAAACTGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagtcgaggcctc                                                                                                                                                                        <  2:49219‑M1/213‑75 (MQ=255)
                                                                                                                    tGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcatcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccac                                                                                                                          >  1:157484‑M1/1‑131 (MQ=255)
                                                                                                                          cagttcttcgggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATATCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATa                                                                                                                     <  1:234742‑M1/131‑1 (MQ=255)
                                                                                                                           agttcttagggcttaccaatatgggccagaattgctgaccagttatcgtgagagtaattcataggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATaa                                                                                                                    <  1:33147‑M1/132‑1 (MQ=255)
                                                                                                                                   gATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgta                                                                                                           <  1:210033‑M1/251‑136 (MQ=255)
                                                                                                                                          ccaatatgggccagaattgctgaccagttatcgtgagagtaattcatcgacacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATCAAATGGTGCCGGAGTCt                                                                                                    <  1:255924‑M1/148‑1 (MQ=255)
                                                                                                                                                     gcACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTTATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccattt                                                                                         <  1:194130‑M1/251‑154 (MQ=255)
                                                                                                                                                            gATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttat                                                                                  >  1:232054‑M1/1‑91 (MQ=255)
                                                                                                                                                                 atatCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGCCAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttaaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacatgtaagctccagatgaccatttttgttattctcc                                                                             >  1:295604‑M1/1‑86 (MQ=255)
                                                                                                                                                                   agttatcgtgagagtaattcatcggaacgttaaatcatatcaggcgtaataccacaaaccttaagttagcgcttatggATATCGCCCAAGCCACATGTTATACATAACCGAACGTTGTAATACAAGGTTTACATCATGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGAAAATCGATCTCTTTACCATGGAATAAAATGGAGCCGGAGTCTTTTTTATAAATAcaaaacaagaatt                                                                           <  2:154635‑M1/173‑7 (MQ=255)
                                                                                                                                                                           ccAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagtt                                                                   <  1:84619‑M1/251‑176 (MQ=255)
                                                                                                                                                                                 cGTCAAAAATCGCTTTTGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttct                                                             <  2:288820‑M1/251‑182 (MQ=255)
                                                                                                                                                                                   attcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATAAAATGGTGCCGGAGTCTTTTTGATAAATACCAAACAGGCATTTTAATAATGTCGAttt                                                           >  1:191880‑M1/63‑251 (MQ=255)
                                                                                                                                                                                    cAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgactggccacctgtaagcaccagatgaccatttttgttattctccacaacgagttagttcttct                                                          >  2:233187‑M1/1‑67 (MQ=255)
                                                                                                                                                                                     tcatcggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAACCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTTCAGCATGGAAATCGATCTCTTTACCCTGGAATAAAATGGTGCCGGAGTCTTTTTGATAAATACCAAACAGGCATTTTAATAATGTCGATTTTc                                                         <  1:160014‑M1/191‑1 (MQ=255)
                                                                                                                                                                                          ggcacgttaaatcatatcaggcgtaataccacaacccttaagttagcgcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATAAAATGGTGCCGGAGTCTTTTTGATAAATACCAAACAGGCATTTTCATAATGTCGATTTTCCTGc                                                     <  1:95622‑M1/195‑1 (MQ=255)
                                                                                                                                                                                               ttGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatcc                                                <  1:229527‑M1/250‑195 (MQ=255)
                                                                                                                                                                                                         atcaggcgtaataccacaacccttaagttagctcttatggATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTTTAATACCAGGTTTAACGCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATAAAATGGTGCCGGAGTCTTTTTGATAAATACCAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCAt                                      >  2:224244‑M1/41‑250 (MQ=255)
                                                                                                                                                                                                                tttGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctggggggga                               <  1:5226‑M1/250‑212 (MQ=255)
                                                                                                                                                                                                                           cGACAACCATGCCTTTGGTGGGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctggggggatactccagcgatg                   >  2:224499‑M1/1‑28 (MQ=255)
                                                                                                                                                                                                                                             tggGATATCGcccataagcgctaacttaagggttgaaccatctgaagaatgcgacgcctcggtgcctcgttaagacgatgcctcgcgttcttcaattgcgttttgtaggctgtcagggatactgtcccacgaatggccacctgtaagctccagatgaccatttttgttattctccacaacgagttagttcttcttttcggatccggcacttctgggggggaaatccagcgatggctggattatgtcgtca  >  2:273339‑M1/1‑10 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
ATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGCTTTGGCGATTTCGATCATCTGCATTTGCGAAACGGATAATGTGCCGACGCGCGCACGCGGATCGATATCAATATCCAGTTCGTCAAAAATCGCTTTGGTTTCGCGGTACATTTTGTCCTGATCGACAAACATGCCTTTGGTGGGATATCGCCCCAGCCACATGTTATCCATCACCGAACGTTGTAATACCAGGTTTAACTCCTGGTGTACCATCGAAATACCATTTTCCAGGGCTTCTTTTGCAGAATGGAAATCGATCTCTTTACCCTGGAATAAAATGGTGCCGGAGTCTTTTTGATAAATACCAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATT  >  NC_000913/2238707‑2239193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ 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.