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)37 (0.860) 32/482 0.4 97.3% noncoding (1/1331 nt) IS2 repeat region
?NC_000913 = 1093485 1 (0.020)coding (610/1359 nt) dgcT putative diguanylate cyclase DgcT
* ? NC_000913 1093481 =1 (0.020)34 (0.790) 33/486 0.3 97.1% coding (606/1359 nt) dgcT putative diguanylate cyclase DgcT
?NC_000913 = 1651548 NA (NA)noncoding (706/706 nt) IS2 repeat region

GAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTCAGTGCGCTACGAGTAACAAAGAACATTGCACATCGCGA  >  NC_000913/1093238‑1093729
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
cggagtccggagagtaacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGtt                                                                                                                                                                                                                                                     <  1:13082‑M1/6‑1 (MQ=255)
       cggagagtaacggaatagcagcgagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgtgcgagcattataacgaatggcagccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgagaacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTg                                                                                                                                                                                                                                              >  2:243787‑M1/239‑251 (MQ=255)
        ggagagtaacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGc                                                                                                                                                                                                                                             >  1:241857‑M1/238‑251 (MQ=255)
           gagtaacggaatagcagagatcttcgtgaaaacgataaagcgtgactacaccagtatcatgcccaatccagacggtttaacggcagcaaagaaccttgcagatgcgttcgagcagtataacgaattgcatccgcatagtgcgctgggttatcgctcgccacgggaatttctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatagaggggcaaatctaGTGGTTATCCGTGCttt                                                                                                                                                                                                                                          >  2:61843‑M1/235‑251 (MQ=255)
               aacggaatagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGt                                                                                                                                                                                                                                      >  1:135979‑M1/231‑251 (MQ=255)
                     atagcagagagcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcatcaaagaaccttgtagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtggtaacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCt                                                                                                                                                                                                                                >  2:423649‑M1/225‑251 (MQ=255)
                           aTTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaacc                                                                                                                                                                                                                        >  1:407582‑M2/1‑221 (MQ=255)
                               gcttcgtgaaaacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAAt                                                                                                                                                                                                                      >  1:143427‑M1/215‑251 (MQ=255)
                                    tAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGAGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcact                                                                                                                                                                                                               <  1:401237‑M2/251‑40 (MQ=255)
                                      gaaaacgataaagcgtgactacatcagtatcatgcccaaaccatacggtttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTTTACTTTAATGGTTTa                                                                                                                                                                                                                >  1:403838‑M1/208‑250 (MQ=255)
                                      aaGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactg                                                                                                                                                                                                                                               >  2:116569‑M2/1‑210 (MQ=255)
                                      aaGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactg                                                                                                                                                                                                                                               <  1:116569‑M2/217‑8 (MQ=255)
                                         aacgataaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCt                                                                                                                                                                                                            >  1:264159‑M1/205‑251 (MQ=255)
                                            gTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaata                                                                                                                                                                                                         <  2:29322‑M2/249‑46 (MQ=255)
                                             tGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGCCTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcccttaaataagt                                                                                                                                                                                                      <  1:378859‑M2/251‑49 (MQ=255)
                                             tGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagt                                                                                                                                                                                                      <  1:34202‑M2/251‑49 (MQ=255)
                                               aaagcgtgactacatcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGat                                                                                                                                                                                                       >  1:222833‑M1/199‑250 (MQ=255)
                                                  gcgtgactacagcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCt                                                                                                                                                                                                   >  2:444105‑M1/196‑251 (MQ=255)
                                                            gctcaatGTATTTTCCAAAAAAATTGATTTGGGGATAAATTTGTTTTTTGTTTTTGGGGGTCCGATAGTTGCTCACATTCATTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagacgggcaccctgaatcgccagacagccaatatcacttaaataagttatagtcgtaattct                                                                                                                                                                                       >  2:199240‑M2/4‑188 (MQ=255)
                                                             tcagtatcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgctaacagatgtctggcaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTATTCTATACTTTAATGGTTTACGCTAGGATATCTGGAATGGAGTg                                                                                                                                                                                        <  2:24518‑M1/67‑1 (MQ=255)
                                                                  atcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgcttggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatatagggtcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGt                                                                                                                                                                                   >  2:88247‑M1/180‑251 (MQ=255)
                                                                   tcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAg                                                                                                                                                                                          >  2:286528‑M1/179‑243 (MQ=255)
                                                                   tcatgcccaaaccagacgggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAg                                                                                                                                                                                          <  1:286528‑M1/65‑1 (MQ=255)
                                                                     tttttcCAAAAAAATGACTTTGTGCTAAAAGTTGTTTTTGATTTTTGTGGCGCCGATAGTTGCTCACATACTTTCTAGTCATTATGATAGCTACTATTTACACATTGCCGAGTTATCCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATTTTACCATAATGATTTTCATGTGGTtagactggcaccctgaatctccagtcatccaatatcacttaaataagtgattggcttaatgctagttttta                                                                                                                                                                               >  2:398‑M2/2‑179 (MQ=255)
                                                                           cAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGTCAGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggcccccttaatctccagacaaccaatatcacttaaataagtgatagtcttaatgctagtttttagactagt                                                                                                                                                                        >  1:215222‑M2/1‑173 (MQ=255)
                                                                               aaaaTGACTTTGTGCATAAGTTTGTTTTTTGTGTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAAGTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcatt                                                                                                                                                                    >  2:24840‑M2/1‑169 (MQ=255)
                                                                                aaaTGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtaattg                                                                                                                                                                   <  2:288910‑M2/251‑84 (MQ=255)
                                                                                 aaTGACTTTGTGCATAAGTTTGTTTTTTGTTTTTTGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACAGTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagaatggcgccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattgg                                                                                                                                                                  >  2:378526‑M2/1‑167 (MQ=255)
                                                                                    gggttaacggcagcaaagatccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATc                                                                                                                                                                                                                                 >  2:36084‑M1/162‑187 (MQ=255)
                                                                                    gggttaacggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATc                                                                                                                                                                                                                                 <  1:36084‑M1/26‑1 (MQ=255)
                                                                                     gACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattgg                                                                                                                                                                  <  1:148586‑M2/248‑85 (MQ=255)
                                                                                           cggcagcaaagaaccttgcagaggcgttcgagcattataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTAc                                                                                                                                                          <  1:409028‑M1/97‑1 (MQ=255)
                                                                                                aTAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgatt                                                                                                                                                     >  1:259247‑M2/1‑154 (MQ=255)
                                                                                                 aaagTTTGTTTTTTGGTTTTGGGGGGCGGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGACAGGTGGTTTGGAAAGCATAATATGTTACCATAAAGATTTTCATGAGGTtagactggccacctgaatctccagacaacaaatatctcttaaataagtgataatcttaatactagtttttagactagtcattggagaaaagatgattg                                                                                                                                                    <  2:19894‑M2/250‑99 (MQ=255)
                                                                                                     tttgttttttgtttTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgt                                                                                                                                                <  2:91153‑M2/251‑103 (MQ=255)
                                                                                                            gcagaggcgttcgggcattataacgaatggcaaccgcatagtgcgcttgtgtagcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGGTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATTATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGttat                                                                                                                                          >  2:280110‑M1/138‑250 (MQ=255)
                                                                                                                               ataacgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCgtg                                                                                                                      <  1:98772‑M1/133‑1 (MQ=255)
                                                                                                                                  acgaatggcatccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTcc                                                                                                                   >  2:238584‑M1/116‑251 (MQ=255)
                                                                                                                                           tACTTTCTAGTCATTATGCGAGCTACAATTTACACATTGCCGAGTTACCCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgttgtcttagtgccggcgaaacgcagacggcgtaccacactgg                                                                                                          <  2:44946‑M2/251‑141 (MQ=255)
                                                                                                                                           atccgcatagtgcgctgggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGtat                                                                                                          <  1:396733‑M1/145‑1 (MQ=255)
                                                                                                                                                       aTTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggcgaacagatgattgatgtcttatggccggagaaacgcagacggcgtaccacacaggaaaagatcgca                                                                                               <  1:9612‑M2/250‑152 (MQ=255)
                                                                                                                                                       aTTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccaaacaggaaaagatcgca                                                                                               <  2:169523‑M2/250‑152 (MQ=255)
                                                                                                                                                           gggttatcgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTg                                                                                          >  1:453813‑M1/91‑251 (MQ=255)
                                                                                                                                                             agagCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttc                                                                                        >  1:111033‑M2/1‑93 (MQ=255)
                                                                                                                                                              gagCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaatacgtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgccgacggcgtaccacacaggaaaagatcgcaattgttcc                                                                                       <  1:327858‑M2/251‑160 (MQ=255)
                                                                                                                                                               agCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttc                                                                                        <  1:351972‑M2/249‑159 (MQ=255)
                                                                                                                                                                  cgctcgccacgggaatatctgcggcagcgggcttgtaatgggttaagtgataactgatttctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATTCTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATATCGTGTCCACTTGGTATATTAGTCGCACTATTGAagtagt                                                                                   >  1:174889‑M1/84‑251 (MQ=255)
                                                                                                                                                                   aCAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcaga                                                                                  >  2:238448‑M2/1‑87 (MQ=255)
                                                                                                                                                                              aTTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAGGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagcaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccg                                                                       <  2:190278‑M2/251‑176 (MQ=255)
                                                                                                                                                                               aatatctgcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTGTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAAcc                                                                      >  1:92111‑M1/71‑251 (MQ=255)
                                                                                                                                                                               aatatctgcggaagagggcttgtaaggggttaagtgataagagatgtctagaaatataggggcaaatccaGTGGTTAACGTTGCTTTCGGAAAATCTATACTTTAATGGATTACGATATGATATCGGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATAACTCTCTGTTATTTAGGAGCAGGTATAGCGTGTCCACTGTGTATATTAGTCGAACTATTGAAGAAGTCATTAAGTTAAcc                                                                      <  2:314311‑M1/181‑1 (MQ=255)
                                                                                                                                                                                 gCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGGGGTtagactggccccctgaatctccagaccaccaatatcacttacataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttaagcagagctttgacccgggg                                                                    <  2:404822‑M2/251‑179 (MQ=255)
                                                                                                                                                                                   cGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggat                                                                  >  2:351972‑M2/1‑71 (MQ=255)
                                                                                                                                                                                      gcggcagcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTAt                                                                 >  1:381640‑M1/64‑249 (MQ=255)
                                                                                                                                                                                       cggcagagggcttgtaatggggtaagtgataacagatgtctggaaagataggggcaaatccaGTGGTTAACCGTGCTTTCGGTTAATCTATACTTTAATTGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCTCTTGGTATATTAGTCGCGCTATTGAAGTAGTCAGTAAGTTAACCGttattgt                                                              >  2:417760‑M1/63‑248 (MQ=255)
                                                                                                                                                                                        tAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacgg                                                             <  1:49129‑M2/251‑186 (MQ=255)
                                                                                                                                                                                           ggcgggcttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACGTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCCCTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAATTTAACCGTTATGGTAata                                                          >  1:113998‑M1/59‑251 (MQ=255)
                                                                                                                                                                                              aTGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcag                                                                                      >  2:58217‑M2/1‑60 (MQ=255)
                                                                                                                                                                                              aTGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcag                                                                                      <  1:58217‑M2/220‑161 (MQ=255)
                                                                                                                                                                                                 cttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGtt                                                                                                                                                                                                                  >  1:375251‑M1/53‑93 (MQ=255)
                                                                                                                                                                                                 cttgtaatgggttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGtt                                                                                                                                                                                                                  <  2:375251‑M1/41‑1 (MQ=255)
                                                                                                                                                                                                           gttaagtgataacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCatatat                                          >  1:195991‑M1/43‑251 (MQ=255)
                                                                                                                                                                                                             gggttTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaagcgccagccaaccaatatcacttaaataagtgatagtcttaatactagtttttagactcgtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccttcgttgcctggcaac                                        >  2:34202‑M2/2‑45 (MQ=255)
                                                                                                                                                                                                                 tgataacagatgtctggacatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCATACTGTAAAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTCAGTg                                    <  2:163658‑M1/215‑1 (MQ=255)
                                                                                                                                                                                                                   ggAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtg                                  <  1:311108‑M2/251‑213 (MQ=255)
                                                                                                                                                                                                                   ggAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTGACCGTTATGGTAATATTTATGTGCCATATATTCAGTgc                                   <  1:161030/250‑1 (MQ=255)
                                                                                                                                                                                                                     aacagatgtctggaaatataggggcaaatccaGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTCAGTGCGCt                                >  1:19559‑M1/33‑251 (MQ=255)
                                                                                                                                                                                                                     aaaGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacgtcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacctggtgt                                 <  1:234110‑M2/250‑214 (MQ=255)
                                                                                                                                                                                                                      aaGCCTCATATGTTACCATAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgt                                 <  2:202047‑M2/249‑214 (MQ=255)
                                                                                                                                                                                                                                       aTAATGATTTTCATGTGGTtagactggccccctgaatctccctacagccaatatcacttaaataagtgacagtcttaatactagtttttagacgaggcattggagaacagatgattgatgtcttagggccggggaaacgcagacggcgtaccacacaggaaaagatcgcatcggttcagca                                                                                    >  2:51337‑M2/1‑19 (MQ=255)
                                                                                                                                                                                                                                       aTAATGATTTTCATGTGGTtagactggccccctgaatctccagacaaccaatatcacttaaataagtgacagtcttaatactagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagca                                                                                    <  1:51337‑M2/181‑163 (MQ=255)
                                                                                                                                                                                                                                                agccatgggTTAACCTTGCTTTCGGTTAATCTGTTCTTTAATGGTTTACGCTATGAGATCTGGAATGGAGTGGCGGTTATTGGTTTTTGTGACGGACTGTACAATATCTCTCTGTTATTTATGAGCAATTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTGAAGTTAATATTTATGTGCCATATATACAGTGCGCTACGAGTAACAATGAACATTGcacagcg     <  2:381276/244‑4 (MQ=255)
                                                                                                                                                                                                                                                   ggggggttagactgtccccctgaatctccagacaaccaatatcacttaaataagtgatagtcttaattctagtttttagactagtcattggagaacagatgattgatgtcttagggccggagaaacgcagacggcgtaccacacaggaaaagatcgcaattgttcagcagagctttgaaccggggatgacggtctccctcgttgcccggcaacatggtgtagcagccagccagttatttctctggcgtaag  >  2:154580‑M2/5‑7 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
GAGTTTGTGCATATTAATATTTCTGGCATTGGTGAATAAAGTTAGTGAAAACACAAAACAGCGCAATTTATTTTCCAAAAAAATGACTTTGTGCATAAGTTTGTTTTTTGTTTTTGGGGGGCCGATAGTTGCTCACATACTTTCTAGTCATTATGAGAGCTACAATTTACACATTGCCGAGTTAACCAATGAAAATGGTCAGGTGGTTTGGAAAGCCTCATATGTTACCATAATGATTTTCATGTGGTTAACCTTGCTTTCGGTTAATCTATACTTTAATGGTTTACGCTATGATATCTGGAATGGAGTGACGGTGATTGCTTTTTGTGCCGTACTGTACAATATCTCTCTGTTATTTATGAGCAGGTATAGCGTGTCCACTTGGTATATTAGTCGCACTATTGAAGTAGTCAGTAAGTTAACCGTTATGGTAATATTTATGTGCCATATATTCAGTGCGCTACGAGTAACAAAGAACATTGCACATCGCGA  >  NC_000913/1093238‑1093729

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.