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

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

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

GGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAAT  >  NC_000913_3_pae_pgi/1978362‑1978649                                                                                                                                             |                                                                                                                                                    ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                             >  1:141937‑M1/1‑141 (MQ=255)         cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATCCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAACAACTCAATTACTTAACATAAATGtatacgtcatactttt                                                                                                                                     >  1:615255‑M1/1‑133 (MQ=255)         cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                     <  1:2087‑M1/149‑17 (MQ=255)              aTCTAAGATTTTTGCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                                <  2:810313‑M1/149‑22 (MQ=255)               tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                               >  1:106349‑M1/1‑127 (MQ=255)               tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                               >  1:805386‑M1/1‑127 (MQ=255)                 tAAGATTTTTCCTAAATCGACGCAACTGAACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACGTAACGTAACTGtataattcatcctgttgttttggg                                                                                                                             >  1:219119‑M1/1‑125 (MQ=255)                       ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                       >  1:54185‑M1/1‑119 (MQ=255)                        tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                      <  2:345593‑M1/149‑32 (MQ=255)                         ttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                                                                     >  1:681415‑M1/1‑117 (MQ=255)                           ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                                                                   >  1:824630‑M1/1‑115 (MQ=255)                            gcaatgcgacgcaACTGTACTCGTCACTACACGCATATACAACGGAGGGGGGCTGGGATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTCCTTAACATAAATGtataagtcatacttttgttttgggggtatttcaat                                                                                                                  <  2:219119‑M1/143‑36 (MQ=255)                            cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                          <  2:497982‑M1/125‑12 (MQ=255)                            cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                          >  1:497982‑M1/1‑114 (MQ=255)                             tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaatc                                                                                                                 >  1:107168‑M1/1‑113 (MQ=255)                                aTCGAAGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGAGCTGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                              <  1:272419‑M1/149‑40 (MQ=255)                                    aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTACCATAAATGtataagccatacttttgttttgggtggatttcaatctgttaaa                                                                                                          >  1:569116‑M1/1‑106 (MQ=255)                                      gCAACTGTATTCGTCACTACACGCAACTACAACGACGGGGGGCTGCGATTTTCAATAATGCTTGAGTCAGATCACACAAAACACTCAATGACGGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                        <  2:812848‑M1/149‑46 (MQ=255)                                        aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                      <  2:824630‑M1/149‑48 (MQ=255)                                           tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                                                                       <  1:276640‑M1/145‑47 (MQ=255)                                           tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                                                                       >  2:276640‑M1/1‑99 (MQ=255)                                             tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttc                                                                                                 >  1:128106‑M1/1‑97 (MQ=255)                                                  gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                            >  1:344408‑M1/1‑92 (MQ=255)                                                  gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                            >  1:789809‑M1/1‑92 (MQ=255)                                                          acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                    <  1:922311‑M1/149‑66 (MQ=255)                                                            gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                  >  2:141614‑M1/1‑82 (MQ=255)                                                                    aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgttaaaaaatttttcgctacgctcgcaatcaaaaatga                                                                          >  2:853119‑M1/1‑74 (MQ=255)                                                                       ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaaaaatgaaac                                                                       >  1:810691‑M1/1‑71 (MQ=255)                                                                       ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCAAACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgagtgtatttcaatctgttaaaaagtttttcgctacgctaacaagcaaaaatgaaac                                                                       >  2:349306‑M1/1‑71 (MQ=255)                                                                        gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaca                                                                      >  1:884961‑M1/1‑70 (MQ=255)                                                                         tgcggggctgCGATGTGAAATAATGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttggtttgggtggatttgaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                     <  2:609144‑M1/146‑81 (MQ=255)                                                                          ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                                                                    >  2:19881‑M1/1‑68 (MQ=255)                                                                           ggacgctgcGATTTTCTATAAGGCGTTATGCAGATCACATAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacatga                                                                   <  2:64298‑M1/145‑83 (MQ=255)                                                                                         tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataaacgaaatgtg                                                     >  1:419717‑M1/1‑53 (MQ=255)                                                                                          cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                                    >  1:600559‑M1/1‑52 (MQ=255)                                                                                           aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                                   >  1:40191‑M1/1‑51 (MQ=255)                                                                                              aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                                <  2:106349‑M1/149‑102 (MQ=255)                                                                                                     gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                         <  1:62678‑M1/149‑109 (MQ=255)                                                                                                           gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcacgcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                                   >  1:352260‑M1/1‑35 (MQ=255)                                                                                                               acacaAAACACTGAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttccgctacgctagcacgcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                               <  1:641112‑M1/149‑119 (MQ=255)                                                                                                               aaacaaaACACTCAATTACTTAACATAAATTtataagtcatagttttgttttgcgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                               <  2:128106‑M1/147‑119 (MQ=255)                                                                                                                 acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                             >  2:45574‑M1/1‑29 (MQ=255)                                                                                                                     aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                                         <  2:488528‑M1/85‑61 (MQ=255)                                                                                                                     aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                                         >  1:488528‑M1/1‑25 (MQ=255)                                                                                                                        acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta                      >  1:502754‑M1/1‑22 (MQ=255)                                                                                                                        acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta                      >  1:502749‑M1/1‑22 (MQ=255)                                                                                                                         cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                     >  1:438938‑M1/1‑21 (MQ=255)                                                                                                                                 cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                 <  1:317817‑M1/129‑117 (MQ=255)                                                                                                                                 cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcacacagtcaaaat                                 >  2:317817‑M1/1‑13 (MQ=255)                                                                                                                                       aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaatactcatctgcta       >  2:227479‑M1/1‑7 (MQ=255)                                                                                                                                       aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaaaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtcaataatcatctgcta       >  2:671583‑M1/1‑7 (MQ=255)                                                                                                                                        tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaactttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccctaataaatcatctcctat      >  2:405700‑M1/1‑6 (MQ=255)                                                                                                                                            tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaat  <  2:805386‑M1/149‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                    GGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAAT  >  NC_000913_3_pae_pgi/1978362‑1978649

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