Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913_3_pae_pgi 1,978,503 Δ776 bp insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_pae_pgi 1978503–1979269 1979278 10–776 52 [0] [0] 52 insB1–insA insB1,insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_pgi = 19785020 (0.000)48 (0.830) 39/264 0.6 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913_3_pae_pgi 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

AGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_pae_pgi/1978355‑1978648
                                                                                                                                                   |                                                                                                                                                  
aGAAATGGGGCTTGACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGt                                                                                                                                                   <  2:147365‑M1/149‑2 (MQ=255)
                acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                                                                   >  1:974914‑M1/1‑132 (MQ=255)
                  ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACCCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACTTAAATGtatcagtcatacttttttt                                                                                                                                 >  2:181630‑M1/1‑130 (MQ=255)
                      cTAAGATTTTTCCTAAAGCGACGCAATTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  2:852470‑M1/149‑24 (MQ=255)
                      cTAAGATCTTTCCTAAACCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  1:181630‑M1/149‑24 (MQ=255)
                            tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               >  1:153865‑M1/1‑120 (MQ=255)
                            tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               <  2:153865‑M1/125‑6 (MQ=255)
                             ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGCGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTACCATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                      >  1:704245‑M1/1‑119 (MQ=255)
                                     aaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                                                                      <  2:753007‑M1/125‑15 (MQ=255)
                                     aaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                                                                      >  1:753007‑M1/1‑111 (MQ=255)
                                       tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                            >  1:277241‑M1/1‑109 (MQ=255)
                                       tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                            >  1:790130‑M1/1‑109 (MQ=255)
                                       tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                            >  2:894079‑M1/1‑109 (MQ=255)
                                                  gTACTCGTCACTACACGCTCATAAAACGAAGTGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                                 <  1:196744‑M1/149‑52 (MQ=255)
                                                        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                   >  2:969362‑M1/1‑92 (MQ=255)
                                                        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                   <  1:969362‑M1/141‑50 (MQ=255)
                                                        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                              <  1:453156‑M1/114‑23 (MQ=255)
                                                        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                              >  2:453156‑M1/1‑92 (MQ=255)
                                                        gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATAACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtttatttcaatctgttaaaaagtttttcgcttc                                                                                           >  2:522279‑M1/1‑92 (MQ=255)
                                                                acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  2:115962‑M1/107‑24 (MQ=255)
                                                                acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             >  1:115962‑M1/1‑84 (MQ=255)
                                                                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                         >  1:620579‑M1/1‑74 (MQ=255)
                                                                             ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgatacgctaacaagcaaaaatgaaac                                                                      >  1:122556‑M1/1‑71 (MQ=255)
                                                                               aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                             >  1:799208‑M1/1‑69 (MQ=255)
                                                                               aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                             <  2:799208‑M1/108‑40 (MQ=255)
                                                                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaataat                                                             >  2:874378‑M1/1‑62 (MQ=255)
                                                                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                            >  2:263798‑M1/1‑61 (MQ=255)
                                                                                        gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                           >  2:765415‑M1/1‑60 (MQ=255)
                                                                                         cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                          >  1:278426‑M1/1‑59 (MQ=255)
                                                                                            ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                       >  1:864345‑M1/1‑56 (MQ=255)
                                                                                              tttacaaatgcgtgaACGAAATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgacacaggaataatcgaaatgg                                                     <  2:511377‑M1/143‑96 (MQ=255)
                                                                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                             >  2:477039‑M1/1‑51 (MQ=255)
                                                                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                             <  1:477039‑M1/90‑40 (MQ=255)
                                                                                                  ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaaaggaataatcgaaatgggatg                                                 >  1:740423‑M1/1‑50 (MQ=255)
                                                                                                    aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaaaggaataatagaaatgggatgtt                                               >  2:153063‑M1/1‑48 (MQ=255)
                                                                                                       gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                       <  1:584112‑M1/58‑14 (MQ=255)
                                                                                                       gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                       >  2:584112‑M1/1‑45 (MQ=255)
                                                                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                      <  2:479158‑M1/149‑111 (MQ=255)
                                                                                                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagaaataatcgaaatgggctattgcgcacagcc                                     >  2:565663‑M1/1‑38 (MQ=255)
                                                                                                              gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagtaataatcgaaatgggatgttgcgcacagtc                                     >  2:530125‑M1/1‑38 (MQ=255)
                                                                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaagtgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                                  >  1:603469‑M1/1‑35 (MQ=255)
                                                                                                                  aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                 >  2:441567‑M1/1‑34 (MQ=255)
                                                                                                                   tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                <  1:263798‑M1/149‑117 (MQ=255)
                                                                                                                           aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                        <  1:643800‑M1/149‑125 (MQ=255)
                                                                                                                                tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcccagtccaaataactcacgcaaaa                   >  2:166112‑M1/1‑20 (MQ=255)
                                                                                                                                  caatacaTAACATCAATGtataaatcatacttttgtttttgctgtatttccatctgttaaaaagtttttcgctacgctagcaagcaaaaatcaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                 <  2:575249‑M1/146‑132 (MQ=255)
                                                                                                                                     tACTAAATATAGATGtataagtcataattttgttttgggtttatttcaatctgttaaaaagtttttcgctaagaaaaaaaacaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  2:122556‑M1/149‑135 (MQ=255)
                                                                                                                                          aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg         >  1:160159‑M1/1‑10 (MQ=255)
                                                                                                                                            cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaattgggatgttgcgcacagtcaaaatacctcacagtaaataatcatctgct       >  2:765949‑M1/1‑8 (MQ=255)
                                                                                                                                            cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct       >  1:254194‑M1/1‑8 (MQ=255)
                                                                                                                                               aaaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctata    <  2:864345‑M1/149‑145 (MQ=255)
                                                                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  <  1:677254‑M1/149‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_pae_pgi/1978355‑1978648

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