Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_pae_tpiA	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_tpiA	1978503–1979269	1979278	10–776	41 [0]	[0] 41	insB1–insA	insB1,insA

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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_pae_tpiA/1978357‑1978644                                                                                                                                                  |                                                                                                                                               aaaTGGTGCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                             <  2:68266‑M1/149‑4 (MQ=255) aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATAGACGCAACTGTACTCGTCACTAACAGCACACACAACGGAGAGGGGCTGCGATTTTCAATAATGCGTGATGCACATCAAACAAAACACCCAATTACTTAACATCAATGtct                                                                                                                                             >  1:375102‑M1/1‑146 (MQ=255)  aaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                                                            >  2:285526‑M1/1‑145 (MQ=255)      ggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                        <  1:122584‑M1/149‑9 (MQ=255)       gTCTTTACACTTATCCAAGAGTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       <  2:158951‑M1/149‑10 (MQ=255)               acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCAGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTGCTTAACATAAATGtataagtcatacttttg                                                                                                                               <  1:374322‑M1/149‑18 (MQ=255)                   aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                           <  1:268485‑M1/149‑22 (MQ=255)                             tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                 <  1:749453‑M1/149‑32 (MQ=255)                                    aaccgacGCAACTGTACTCGGCACCACACGGAAATACAACGGAGGGCGGCGGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                                                                          <  1:175664‑M1/146‑39 (MQ=255)                                         aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     >  1:394158‑M1/1‑106 (MQ=255)                                                tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                              >  1:313591‑M1/1‑99 (MQ=255)                                                     tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                         >  1:96441‑M1/1‑94 (MQ=255)                                                               acGCGCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataattcatacttttgttttgggtgtatttctatctgttaaacagttattagctagtctagcaa                                                                               >  1:272881‑M1/1‑84 (MQ=255)                                                               acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                               >  1:387892‑M1/1‑84 (MQ=255)                                                                       aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                       >  1:744070‑M1/1‑76 (MQ=255)                                                                         aaCGGAGGGGCGCTGCTATTTTAAATAAGGAGTGAGGCAGATCACAAAAAACACTCAAGTACTTAACAAAAATGtataagtcattcttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                     <  2:444082‑M1/149‑76 (MQ=255)                                                                           cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                         <  1:102184‑M1/111‑40 (MQ=255)                                                                           cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                         >  2:102184‑M1/1‑72 (MQ=255)                                                                               ggggggCTGCGATTTTCAATAATCCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaacaatgaaacagg                                                               <  1:521499‑M1/149‑82 (MQ=255)                                                                                gggggCGGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                              <  1:479760‑M1/149‑83 (MQ=255)                                                                                    gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                          <  2:96441‑M1/149‑87 (MQ=255)                                                                                    gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                          <  2:387892‑M1/149‑87 (MQ=255)                                                                                    gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgataaaaagtttttcgctacgctagcaaccaaaaatgaaacaggaataa                                                          >  2:728451‑M1/1‑63 (MQ=255)                                                                                     cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                         <  1:568926‑M1/149‑88 (MQ=255)                                                                                         gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaaaagaaataatagaa                                                     >  2:344863‑M1/1‑58 (MQ=255)                                                                                          aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                    >  1:220642‑M1/1‑57 (MQ=255)                                                                                                   aatGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                           <  1:640004‑M1/149‑102 (MQ=255)                                                                                                       cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       >  1:583726‑M1/1‑44 (MQ=255)                                                                                                            tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                  >  1:585281‑M1/1‑39 (MQ=255)                                                                                                               aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                               <  1:715579‑M1/149‑114 (MQ=255)                                                                                                                    acaaaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                          <  2:356543‑M1/149‑119 (MQ=255)                                                                                                                      gcaaAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcccacagtcaaaataact                        >  2:462709‑M1/2‑29 (MQ=255)                                                                                                                        aaaaCAATCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                      <  2:641394‑M1/149‑123 (MQ=255)                                                                                                                             acTCAATTACTTAACATAAATGtataagtgatacttttgttttggctgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtgatgttgcgcacagtacaaatacctccccgta                 >  2:773472‑M1/1‑22 (MQ=255)                                                                                                                             acTCAATTACTTAACATAAATGtataagtcatacttttgttttggatgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtaatgttgcgcacagtcaaaataattcaccgta                 >  2:309306‑M1/1‑22 (MQ=255)                                                                                                                                       ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagcatttagcttctgtagcgtgctaaaatgaaacaagagtaattgaaatgggttgttgagcacagtcaaaataactcaccgtaaataatcatc       >  2:775800‑M1/1‑12 (MQ=255)                                                                                                                                        tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct      <  2:744070‑M1/149‑139 (MQ=255)                                                                                                                                         tatataaATGtataagtcatacttttgttttgggtttatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaacgaaacaagtatgttctagatgggattttgcgcacagtcaaaataactcaccgtaaataatcatctg     <  2:393441‑M1/146‑140 (MQ=255)                                                                                                                                         aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgatacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaactaatcatctg     >  1:22857‑M1/1‑10 (MQ=255)                                                                                                                                         aaCATAAATGtataagtcatactatcgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg     <  1:644581‑M1/149‑140 (MQ=255)                                                                                                                                            aTAAATGtataagtcatactgttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  <  2:411408‑M1/149‑143 (MQ=255)                                                                                                                                                  |                                                                                                                                               AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_pae_tpiA/1978357‑1978644

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