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	90 [0]	[0] 89	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)	69 (1.570)	24/112	0.2	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:3871892‑M2/64‑65 (MQ=255) gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                               >  1:2059998‑M2/64‑65 (MQ=255)    gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                            <  1:2867796‑M2/5‑1 (MQ=255)     cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                           <  1:941380‑M2/6‑1 (MQ=255)     cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                           <  1:3584716‑M2/6‑1 (MQ=255)      gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                          >  1:658308‑M2/59‑65 (MQ=255)      gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                          <  1:2888520‑M2/7‑1 (MQ=255)      gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                          <  1:1715844‑M2/7‑1 (MQ=255)            caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                    <  1:1984203‑M2/13‑1 (MQ=255)            caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                    <  1:1433760‑M2/13‑1 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                <  1:2022796‑M2/17‑1 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                <  1:573923‑M2/17‑1 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                >  1:3192392‑M2/49‑65 (MQ=255)                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                >  1:3093513‑M2/49‑65 (MQ=255)                      tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                          <  1:3251279‑M2/23‑1 (MQ=255)                      tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                          <  1:928369‑M2/23‑1 (MQ=255)                       gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                         <  1:700233‑M2/24‑1 (MQ=255)                       gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                         <  1:3171833‑M2/24‑1 (MQ=255)                        atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                        <  1:1213094‑M2/25‑1 (MQ=255)                        atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                        <  1:213833‑M2/25‑1 (MQ=255)                         tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                       >  1:2022102‑M2/40‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:2238571‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:1645845‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      <  1:709530‑M2/27‑1 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      <  1:919861‑M2/27‑1 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:3264394‑M2/39‑65 (MQ=255)                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                      >  1:112116‑M2/39‑65 (MQ=255)                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                    <  1:1987355‑M2/29‑1 (MQ=255)                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                    <  1:56465‑M2/29‑1 (MQ=255)                               tcacacaaaacactcgattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:1630414‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:1954267‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:55393‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:200848‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:2942215‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:2921484‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:616310‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:1775429‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:1748426‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:3997077‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:140470‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:3681100‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:2648246‑M2/34‑65 (MQ=255)                               tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                 >  1:2724940‑M2/34‑65 (MQ=255)                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                             >  1:447955‑M2/30‑65 (MQ=255)                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                             >  1:3740237‑M2/30‑65 (MQ=255)                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                             >  1:242649‑M2/30‑65 (MQ=255)                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                             >  1:1346639‑M2/30‑65 (MQ=255)                                    caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                            >  1:2685420‑M2/29‑65 (MQ=255)                                    caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                            >  1:2930052‑M2/29‑65 (MQ=255)                                    caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                            >  1:1733‑M2/29‑65 (MQ=255)                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                          <  1:2282822‑M2/39‑1 (MQ=255)                                        acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                        >  1:858962‑M2/25‑65 (MQ=255)                                        acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                        >  1:3733253‑M2/25‑65 (MQ=255)                                               attacttaacataaatgTATAAGTCATGCTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                 >  1:3933412‑M2/18‑65 (MQ=255)                                               attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                 >  1:1755401‑M2/18‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAACGttttt              >  1:1570877‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:137666‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:2965880‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:1930555‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:369939‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:882887‑M2/15‑65 (MQ=255)                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt              >  1:3534203‑M2/15‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:1066406‑M2/13‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:3366421‑M2/13‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:1188974‑M2/13‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:1439219‑M2/13‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:2910834‑M2/13‑65 (MQ=255)                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg            >  1:1938105‑M2/13‑65 (MQ=255)                                                      aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct          <  1:2498562‑M2/55‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:1912630‑M2/57‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:1339905‑M2/57‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:1807244‑M2/57‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:35810‑M2/57‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:2056114‑M2/57‑1 (MQ=255)                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac        <  1:2355733‑M2/57‑1 (MQ=255)                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg       >  1:394378‑M2/8‑65 (MQ=255)                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg       >  1:3731364‑M2/8‑65 (MQ=255)                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg       >  1:1120763‑M2/8‑65 (MQ=255)                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg       >  1:1041733‑M2/8‑65 (MQ=255)                                                          taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc      >  1:604334‑M2/7‑65 (MQ=255)                                                          taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc      >  1:1525915‑M2/7‑65 (MQ=255)                                                              tgTATAAGTCATGCTTTTGTTTTGGGGGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:508066‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:2444938‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:1782095‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:3252079‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:3280832‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:1183038‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:971901‑M2/3‑65 (MQ=255)                                                              tgTATAAGTCATACTTTTGTTTTGGGCGGATTTCACTCTGTTAAAAAGTTTTTCGCTACGCTAGc  >  1:1570851‑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 < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑