Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,066,152 Δ12,778 bp insH6[yoeF] insH6, yoeA, insD3, insC3, yoeA, isrC, flu, yeeR, yeeS, yeeT, cbeA, cbtA, yeeW, [yoeF]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2066152 2078929 12778 49 [0] [0] 49 insH6–[yoeF] insH6,yoeA,insD3,insC3,yoeA,isrC,flu,yeeR,yeeS,yeeT,cbeA,cbtA,yeeW,[yoeF]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 20661510 (0.000)49 (0.620) 44/492 0.4 100% intergenic (‑387/+154) cobU/insH6 cobinamide‑P guanylyltransferase/cobinamide kinase/IS5 family transposase and trans‑activator
?NC_000913 2078930 = 0 (0.000)coding (356/357 nt) yoeF putative uncharacterized protein YoeF

GGCAACAAGCAGGAATCTCAGCGCTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGC  >  NC_000913/2078687‑2079174
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
gatggcattacggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTc                                                                                                                                                                                                                                                 >  2:83070‑M2/244‑249 (MQ=255)
gatggcattacggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTc                                                                                                                                                                                                                                                 <  1:83070‑M2/6‑1 (MQ=255)
gatggcattacggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCg                                                                                                                                                                                                                                                >  1:77884‑M2/244‑250 (MQ=255)
 atggcattacggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGAt                                                                                                                                                                                                                                              >  2:100215‑M2/243‑251 (MQ=255)
           ggaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCgg                                                                                                                                                                                                                                    >  1:627899‑M2/233‑251 (MQ=255)
            gaatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtaagcgccagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgcctatAATTTCGATTGGTTCTCgg                                                                                                                                                                                                                                    >  2:375992‑M2/232‑250 (MQ=255)
             aatagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGAt                                                                                                                                                                                                                                  <  1:6820‑M2/21‑1 (MQ=255)
              atagccagactttctcgtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATg                                                                                                                                                                                                                                 <  1:100215‑M2/22‑1 (MQ=255)
                              gtgttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTg                                                                                                                                                                                                                  >  2:438923‑M2/214‑250 (MQ=255)
                                gttaacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcatgcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATa                                                                                                                                                                                                               <  2:108932‑M2/40‑1 (MQ=255)
                                  taacgtccgacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATac                                                                                                                                                                                                              >  1:370329‑M2/210‑250 (MQ=255)
                                      gtcagacaaaatgcgtcattgacagcggtttctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATacac                                                                                                                                                                                                            >  2:84214‑M2/206‑248 (MQ=255)
                                                   cgtcattgacagcggttgctgtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatcggcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATcaatga                                                                                                                                                                                            <  2:481112‑M2/59‑1 (MQ=255)
                                                                       gtggattgacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTGGCCggg                                                                                                                                                                          <  2:235027‑M2/77‑1 (MQ=255)
                                                                              gacacatttcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTc                                                                                                                                                                  <  2:474023‑M2/85‑1 (MQ=255)
                                                                                      tcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGgtta                                                                                                                                                                                                                                        <  2:356762‑M2/15‑2 (MQ=255)
                                                                                      tcattcagataatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTtc                                                                                                                                                                                                                                        >  1:356762‑M2/158‑172 (MQ=255)
                                                                                             gataatgaattaatgctactgacggaacaaccagcaaaagagcattgtgttctgaaatccaggctattgattcaacgtcagcgccagatatgcgctgtaaaagaacttgtacaacaacgtaatagaataatccggaatggtggcgactatAATTTCGATTGCTTCTCGGATGAGGATGATGACGGCGATACACAGATAAGAATCAATGCGCAATGAGGGTTGCCGGGCAACACTCATTGAATAAAACggg                                                                                                                                                   >  1:86608‑M2/151‑250 (MQ=255)
                                                                                                aatgaattaatgctactgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAAtt                                                                                                                                               <  2:260758‑M2/104‑1 (MQ=255)
                                                                                                            ctactgccggaacaaccagcgaaagggaattgtgttcggaaatccaggctattgaatcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatacggaatggtggagactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACttcgcttcgc                                                                                                                                   >  2:220915‑M2/136‑251 (MQ=255)
                                                                                                              actgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGcc                                                                                                                                  <  1:84214‑M2/117‑1 (MQ=255)
                                                                                                              actgccggaacaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCg                                                                                                                                 >  1:484769‑M2/134‑251 (MQ=255)
                                                                                                                    ggaacaaccatcaaacgggcatcgtgttcttaaatccaggctattgattcaacgtcagagacagatctgcgctgtaaaagcacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGCCGGTGATACACAGATAAGAATCAATGAGAAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCCCTTTGCCGTTTTTc                                                                                                                           >  2:29652‑M2/128‑251 (MQ=255)
                                                                                                                       acaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCAtt                                                                                                                        <  1:532939‑M2/127‑1 (MQ=255)
                                                                                                                       acaaccagcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCAtt                                                                                                                        >  1:424639‑M2/125‑251 (MQ=255)
                                                                                                                             agcaaacgggcattgtgtgctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtcaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAAGTAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGccc                                                                                                                  <  2:241495‑M2/133‑1 (MQ=255)
                                                                                                                              gcaaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagctctgcgctgtaatagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCCACCCTCCTTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGccc                                                                                                                  >  1:451277‑M2/118‑250 (MQ=255)
                                                                                                                                aaacgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCAt                                                                                                                >  2:324258‑M2/116‑250 (MQ=255)
                                                                                                                                  acgggcattgtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGc                                                                                                             >  2:286599‑M2/114‑251 (MQ=255)
                                                                                                                                           gtgttctgaaatccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCt                                                                                                     <  2:484769‑M2/146‑1 (MQ=255)
                                                                                                                                                     atccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtcatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGCGat                                                                                           >  2:418183‑M2/95‑250 (MQ=255)
                                                                                                                                                      tccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGatat                                                                                         >  1:476477‑M2/94‑251 (MQ=255)
                                                                                                                                                       ccaggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtataacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATc                                                                                        <  2:442329‑M2/159‑1 (MQ=255)
                                                                                                                                                        caggctattgattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATcg                                                                                       >  1:295887‑M2/92‑251 (MQ=255)
                                                                                                                                                                 gattcaacgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTTCTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTc                                                                               <  1:462341‑M2/168‑1 (MQ=255)
                                                                                                                                                                        cgtcagcgacagatctgcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGgaga                                                                       <  1:11503‑M2/176‑1 (MQ=255)
                                                                                                                                                                                        gcgctgtaaaagaacttgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCAtttt                                                           <  1:464636‑M2/188‑1 (MQ=255)
                                                                                                                                                                                                     actgttaaaacaacgtaatagaattatccggaatggtggcgacgatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCATTTTTTTGGAGATATCGCGGCGTCCTTTGGAGAGCTCGGCATTCTTGATGATGTAGTCGTCa                                          <  2:461993‑M2/205‑1 (MQ=255)
                                                                                                                                                                                                        tgtaaaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGTTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCaa                                         <  2:504650‑M2/206‑1 (MQ=255)
                                                                                                                                                                                                            aaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGAt                                   <  2:34424‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                            aaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGAt                                   >  1:269297‑M2/40‑251 (MQ=255)
                                                                                                                                                                                                            aaacaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGAt                                   <  1:233832‑M2/212‑1 (MQ=255)
                                                                                                                                                                                                              acaacgtaatagaattatccggaatggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCt                                 >  1:608335‑M2/38‑251 (MQ=255)
                                                                                                                                                                                                                                ccggaatggtggcgactatAATTTCGAGTGGTGCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGTGTTGCCGGGCAACCCTTATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAAGCGTCTACGCGATCTTCATAGTCACCTTTCAt                <  1:477277‑M2/231‑1 (MQ=255)
                                                                                                                                                                                                                                     atggtggcgactatAATTTCGATTGGTTCTCGGATGAGTATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGAGATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACTCGATCTTCATAGTCGCCTTTCATGCTGGc          <  1:492451‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                     atggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGc          <  1:492469‑M2/237‑1 (MQ=255)
                                                                                                                                                                                                                                       ggtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCg         >  1:408272‑M2/13‑250 (MQ=255)
                                                                                                                                                                                                                                        gtggcgactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCgat       >  2:148439‑M2/12‑251 (MQ=255)
                                                                                                                                                                                                                                             gactatAATTTCGATTGGTTCTCGGATGAGGATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGc  >  1:213491‑M2/7‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
GGCAACAAGCAGGAATCTCAGCGCTGAAGCAGGAATTAGAGGTAAAAACGCCCTACCGGGCCATGAATCACCCGGTCATTGGAGTAGTAACCAAAGCAGATTTAGCCAGCATGGAACAAATATCTTTGGTGAAAAGTTGGTTACGGGAAGCTGGAGCGCACAACGTGTTAGTAACCAGTGCAGTTAATAACAATGGAGTTACAGAACTTTTTGCCTTGCTGCATACAGAAGAAGGCTGTTGTTAATTTCGATTGGTTCTCGGATGAGAATGATGACGGTGATACACAGATAAGAATCAATGAGCAATGAGGGTTGCCGGGCAACCCTCATTGAATAAAACGGGAATTACTTCGCTTCGCCGTTTTTCATTTCGCCCATAGCTTTCAGCTTTTTGGAGATATCGCGGCGTTCTTTGGAGAGCTCGGCATTTTTGATGATGTAATCGTCAACGCGATCTTCATAGTCACCTTTCATGCTGGCGATGATGC  >  NC_000913/2078687‑2079174

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