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–1979270	1979278	9–776	64 [0]	[0] 64	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)	52 (1.430)	23/112	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

TGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGC  >  NC_000913/1979216‑1979341                                                                |                                                               gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                               >  1:1269699‑M2/64‑65 (MQ=255) gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                               >  1:577241‑M2/64‑65 (MQ=255) gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                               >  1:527358‑M2/64‑65 (MQ=255)    gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                            >  1:472454‑M2/61‑65 (MQ=255)     cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                           <  1:1118923‑M2/6‑1 (MQ=255)     cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                           <  1:469409‑M2/6‑1 (MQ=255)            caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                    <  1:3192299‑M2/13‑1 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                >  1:257211‑M2/49‑65 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                >  1:83391‑M2/49‑65 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                >  1:1648601‑M2/49‑65 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                >  1:2416926‑M2/49‑65 (MQ=255)                 atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                               <  1:2977000‑M2/18‑1 (MQ=255)                     gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                           >  1:1274651‑M2/44‑65 (MQ=255)                     gtgatgcaaatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                           >  1:2672615‑M2/44‑65 (MQ=255)                        atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                        <  1:2223099‑M2/25‑1 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:252042‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:94552‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:401750‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:3237869‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:2061013‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:3199390‑M2/39‑65 (MQ=255)                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                    <  1:2289318‑M2/29‑1 (MQ=255)                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                    <  1:2965540‑M2/29‑1 (MQ=255)                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                    <  1:2529019‑M2/29‑1 (MQ=255)                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                    <  1:1034958‑M2/29‑1 (MQ=255)                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                    <  1:1045311‑M2/29‑1 (MQ=255)                             gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                   <  1:266512‑M2/30‑1 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:2136569‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:2056489‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:582923‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:960147‑M2/34‑65 (MQ=255)                                 acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                               >  1:76994‑M2/32‑65 (MQ=255)                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                             >  1:2912398‑M2/30‑65 (MQ=255)                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                             >  1:2688732‑M2/30‑65 (MQ=255)                                    caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                            >  1:1708644‑M2/29‑65 (MQ=255)                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                          <  1:2710506‑M2/39‑1 (MQ=255)                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                          <  1:980684‑M2/39‑1 (MQ=255)                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                          <  1:1040029‑M2/39‑1 (MQ=255)                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                          >  1:2212676‑M2/27‑65 (MQ=255)                                        acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                        >  1:3187115‑M2/25‑65 (MQ=255)                                         cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                       >  1:1464046‑M2/24‑65 (MQ=255)                                         cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                       >  1:669654‑M2/24‑65 (MQ=255)                                               attacttaacataaatgTATAAGTCATCCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                 >  1:1947133‑M2/18‑65 (MQ=255)                                               attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                 >  1:2274572‑M2/18‑65 (MQ=255)                                               attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                 >  1:1761563‑M2/18‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:2710704‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:767751‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:649703‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:2737449‑M2/15‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            <  1:1703783‑M2/53‑1 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:1520166‑M2/13‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:1254969‑M2/13‑65 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:193457‑M2/57‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:555288‑M2/57‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:1415410‑M2/57‑1 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:514981‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:1956753‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:1666670‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:2063952‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:2279711‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:2349000‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:834807‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:242576‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGAGCATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:1921893‑M2/3‑65 (MQ=255)                                                                |                                                               TGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGC  >  NC_000913/1979216‑1979341

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