Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,640,057 Δ7,591 bp [rrrQ][dicC] 18 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1640057 1647647 7591 20 [0] [0] 20 [rrrQ]–[dicC] [rrrQ],[ydfR],essQ,ynfU,cspB,cspF,ynfS,ydfT,ydfU,rem,hokD,relE,relB,ydfV,flxA,ydfW,ydfX,[dicC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 16400560 (0.000)20 (0.650) 20/490 0.5 100% coding (2/534 nt)
coding (310/312 nt)
rrrQ
ydfR
putative lysozyme
protein YdfR
?NC_000913 1647648 = 0 (0.000)coding (203/231 nt) dicC putative DNA‑binding transcriptional regulator DicC

GGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTAC  >  NC_000913/1639832‑1640305
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
ggTCGCATTTTTCCTTCGACAGTTTCATATTGGAAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTCCCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgata                                                                                                                                                                                                                                    <  2:230381‑M1/249‑24 (MQ=255)
              ttCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGGCAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttg                                                                                                                                                                                                                    >  1:193184‑M1/1‑212 (MQ=255)
                                                                      aCAGATGGTCCAGATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagattcctcctgtagacgcatcgcgcgaccttc                                                                                                                                                            <  1:263063‑M1/251‑96 (MQ=255)
                                                                             gTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTGCGGCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCTGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgagattcatcataaactttgggatcatactgaagctccccgccataggcctcctgtagacgcatcgcgcgaccttcgggaact                                                                                                                                                     <  1:189856‑M1/251‑103 (MQ=255)
                                                                                 atatGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaat                                                                                                                                                 >  2:111692‑M1/1‑145 (MQ=255)
                                                                                     gCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatccc                                                                                                                                              <  1:284877‑M1/250‑110 (MQ=255)
                                                                                         gAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgaagacgcatcgcgcgaccttcgggaactaaatcccctttc                                                                                                                                         >  2:10289‑M1/1‑137 (MQ=255)
                                                                                                           aTTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagcc                                                                                                                       >  1:216645‑M1/1‑119 (MQ=255)
                                                                                                                  ggTTACCTTCTTTTTAGTCCTGATACTGGTAAAGTATCTGAGGAGCAGATGCGCTAGCACCAATCAGCGCCAGAACGGCAGAAGACAGGCCGTATCTGATATTTTGCGTTCCcccgcccgcttcgtcttacgatcttcataataaactttgggatcatactgaagctccccgccagctgcctcctgtagacgcatcgcgcgaccttcgggaactaactcccctttcaagctataacgcgaagccaaacga                                                                                                                 >  2:263005‑M1/1‑112 (MQ=255)
                                                                                                                                      aGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCGTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagt                                                                                            <  1:192340‑M1/251‑160 (MQ=255)
                                                                                                                                             ggTCAAGTATCTGAGGAGCAGATGCGCATGCACCAAGCAGCGCCAGAACGGCAGCCGACAGGCGATATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttactatattcatcataaactttgggatcatactgaagctccccgccagatgcctcgtgtagactgatcgcgcgaccgtcgggagctaaatcccctgtccggctataaagcggagccaaacgaatacctgctgcttgtgcaagttttgtt                                                                                      <  1:285817‑M1/250‑166 (MQ=255)
                                                                                                                                                cAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttg                                                                                  >  2:95840‑M1/1‑82 (MQ=255)
                                                                                                                                                    tATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaacc                                                                              >  2:209468‑M1/1‑78 (MQ=255)
                                                                                                                                                              gCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagggttggttttgaaccgaaatacaaa                                                                    <  1:111692‑M1/251‑184 (MQ=255)
                                                                                                                                                                                  cAGCGCCTGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCCcccgcccgcttcgtcttacgatattcatcataaactttgggatcatgctgaagctccccgccagatgcctcctgtagacgcatcgcgcgacctgcgggaacgaaatcccctttgcagctataaagcgaaggcaaacgaatacctgctgcttgtgcgattgtggtttttgaaccgaaatacaaaagagcgtcagttttaagcat                                                <  1:19303‑M1/251‑204 (MQ=255)
                                                                                                                                                                                   aGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatt                                               <  1:209468‑M1/251‑205 (MQ=255)
                                                                                                                                                                                     cgcCCGAACGGCAGCCGACGGGCCGTATCGGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaattcctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcattta                                             <  2:241067‑M1/251‑207 (MQ=255)
                                                                                                                                                                                                               aTCTGATCTTTTGTGTTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgacgctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcataact                   >  2:48979‑M1/1‑19 (MQ=255)
                                                                                                                                                                                                                   gatcttttgtgtTCAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcataactaac                <  1:71380‑M1/246‑236 (MQ=255)
                                                                                                                                                                                                                                cAcccgcccgcttcgtcttacgatattcatcataaactttgggatcatactgaagctccccgccagatgcctcctgtagacgcatcgcgcgaccttcgggaactaaatcccctttccagctataaagcgaagccaaacgaatacctgctgcttgtgcaagttttgtttttgaaccgaaatacaaaagagcgtcagttttaagcatttaaaacacctttattgttagtcataactaacaagatagatgttaa  <  1:270435‑M1/251‑250 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                         
GGTCGCATTTTTCCTTCGACAGTTTCATATTGGGAAAAACGGTTTTTCCATCCACCACTGTGGCACCCCGACAGATGGTCCATATGCCAGAACCATCGCGGTATGCCATTGTGTGGTTACCTTCTTTTTCGTCCAGAAACTGGTCAAGTATCTGAGGAGCAGATGCGCCAGCACCAATCAGCGCCAGAACGGCAGCCGACAGGCCGTATCTGATTTTTGTGTTCATAGATATTTATGATGAGGACGCTCGTGCTTATTGGCAGGATTTTCAATCTTAAAGGAGTACTGATGCTGCAGATAAGACTCAACTTTTTCTGACAATTTTTCTGCTACTTCCAGGAAGACTTGCCGGACGCTCCTTCTGGCTGCTGCCTCATAAAACTCCAGCGCAGCTCCTTCAACACGGTCCATGGCGACATCCAGGTCAAAAATTTCACCGTCAAAGCGTTCTTTGTCCTGTAAGGCTACAGTTAC  >  NC_000913/1639832‑1640305

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.