Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,997,112 Δ26,567 bp [tcyN][fliR] 32 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1997112 2023678 26567 29 [0] [0] 27 [tcyN]–[fliR] [tcyN],tcyL,dcyD,tcyJ,fliZ,fliA,fliC,fliD,fliS,fliT,amyA,yedD,yedE,yedF,yedK,yedL,yedN,intG,fliE,fliF,fliG,fliH,fliI,fliJ,fliK,fliL,fliM,fliN,fliO,fliP,fliQ,[fliR]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19971110 (0.000)27 (0.890) 25/98 0.1 100% coding (704/753 nt) tcyN cystine ABC transporter ATP binding subunit
?NC_000913 2023679 = 0 (0.000)intergenic (+1/‑289) fliR/rcsA flagellar biosynthesis protein FliR/DNA‑binding transcriptional activator RcsA

TTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAAT  >  NC_000913/2023632‑2023725
                                               |                                              
ttgcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggtt                                               <  1:1887690‑M2/2‑1 (MQ=255)
   cagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGt                                           <  1:1548074‑M2/6‑1 (MQ=255)
       agaaacttctcgaggaactggcgggtgcgagggtgctcggTTCCGTa                                          >  1:2488490‑M2/41‑47 (MQ=255)
       agaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACg                                       <  1:1776982‑M2/10‑1 (MQ=255)
          aacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGttt                                    <  1:373947‑M2/13‑1 (MQ=255)
          aacttctcgaggaactggcgggtgcgaggctgctcggTTCCCTAACGtt                                     >  1:571883‑M2/38‑49 (MQ=255)
              tctcgaggaactggcgggtgcgaggctgctcggGTCCGTAACGTTTATCa                                >  1:1647219‑M2/34‑50 (MQ=255)
                   aggaactggcgggtgcggggctgctccgTTCCGTAACGTTTATCATGTTa                           >  1:183935‑M2/29‑50 (MQ=255)
                    ggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTAt                          <  1:2404503‑M2/23‑1 (MQ=255)
                     gaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATc                         <  1:878902‑M2/24‑1 (MQ=255)
                      aactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATcc                        <  1:2280440‑M2/25‑1 (MQ=255)
                       actggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                       <  1:716022‑M2/26‑1 (MQ=255)
                        ctggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTa                      <  1:572172‑M2/27‑1 (MQ=255)
                            cgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGa                  <  1:814804‑M2/31‑1 (MQ=255)
                            cgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGa                  >  1:1510223‑M2/20‑50 (MQ=255)
                              ggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                 >  1:512741‑M2/18‑49 (MQ=255)
                                tgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTAt              <  1:2690129‑M2/35‑1 (MQ=255)
                                    aggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCg           >  1:1532812‑M2/12‑49 (MQ=255)
                                    aggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGa          <  1:1259788‑M2/39‑1 (MQ=255)
                                      gctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa        <  1:550636‑M2/41‑1 (MQ=255)
                                      gctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaa         >  1:1493687‑M2/10‑49 (MQ=255)
                                       ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaa       <  1:1832385‑M2/42‑1 (MQ=255)
                                       ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaa       <  1:226919‑M2/42‑1 (MQ=255)
                                       ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa        >  1:2666333‑M2/9‑49 (MQ=255)
                                       ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa        <  1:2235784‑M2/41‑1 (MQ=255)
                                          ctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaata    <  1:2674844‑M2/45‑1 (MQ=255)
                                             ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataat  <  1:589993‑M2/47‑1 (MQ=255)
                                               |                                              
TTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAAT  >  NC_000913/2023632‑2023725

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: