Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,492,211 IS150 (–) +3 bp intergenic (+36/‑257) sokB → / → trg putative small regulatory RNA antitoxin SokB/methyl‑accepting chemotaxis protein Trg

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1492211 =0 (0.000)25 (0.440) 24/488 0.9 100% intergenic (+36/‑259) sokB/trg putative small regulatory RNA antitoxin SokB/methyl‑accepting chemotaxis protein Trg
?NC_000913 3720633 = NA (NA)noncoding (1/1443 nt) IS150 repeat region
* ? NC_000913 = 14922131 (0.020)21 (0.360) 20/492 1.1 95.4% intergenic (+38/‑257) sokB/trg putative small regulatory RNA antitoxin SokB/methyl‑accepting chemotaxis protein Trg
?NC_000913 = 3722075 NA (NA)noncoding (1443/1443 nt) IS150 repeat region

TGTCCAGAAAAGGCAATTGTCTGGATAATGCTGTGGTGGAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGG  >  NC_000913/3721829‑3722322
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
tGTCCAGAAAAGGCAATTGTCTGGATAATGCTGTGGTGGAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                                                                                                                                                                                                                                     >  1:321283/1‑251 (MQ=255)
    cAGAAAAGGCAATTGTCTGGATAATGCTGTGGTGGAGTGATTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTg                                                                                                                                                                                                                                                  >  2:134275/1‑250 (MQ=255)
           ggCAATTGTCTGGATAATGCTGTGGTGGAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTtat                                                                                                                                                                                                                                            >  2:407692/1‑249 (MQ=255)
             cAATTGTCTGGATAATGCTGTGGTGGAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                                                                                                                                                                                                                                        >  1:514396/1‑251 (MQ=255)
               aTTGTCTGGATAATGCTGTGGTGGAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgtttt                                                                                                                                                                                                                                      >  1:329419‑M2/1‑232 (MQ=255)
                      ggATAATGCTGTGGTGGAGTGTTTCTGTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAGTTATAGTCAGGTggg                                                                                                                                                                                                                                 <  1:85168/249‑1 (MQ=255)
                          aaTGCTGTGGTGGAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccg                                                                                                                                                                                                                           <  1:232340‑M2/251‑31 (MQ=255)
                                   gtggAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcaga                                                                                                                                                                                                                  >  1:256054‑M2/1‑212 (MQ=255)
                                            gtCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACATACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTTTCTGCCGTTCgg                                                                                                                                                                                                          <  2:563781/249‑1 (MQ=255)
                                                ttGGAACCGTAAAGGCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggcc                                                                                                                                                                                                      >  2:492949‑M2/1‑199 (MQ=255)
                                                                          tATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgagg                                                                                                                                                                            >  1:500087‑M2/1‑173 (MQ=255)
                                                                           aTCTTGATGAGTTCAGTAATATAAGCGAACTTAAGGATGCTGTTACTGAATATATTGAATACTACAACATCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCTAATTGAATATCGTAATCAGACCTATATGCATCGTGGTTAACTGTCCACCTTTTTGGGGTTAGTACAAACTTTGATTTATAGGCAGGGGGTGCTTTGCTCTGTCGGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCAcc                                                                                                                                                                           <  2:106795/250‑1 (MQ=255)
                                                                                gatgaGTTCAGTAATATAAGCGAACTGAAGGAAGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggctccgaccgatatttctatcaaccgaggccaacg                                                                                                                                                                      >  1:119883‑M2/1‑167 (MQ=255)
                                                                                gatgaCTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGCAGCATTAGCCTGAAATTAAAAGGACTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGCTTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGCGACAAACAGTCTCAAGCACCCGTGGc                                                                                                                                                                     >  2:347876/1‑251 (MQ=255)
                                                                                 atgaGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacga                                                                                                                                                                     <  1:384493‑M2/250‑85 (MQ=255)
                                                                                   gaGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacgaatg                                                                                                                                                                  <  2:309082‑M2/251‑88 (MQ=255)
                                                                                            aaTATAAGCGAACTGAAGGGTGCTGTTACGGAATATATTGAATACTACAACAGAAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGCACAgcctgagactttttgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatattgttatcaaccgaggccaacgaatgaacatcga                                                                                                                                                          <  2:487851‑M2/250‑96 (MQ=255)
                                                                                                   gCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacgaatgaacctagcaattc                                                                                                                                                     <  1:595855‑M2/248‑101 (MQ=255)
                                                                                                         tGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccg                                                                                                                                                                                                                           <  1:183424‑M2/172‑31 (MQ=255)
                                                                                                         tGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccg                                                                                                                                                                                                                           >  2:183424‑M2/1‑142 (MQ=255)
                                                                                                          gAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAgtatgttt                                                                                                                                             <  1:134275/249‑1 (MQ=255)
                                                                                                           aaGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTttctt                                                                                                                                          <  1:287644/251‑1 (MQ=255)
                                                                                                           aaGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGGCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTttctt                                                                                                                                          <  1:41442/251‑1 (MQ=255)
                                                                                                                  cTGTTAAGGAATATATTGAATACTCCAACAGCAGATGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATAACTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTg                                                                                                                                   >  2:65776/1‑251 (MQ=255)
                                                                                                                      tACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacgaatgaacctagcaattcagatcttagcctcttacccgcc                                                                                                                               <  1:526404‑M2/251‑123 (MQ=255)
                                                                                                                        cGGAATATATTGAATACTACAACAGCAGAAGACTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGACTATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGAtata                                                                                                                              >  2:72457/1‑250 (MQ=255)
                                                                                                                              atatTGAATACTAAAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTACCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCGGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTTATATAATccga                                                                                                                        <  2:321283/250‑3 (MQ=255)
                                                                                                                                       actacAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgac                                                                                                                                                                                                 >  1:480225‑M2/1‑112 (MQ=255)
                                                                                                                                                    aagaaTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTATCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGAGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACTAACAGTCTCAAGCACCCGTGGCTATTCTAACTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAatat                                                                                                 <  2:183644/251‑1 (MQ=255)
                                                                                                                                                       aaTTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatgtctatcaaccgaggccaacgaatgaacctagcaattcagatcttagcctcttacccgccatccggcaaggagaaaggctatgaagcacaacc                                                                                              >  2:303767‑M2/1‑96 (MQ=255)
                                                                                                                                                          tagccgcaaattaaaagGTCTGAATCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCATTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagtctatggccccgaccgatatttcga                                                                                                                                                                                      >  2:25753‑M2/8‑93 (MQ=255)
                                                                                                                                                          tAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttcta                                                                                                                                                                                      <  1:25753‑M2/160‑68 (MQ=255)
                                                                                                                                                                          ggTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacgaatgaacctagcaattcagatcttagcctcttacccgccatccggcaaggagaaaggctatgaagcacaaccctccggtggtgtgtctgct                                                                           >  2:139046‑M2/1‑77 (MQ=255)
                                                                                                                                                                                       ttGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacgaatgaacctagcaattcagatcttagcctcttacccgccatccggcaaggagaaaggctatgaagcacaaccctctggtggtgtgtctgctcattatctgcatt                                                              <  1:568796‑M2/251‑188 (MQ=255)
                                                                                                                                                                                        agAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAATTCTCATAAAATTAATATAAGAGATAATCTTTATTCAGCAGAAGATTATTAAAg                                                             <  2:555874/250‑1 (MQ=255)
                                                                                                                                                                                         gAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAg                                                             <  1:511540/250‑1 (MQ=255)
                                                                                                                                                                                             atCGGAGTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTg                                                         <  1:113864/250‑1 (MQ=255)
                                                                                                                                                                                             atCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATTCCTGAGACAAACAGTCTCAAGCACCCGGGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTg                                                         <  1:528913/250‑1 (MQ=255)
                                                                                                                                                                                                ggAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGACAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGTTTGCTGt                                                     >  2:119694/1‑251 (MQ=255)
                                                                                                                                                                                                 gAATCAGACCTATATTCCGCGTGTTGAACTTTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacgaatgaacctagcaattcagatcttagcctcttacccgccatccggcaaggagaaaggctatgaagcacaaccctctggtggtgtgtctgctcattatcttcattacgattctga                                                    >  2:568796‑M2/1‑54 (MQ=255)
                                                                                                                                                                                                      agatcTATATGCCGCGTGTTTAACTGTCGAACTTTTTGGTGTCTGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTTTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAAACCTAAAAAAATCTCATAAAATTAATTTCTGAGATAATCTTTATTCACCCGAAGACTCTTCAAGTTGGCTGTattatt                                               <  2:19809/247‑1 (MQ=255)
                                                                                                                                                                                                         ccTATATGCCTCGTGTTTAACTTTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCGTTAGTCAGCAGAAGATTATTAAAGGTTGCTGTAGTATTTAg                                            <  2:583478/251‑1 (MQ=255)
                                                                                                                                                                                                                       gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGcc                              <  1:435968/251‑1 (MQ=255)
                                                                                                                                                                                                                          tAACTGTCCAACTTTTTGGGGTCAGTACAgcctgagactttctgttttcaggctaaccgggaatcagaaagacgaaggccccgaccgatatttctatcaaccgaggccaacgaatgaacctagcaattcagatcttagcctcttacccgccatccggcaaggagaaaggctatgaagcacaaccctctggtggtgtgtctgctcattatctgcattacgattctgacattcacactcctgacccgacaaac                           <  1:472869‑M2/251‑223 (MQ=255)
                                                                                                                                                                                                                                                   tACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGACAGGCTCTTTAATAACgg  >  2:222843/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
TGTCCAGAAAAGGCAATTGTCTGGATAATGCTGTGGTGGAGTGTTTCTTTGGAACCTTAAAGTCGGAGTGTTTTTATCTTGATGAGTTCAGTAATATAAGCGAACTGAAGGATGCTGTTACGGAATATATTGAATACTACAACAGCAGAAGAATTAGCCTGAAATTAAAAGGTCTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGCACCCGTGGCTATTCTAGCTTAATAAGTTTGTTTCTTCTCCTTGATATAATCCTAAAAAAATCTCATAAAATTAATATATGAGATAATCTTTATTCAGCAGAAGATTATTAAAGGTTGCTGTATTATTTAGCGATAAAAAAAGCCTGCCAGATGGCAGGCTATTTAATAACGG  >  NC_000913/3721829‑3722322

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