Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp		insB5–insA5	insB5, insA5

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	24 [0]	[1] 25	insB5–insA5	insB5,insA5

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	24 (0.430)	24/474	1.1	100%	intergenic (‑305/+16)	flhD/insB5	DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA5/uspC	IS1 family protein InsA/universal stress protein C

TCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAAT  >  NC_000913/1979043‑1979528                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          ccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacgtaacataaattTATAAGTCATACttt                                                                                                                                                                                                                                             <  2:385821‑M2/15‑1 (MQ=255)                      attagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacgtaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                                                                                                                        <  1:215513‑M2/36‑1 (MQ=255)                           aatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                                                                                                                                    >  1:119628‑M2/211‑250 (MQ=255)                                            gcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcccaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                                                                                                                                   <  2:130455‑M2/57‑1 (MQ=255)                                                 aagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggagggggtctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                                                                                                                              <  1:303966‑M2/62‑1 (MQ=255)                                                    cagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgccactgtactcgtcactacacgcccatccaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                                                                                                                                                                            >  1:95784‑M2/186‑249 (MQ=255)                                                         actcctttttgggttaagatggccttacaaatagaaatgggtcgttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcacgacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacaatcaattacttaacataaatgTATAAGTCATACTTGTGGTTTGGGTGTATTTCAATCTGTGAAAAAGTGGTTCGCTACGCTAGCAAGCaaaa                                                                                                                                                                                     <  2:256723‑M2/71‑1 (MQ=255)                                                                                               gggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                                                                                                                               <  1:450757‑M2/109‑1 (MQ=255)                                                                                                                  agatttttcctaaatcgacgcaactgtactcgtcacgacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                                                                                                                            <  1:385462‑M2/128‑1 (MQ=255)                                                                                                                        ttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataat                                                                                                                       <  2:129336‑M2/133‑1 (MQ=255)                                                                                                                                                       acacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCaa                                                                                         <  1:322733‑M2/163‑1 (MQ=255)                                                                                                                                                        cacgcacatacaacggaggggggctgcgattttcactactgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGGTGTTGCGCACAGTCAGAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATAc                                                                                      >  2:118148‑M2/86‑251 (MQ=255)                                                                                                                                                                                     ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACtt                                                          >  1:97793‑M2/57‑250 (MQ=255)                                                                                                                                                                                      ttgcaaaaatgcttgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTGCAATCTGATAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACGGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAAGAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTa                                                         <  1:567629‑M2/195‑1 (MQ=255)                                                                                                                                                                                        tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTAcc                                                       <  1:482421‑M2/197‑1 (MQ=255)                                                                                                                                                                                             aatgcgtgatgcagatcacacaaaacactcaattacttaacatacatgAATAAGTCATACTTTTGTTTTGGGAGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGCATACTCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAACTAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGcaca                                                  >  1:32516‑M2/49‑250 (MQ=255)                                                                                                                                                                                                 cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCAcaacaa                                              >  2:535192‑M2/45‑250 (MQ=255)                                                                                                                                                                                                       gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAATTTTGTTTTGGGTGTATTTCAGTCTGTTACACAGTTTTTCGCTACGCTAGCAAGCAACAATGAAACAGGTATACTCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGCTAAGCTATTTTTAAAAAGACACTTACCGCACAACACACTAAtt                                       >  1:580555‑M2/39‑251 (MQ=255)                                                                                                                                                                                                          gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTaaa                                    <  2:97793‑M2/216‑1 (MQ=255)                                                                                                                                                                                                                  aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACaaaga                            <  1:16839‑M2/224‑1 (MQ=255)                                                                                                                                                                                                                         tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGca                      >  2:148224‑M2/21‑250 (MQ=255)                                                                                                                                                                                                                                cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCa               >  1:370668‑M2/14‑250 (MQ=255)                                                                                                                                                                                                                                 ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACaaag                             >  1:464736‑M2/13‑235 (MQ=255)                                                                                                                                                                                                                                 ttaacataaatgTATAAGGCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACaaag                             <  2:464736‑M2/223‑1 (MQ=255)                                                                                                                                                                                                                                            gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACttaatgttaat  <  2:32273/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          TCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAAT  >  NC_000913/1979043‑1979528

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