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

CGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACG  >  NC_000913/3685308‑3685602
                                                                                                                                                  |                                                                                                                                                    
gcgttaagttcgcgagcaatattgcgtaacggatgcaaaatcaggcggttaatgcaccagctgatggcgacggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTcg                                                                                                                                                   <  1:333572‑M2/4‑1 (MQ=255)
         tcgcgagcaatattgcgtaacggatgcaaaatcaggcggttaatgcaccagctgatggcgacggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACagag                                                                                                                                          <  2:299237‑M2/13‑1 (MQ=255)
            cgagcaatattgcgtaacggatgcaaaatcaggcggttaatgcaccagctgatggcgacggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAgagcgag                                                                                                                                      <  2:188521‑M2/17‑1 (MQ=255)
                    attgcgtaacggatgcaaaatcaggcggttaatgcaccagctgatggtgacgggcaggataagcgacaaaaggaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCtt                                                                                                                              >  1:125927‑M2/127‑151 (MQ=255)
                       gcgtaacggatgcaaaatcaggcggttaatgcaccagctgatggcgacggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTcc                                                                                                                           >  2:262901‑M2/124‑151 (MQ=255)
                                                      caccagctgatggcgacggtcacgataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTctgct                                                                                            <  1:75921‑M2/59‑1 (MQ=255)
                                                              gatggcgacggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAg                                                                                     <  2:388953‑M2/66‑1 (MQ=255)
                                                                     acggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAAc                                                                              >  1:167580‑M2/78‑150 (MQ=255)
                                                                      cggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAg                                                                            <  1:57094‑M2/75‑1 (MQ=255)
                                                                      cggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAg                                                                            <  1:20587‑M2/75‑1 (MQ=255)
                                                                       ggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGc                                                                           <  1:296418‑M2/76‑1 (MQ=255)
                                                                       ggtcaggataagcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGc                                                                           <  2:84003‑M2/76‑1 (MQ=255)
                                                                                  gcgacaaaagtaagtaaatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCa                                                                <  2:235758‑M2/87‑1 (MQ=255)
                                                                                        aaagtaagtaaatggtcactaacgttgagcgggtgctaatcacgaccttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAAt                                                          <  1:330246‑M2/93‑1 (MQ=255)
                                                                                                  aatggtcactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAac                                                <  2:293156‑M2/103‑1 (MQ=255)
                                                                                                        cactaacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAgg                                          >  1:198440‑M2/43‑151 (MQ=255)
                                                                                                            aacgttgagagggtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACTCCAGGAc                                        <  1:291048‑M2/111‑1 (MQ=255)
                                                                                                                        gtgctcatcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGccacca                          <  1:113449‑M2/125‑1 (MQ=255)
                                                                                                                             catcacgaacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGcc                     >  1:351240‑M2/22‑151 (MQ=255)
                                                                                                                                    aacttatacatacgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGAGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAAAGAACACCAGGACAACGGCTGCAACCATTGCCATCt                 >  1:82651‑M2/15‑148 (MQ=255)
                                                                                                                                               acgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAAc   <  2:212464‑M2/148‑1 (MQ=255)
                                                                                                                                                cgTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACg  <  2:337081‑M2/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTCTGGAATAAAACTTTTACGCAGCGCCTGGAACTGGTTAGGCAGCACTACATCGGCACGGCTGACCACGCCTGCCGGTTTAATGCTGGCAAGAATAGCTTCCGCTTCGGGAATATCGCCTTTCAAAATAGCTGAAGATAAGGGTTCGCGGACAGAGCGAGCAATGCTTTCCAGTTGCGTAGCCGTGTTATAGCGATTCTGCTGGACCAGATGGAACAGCAAAACGGTGCAAAAAATAAAAACGAACACCAGGACAACGGCTGCCACCATTGCCATCTGCTTGATTGTTAACG  >  NC_000913/3685308‑3685602

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.