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 44 [0] [0] 44 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)44 (1.480) 32/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

GACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGT  >  NC_000913/1962518‑1962613
                                               |                                                
gACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcg                                                <  1:1186659‑M1/50‑3 (MQ=255)
   aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                             <  1:1118940‑M1/50‑6 (MQ=255)
    tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                            <  1:463101‑M1/50‑7 (MQ=255)
    tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                            <  1:2613259‑M1/50‑7 (MQ=255)
    tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                            <  1:2593871‑M1/50‑7 (MQ=255)
    tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                            <  1:1723816‑M1/50‑7 (MQ=255)
     cGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacg                                           >  1:891163‑M1/1‑43 (MQ=255)
         ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                       <  1:1023175‑M1/50‑12 (MQ=255)
         ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                       <  1:507581‑M1/50‑12 (MQ=255)
          cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                      <  1:2276801‑M1/50‑13 (MQ=255)
          cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                       >  1:1961909‑M1/1‑38 (MQ=255)
              ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                   >  1:565541‑M1/1‑34 (MQ=255)
              ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                   >  1:670047‑M1/1‑34 (MQ=255)
               tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                  <  1:2505840‑M1/49‑17 (MQ=255)
                cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                 >  1:220484‑M1/1‑32 (MQ=255)
                 aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaa                                <  1:827444‑M1/49‑19 (MQ=255)
                   tGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaa                             >  1:466251‑M1/1‑29 (MQ=255)
                     aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                           <  1:2222362‑M1/50‑24 (MQ=255)
                     aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                            <  1:2452560‑M1/49‑23 (MQ=255)
                      gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                           >  1:399422‑M1/1‑26 (MQ=255)
                       ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                          >  1:1437363‑M1/1‑25 (MQ=255)
                       ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                          >  1:627496‑M1/1‑25 (MQ=255)
                          ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaaca                      <  1:948295‑M1/50‑29 (MQ=255)
                            aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                     >  1:1533281‑M1/1‑20 (MQ=255)
                             aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagga                   <  1:686332‑M1/50‑32 (MQ=255)
                             aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                    >  1:2679344‑M1/1‑19 (MQ=255)
                             aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                    >  1:432989‑M1/1‑19 (MQ=255)
                             aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                    >  1:41809‑M1/1‑19 (MQ=255)
                              aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                  <  1:2153649‑M1/50‑33 (MQ=255)
                               cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                 <  1:1527865‑M1/50‑34 (MQ=255)
                                  aCTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat              >  1:1420286‑M1/1‑14 (MQ=255)
                                  aCTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat              >  1:2115614‑M1/1‑14 (MQ=255)
                                   cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat              <  1:73133‑M1/49‑37 (MQ=255)
                                   cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat              <  1:106‑M1/49‑37 (MQ=255)
                                     cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg            >  1:2189809‑M1/1‑11 (MQ=255)
                                     cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg            >  1:120991‑M1/1‑11 (MQ=255)
                                      gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaa          <  1:2519725‑M1/50‑41 (MQ=255)
                                      gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga           >  1:2391689‑M1/1‑10 (MQ=255)
                                      gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg            <  1:496061‑M1/48‑39 (MQ=255)
                                         aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       <  1:354665‑M1/50‑44 (MQ=255)
                                         aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg       <  1:164682‑M1/50‑44 (MQ=255)
                                           cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     <  1:314876‑M1/50‑46 (MQ=255)
                                              gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  <  1:67185‑M1/50‑49 (MQ=255)
                                              gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  >  1:244776‑M1/1‑2 (MQ=255)
                                               |                                                
GACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGT  >  NC_000913/1962518‑1962613

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