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 74 [0] [0] 74 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)73 (1.290) 47/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

AGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962517‑1962614
                                                |                                                 
aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:1297018‑M1/50‑2 (MQ=255)
    aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:89741‑M1/1‑45 (MQ=255)
    aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:1990104‑M1/1‑45 (MQ=255)
    aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:2655590‑M1/1‑45 (MQ=255)
    aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:2974084‑M1/1‑45 (MQ=255)
     tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:214888‑M1/50‑7 (MQ=255)
       gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:3491913‑M1/1‑42 (MQ=255)
       gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacg                                            <  1:2183929‑M1/49‑8 (MQ=255)
        aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgct                                          <  1:299706‑M1/50‑10 (MQ=255)
        aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgct                                          <  1:1795366‑M1/50‑10 (MQ=255)
        aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgct                                          <  1:2443453‑M1/50‑10 (MQ=255)
        aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgct                                          <  1:57720‑M1/50‑10 (MQ=255)
         gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                         <  1:365084‑M1/50‑11 (MQ=255)
         gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                         <  1:2369718‑M1/50‑11 (MQ=255)
         gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                         <  1:2249528‑M1/50‑11 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    >  1:3815626‑M1/1‑35 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaa                                     <  1:1346225‑M1/49‑15 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaa                                     <  1:1806190‑M1/49‑15 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaa                                     <  1:2384995‑M1/49‑15 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagca                                      <  1:1520669‑M1/48‑14 (MQ=255)
               ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    <  1:4107302‑M1/49‑16 (MQ=255)
               ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    <  1:2660081‑M1/49‑16 (MQ=255)
                tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                   >  1:2350841‑M1/1‑33 (MQ=255)
                 cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                  <  1:1072552‑M1/49‑18 (MQ=255)
                  aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaa                                 <  1:2540152‑M1/49‑19 (MQ=255)
                  aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaa                                 >  1:3282943‑M1/1‑31 (MQ=255)
                  aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaa                                 <  1:1741380‑M1/49‑19 (MQ=255)
                   cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaa                               <  1:269587‑M1/50‑21 (MQ=255)
                   cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                                >  1:598774‑M1/1‑30 (MQ=255)
                      aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:1901132‑M1/50‑24 (MQ=255)
                      aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:1190422‑M1/50‑24 (MQ=255)
                      aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:852425‑M1/50‑24 (MQ=255)
                      aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:1181689‑M1/50‑24 (MQ=255)
                      aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:833934‑M1/49‑23 (MQ=255)
                      aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:3534042‑M1/1‑27 (MQ=255)
                      aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:3723268‑M1/1‑27 (MQ=255)
                       gTTCTTAAACGTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           <  1:829283‑M1/50‑25 (MQ=255)
                        ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaa                          <  1:1829669‑M1/50‑26 (MQ=255)
                        ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaa                          <  1:2398564‑M1/50‑26 (MQ=255)
                        ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           >  1:728073‑M1/1‑25 (MQ=255)
                         tattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:2085470‑M1/1‑24 (MQ=255)
                          attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         <  1:449176‑M1/49‑27 (MQ=255)
                          attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:2911963‑M1/1‑23 (MQ=255)
                           ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:576498‑M1/1‑22 (MQ=255)
                             aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      <  1:2003026‑M1/49‑30 (MQ=255)
                             aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      <  1:3710355‑M1/49‑30 (MQ=255)
                                cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  <  1:929511‑M1/50‑34 (MQ=255)
                                cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   >  1:231448‑M1/1‑17 (MQ=255)
                                 aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:810606‑M1/1‑16 (MQ=255)
                                  tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:699684‑M1/50‑36 (MQ=255)
                                    cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              <  1:597914‑M1/50‑38 (MQ=255)
                                    cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              <  1:1118876‑M1/50‑38 (MQ=255)
                                    cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              <  1:1896526‑M1/50‑38 (MQ=255)
                                    cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat               <  1:2257130‑M1/49‑37 (MQ=255)
                                      cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg             >  1:897555‑M1/1‑11 (MQ=255)
                                       gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaa           <  1:1964992‑M1/50‑41 (MQ=255)
                                       gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaa           <  1:2492677‑M1/50‑41 (MQ=255)
                                       gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga            >  1:2701114‑M1/1‑10 (MQ=255)
                                       gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga            >  1:3841840‑M1/1‑10 (MQ=255)
                                        cgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa          <  1:4243240‑M1/50‑42 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:1407066‑M1/1‑8 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:1712272‑M1/1‑8 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:414398‑M1/1‑8 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         <  1:2645832‑M1/50‑43 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:1200744‑M1/1‑8 (MQ=255)
                                           gcgcGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg       <  1:4326172‑M1/50‑45 (MQ=255)
                                            cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg      <  1:227848‑M1/50‑46 (MQ=255)
                                            cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg      <  1:476104‑M1/50‑46 (MQ=255)
                                            cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg       >  1:2712475‑M1/1‑5 (MQ=255)
                                             gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     >  1:1002051‑M1/1‑4 (MQ=255)
                                             gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     <  1:1337696‑M1/50‑47 (MQ=255)
                                              cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    <  1:1955910‑M1/50‑48 (MQ=255)
                                              cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    >  1:328716‑M1/1‑3 (MQ=255)
                                                tcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt  <  1:2701281‑M1/50‑50 (MQ=255)
                                                |                                                 
AGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962517‑1962614

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