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

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46208272 (0.040)19 (0.550)
+AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT
18/354 0.8 92.9% intergenic (+1/‑56) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase
?NC_000913 4620871 = 2 (0.040)intergenic (+45/‑12) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

CGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAAT  >  NC_000913/4620584‑4621055
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
cgTCTTCCGGAGATGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCGCGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggt                                                                                                                                                                                                                               <  1:3585‑M1/251‑8 (MQ=255)
      ccGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggct                                                                                                                                                                                                                         >  1:473710‑M1/1‑238 (MQ=255)
       cGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggctt                                                                                                                                                                                                                        <  2:291553‑M1/251‑15 (MQ=255)
                              gATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggcccccgc                                                                                                                                                                                                 <  2:163271‑M1/251‑38 (MQ=255)
                                      cATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACaaaaa                                                                                                                                                                                         <  2:145314/251‑1 (MQ=255)
                                       aTCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgcccccaccccggccccctcccacagggg                                                                                                                                                                                        >  2:135154‑M1/1‑205 (MQ=255)
                                            gttccTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagaggg                                                                                                                                                                                   <  1:493516‑M1/249‑52 (MQ=255)
                                                                                          gaacacacgcGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacccatt                                                                                                                                      <  2:241578‑M1/247‑97 (MQ=255)
                                                                                               cacGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacaggttgataaggataaaacaatggctaccccacacattaatgca                                                                                                                                <  1:437673‑M1/251‑103 (MQ=255)
                                                                                                  gcgTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaa                                                                                                                             <  1:325779‑M1/251‑106 (MQ=255)
                                                                                                            aTTCCGGGACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacacggggaggggggaaacatgttggtaaggctaaaacaattgctaccccaacacttcattcagaaatgggggatt                                                                                                                   >  2:206605‑M1/1‑136 (MQ=255)
                                                                                                                       ggtttacggcccgaaaGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacg                                                                                                          <  1:390609‑M1/242‑125 (MQ=255)
                                                                                                                          ttacggcccgAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagt                                                                                                      <  1:484740‑M1/246‑129 (MQ=255)
                                                                                                                               gCCCGAAAGTAAAACCGGGCTCAGTGGGTCACCGTGAAACCCTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCCGATATGGAATATGGCAAAGCCATGTTCTGATaattggtgcgtcttgatgccctcaccccggccctctcccacagggatagggtgaaaacatgttgataaggataaatcaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgat                                                                                                <  1:215489‑M1/251‑135 (MQ=255)
                                                                                                                                   gAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccaacgggagagggtgaaaacatgttgataaggataaaacaatggccacccccaaacttaatgacaaaaagggagatttcgatgacgtagttttgaagccc                                                                                            >  2:282736‑M1/1‑113 (MQ=255)
                                                                                                                                       gTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcg                                                                                        <  2:474111‑M1/251‑143 (MQ=255)
                                                                                                                                        tAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctccccccggccctctcccacagggggagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcg                                                                                        >  1:191584‑M1/1‑108 (MQ=255)
                                                                                                                                         aaaaCCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcacccccgccctctcccacagggggggggtgaaaacatgttgataagggtaaaacaatggctaccccacaaattaatgcagaaatgggcgattttgctttcgtaggtttgatgcccggcagc                                                                                      >  1:162453‑M1/1‑107 (MQ=255)
                                                                                                                                         aaaaCCGGGCTCACGGGGTAACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatggtgataaggataaaacaatggcttccccacacattaatgcagaaatgggcgatttcgctgacgtagttgtgatgccaggcgac                                                                                      <  2:61734‑M1/251‑145 (MQ=255)
                                                                                                                                            aCCGGGCTCCCTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccccagggagggggtgaaaacatgtttataaagataaaacaatggctaccccaccccttaatgcagaaatgggcgatttccctcacgtagttttgatgcccagcgacccc                                                                                   >  2:59115‑M1/1‑104 (MQ=255)
                                                                                                                                                                                 attggtcagactcTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaa                                               <  1:401039‑M1/244‑184 (MQ=255)
                                                                                                                                                                                   tggtcagacTCGGGCAAAATAGTTTGGTACTTCGGATATGGAATAGGGCAAAGCCATGTTCTGATaatgggggcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatgtctaccccacacattaatgtagaaatgggcgatttcgctgtcgtagttttgatgctaggcgacacgctgcgtgcgaagtatattgctgaaactttcccggaaga                                             <  2:52598‑M1/246‑186 (MQ=255)
                                                                                                                                                                                      gcaGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgattccgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcc                                         <  1:135154‑M1/250‑190 (MQ=255)
                                                                                                                                                                                       cAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatggcctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatgtctaccccacacattaaggccgaaatgggagattttgcttgcgtagttttgatgcccggcgacccgctgcgtgcgaagtgtgtttctgaaacttttctttaaaatgccc                                        >  2:277667‑M1/1‑61 (MQ=255)
                                                                                                                                                                                                                          ggAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgagaaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttc     <  2:409257‑M1/251‑226 (MQ=255)
                                                                                                                                                                                                                           gAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttggtaagggtaaaaaaatggctaccccacacattaatgcagaaatgggggatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcggaggatattgctgaaaatttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttc     >  1:422379‑M1/1‑25 (MQ=255)
                                                                                                                                                                                                                             atatGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACACCACACATTAACGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAAt  <  1:140268/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                    
CGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAAT  >  NC_000913/4620584‑4621055

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