Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,962,566 Δ16,764 bp flhEinsA‑5 17 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1962566 1979329 16764 98 [0] [0] 96 flhE–insA‑5 flhE,flhA,flhB,cheZ,cheY,cheB,cheR,tap,tar,cheW,cheA,motB,motA,flhC,flhD,insB‑5,insA‑5

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19625650 (0.000)94 (0.980) 57/96 0.0 100% intergenic (+105/+15) yecT/flhE protein YecT/flagellar protein
?NC_000913 1979330 = 0 (0.000)intergenic (‑115/‑423) insA‑5/uspC IS1 protein InsA/universal stress protein C

AAGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962516‑1962614
                                                 |                                                 
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   <  1:593583‑M1/50‑1 (MQ=255)
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   >  1:577285‑M1/1‑50 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGGGCGCGTc                                                  <  1:373052‑M1/50‑2 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  >  1:1374074‑M1/1‑49 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:2336379‑M1/50‑2 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:1624705‑M1/50‑2 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:6401832‑M1/50‑2 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:3366773‑M1/50‑2 (MQ=255)
  gACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcg                                                 >  1:6346766‑M1/1‑48 (MQ=255)
   aCATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgc                                                >  1:7060043‑M1/1‑47 (MQ=255)
     aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:2060841‑M1/1‑45 (MQ=255)
     aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              <  1:2161147‑M1/50‑6 (MQ=255)
     aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              <  1:1845624‑M1/50‑6 (MQ=255)
      tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:1868635‑M1/50‑7 (MQ=255)
      tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             >  1:6779467‑M1/1‑44 (MQ=255)
      tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:2370265‑M1/50‑7 (MQ=255)
        gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:1668972‑M1/1‑42 (MQ=255)
        gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:5587448‑M1/1‑42 (MQ=255)
        gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:3132635‑M1/1‑42 (MQ=255)
        gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacg                                            <  1:2920289‑M1/49‑8 (MQ=255)
         aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgct                                          <  1:690294‑M1/50‑10 (MQ=255)
         aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgct                                          <  1:4931795‑M1/50‑10 (MQ=255)
          gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                         <  1:2442059‑M1/50‑11 (MQ=255)
           ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGTGTcgctacgctag                                        >  1:1219817‑M1/1‑39 (MQ=255)
            cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                        >  1:4096608‑M1/1‑38 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    >  1:3206470‑M1/1‑35 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaa                                     <  1:5607085‑M1/49‑15 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagca                                      <  1:2061757‑M1/48‑14 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagca                                      <  1:5876067‑M1/48‑14 (MQ=255)
                ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    >  1:5168482‑M1/1‑34 (MQ=255)
                  cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                  >  1:301715‑M1/1‑32 (MQ=255)
                  cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                  >  1:5289358‑M1/1‑32 (MQ=255)
                   aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                                <  1:1418252‑M1/50‑20 (MQ=255)
                   aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                                <  1:1932164‑M1/50‑20 (MQ=255)
                   aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaa                                 <  1:6014828‑M1/49‑19 (MQ=255)
                    cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaa                               <  1:3991067‑M1/50‑21 (MQ=255)
                    cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                                >  1:2645715‑M1/1‑30 (MQ=255)
                    cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                                >  1:3915401‑M1/1‑30 (MQ=255)
                     tGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaa                              >  1:7110918‑M1/1‑29 (MQ=255)
                     tGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaa                              >  1:2124606‑M1/1‑29 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:2478515‑M1/49‑23 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:7010381‑M1/49‑23 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:7361746‑M1/49‑23 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:6953916‑M1/49‑23 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:2068738‑M1/1‑27 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaa                              <  1:6827739‑M1/48‑22 (MQ=255)
                        gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           <  1:4127694‑M1/50‑25 (MQ=255)
                         ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           >  1:807908‑M1/1‑25 (MQ=255)
                         ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           >  1:4644293‑M1/1‑25 (MQ=255)
                           attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:4326664‑M1/1‑23 (MQ=255)
                            ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaaca                       <  1:1198694‑M1/50‑29 (MQ=255)
                            ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaac                        <  1:860136‑M1/49‑28 (MQ=255)
                            ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaac                        <  1:6436586‑M1/49‑28 (MQ=255)
                            ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaac                        <  1:2771360‑M1/49‑28 (MQ=255)
                             taAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:7131187‑M1/1‑21 (MQ=255)
                              aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:4742572‑M1/1‑20 (MQ=255)
                              aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:474886‑M1/1‑20 (MQ=255)
                              aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:906882‑M1/1‑20 (MQ=255)
                               aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                     <  1:4496526‑M1/49‑31 (MQ=255)
                                aCATACTCGCGAGCGCGTcgctacgttagcaagcaaaaacgaaacaggaa                   >  1:202755‑M1/1‑18 (MQ=255)
                                aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   <  1:741463‑M1/50‑33 (MQ=255)
                                 cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   >  1:2122683‑M1/1‑17 (MQ=255)
                                 cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   >  1:3887470‑M1/1‑17 (MQ=255)
                                  aTACTCGGGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaata                 <  1:3041844‑M1/50‑35 (MQ=255)
                                  aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaata                 <  1:6624847‑M1/50‑35 (MQ=255)
                                  aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaata                 <  1:6557748‑M1/50‑35 (MQ=255)
                                  aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaata                 <  1:3901242‑M1/50‑35 (MQ=255)
                                  aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaata                 <  1:6067002‑M1/50‑35 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagccagcaaaaatgaaacaggaataa                >  1:139861‑M1/1‑15 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcagaaatgaaacaggaataa                <  1:5019063‑M1/50‑36 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:6203657‑M1/50‑36 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                >  1:6304605‑M1/1‑15 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:449405‑M1/50‑36 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:3810437‑M1/50‑36 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                >  1:2189573‑M1/1‑15 (MQ=255)
                                    aCTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat               >  1:1802376‑M1/1‑14 (MQ=255)
                                     cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              >  1:4565759‑M1/1‑13 (MQ=255)
                                     cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              <  1:3211691‑M1/50‑38 (MQ=255)
                                       cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg             >  1:4143649‑M1/1‑11 (MQ=255)
                                        gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga            >  1:968932‑M1/1‑10 (MQ=255)
                                        gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga            >  1:4550225‑M1/1‑10 (MQ=255)
                                         cgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa          >  1:6192094‑M1/1‑9 (MQ=255)
                                         cgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa          <  1:7462140‑M1/50‑42 (MQ=255)
                                          gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:6939619‑M1/1‑8 (MQ=255)
                                          gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:2139115‑M1/1‑8 (MQ=255)
                                           aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        <  1:3505668‑M1/50‑44 (MQ=255)
                                             cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg      <  1:3969505‑M1/50‑46 (MQ=255)
                                             cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg       >  1:60948‑M1/1‑5 (MQ=255)
                                              gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     >  1:4142222‑M1/1‑4 (MQ=255)
                                              gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     <  1:871998‑M1/50‑47 (MQ=255)
                                               cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    >  1:2445119‑M1/1‑3 (MQ=255)
                                               cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    >  1:564022‑M1/1‑3 (MQ=255)
                                               cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    <  1:5683610‑M1/50‑48 (MQ=255)
                                                gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg   <  1:5654956‑M1/50‑49 (MQ=255)
                                                gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg   <  1:2904828‑M1/50‑49 (MQ=255)
                                                gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg   <  1:4941205‑M1/50‑49 (MQ=255)
                                                 tcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt  <  1:6571059‑M1/50‑50 (MQ=255)
                                                 tcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt  <  1:516669‑M1/50‑50 (MQ=255)
                                                 |                                                 
AAGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962516‑1962614

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