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	85 [0]	[0] 80	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)	68 (1.450)	28/118	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

CCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAG  >  NC_000913/1979212‑1979344                                                                    |                                                                  gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                   <  1:1760820‑M2/1‑1 (MQ=255) gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                   <  1:2137920‑M2/1‑1 (MQ=255) gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgt                                                                   <  1:3641421‑M2/1‑1 (MQ=255)  ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                  >  1:3564491‑M2/67‑68 (MQ=255)  ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                  >  1:506981‑M2/67‑68 (MQ=255)  ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                  >  1:2017766‑M2/67‑68 (MQ=255)  ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                  >  1:2444338‑M2/67‑68 (MQ=255)  ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                  >  1:2931359‑M2/67‑68 (MQ=255)  ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                  >  1:933196‑M2/67‑68 (MQ=255)    ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                >  1:1786416‑M2/65‑68 (MQ=255)     gctgcgattttccataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                               >  1:2554660‑M2/64‑68 (MQ=255)     gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                               >  1:138734‑M2/64‑68 (MQ=255)     gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                               >  1:500266‑M2/64‑68 (MQ=255)      ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                              <  1:1613140‑M2/6‑1 (MQ=255)      ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                              <  1:3395938‑M2/6‑1 (MQ=255)       tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                             >  1:2345969‑M2/62‑68 (MQ=255)             tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                       <  1:1405816‑M2/13‑1 (MQ=255)                 aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                   <  1:1682138‑M2/17‑1 (MQ=255)                    aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                >  1:1940464‑M2/49‑68 (MQ=255)                        cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                            <  1:1654549‑M2/24‑1 (MQ=255)                         gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                           <  1:740164‑M2/25‑1 (MQ=255)                         gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                           <  1:1320655‑M2/25‑1 (MQ=255)                         gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                           <  1:956040‑M2/25‑1 (MQ=255)                           gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCNTACTTTAGTTTTGGgtgt                                         <  1:3093819‑M2/27‑1 (MQ=255)                             tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       <  1:941600‑M2/29‑1 (MQ=255)                             tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       <  1:553504‑M2/29‑1 (MQ=255)                             tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       <  1:3171131‑M2/29‑1 (MQ=255)                             tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                       <  1:2215847‑M2/29‑1 (MQ=255)                              gcagatcacacaaaacactcaattacttatcataaatgTATAAGTCATACTTTTGCTTTGGGTGTAtt                                      >  1:3591117‑M2/39‑68 (MQ=255)                              gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                      >  1:2680484‑M2/39‑68 (MQ=255)                              gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                      >  1:3762938‑M2/39‑68 (MQ=255)                              gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                      >  1:1958089‑M2/39‑68 (MQ=255)                                agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                    >  1:3337903‑M2/37‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:1653909‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:3758831‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:1646669‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:1644574‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:2523531‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:1729542‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:1804633‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:3251915‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                 >  1:3693614‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTACTTCAAt                                 >  1:2814318‑M2/34‑68 (MQ=255)                                   tcacacaaaacactcaattacttaacataaatgTATAAGACATACTTCTGTCCTGGCTGTATTtcact                                 >  1:3578388‑M2/34‑66 (MQ=255)                                     acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                               >  1:3721411‑M2/32‑68 (MQ=255)                                     acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                               >  1:2932412‑M2/32‑68 (MQ=255)                                       acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                             >  1:322320‑M2/30‑68 (MQ=255)                                       acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                             <  1:1459957‑M2/39‑1 (MQ=255)                                       acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                             >  1:2739543‑M2/30‑68 (MQ=255)                                       acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                             >  1:1751830‑M2/30‑68 (MQ=255)                                       acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                             >  1:3491295‑M2/30‑68 (MQ=255)                                        caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                            >  1:65954‑M2/29‑68 (MQ=255)                                            acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                        >  1:1101676‑M2/25‑68 (MQ=255)                                            acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                        >  1:3771944‑M2/25‑68 (MQ=255)                                            acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                        >  1:3042069‑M2/25‑68 (MQ=255)                                            acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                        >  1:1348961‑M2/25‑68 (MQ=255)                                             cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCGATCTGTTaaaaa                       >  1:3593599‑M2/24‑68 (MQ=255)                                               ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                     >  1:610374‑M2/22‑68 (MQ=255)                                               ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                     >  1:1962418‑M2/22‑68 (MQ=255)                                                 caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAaaagcct                   >  1:1778232‑M2/20‑65 (MQ=255)                                                   attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                 >  1:3854686‑M2/18‑68 (MQ=255)                                                      acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc              >  1:912632‑M2/15‑68 (MQ=255)                                                      acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc              >  1:70598‑M2/15‑68 (MQ=255)                                                      acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc              >  1:1751954‑M2/15‑68 (MQ=255)                                                       cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct             <  1:1437735‑M2/55‑1 (MQ=255)                                                        ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta            >  1:1878599‑M2/13‑68 (MQ=255)                                                        ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta            >  1:58138‑M2/13‑68 (MQ=255)                                                         taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac           <  1:721829‑M2/57‑1 (MQ=255)                                                             ataaatgTATAAGTCATANTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta       >  1:3098100‑M2/8‑68 (MQ=255)                                                             ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta       >  1:205100‑M2/8‑68 (MQ=255)                                                             ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta       >  1:2990215‑M2/8‑68 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAACGTTTTTCGCTACGctagcacg  >  1:3115628‑M2/3‑66 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCtagcaac  >  1:81951‑M2/3‑67 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  1:971013‑M2/3‑68 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  1:746685‑M2/3‑68 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  1:473136‑M2/3‑68 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  1:3631358‑M2/3‑68 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  1:1915260‑M2/3‑68 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  1:1538010‑M2/3‑68 (MQ=255)                                                                  tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag  >  1:1158194‑M2/3‑68 (MQ=255)                                                                    |                                                                  CCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAG  >  NC_000913/1979212‑1979344

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