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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979270	1979278	9–776	56 [0]	[0] 57	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)	53 (1.140)	44/268	0.1	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

GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979137‑1979421                                                                                                                                               |                                                                                                                                               gggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacaccaaccactcaattccttaacataaatgt                                                                                                                                                >  1:792233‑M2/143‑143 (MQ=255)      tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacatcaatgTATAAg                                                                                                                                           >  1:770143‑M2/138‑143 (MQ=255)      tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                           >  2:631420‑M2/138‑143 (MQ=255)      tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                           >  1:180635‑M2/138‑143 (MQ=255)           acttatctaagatttttcctaaatcgacgcaactgtactcgtcactaccctcacatacaacggaggggggctgcgattttcaaaactgcgccatacagatcgcacaagacactcaattacttaacataaatgCATAactccac                                                                                                                                      >  2:592349‑M2/133‑140 (MQ=255)            cttatctaagaattttcctaaagcgacgcaactgtacgcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                     <  1:670223‑M2/12‑1 (MQ=255)                tctaagtttcttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggtgctgcgattttcaataattcgtgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                 <  2:479111‑M2/16‑1 (MQ=255)                  taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActtttga                                                                                                                               >  2:37508‑M2/126‑142 (MQ=255)                   aagtcttttcgcaaatcgacgcaacgggagtcgtcaccacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacgctcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                                                              <  2:476286‑M2/19‑1 (MQ=255)                     gatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                            <  1:52947‑M2/21‑1 (MQ=255)                      attttgtctaagtcgacgcaactgtactcgtcactacatgcacatacaactgagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                           <  2:247399‑M2/22‑1 (MQ=255)                       tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  2:437300‑M2/23‑1 (MQ=255)                         tttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataactgTATAATTCATACTTTTGTTTtggga                                                                                                                        >  1:752136‑M2/119‑142 (MQ=255)                             ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                    >  1:830833‑M2/115‑143 (MQ=255)                              taaatcgacgcaactgtactcgtcactacccgcacatacgacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacacgcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                   >  1:38738‑M2/114‑143 (MQ=255)                              taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacaaaaggatgagaagtacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                   >  1:659607‑M2/114‑143 (MQ=255)                               aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                  >  1:340552‑M2/113‑143 (MQ=255)                               aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                  <  2:827375‑M2/31‑1 (MQ=255)                                 atcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                                >  1:890834‑M2/111‑143 (MQ=255)                                   cgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                              <  1:156527‑M2/35‑1 (MQ=255)                                    gacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                             >  1:624808‑M2/108‑143 (MQ=255)                                     acgcaactgtactcgtcactacacgcacataccacggcgcgaccctgcgattttcaataatacgtgatgcagatcacacaaaacactcaattacttaacataactgTATAAGTCCTACTTCTGCTTTCGGTCTATTTCAATCt                                                                                                            >  2:404781‑M2/107‑143 (MQ=255)                                        caactgtacccggcactacacgcacacacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                         <  1:93863‑M2/40‑1 (MQ=255)                                         aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaaccctcaattccttaacataaatgTATAAGTCATACTTTTGTTTTTGTTCTATTTCCATCTcttc                                                                                                        >  1:454541‑M2/103‑142 (MQ=255)                                            tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                                                                                     >  1:304315‑M2/100‑143 (MQ=255)                                             gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                    <  2:1510‑M2/45‑1 (MQ=255)                                                 tcgtcacgacacgaacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                <  1:86543‑M2/49‑1 (MQ=255)                                                        tccacgcacatagcactgagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                         <  2:771060‑M2/56‑1 (MQ=255)                                                        tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattccttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTtcgctc                                                                                         >  2:269955‑M2/88‑142 (MQ=255)                                                        tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                         >  1:65340‑M2/88‑143 (MQ=255)                                                           acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                      <  2:431983‑M2/59‑1 (MQ=255)                                                               acatacagcggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacgctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                  <  2:289857‑M2/63‑1 (MQ=255)                                                                   acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                                                              <  2:850917‑M2/67‑1 (MQ=255)                                                                     aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                            >  1:526829‑M2/75‑143 (MQ=255)                                                                         gaggggggctgtgattttcaataatgcgtgatgcagatcacacaaaccactccattacttatcataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                        <  2:121088‑M2/73‑1 (MQ=255)                                                                         gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaaccctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                        <  1:512862‑M2/73‑1 (MQ=255)                                                                                    cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                             >  2:130510‑M2/60‑143 (MQ=255)                                                                                    cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                             >  1:158780‑M2/60‑143 (MQ=255)                                                                                      attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                                                           >  1:274349‑M2/58‑143 (MQ=255)                                                                                       ttttcaataatgcgtgatgcagatcacactaatcactccagtacctaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                          <  1:165425‑M2/87‑1 (MQ=255)                                                                                          tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa                                                       <  1:695457‑M2/90‑1 (MQ=255)                                                                                           caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa                                                      >  1:346399‑M2/53‑143 (MQ=255)                                                                                                 ttggtgatggagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                                <  1:727726‑M2/97‑1 (MQ=255)                                                                                                 tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAACAATGAAACAGGAATAATCGAAATGGGa                                                >  2:439666‑M2/47‑143 (MQ=255)                                                                                                             atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                    >  1:152866‑M2/35‑143 (MQ=255)                                                                                                              ccacacaaaacactcaattacttaacataaacgTATAAGTCATACTTTTGTTTGGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                   <  1:244337‑M2/110‑1 (MQ=255)                                                                                                                      aacactcaattacttagcataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                           >  2:744372‑M2/26‑143 (MQ=255)                                                                                                                             aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGAGCACAGTCAAAATAACTCACCg                    >  2:3935‑M2/19‑143 (MQ=255)                                                                                                                                   ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata              >  1:437141‑M2/13‑143 (MQ=255)                                                                                                                                   ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATATTGCGCACAGTCAAAATAAGTCACCGTCaata              >  1:461129‑M2/13‑143 (MQ=255)                                                                                                                                   ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAACCAAAATGGGATCGCGCACACAGTCAAAATAACTCACCGTAaata              >  2:434275‑M2/13‑143 (MQ=255)                                                                                                                                     aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataat            >  1:15189‑M2/11‑143 (MQ=255)                                                                                                                                      acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc           >  1:185482‑M2/10‑143 (MQ=255)                                                                                                                                        ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGTGATGTTGCGCACAGTCAAAATAACTCACCGTCCATAatcat         >  2:152189‑M2/8‑143 (MQ=255)                                                                                                                                           aatgTATAAGTCATCCTTGTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATTTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg      <  1:305980‑M2/139‑1 (MQ=255)                                                                                                                                             tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCt    <  1:137680‑M2/141‑1 (MQ=255)                                                                                                                                               tataAGTCATACTTTTGTTTTGAGTGTAAAAAAAAAAAAAAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtat  <  2:659607/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979137‑1979421

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