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

GACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTAT  >  NC_000913/4620578‑4621076
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
atcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgatcataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAATCCATGTTCTGata                                                                                                                                                                                                                                                          <  1:193515‑M2/80‑2 (MQ=255)
atcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGAt                                                                                                                                                                                                                                                           >  2:377701‑M2/172‑250 (MQ=255)
  cagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCtgataat                                                                                                                                                                                                                                                        <  1:62792‑M2/82‑4 (MQ=255)
             atacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAaagccatgttctgataatgggggcggctt                                                                                                                                                                                                                                             >  1:105492‑M2/159‑237 (MQ=255)
             atacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAagccatgttctgataatgggtgcggct                                                                                                                                                                                                                                              <  2:167743‑M2/92‑14 (MQ=255)
                      gaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGaatatggcaaagccatgttctgataatgggtgcggcttgatgccctc                                                                                                                                                                                                                                    >  2:57547‑M2/150‑228 (MQ=255)
                                                                                                                                                           agttggagtcattaccaccttcgcggatattggTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcg                                                                                               <  2:158780‑M1/235‑157 (MQ=255)
                                                                                                                                                               ggagtcattaccaccttcgcggataTTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggggggggtgaaaacctgttgataaggataaaacaatggctacaccacacattaatgcaggaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagt                                                                                           >  2:601573‑M1/13‑91 (MQ=255)
                                                                                                                                                                 agtcattaccaccttcgcggaTATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccccctaccccagggggagggggaaaacatgttgataaggataaaacaaaggctaccccacacattaatgcagaaatgggcggtttcgctgacgtaggttttattccaagcgacccgctgcgtgcgaagtat                                                                                         >  1:567934‑M1/11‑89 (MQ=255)
                                                                                                                                                                    cattaccaccttcgcGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccccagggggagggggaaaacaagttgataaggataaaacaatggctaccccacccattaatgcagaaaagggcgatttcgctgacgtagttttgctgccaggcgaccccctgcgggcgaaggatatt                                                                                      >  1:183732‑M1/8‑86 (MQ=255)
                                                                                                                                                                     attaccaccttcgCGGATAGTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctgccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattg                                                                                     <  2:314312‑M1/245‑167 (MQ=255)
                                                                                                                                                                         ccaccTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctga                                                                                 <  2:591464‑M1/249‑171 (MQ=255)
                                                                                                                                                                         ccaccTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctga                                                                                 <  1:462901‑M1/249‑171 (MQ=255)
                                                                                                                                                                                cgcgGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATAGGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacgatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtaggtttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaacttt                                                                           <  2:432537‑M1/250‑177 (MQ=255)
                                                                                                                                                                                   ggATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttcctt                                                                       <  2:251231‑M1/251‑181 (MQ=255)
                                                                                                                                                                                             cAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccccagggagagggggaaaacatgttgataaggataaaacaaatgctaccccccacattaatgcagaaaagggcgatttcgctgacgtagtttttatgcccg                                                                                                                >  1:355154‑M1/1‑61 (MQ=255)
                                                                                                                                                                                             cAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccag                                                                                                                <  2:355154‑M1/200‑140 (MQ=255)
                                                                                                                                                                                             cAGACTCTGGCAAAATATTTTGGTACGTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggattgatgccctcacaccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacagtaatgcagaaatgggcgatttcgctgacgtagttgtgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgccc                                                             <  2:464069‑M1/251‑191 (MQ=255)
                                                                                                                                                                                                     ggCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagt                                                      <  1:57547‑M1/250‑198 (MQ=255)
                                                                                                                                                                                                             atTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacg                                              <  1:400342‑M1/250‑206 (MQ=255)
                                                                                                                                                                                                                           gATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagaggcccgtgaagtgaacaacgttcgcggtatgctgg                               <  2:105492‑M1/251‑221 (MQ=255)
                                                                                                                                                                                                                            atatGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctggg                              <  1:420090‑M1/251‑222 (MQ=255)
                                                                                                                                                                                                                            atatGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctggg                              <  1:510313‑M1/251‑222 (MQ=255)
                                                                                                                                                                                                                                ggAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggccaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcggagtatattgctgaaactttccttgaagatgcccgtgcagtgaaccacgttcgcggtattctgggcttt                          >  2:151755‑M1/1‑26 (MQ=255)
                                                                                                                                                                                                                                    tatGGCAAAGCCATGTTCTGATaatgggtgcggcttgatggcctcaccccgggcctctcccacagggagagggggaaaacaagttgataaggataaaacaatggctaccccccacattaaagccgaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatatttcttaaaatttccttgaagatgcccgtgaagtggacaacgttcgcgggatgctgggcttcaccc                      >  1:571319‑M1/1‑22 (MQ=255)
                                                                                                                                                                                                                                        gCAAAGCCATGTTCTGATaatgggtgcggctttatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtac                  <  1:377701‑M1/251‑234 (MQ=255)
                                                                                                                                                                                                                                                       gATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaa   <  1:435167‑M1/251‑249 (MQ=255)
                                                                                                                                                                                                                                                        aTaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaa  <  1:601573‑M1/251‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                         
GACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTAT  >  NC_000913/4620578‑4621076

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