Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp	100%		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503	1979278	776	45 [0]	[0] 46	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	41 (0.900)	35/260	NT	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913/1979133‑1979427                                                                                                                                                   |                                                                                                                                                     aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    <  1:512660‑M2/3‑1 (MQ=255)     gggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcaggtcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                                                                                <  2:124573‑M2/7‑1 (MQ=255)         ctttacacttatcgaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                            <  1:509180‑M2/11‑1 (MQ=255)          tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                           >  1:337939‑M2/138‑149 (MQ=255)           ttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacgacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                                                                          <  1:580267‑M2/13‑1 (MQ=255)             acacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttt                                                                                                                                        <  2:26491‑M2/15‑1 (MQ=255)             acacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttt                                                                                                                                        <  1:659234‑M2/15‑1 (MQ=255)               acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                                                                      >  1:457104‑M2/133‑149 (MQ=255)                       aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattaattaacataaatgTATAAGTCATACTTTTGTTTTgtgt                                                                                                                              >  1:10451‑M2/125‑149 (MQ=255)                           tttttcctaaatcgacgcaactgtactcgtcactacacgtacatacaacggaggggggctgcgatttttaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCTTATTTTTGTTTTGGGTGTAt                                                                                                                          >  2:427882‑M2/121‑149 (MQ=255)                                  taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                   >  2:711149‑M2/114‑149 (MQ=255)                                  taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                   >  2:465478‑M2/114‑149 (MQ=255)                                  taaatcgacgcaactgtactcgtcaatgaaagataatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATACTTTTGTTTTCGGTGTATTTGAATc                                                                                                                   >  2:825646‑M2/114‑149 (MQ=255)                                          cgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                                                                                           <  1:653736‑M2/44‑1 (MQ=255)                                                tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                     >  1:533787‑M2/100‑149 (MQ=255)                                                 gtactcctcacgacacctccatacaaccgaggggcggtcccattttcaataatgcgtcattcagatcacacaaaaccctcaattacttaacataaatgTATAATTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAgttttc                                                                                                    <  2:254188‑M2/51‑2 (MQ=255)                                                  tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                                   <  2:337939‑M2/52‑1 (MQ=255)                                                   actcgtcaccacacgcacatacaacggagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                                  <  1:235915‑M2/53‑1 (MQ=255)                                                     gcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                                >  2:632263‑M2/95‑149 (MQ=255)                                                        tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                             >  1:345356‑M2/92‑149 (MQ=255)                                                         cactacacgaacatagaaaggagggaggctgagattttcaataatgcgtgatgcagatcacacaaaacactcaatgacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTAAAAACAAAAAAAAAGTTTTTcgctacgc                                                                                            <  2:665697‑M2/59‑1 (MQ=255)                                                               acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacacccaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTcgcaa                                                                                      >  1:93685‑M2/85‑149 (MQ=255)                                                                     atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                                >  2:33781‑M2/79‑149 (MQ=255)                                                                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                               <  1:613405‑M2/72‑1 (MQ=255)                                                                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaaaaatcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTTCGCTAGCAAGCaaaaa                                                                               >  2:465756‑M2/78‑149 (MQ=255)                                                                         aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaattTATAAGTCATACGTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                            <  2:290563‑M2/75‑1 (MQ=255)                                                                         aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTATCAAGCAAAAATGa                                                                            >  2:275739‑M2/75‑149 (MQ=255)                                                                                   ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                                  >  1:145026‑M2/65‑149 (MQ=255)                                                                                    gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                                                                 >  1:297346‑M2/64‑149 (MQ=255)                                                                                      tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGgaataatg                                                               >  2:252795‑M2/62‑148 (MQ=255)                                                                                               caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                                      >  1:826622‑M2/53‑149 (MQ=255)                                                                                                aataatgcgtgatgcagatcacacaaaacactcgattacttaacataaatgTATAAGTCATACTTTTGTCATGGGTGTATTTCGATCTGTTAAAAAATTTTTCGCTACGCTACCAAGCAAAAGTGAAAAAGGAATAAACGAAATGGGAt                                                     >  2:244867‑M2/52‑149 (MQ=255)                                                                                                aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                                     <  2:488511‑M2/98‑1 (MQ=255)                                                                                                     tccgggatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTACAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                                <  2:467435‑M2/103‑1 (MQ=255)                                                                                                          gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           >  1:216041‑M2/42‑149 (MQ=255)                                                                                                             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                                        <  2:379128‑M2/111‑1 (MQ=255)                                                                                                             gcagatcacacaaaacactcaattacttaacataaataTATGAGTCATACTTTTGTTTTGGGTGTATTTCAATCTATTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAAACGAAATGGGATGTTGCGCACAGTc                                        >  1:418805‑M2/39‑149 (MQ=255)                                                                                                                    acacaaaacactcaattccttaacattactgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTATTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                                 <  1:587511‑M2/118‑1 (MQ=255)                                                                                                                    acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCCAAATaa                                 >  1:147893‑M2/32‑149 (MQ=255)                                                                                                                              ctcaattactaaacataaatcTATAAGTCATACTGGTGTTTGGGGTTTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                       <  1:404529‑M2/128‑1 (MQ=255)                                                                                                                               tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                      >  1:516722‑M2/21‑149 (MQ=255)                                                                                                                                   ttatttaacataaatgGAGAAGTCATACTTTTGTTTGGGGTGTATTTCAATCTGTTAAAAAGTGTCTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCgtaaatcaa                  <  2:35129‑M2/133‑4 (MQ=255)                                                                                                                                       ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc              >  2:240584‑M2/13‑149 (MQ=255)                                                                                                                                          acataaatgTATCAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGc           <  2:722244‑M2/140‑1 (MQ=255)                                                                                                                                                atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa     >  1:844951‑M2/4‑149 (MQ=255)                                                                                                                                                   tatgAGTCATACTTTGTTTTGGGGTGTAGTTCAATCTGTAAAAAAGTATGTCGCTACGCGAGCAAGCAAAAATGAAACAGGAACATTCGAAATGAGATGTTGCGCACAGTCAAAATAACGCACCGTAAATAATCATCTGCTataaataa  <  2:26957/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913/1979133‑1979427

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: ‑