Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 1,093,481 IS2 (+) +5 bp coding (606‑610/1359 nt) dgcT → putative diguanylate cyclase DgcT

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 381260 =NA (NA)28 (0.810) 26/486 0.4 100% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 1093485 0 (0.000)coding (610/1359 nt) dgcT putative diguanylate cyclase DgcT
* ? NC_000913 1093481 =1 (0.030)20 (0.570) 20/490 0.7 95.1% coding (606/1359 nt) dgcT putative diguanylate cyclase DgcT
?NC_000913 = 1651548 NA (NA)noncoding (706/706 nt) IS2 repeat region

ATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTCAGTGCGCTACGAGTAACAAAGAACATTGCACATCGCGAT  >  NC_000913/1093236‑1093730
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
aTGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGtt                                                                                                                                                                                                                                                      <  1:29692/251‑1 (MQ=255)
   agtccggagagtaacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAAc                                                                                                                                                                                                                                                   >  1:138295‑M1/243‑251 (MQ=255)
    gtccggagagtaacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAAcc                                                                                                                                                                                                                                                  >  2:79625‑M1/242‑251 (MQ=255)
            agtaacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTc                                                                                                                                                                                                                                          >  1:269313‑M1/234‑251 (MQ=255)
               ttAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctg                                                                                                                                                                                                                                        <  1:352839‑M2/250‑16 (MQ=255)
                 aaTATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaat                                                                                                                                                                                                                                     <  2:79234‑M2/251‑19 (MQ=255)
                           gagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACtt                                                                                                                                                                                                                           >  2:179922‑M1/219‑251 (MQ=255)
                              agcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTaa                                                                                                                                                                                                                        >  2:147534‑M1/216‑251 (MQ=255)
                                     aTAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTCCCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcac                                                                                                                                                                                                                 >  2:232168‑M2/1‑213 (MQ=255)
                                       aaaGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcactt                                                                                                                                                                                                               >  1:120140‑M2/1‑211 (MQ=255)
                                       aaaGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcactt                                                                                                                                                                                                               >  2:298776‑M2/1‑211 (MQ=255)
                                           ttAGTGAAAACACACAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCAGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactgcccccctgaatctcaagacacccaatatcacttaaat                                                                                                                                                                                                           <  2:343018‑M2/251‑45 (MQ=255)
                                            gataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATg                                                                                                                                                                                                          <  2:78147‑M1/50‑1 (MQ=255)
                                                 agcgtgactacagcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAAGCTATACTTTAATGGTTTACGCTATGATATc                                                                                                                                                                                                     >  2:169705‑M1/197‑251 (MQ=255)
                                                                gcAATTTATTTTCCAAAAAAATGACTTTGGGCATATGTTTGTTTTTTGATTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATACTGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatact                                                                                                                                                                                        <  1:311909‑M2/249‑64 (MQ=255)
                                                                  atcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACgg                                                                                                                                                                                     >  1:168165‑M1/180‑250 (MQ=255)
                                                                         ccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCt                                                                                                                                                                             <  1:96279‑M1/79‑1 (MQ=255)
                                                                                                              tgtttttgGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagaaaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggcc                                                                                                                                        <  2:359606‑M2/251‑112 (MQ=255)
                                                                                                                 ttttgGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccgg                                                                                                                                      <  1:190087‑M2/250‑114 (MQ=255)
                                                                                                                   cgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGc                                                                                                                                   >  1:302375‑M1/131‑251 (MQ=255)
                                                                                                                       cgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAgg                                                                                                                                >  2:42941‑M1/127‑250 (MQ=255)
                                                                                                                       cgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGt                                                                                                                               >  2:136958‑M1/127‑251 (MQ=255)
                                                                                                                         agcattataacgaatggtatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGtat                                                                                                                             >  1:81502‑M1/125‑251 (MQ=255)
                                                                                                                          gcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGtat                                                                                                                             <  2:215366‑M1/127‑1 (MQ=255)
                                                                                                                           ccGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgc                                                                                                                              <  1:89991‑M2/248‑122 (MQ=255)
                                                                                                                           ccGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTGCATGTGGTtagactggccccctgaatctccagacaaccaatatcccttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcaga                                                                                                                           >  2:8696‑M2/1‑127 (MQ=255)
                                                                                                                             gATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAACGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacg                                                                                                                         <  2:41778‑M2/251‑127 (MQ=255)
                                                                                                                                aGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcg                                                                                                                      <  2:48674‑M2/251‑130 (MQ=255)
                                                                                                                                 gTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccgtagaaacgcagacggcgt                                                                                                                     >  2:151288‑M2/1‑121 (MQ=255)
                                                                                                                                           tACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacagg                                                                                                           <  2:24320‑M2/251‑141 (MQ=255)
                                                                                                                                                            gagaGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgtt                                                                                          <  1:324975‑M2/251‑158 (MQ=255)
                                                                                                                                                             gttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGaa                                                                                         >  2:338261‑M1/89‑251 (MQ=255)
                                                                                                                                                                                   cGAGTTAACCAATGACAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggat                                                                   <  1:227197‑M2/251‑181 (MQ=255)
                                                                                                                                                                                     aGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatga                                                                 >  1:148986‑M2/1‑69 (MQ=255)
                                                                                                                                                                                      gTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgac                                                                <  2:193823‑M2/251‑184 (MQ=255)
                                                                                                                                                                                         aaCCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCAGATTTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcagtggataacagatgattgatgtcttagggccggagaaacgcatacggcgttccacacatgaaaagatcgcaatttttcagcagatctttgaaccggggatgacgtt                                                             >  2:305745‑M2/1‑65 (MQ=255)
                                                                                                                                                                                                              aagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTc                                         >  1:43110‑M1/40‑250 (MQ=255)
                                                                                                                                                                                                                        agatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAAGATTTATGTGCCATATATTCAGTGCGCTACg                              >  1:320589‑M1/30‑251 (MQ=255)
                                                                                                                                                                                                                          atgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGttatttat                                                                                                                                       <  1:140180‑M1/117‑1 (MQ=255)
                                                                                                                                                                                                                          atgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGttatttat                                                                                                                                       >  2:140180‑M1/28‑144 (MQ=255)
                                                                                                                                                                                                                            cATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaagtgttcagcagagctgtgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagcc                           >  1:41778‑M2/1‑30 (MQ=255)
                                                                                                                                                                                                                                 gTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagccagcc                       <  1:624‑M2/249‑225 (MQ=255)
                                                                                                                                                                                                                                      cATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagccagccagttatt                <  1:64245‑M2/251‑232 (MQ=255)
                                                                                                                                                                                                                                            gATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagccagccaggtatttctct           >  2:54506‑M2/1‑14 (MQ=255)
                                                                                                                                                                                                                                              aaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTCAGTGCGCTACGAGTAACAAAGAACATTGCACAt        <  2:86888‑M1/244‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTGTGGTtagactggccccctgactctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagccagccagttatttctctggcgta     <  2:216283‑M2/249‑243 (MQ=255)
                                                                                                                                                                                                                                                    tgtgGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagccagccagttatttctctggcgtaagc  <  1:70557‑M2/251‑246 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         
ATGAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTCAGTGCGCTACGAGTAACAAAGAACATTGCACATCGCGAT  >  NC_000913/1093236‑1093730

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