Predicted mutation
evidence seq id position mutation annotation gene description
MC JC 1 1,010,036 Δ641 bp coding (293‑933/1236 nt) unknown → Type III restriction‑modification system restriction subunit (EC 3.1.21.5)

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ 1 1010036 1010676 641 41 [0] [0] 44 unknown Type III restriction‑modification system restriction subunit (EC 3.1.21.5)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? 1 = 10100350 (0.000)41 (0.490) 28/426 1.8 100% coding (292/1236 nt) unknown Type III restriction‑modification system restriction subunit (EC 3.1.21.5)
?1 1010677 = 0 (0.000)coding (934/1236 nt) unknown Type III restriction‑modification system restriction subunit (EC 3.1.21.5)

GAAAAATATCGCATTGCGTTTGATTAATCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTATGCGTGAAGGCAAATATGAAGAGTTAAAAGCCTTGTGGGAGCTGATTCATCACAAAGCGATTTTGCAGTATAAAATTGATTCCGAAGCCGACTTTTTAGCCTTATTTGTCGTTTATTTGCAGGAAAATGCAGATAAATTCAAACAAGCAGGCATTCGCACAAAAGTGTCGGAAACCTACATAAAAAATAACGTGATGCTGACCCGAGACAAGGAAAATCTAGACGATACGCCGATTCGTTTTAACA  >  1/1009869‑1010281
                                                                                                                                                                      |                                                                                                                                                                                                                                                      
tagtattaccaataaattggcggataatcagctgattgatgataacgatCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatc                                                                                                                                                                    <  1:272779‑M1/227‑85 (MQ=255)
tagtattaccaataaattggcggataatcagctgattgatgataacgatCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatat                                                                                                                                                                     <  1:490172‑M1/226‑84 (MQ=255)
tagtattaccaataaattggcggataatcagctgattgatgataacgatCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcacgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatc                                                                                                                                                                    <  1:371442‑M1/227‑85 (MQ=255)
        ccaataaattggcggataatcagctgattgatgATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaa                                                                                                                                                             >  1:403156‑M1/17‑159 (MQ=255)
             aaattggcggataatcagctgatTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaa                                                                                                                                                       >  1:186513‑M1/12‑154 (MQ=255)
              aattggcggataatcagctgaTTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaac                                                                                                                                                      <  1:735720‑M1/241‑99 (MQ=255)
                      gataaTCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgaga                                                                                                                                              >  1:7521‑M1/3‑145 (MQ=255)
                      gataaTCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaa                                                                                                                                                                                                                                          >  1:595683‑M1/3‑145 (MQ=255)
                                                            cagTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatg                                                                                                          <  1:531151‑M1/249‑143 (MQ=255)
                                                            cagTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggat                                                                                                           <  1:663044‑M1/248‑142 (MQ=255)
                                                                       aTCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataactt                                                                                             >  1:37742‑M1/1‑96 (MQ=255)
                                                                        tCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataactta                                                                                            <  1:901634‑M1/251‑157 (MQ=255)
                                                                             aGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaa                                                                                       >  1:364251‑M1/1‑90 (MQ=255)
                                                                                tGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacgtattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaagg                                                                                    <  1:71499‑M1/251‑165 (MQ=255)
                                                                                             cccGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaa                                                                       >  1:813174‑M1/1‑74 (MQ=255)
                                                                                                     ccTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaaca                                                               >  1:786107‑M1/1‑66 (MQ=255)
                                                                                                                gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtt                                                    >  1:117197‑M1/1‑55 (MQ=255)
                                                                                                                gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtt                                                    >  1:176467‑M1/1‑55 (MQ=255)
                                                                                                                gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaa                                                                                        >  1:185152‑M1/1‑55 (MQ=255)
                                                                                                                gTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggagg                                                                                                                                                                                                   >  1:640665‑M1/1‑55 (MQ=255)
                                                                                                                     gCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatc                                                                    >  1:626861‑M1/1‑50 (MQ=255)
                                                                                                                         aaaGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgccc                                           <  1:733705‑M1/251‑206 (MQ=255)
                                                                                                                            ggCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatat                                                                                                                                                                     >  1:900655‑M1/1‑43 (MQ=255)
                                                                                                                                     aCTGATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaag                               <  1:162789‑M1/251‑218 (MQ=255)
                                                                                                                                        gATGCGAAAACCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaa                                                                                       >  1:620831‑M1/1‑31 (MQ=255)
                                                                                                                                              aaaaCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaac                                                                >  1:629056‑M1/1‑25 (MQ=255)
                                                                                                                                                aaCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagt                     >  1:75313‑M1/1‑23 (MQ=255)
                                                                                                                                                aaCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagt                     >  1:571232‑M1/1‑23 (MQ=255)
                                                                                                                                                 aCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcg                   >  1:605010‑M1/1‑22 (MQ=255)
                                                                                                                                                 aCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcg                   >  1:230783‑M1/1‑22 (MQ=255)
                                                                                                                                                 aCCAAAGGCGATGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcg                   >  1:94995‑M1/1‑22 (MQ=255)
                                                                                                                                                           aTGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcacaagaggtgaaagtcgtgtttaaaa          >  1:858713‑M1/1‑12 (MQ=255)
                                                                                                                                                           aTGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataactt                                                                                             <  1:410335‑M1/167‑156 (MQ=255)
                                                                                                                                                           aTGTTATTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttc                                                                                                      <  1:355201‑M1/158‑147 (MQ=255)
                                                                                                                                                                aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcgtgtttaaaactcagt    >  1:865805‑M1/1‑7 (MQ=255)
                                                                                                                                                                aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtt                                                    <  1:212968‑M1/203‑197 (MQ=255)
                                                                                                                                                                aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgc                                                            <  1:612904‑M1/195‑189 (MQ=255)
                                                                                                                                                                aTTCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaat                                                                                                                                      <  1:896995‑M1/121‑115 (MQ=255)
                                                                                                                                                                 ttCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaagacgaagaacgcaaaatcaaacac                                                              >  1:446118‑M1/1‑6 (MQ=255)
                                                                                                                                                                  tCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcgtgtttaaaactcagttt  >  1:492655‑M1/1‑5 (MQ=255)
                                                                                                                                                                  tCCTAaccgtgtttactaaaatcccgaaaaattctatcaaaatccccgtggcaggaggcgggacttattcccctgattttgcctatatcgtcaaaaccaaaacaggcgagattttgaattttgtgattgaagcaaaaggcgtggatggttcggataacttacgaaaaggcgaagaacgcaaaatcaaacacgccgaacgtttatttgcccaaatctcaaaagaggtgaaagtcgtgtttaaaactcagttt  <  1:485454‑M1/251‑247 (MQ=255)
                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GAAAAATATCGCATTGCGTTTGATTAATCAGCTGATTGATGATAACGATCATTTCGTCAGCAGTGCCGCTTATCAACAGGTGAAACAGCTGTTCCCGAATGCCTTTGCCGCAGTCTTGCATAAAGGCAAAATCACTGATGCGAAAACCAAAGGCGATGTTATTCCTATGCGTGAAGGCAAATATGAAGAGTTAAAAGCCTTGTGGGAGCTGATTCATCACAAAGCGATTTTGCAGTATAAAATTGATTCCGAAGCCGACTTTTTAGCCTTATTTGTCGTTTATTTGCAGGAAAATGCAGATAAATTCAAACAAGCAGGCATTCGCACAAAAGTGTCGGAAACCTACATAAAAAATAACGTGATGCTGACCCGAGACAAGGAAAATCTAGACGATACGCCGATTCGTTTTAACA  >  1/1009869‑1010281

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