Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,620,828 43 bp→AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT intergenic (+2/‑13) deoB → / → deoD phosphopentomutase/purine nucleoside phosphorylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4620828 4620870 43 28 [0] [0] 41 deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46208270 (0.000)22 (0.480)
+AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT
17/358 1.1 100% intergenic (+1/‑56) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase
?NC_000913 4620871 = 0 (0.000)intergenic (+45/‑12) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

GGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGA  >  NC_000913/4620710‑4621121
                                                                                                                                                                 |                                                                                                                                                                                                                                                          
agttggagtcattaccaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCtgcgtgcg                                                                                                                                                                   <  2:158780‑M2/90‑1 (MQ=255)
    ggagtcattaccaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggggggggtgaaaacctgttgatAAGGATAAAACAATGGCTACACCACACATTAATGCAGGAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGt                                                                                                                                                               >  2:601573‑M2/158‑251 (MQ=255)
      agtcattaccaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccccctaccccagggggagggggaaaacatgttgatAAGGATAAAACAAAGGCTACCCCACACATTAATGCAGAAATGGGCGGTTTCGCTGACGTAGGTTTTATTCCAAGCGACCCGCTGCGTGCGAAGtat                                                                                                                                                             >  1:567934‑M2/156‑251 (MQ=255)
         cattaccaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccccagggggagggggaaaacaagttgatAAGGATAAAACAATGGCTACCCCACCCATTAATGCAGAAAAGGGCGATTTCGCTGACGTAGTTTTGCTGCCAGGCGACCCCCTGCGGGCGAAGGATAtt                                                                                                                                                          >  1:183732‑M2/153‑251 (MQ=255)
          attaccaccttcgcggatagtggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTGCCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTg                                                                                                                                                         <  2:314312‑M2/100‑1 (MQ=255)
              ccaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGa                                                                                                                                                     <  2:591464‑M2/104‑1 (MQ=255)
              ccaccttcgcggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGa                                                                                                                                                     <  1:462901‑M2/104‑1 (MQ=255)
                     cgcggatattggtcagactctggcaaaatattttggtacttctgatagggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACGATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGGTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACttt                                                                                                                                               <  2:432537‑M2/110‑1 (MQ=255)
                        ggatattggtcagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCtt                                                                                                                                           <  2:251231‑M2/114‑1 (MQ=255)
                                  cagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccccagggagagggggaaaacatgttgatAAGGATAAAACAAATGCTACCCCCCACATTAATGCAGAAAAGGGCGATTTCGCTGACGTAGTTTTTAtgcccg                                                                                                                                                                                    >  1:355154‑M2/128‑198 (MQ=255)
                                  cagactctggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAg                                                                                                                                                                                    <  2:355154‑M2/73‑1 (MQ=255)
                                  cagactctggcaaaatattttggtacgtctgatatggaatatggcaaagccatgttctgataatgggtgcggattgatgccctcacaccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACAGTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTGTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGccc                                                                                                                                 <  2:464069‑M2/124‑1 (MQ=255)
                                          ggcaaaatattttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCgtgaagt                                                                                                                          <  1:57547‑M2/131‑1 (MQ=255)
                                                  attttggtacttctgatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACg                                                                                                                  <  1:400342‑M2/139‑1 (MQ=255)
                                                                gatatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGAGGCCCGTGAAGTGAACAACGTTCGCGGTATGCTgg                                                                                                   <  2:105492‑M2/154‑1 (MQ=255)
                                                                 atatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTggg                                                                                                  <  1:510313‑M2/155‑1 (MQ=255)
                                                                 atatggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTggg                                                                                                  <  1:420090‑M2/155‑1 (MQ=255)
                                                                     ggaatatggcaaagccatgttctgataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCCACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGGAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGCAGTGAACCACGTTCGCGGTATTCTGGGcttt                                                                                              >  2:151755‑M2/93‑250 (MQ=255)
                                                                         tatggcaaagccatgttctgataatgggtgcggcttgatggcctcaccccgggcctctcccacagggagagggggaaaacaagttgatAAGGATAAAACAATGGCTACCCCCCACATTAAAGCCGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTTCTTAAAATTTCCTTGAAGATGCCCGTGAAGTGGACAACGTTCGCGGGATGCTGGGCTTCaccc                                                                                          >  1:571319‑M2/89‑250 (MQ=255)
                                                                             gcaaagccatgttctgataatgggtgcggctttatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTAc                                                                                      <  1:377701‑M2/167‑1 (MQ=255)
                                                                                            gataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCaa                                                                       <  1:435167‑M2/182‑1 (MQ=255)
                                                                                             ataatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCaaa                                                                      <  1:601573‑M2/183‑1 (MQ=255)
                                                                                                   ggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTcc                                                                <  1:560596‑M2/189‑1 (MQ=255)
                                                                                                        ggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTa                                                                                   <  1:464973‑M2/170‑1 (MQ=255)
                                                                                                             gatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGGAATGGGTc                                                      <  1:204477‑M2/199‑1 (MQ=255)
                                                                                                                gccctcaccccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCTATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTGGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTTGGCTTCACCGGTTCTTACAAAGGCCGCATAATTTCCGTAATGGGTCact                                                   <  2:49018‑M2/202‑2 (MQ=255)
                                                                                                                         ccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTTCTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGTGTa                                          <  2:412230‑M2/211‑1 (MQ=255)
                                                                                                                         ccggccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTa                                          <  1:250037‑M2/211‑1 (MQ=255)
                                                                                                                             ccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATccc                                      <  1:413704‑M2/215‑1 (MQ=255)
                                                                                                                             ccctctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCGACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACTTAGTTTTGATGCCAGGCGACCCTCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACATCGTTCGCGGTATGCTGTGCTTCTCCGGTTCTTACAATGGCCGCAAAATTTCCTTAagggttcacggtcggtggatccc                                      <  2:256947‑M2/215‑12 (MQ=255)
                                                                                                                                tctcccacagggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTc                                   >  1:444164‑M2/34‑251 (MQ=255)
                                                                                                                                          ggagagggtgaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                         <  1:337957‑M2/228‑1 (MQ=255)
                                                                                                                                                   gaaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGCTGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCaa                 >  2:493441‑M2/15‑250 (MQ=255)
                                                                                                                                                    aaaacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAg               <  1:108832‑M2/238‑1 (MQ=255)
                                                                                                                                                      aacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGa              >  1:285538‑M2/12‑250 (MQ=255)
                                                                                                                                                      aacatgttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGa              >  1:314858‑M2/12‑250 (MQ=255)
                                                                                                                                                           gttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGa         <  1:222425‑M2/244‑1 (MQ=255)
                                                                                                                                                            ttgatAAGGATAAAACCATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGCTTTGATGCACGGCGAACCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGAGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCACAATTTCCGTAATGGGTCACGGTATGGGAATCCCGTCCTGCTCCATCTACACCAAAGAACTGATc       >  1:234544‑M2/6‑251 (MQ=255)
                                                                                                                                                            ttgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGACCTGATc       >  1:399688‑M2/6‑251 (MQ=255)
                                                                                                                                                             tgatAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATc       >  1:247480‑M2/5‑250 (MQ=255)
                                                                                                                                                                 aaGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGa  <  2:603585/251‑1 (MQ=255)
                                                                                                                                                                 |                                                                                                                                                                                                                                                          
GGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGA  >  NC_000913/4620710‑4621121

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