Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NZ_CP009273 2,399,465 Δ7 bp coding (650‑656/939 nt) lrhA ← transcriptional regulator LrhA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NZ_CP009273 2399465 2399471 7 26 [0] [0] 26 lrhA transcriptional regulator LrhA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 23994640 (0.000)26 (1.160) 14/164 0.4 100% coding (657/939 nt) lrhA transcriptional regulator LrhA
?NZ_CP009273 2399472 = 0 (0.000)coding (649/939 nt) lrhA transcriptional regulator LrhA

GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCA  >  NZ_CP009273/2399383‑2399549
                                                                                 |                                                                                     
gTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagt                                                                               >  2:77537‑M1/1‑82 (MQ=255)
     ggCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgcca                                                                          <  1:553615‑M1/90‑14 (MQ=255)
         catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatgga                                                                      >  2:538562‑M1/1‑73 (MQ=255)
         catcatCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatgga                                                                      >  2:366881‑M1/1‑73 (MQ=255)
              tctcAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatc                                                                 <  2:141854‑M1/90‑23 (MQ=255)
                               gTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcg                                                <  2:441794‑M1/90‑40 (MQ=255)
                                 cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                              >  1:437364‑M1/1‑49 (MQ=255)
                                 cACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgcc                                              >  2:387391‑M1/1‑49 (MQ=255)
                                             ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                  >  2:681164‑M1/1‑37 (MQ=255)
                                             ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                  >  2:609236‑M1/1‑37 (MQ=255)
                                             ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                  >  1:597477‑M1/1‑37 (MQ=255)
                                             ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                  >  2:109656‑M1/1‑37 (MQ=255)
                                             ggCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcg                                  >  2:252742‑M1/1‑37 (MQ=255)
                                                   cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                            <  2:227112‑M1/90‑60 (MQ=255)
                                                   cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                            <  1:109656‑M1/90‑60 (MQ=255)
                                                   cACTGCCGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaac                            <  1:604238‑M1/90‑60 (MQ=255)
                                                        ccGCACGAACGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggct                       <  1:33226‑M1/90‑65 (MQ=255)
                                                                ccgGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcat               <  2:591132‑M1/89‑73 (MQ=255)
                                                                aCGGCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcat               <  2:597477‑M1/90‑73 (MQ=255)
                                                                   gCCGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatcca            <  1:563016‑M1/90‑76 (MQ=255)
                                                                     cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc          <  1:333197‑M1/90‑78 (MQ=255)
                                                                     cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc          <  2:554277‑M1/90‑78 (MQ=255)
                                                                     cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc          <  2:604565‑M1/90‑78 (MQ=255)
                                                                     cGGAAGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagc          <  1:206498‑M1/90‑78 (MQ=255)
                                                                        aaGCGTCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaat       <  1:252742‑M1/90‑81 (MQ=255)
                                                                             tCGAGaagccagtcgccatggaatatctgctttattaagcgtcgccagcaccatatcgcgaaacgggctaggatcatccagcaatacaag  >  2:215992‑M1/1‑5 (MQ=255)
                                                                                 |                                                                                     
GTCCGGGCTCATCATCTCAACCGGCCTTGCCGTCACGCCAAGACCGGCTTTCACTGCCGCACGAACGGCCGGAAGCGTCGAGGCGACATAAGCCAGTCGCCATGGAATATCTGCTTTATTAAGCGTCGCCAGCACCATATCGCGAAACGGGCTAGGATCATCCAGCA  >  NZ_CP009273/2399383‑2399549

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ 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.