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	61 [0]	[0] 62	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)	59 (0.780)	48/262	0.7	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

GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_hsa_tpiA/1978356‑1978641                                                                                                                                                   |                                                                                                                                            gAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                                                           >  2:155326‑M1/1‑147 (MQ=255)  aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtct                                                                                                                                          >  1:681291‑M1/1‑146 (MQ=255)      gggTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                                                                      >  1:220267‑M1/1‑142 (MQ=255)      ggctcttTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                                                                      <  1:54561‑M1/146‑8 (MQ=255)        gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACACCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                    >  1:179082‑M1/1‑140 (MQ=255)                  ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATCAATTACTTAACCTACATGtataagtcatacttttgtt                                                                                                                          >  2:407281‑M1/1‑130 (MQ=255)                  ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTCCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                          >  2:839270‑M1/1‑130 (MQ=255)                     tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACACAACACTCAATTACTTAACATAAATGtataattcatacttttgttttg                                                                                                                       >  1:168653‑M1/1‑127 (MQ=255)                     tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                       <  1:1080577‑M1/149‑23 (MQ=255)                      cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                      <  1:921748‑M1/149‑24 (MQ=255)                      cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                      <  1:537770‑M1/149‑24 (MQ=255)                       tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                     >  2:215237‑M1/1‑125 (MQ=255)                        aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                    >  2:957539‑M1/1‑124 (MQ=255)                                  ctaagctgacggagcTGTACTCGTCACTACAGGCGCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTGCTTTACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                          <  2:566291‑M1/142‑36 (MQ=255)                                   tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatgcttttgttttgggtgtatttcaatc                                                                                                         >  2:1111441‑M1/1‑113 (MQ=255)                                   tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttctgggtgtatttcaatc                                                                                                         >  2:857248‑M1/1‑113 (MQ=255)                                    aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCGATAATGCGTGATGCAGATCACACAAAACACTCAATTAATTAACATAAATGtataagtcatacttttgttttgggtgtttttcaatct                                                                                                        >  1:465618‑M1/1‑112 (MQ=255)                                        cGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                                                                    <  1:1111441‑M1/149‑42 (MQ=255)                                            gCAACTGTACTCGTCACTAGACGCACATACAACGGAGGGGAGCTGCGGTTTTCAATAGTGCGTGATGCAGATCACACAAAACACTCTATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                <  2:81939‑M1/149‑46 (MQ=255)                                            gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                <  1:233593‑M1/149‑46 (MQ=255)                                             cAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                                                               >  2:684044‑M1/1‑103 (MQ=255)                                               aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                             <  2:395393‑M1/149‑49 (MQ=255)                                                  gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                          <  2:218261‑M1/149‑52 (MQ=255)                                                     cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                                                                       <  2:769357‑M1/149‑55 (MQ=255)                                                       cGTCACTACACACACCGACAATGGATGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                     <  2:220267‑M1/149‑57 (MQ=255)                                                             tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                               <  1:957539‑M1/149‑63 (MQ=255)                                                                  gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctcgcaagc                                                                          >  1:434576‑M1/1‑82 (MQ=255)                                                                  gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                          >  1:578191‑M1/1‑82 (MQ=255)                                                                          aatggaggGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                  <  1:1313337‑M1/146‑76 (MQ=255)                                                                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagccagcacaaatga                                                                  >  2:529215‑M1/1‑74 (MQ=255)                                                                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcataattttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaaaataa                                                                  >  2:530555‑M1/1‑74 (MQ=255)                                                                           aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                 >  2:449297‑M1/1‑73 (MQ=255)                                                                                  ggggCTGCGATTTTGAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTTACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                          <  2:138449‑M1/149‑84 (MQ=255)                                                                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                      <  1:155326‑M1/149‑88 (MQ=255)                                                                                      cTGCGATTTTCAACAATGCGTTATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                      <  1:801242‑M1/149‑88 (MQ=255)                                                                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                     >  1:1192211‑M1/1‑61 (MQ=255)                                                                                              ttCAATAATGCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                              <  1:740318‑M1/149‑96 (MQ=255)                                                                                                   taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                         >  2:620247‑M1/1‑49 (MQ=255)                                                                                                   taatGCGTGATGCAGATCACACAAAAAACTCAATTACTTAACATAAATGtataagtcatacttttgatttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaataaaaccggaaaaatcaaaatggaatgg                                         >  2:701102‑M1/1‑49 (MQ=255)                                                                                                        cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                    <  2:578191‑M1/149‑106 (MQ=255)                                                                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                 <  1:351534‑M1/149‑109 (MQ=255)                                                                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                             >  1:276262‑M1/1‑39 (MQ=255)                                                                                                             tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                             <  2:276262‑M1/55‑17 (MQ=255)                                                                                                               cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                             >  2:1234865‑M1/1‑37 (MQ=255)                                                                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             >  2:1224090‑M1/1‑35 (MQ=255)                                                                                                                 gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                             <  1:1224090‑M1/115‑81 (MQ=255)                                                                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatatttttcttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                       <  1:474035‑M1/149‑119 (MQ=255)                                                                                                                            acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                <  1:449297‑M1/149‑126 (MQ=255)                                                                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta              >  1:897025‑M1/1‑22 (MQ=255)                                                                                                                                tCAATTACTTAACATAAATGtataagtcatacttttgttttggtgatatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaccaggaataagcgaaataggctgttgcgcacattcaaaataactcaccgtaaa            >  2:164001‑M1/1‑20 (MQ=255)                                                                                                                                tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcactcggttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatagaaatgggatgttgcgcacagtcaacataactcaacgtaaa            >  1:997951‑M1/1‑20 (MQ=255)                                                                                                                                tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgaacagtcaaaataactaaccgtaaa            >  2:1203987‑M1/1‑20 (MQ=255)                                                                                                                                 cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaat           <  1:684044‑M1/149‑131 (MQ=255)                                                                                                                                    ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataat        <  2:1162618‑M1/149‑134 (MQ=255)                                                                                                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                    >  2:1196153‑M1/1‑12 (MQ=255)                                                                                                                                        ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                    <  1:1196153‑M1/69‑58 (MQ=255)                                                                                                                                          aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatcttttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  >  1:925260‑M1/1‑10 (MQ=255)                                                                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                     <  1:919113‑M1/65‑57 (MQ=255)                                                                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                     >  2:919113‑M1/1‑9 (MQ=255)                                                                                                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                      >  2:597041‑M1/1‑7 (MQ=255)                                                                                                                                             aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                                                                      <  1:597041‑M1/78‑72 (MQ=255)                                                                                                                                                   |                                                                                                                                            GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_hsa_tpiA/1978356‑1978641

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