Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,685,181 Δ273 bp coding (166‑438/1989 nt) yhjK ← cyclic‑di‑GMP phosphodiesterase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3685181 3685453 273 22 [0] [0] 22 yhjK cyclic‑di‑GMP phosphodiesterase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 36851800 (0.000)22 (0.890) 20/296 0.4 100% coding (439/1989 nt) yhjK cyclic‑di‑GMP phosphodiesterase
?NC_000913 3685454 = 0 (0.000)coding (165/1989 nt) yhjK cyclic‑di‑GMP phosphodiesterase

GCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAGGCGGTTAATGCACCAGCTGATGGCGACGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGGAAGGAATCCGCCTGTAGTACCAGATAGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTAC  >  NC_000913/3685035‑3685329
                                                                                                                                                 |                                                                                                                                                     
gCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAGGCGGTTAATGCACCAGCTGATGGCGACGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcg                                                                                                                                                   <  1:333572‑M1/150‑5 (MQ=255)
         tCGCGAGCAATATTGCGTAACGGATGCAAAATCAGGCGGTTAATGCACCAGCTGATGGCGACGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagag                                                                                                                                          <  2:299237‑M1/150‑14 (MQ=255)
            cgAGCAATATTGCGTAACGGATGCAAAATCAGGCGGTTAATGCACCAGCTGATGGCGACGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgag                                                                                                                                      <  2:188521‑M1/151‑18 (MQ=255)
                    atTGCGTAACGGATGCAAAATCAGGCGGTTAATGCACCAGCTGATGGTGACGGGCAGGATAAGCGACAAAAGGAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctt                                                                                                                              >  1:125927‑M1/1‑126 (MQ=255)
                       gCGTAACGGATGCAAAATCAGGCGGTTAATGCACCAGCTGATGGCGACGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttcc                                                                                                                           >  2:262901‑M1/1‑123 (MQ=255)
                                                      cACCAGCTGATGGCGACGGTCACGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgct                                                                                            <  1:75921‑M1/151‑60 (MQ=255)
                                                              gATGGCGACGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccag                                                                                     <  2:388953‑M1/150‑67 (MQ=255)
                                                                     aCGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaac                                                                              >  1:167580‑M1/1‑77 (MQ=255)
                                                                      cGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacag                                                                            <  1:57094‑M1/151‑76 (MQ=255)
                                                                      cGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacag                                                                            <  1:20587‑M1/151‑76 (MQ=255)
                                                                       ggTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagc                                                                           <  1:296418‑M1/151‑77 (MQ=255)
                                                                       ggTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagc                                                                           <  2:84003‑M1/151‑77 (MQ=255)
                                                                                  gCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgca                                                                <  2:235758‑M1/151‑88 (MQ=255)
                                                                                        aaaGTAAGTAAATGGTCACTAACGTTGAGCGGGTGCTAATCACGACCTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaat                                                          <  1:330246‑M1/151‑94 (MQ=255)
                                                                                                  aaTGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaacgaac                                                <  2:293156‑M1/151‑104 (MQ=255)
                                                                                                        cACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaacgaacaccagg                                          >  1:198440‑M1/1‑42 (MQ=255)
                                                                                                            aaCGTTGAGAGGGTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaacgaactccaggac                                        <  1:291048‑M1/149‑112 (MQ=255)
                                                                                                                        gTGCTCATCACGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaacgaacaccaggacaacggctgccacca                          <  1:113449‑M1/151‑126 (MQ=255)
                                                                                                                             catcaCGAACTTATACATACGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaacgaacaccaggacaacggctgccaccattgcc                     >  1:351240‑M1/1‑21 (MQ=255)
                                                                                                                                    aaCTTATACATACGttcgcggacagagcgagcaatgctttccagttgagtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaaagaacaccaggacaacggctgcaaccattgccatct                 >  1:82651‑M1/1‑14 (MQ=255)
                                                                                                                                               acGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaacgaacaccaggacaacggctgccaccattgccatctgcttgattgttaac   <  2:212464‑M1/151‑149 (MQ=255)
                                                                                                                                                cGttcgcggacagagcgagcaatgctttccagttgcgtagccgtgttatagcgattctgctggaccagatggaacagcaaaacggtgcaaaaaataaaaacgaacaccaggacaacggctgccaccattgccatctgcttgattgttaacg  <  2:337081‑M1/151‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAGGCGGTTAATGCACCAGCTGATGGCGACGGTCAGGATAAGCGACAAAAGTAAGTAAATGGTCACTAACGTTGAGAGGGTGCTCATCACGAACTTATACATACGGAAGGAATCCGCCTGTAGTACCAGATAGGCAATTGGCTGCGGGTTTGCCGGACGTTCGAGCGAGTAAACGCCCAGCGAGATTTGAACCGGTAGCTCAAACAGGCGAGTAACCATTACCGGCACCGGGCGCTCTGGAATAAAACTTTTAC  >  NC_000913/3685035‑3685329

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