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

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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCA  >  NC_000913_3_hsa_tpiA/1978357‑1978647                                                                                                                                                   |                                                                                                                                                  aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtat                                                                                                                                                >  1:492412‑M1/1‑146 (MQ=255) aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtat                                                                                                                                                >  1:45363‑M1/1‑146 (MQ=255)   aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataa                                                                                                                                              >  1:960169‑M1/1‑144 (MQ=255)     gggTCTTTACACTTAAATAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATATACAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagt                                                                                                                                            >  2:663521‑M1/1‑142 (MQ=255)       gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtca                                                                                                                                          >  2:499396‑M1/1‑140 (MQ=255)              cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTCCACGCACATACAACGGAGGGGGGCTGCGCTTTTCAATAATGCCTGATGCAGATCACACAAAACTCTCCAT‑TATTTAATTTATATGtattagtcatactttt                                                                                                                                   <  2:893735‑M1/149‑17 (MQ=255)                     cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                            <  1:675694‑M1/149‑24 (MQ=255)                      tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                           >  2:661988‑M1/1‑125 (MQ=255)                         gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgt                                                                                                                        <  2:827395‑M1/149‑28 (MQ=255)                           tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                      <  1:99987‑M1/149‑30 (MQ=255)                                  tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                               >  2:400803‑M1/1‑113 (MQ=255)                                  tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                               >  2:771885‑M1/1‑113 (MQ=255)                                     aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                            >  1:19410‑M1/1‑110 (MQ=255)                                     aTCGACGCAACTGGACTCGTCAATACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                            <  2:106751‑M1/149‑40 (MQ=255)                                       cGACGCACCTGTACTCGTCACTACACGCGCACACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                                                                          <  1:669623‑M1/149‑42 (MQ=255)                                       cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                                                                          >  2:150893‑M1/1‑108 (MQ=255)                                           gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                      <  1:731499‑M1/149‑46 (MQ=255)                                           gCAAATGTAATCGTGACTAAACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCATATCAAACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                      <  2:61477‑M1/149‑46 (MQ=255)                                             aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                    >  1:538093‑M1/1‑102 (MQ=255)                                             aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttggctgtatttcaatatgttaaaaagt                                                                                                    >  1:167536‑M1/1‑102 (MQ=255)                                               cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                  >  2:88904‑M1/1‑100 (MQ=255)                                               cTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                  <  1:206712‑M1/149‑50 (MQ=255)                                                    cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttcaaaagtttttcgc                                                                                             >  1:954875‑M1/1‑95 (MQ=255)                                                     tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatactttttttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                            >  2:224441‑M1/1‑94 (MQ=255)                                                              cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtattttaatctgttaaaaagtttttcgctacgctagca                                                                                   >  1:585764‑M1/1‑85 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                >  1:739748‑M1/1‑82 (MQ=255)                                                                    caTACAACGGAGGGGGGCTGCGCTTTTACATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcataattttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                                                                             <  1:368892‑M1/149‑71 (MQ=255)                                                                         aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATACAGATCACACAAAA‑ATTTATTGTTCGTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                        <  1:1039870‑M1/149‑76 (MQ=255)                                                                           cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                      <  1:88904‑M1/149‑78 (MQ=255)                                                                           cGGAGGGGGGCTACGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                      <  1:658397‑M1/149‑78 (MQ=255)                                                                                gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagt                                                                                                                                            <  2:475845‑M1/74‑8 (MQ=255)                                                                                gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagt                                                                                                                                            >  1:475845‑M1/1‑67 (MQ=255)                                                                                gggggCTGCGATTTTCAATAATGCGTGAGGCATATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatactttttttttgggtgtatttcaatctgttaaaaagtttttcgctacactagcaagcaaaaatgaaacagga                                                                 >  2:10044‑M1/1‑67 (MQ=255)                                                                                        cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                         >  1:303144‑M1/1‑59 (MQ=255)                                                                                          aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                       >  2:95515‑M1/1‑57 (MQ=255)                                                                                             ttCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TCCTTAACATAAATGtataagtcatacttttggtttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaaaggaataatggaaatgg                                                    >  1:788027‑M1/1‑54 (MQ=255)                                                                                             ttCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                  <  2:649034‑M1/103‑50 (MQ=255)                                                                                             ttCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                  >  1:649034‑M1/1‑54 (MQ=255)                                                                                                      gagtgATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcctacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcg                                           <  2:344848‑M1/147‑105 (MQ=255)                                                                                                       cGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                          >  2:541911‑M1/1‑44 (MQ=255)                                                                                                               aaatcACACAAAACACTCAAT‑TACTTAACATAAATGtagaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                                  <  2:244987‑M1/147‑114 (MQ=255)                                                                                                               aGATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaa                                  >  1:499881‑M1/1‑36 (MQ=255)                                                                                                                gATCACACAAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                                 >  1:727200‑M1/1‑35 (MQ=255)                                                                                                                    acacaAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                             <  1:150893‑M1/149‑119 (MQ=255)                                                                                                                    acacaAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcgacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                             <  2:516224‑M1/149‑119 (MQ=255)                                                                                                                    acacaAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgtttaaggtgtatttcaatctgttaaaaagtttttcgctacgctaacaagccaaaatgtaaccggaataatcgacatgggatgttgcgaacagtcaaaaaaa                             >  1:938772‑M1/1‑31 (MQ=255)                                                                                                                      acaAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                           <  1:202634‑M1/149‑121 (MQ=255)                                                                                                                      acaAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                           >  1:80708‑M1/1‑29 (MQ=255)                                                                                                                      acaAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                             >  2:310957‑M1/1‑29 (MQ=255)                                                                                                                      acaAAACACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                             <  1:310957‑M1/131‑103 (MQ=255)                                                                                                                        aaaaCACTCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                         <  1:353917‑M1/149‑123 (MQ=255)                                                                                                                               tCAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                  <  1:771885‑M1/149‑130 (MQ=255)                                                                                                                                gtgatggcttgacatAAATGtataaggcatacttttggtttgggtgtatttcaaactgttaaaaagtttttcgctacgctagcaagcaaaaatcaaacaagaacaatcgaactggtatgttgcgcatagtcaaaacaacccaccgtaaa                  <  2:1048173‑M1/142‑130 (MQ=255)                                                                                                                                cAAT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat                 <  2:45363‑M1/149‑131 (MQ=255)                                                                                                                                  aT‑TACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa               <  1:661988‑M1/149‑133 (MQ=255)                                                                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                            >  2:648911‑M1/1‑13 (MQ=255)                                                                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                            <  1:648911‑M1/100‑88 (MQ=255)                                                                                                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                    <  2:12614‑M1/107‑96 (MQ=255)                                                                                                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                                    >  1:12614‑M1/1‑12 (MQ=255)                                                                                                                                        ttAACATAAATGtataagtcatactgttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcgcagtcaaaataactcaccataaataatcatc          >  1:427459‑M1/1‑12 (MQ=255)                                                                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc       >  2:1038074‑M1/1‑9 (MQ=255)                                                                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                            >  2:371612‑M1/1‑9 (MQ=255)                                                                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                            <  1:371612‑M1/128‑120 (MQ=255)                                                                                                                                                aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataa  >  2:517592‑M1/1‑4 (MQ=255)                                                                                                                                                   |                                                                                                                                                  AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAAT‑TACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCA  >  NC_000913_3_hsa_tpiA/1978357‑1978647

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑