Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB‑5–insA‑5	insB‑5, insA‑5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979269	1979278	10–776	39 [0]	[1] 41	insB‑5–insA‑5	insB‑5,insA‑5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	38 (0.930)	35/264	0.1	100%	intergenic (‑305/+16)	flhD/insB‑5	DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA‑5/uspC	IS1 protein InsA/universal stress protein C

CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913/1978365‑1978642                                                                                                                                          |                                                                                                                                             cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTCCTCGTCACTACACGCACATACAACGGAGTGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                         <  2:509434‑M1/143‑6 (MQ=255)            tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagacaaaccata                                                                                                                              >  2:183361‑M1/1‑127 (MQ=255)              tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGCTCCAACAAATCACTCAATTACTCAACATAAATGtttaagtcatacttttgt                                                                                                                            >  2:276589‑M1/1‑125 (MQ=255)                aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                                                          <  1:80039‑M1/143‑21 (MQ=255)                  aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                        <  1:768050‑M1/143‑23 (MQ=255)                   tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                       <  1:93871‑M1/143‑24 (MQ=255)                   tttttCCTAAATCCACGCAACCGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                       <  1:703599‑M1/143‑24 (MQ=255)                    tctgcctaaaGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                      <  1:585133‑M1/139‑25 (MQ=255)                     gttcCTAAAGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                     <  2:411700‑M1/142‑26 (MQ=255)                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACGTAAATGtataagtcatacttctggtttgggtgtatt                                                                                                                >  2:419486‑M1/1‑113 (MQ=255)                          tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                >  1:524131‑M1/1‑113 (MQ=255)                               cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                           <  1:295247‑M1/143‑36 (MQ=255)                                 aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggggttttttcactct                                                                                                         >  1:293345‑M1/1‑106 (MQ=255)                                   gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAGACAAAACACTCAATTACTTAACATAAATGtataagtcatgcttttgttttgggtgtatttcaatctgt                                                                                                       >  1:260879‑M1/1‑104 (MQ=255)                                    ctactGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                      <  1:663049‑M1/141‑41 (MQ=255)                                        tGTACTCGTCACTATATTGACATACAACGGAGGTGGGCTGCGATTTTCAATAATGCGTGATGCAGATTACACTTATCATTCTATTACTTAACATATATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                  <  2:200071‑M1/143‑45 (MQ=255)                                        tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttatgggtgtatttcaatctgttaaaa                                                                                                  >  1:744583‑M1/1‑99 (MQ=255)                                             tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                             >  1:45982‑M1/1‑94 (MQ=255)                                                            caTACAACGGAGGGGTGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                                                              <  1:188577‑M1/143‑65 (MQ=255)                                                                 aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatcctcttgttttggatgtatttcaatctgttaaaaagtttttcgctacactcgcacacaa                                                                         >  2:499659‑M1/1‑74 (MQ=255)                                                                 aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                         >  2:470623‑M1/1‑74 (MQ=255)                                                                    ggcgtggggctgcgATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                      <  2:117220‑M1/138‑73 (MQ=255)                                                                     gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                     <  1:744066‑M1/143‑74 (MQ=255)                                                                      gggggggcTGCGTGTTTCAATGATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                                                                    <  2:329682‑M1/142‑75 (MQ=255)                                                                                 gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcactacgctagcaagcaaaaatgaaacaggaata                                                         >  2:444441‑M1/1‑58 (MQ=255)                                                                                   ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                       <  2:37561‑M1/143‑88 (MQ=255)                                                                                        aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                  >  2:675888‑M1/1‑51 (MQ=255)                                                                                         ataatGCGTGATGCAGATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                 <  2:458960‑M1/143‑94 (MQ=255)                                                                                           aaggcggGAGGCGGATCACAGAAAACACGCAATTACGTAACATAAATGgatcagtcatacttgtgttttgggtgtatttcaagctggtaaaaagtttgtcgcgacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                               <  1:243120‑M1/140‑96 (MQ=255)                                                                                                 tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                         >  1:461704‑M1/1‑42 (MQ=255)                                                                                                  gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                        >  2:28454‑M1/1‑41 (MQ=255)                                                                                                   aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                       >  2:139842‑M1/1‑40 (MQ=255)                                                                                                     gCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                     >  2:573148‑M1/1‑38 (MQ=255)                                                                                                          tCAAACAAAACACTCAATTACGTAACATAAATGtataagtcatacttttgttgtgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                <  1:419486‑M1/143‑111 (MQ=255)                                                                                                             cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                             >  1:209737‑M1/1‑30 (MQ=255)                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttgtgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaggcaaaaatgaaacaggaataatcgaaatgggatgtggcgcacagtgaaa                            >  2:77606‑M1/1‑29 (MQ=255)                                                                                                                aaaaCACTCAATTACATAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                          <  1:618550‑M1/143‑117 (MQ=255)                                                                                                                  aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                        >  1:277518‑M1/1‑25 (MQ=255)                                                                                                                                        ctgtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc  <  2:509261‑M1/142‑141 (MQ=255)                                                                                                                                          |                                                                                                                                             CTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913/1978365‑1978642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑