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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_pae_tpiA	1978503	1979278	776	53 [0]	[1] 55	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)	47 (1.280)	37/260	NT	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                                                                                                                                                  |                                                                                                                                               aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGCGGGGGGCTGCGATTTTCAATAATGCGTGATGCCGATCCCGCAAAACACTCATTTACTTACCATAAATGtat                                                                                                                                             >  2:384400‑M1/1‑146 (MQ=255)   aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTGAACATAAATGtataa                                                                                                                                           <  1:535799‑M1/149‑6 (MQ=255)      gctctTGACTCTTATCTAAGGTTTTTCCTAAATCGACGCATCTGTACTCGTCACTATACGCACATACAACGGAGGGGGGCTGCTATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                                                                        <  2:451621‑M1/147‑9 (MQ=255)        tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAAAACTCAATTACTTCACATAAATGtataagtcat                                                                                                                                      >  1:416451‑M1/1‑139 (MQ=255)          tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                    <  1:658967‑M1/149‑13 (MQ=255)                   ttcTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                           <  1:635422‑M1/148‑22 (MQ=255)                   ctgtaagATTGTTCCTAAATTGACTCAAGTGTATTTGTCACTACCCGTAGATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                           <  2:103838‑M1/146‑22 (MQ=255)                    tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGAACTACAATAATGCGTTATGCAGATCACACAAAACACTCAATTACATAACAAAAATGtataagtcatactattgttttg                                                                                                                          <  2:671226‑M1/149‑23 (MQ=255)                    tCTAAGATTGTGCCTAAATCGACGCAACTGTAGCCGTCACTACACGCACATACAACGGAGGGGGGCTGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                          <  2:221670‑M1/149‑23 (MQ=255)                       aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTAAATAATGCGTGATGCAGATCACACCAAGAACTCAAATACTTACCATAAATGtatcaaatctacttttgttgttggt                                                                                                                       >  1:249123‑M1/1‑124 (MQ=255)                        aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgttttcggtg                                                                                                                      >  2:317290‑M1/1‑123 (MQ=255)                            ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatt                                                                                                                  >  1:162534‑M1/1‑119 (MQ=255)                                 cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                             >  1:561462‑M1/1‑114 (MQ=255)                                 cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGACTGCGAATTTCAATAAAGCGTGAAGCAGAACACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttcgttgtgggtgtatttcaat                                                                                                             >  2:573578‑M1/1‑114 (MQ=255)                                  tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTAATTAACATAAATGtataagtcatacttttgttttgggtgtatttcactc                                                                                                            >  2:127793‑M1/1‑113 (MQ=255)                                      tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                        >  1:122911‑M1/1‑109 (MQ=255)                                         aCGCAACTGTACTCGTCACTACACGCATATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                     <  1:616076‑M1/149‑44 (MQ=255)                                         aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACCTAAATGtataagtcatacttctgttttggatttcttccaatctcttcaa                                                                                                     >  2:583399‑M1/1‑106 (MQ=255)                                           gtaactGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:243590‑M1/147‑46 (MQ=255)                                             aaCTGTACTCGTCACTACACGCACATACAACGGCGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATCACTTAACATAAATGtataagaaagacttttgatttgagtttatatcaatctggtaaaaagt                                                                                                 >  2:507954‑M1/1‑102 (MQ=255)                                                cgaactcGGCGCTCCACGCCCATACATCGGAGGTGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAATCCCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtggatttcaatctgttaaaaagtttt                                                                                              <  2:441604‑M1/146‑51 (MQ=255)                                                 gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                             <  1:447332‑M1/149‑52 (MQ=255)                                                     tCGTCACTACACGCACCTACAACGGAGAGGGATTGGAGAACGCATTAATCCGTGATGCAGATGACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcggt                                                                                         >  1:729171‑M1/1‑94 (MQ=255)                                                               acGCACATACAACGGAGGGGGGCTGCTATTTTCAATAATGCGTCATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                               <  1:253293‑M1/149‑66 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                         <  2:100312‑M1/137‑56 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgttaaaaagtttttcgct                                                                                         >  1:100312‑M1/1‑82 (MQ=255)                                                                      tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                        >  2:230311‑M1/1‑77 (MQ=255)                                                                               ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                                                               <  1:83814‑M1/149‑82 (MQ=255)                                                                                    gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttacctccgctagcgagcaaaaattaaaaacaaataa                                                          >  2:591819‑M1/1‑63 (MQ=255)                                                                                      tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                        >  1:176170‑M1/1‑61 (MQ=255)                                                                                        cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatccttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagaaagcaaaaatgaaacagaaagaataga                                                      >  1:180500‑M1/1‑59 (MQ=255)                                                                                        cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                      >  1:270916‑M1/1‑59 (MQ=255)                                                                                         gATTTTCAATAATGCGTGATGCAGATCACACCAAACCCTCAACTATTTTGCTTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaa                                                     <  1:318488‑M1/149‑92 (MQ=255)                                                                                             ttCAATAAAGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                 <  2:352572‑M1/149‑96 (MQ=255)                                                                                                aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaccaaaaatgaaacaggaataatcgaaatgggat                                              >  2:511782‑M1/1‑51 (MQ=255)                                                                                                         tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     >  1:599950‑M1/1‑42 (MQ=255)                                                                                                         tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     >  1:599944‑M1/1‑42 (MQ=255)                                                                                                                gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                              >  1:369519‑M1/1‑35 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTAGTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                             <  2:455629‑M1/149‑116 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                             >  1:468688‑M1/1‑34 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaaagggatgttgcgcacagtcaaaa                             >  1:378539‑M1/1‑34 (MQ=255)                                                                                                                  tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                            >  1:373704‑M1/1‑33 (MQ=255)                                                                                                                    acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                          >  2:441959‑M1/1‑31 (MQ=255)                                                                                                                               tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggcataatcgaaatgggatgttgcgcacagtcaaaataaataaccgtaaa               >  2:597751‑M1/1‑20 (MQ=255)                                                                                                                               tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgggcacagtcaaaataactcaccgtaaa               >  1:104075‑M1/1‑20 (MQ=255)                                                                                                                               tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactaaccgtaaa               >  2:597810‑M1/1‑20 (MQ=255)                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgggcacagtcaaaataactcaccgtaaata             >  1:133784‑M1/1‑18 (MQ=255)                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata             >  1:441374‑M1/1‑18 (MQ=255)                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaaagaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata             >  1:349333‑M1/1‑18 (MQ=255)                                                                                                                                        tAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct      <  1:47377‑M1/149‑139 (MQ=255)                                                                                                                                         aaaataaTTGtataaggcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgagaacacgtaaaataactcaccgtaaataatcatctg     <  2:157092‑M1/146‑140 (MQ=255)                                                                                                                                         aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg     >  1:80478‑M1/1‑10 (MQ=255)                                                                                                                                            aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtccaaataactcaccgtaaataatcatctgcta  >  2:261212‑M1/1‑7 (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 < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑