Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,398,017 312 bp→TTAGTTCCTATTCCGAAGTTCCTATTCTCTAGAAAGTATAGGAACTTCGAAGCAGCTCCAGCCTACAC coding (51‑362/594 nt) yhdE ← nucleoside triphosphate pyrophosphatase YhdE

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3398017 3398328 312 7 [0] [0] 10 yhdE nucleoside triphosphate pyrophosphatase YhdE

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 33980160 (0.000)5 (0.780)
+TTAGTTCCTATTCCGAAGTTCCTATTCTCTAGAAAGTATAGGAACTTCGAAGCAGCTCCAGCCTACAC
5/168 0.2 100% coding (363/594 nt) yhdE nucleoside triphosphate pyrophosphatase YhdE
?NC_000913 3398329 = 0 (0.000)coding (50/594 nt) yhdE nucleoside triphosphate pyrophosphatase YhdE

TCAGGATAACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGC  >  NC_000913/3398153‑3398573
                                                                                                                                                                                |                                                                                                                                                                                                                                                    
ccacccagcccctgaataccgtatgcacctgctttatctaacggttcatcgctggcgacatagcccgcgatgtcttcgtctgttaacgttctgaaagtcacagcggtgttagttcctattccgaagttcctattctctagaaagtataggaacttcgaagcagctccagcctacacAGTAACTCCTGACGACGCGGAGAACCGGAAg                                                                                                                                                                                                                        <  1:133427‑M2/31‑1 (MQ=255)
  acccagcccctgaataccgtatgcacctgctttatctaacggttcatcgctggcgacatagcccgcgatgtcttcgtctgttaacgttctgaaagtcacatcggtgttagttcctattccgaagttcctattctctagaaagtataggaacttcgaagcagctccagcctacacAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAActgc                                                                                                                                                                          <  1:119481‑M2/77‑1 (MQ=255)
              aataccgtatgcacctgctttatctaacggttcatcgctggcgacatagcccgcgatgtcttcgtctgttaacgttctgaaagtcacatcggtgttagttcctattccgaagttcctattctctagaaagtataggaacttcgaagcagctccagcctacacAGTAACTCCTGACGACGCGGAGAACCGTAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTg                                                                                                                                                              <  1:147426‑M2/89‑1 (MQ=255)
                                                               cccgcgatgtcttcgtctgttaacgttctgaaagtcacatcggtgttagttcctattccgaagttcctattctctagaaagtataggaacttcgaagcagctccatcctacacAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAg                                                                                                                                                                                                         <  1:61130‑M2/46‑1 (MQ=255)
                                                               cccgcgatgtcttcgtctgttaacgtnctgaaanncanatcggtgttagttcctattncgnnntnnntnttctcnagnnagtataggancttcgaagcagctccagcctacacAGTAACTCCTGACGACGCGGAGAACCGGANGCTAAATACNNAGAAg                                                                                                                                                                                                         >  2:61130‑M2/114‑159 (MQ=255)
                                                                                                                                                               gcagctccagcctacacAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCa                                                                                                                                                       >  1:5279‑M2/18‑113 (MQ=255)
                                                                                                                                                                nagctccagcctacacAGTAACTCCTGACGACGCGGAGGACCGGAAGCTAAATACAGAGAAGTCATAGAGACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGACCAGAGCACACCATTGACTACACTACTCCAGAACACTTCCGGACTGAAag                                                                                     >  1:48477‑M2/17‑178 (MQ=255)
                                                                                                                                                                  gctccagcctacacAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATg                    <  1:2057‑M2/227‑1 (MQ=255)
                                                                                                                                                                  gctccagcctacacAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAgcg          >  1:156705‑M2/15‑251 (MQ=255)
                                                                                                                                                                            acacAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCtgctgc  <  1:85505‑M2/245‑1 (MQ=255)
                                                                                                                                                                                |                                                                                                                                                                                                                                                    
TCAGGATAACGATAGTATCCGCACCCAGCACCGGGAGATCCTTCGCCGTTTGCGCGACACCTGCCCGTGCTTTCTCGCGCGCCAGACGCACAACATACTGCTGCGCGCTCTCCTGCGGCTGACGCTGCTCCTCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGC  >  NC_000913/3398153‑3398573

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.