Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_vch_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_vch_tpiA	1978503–1979270	1979278	9–776	37 [0]	[0] 36	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_tpiA	= 1978502	0 (0.000)	35 (1.170)	30/272	0.1	100%	noncoding (1978502/4641682 nt)	repeat_region	repeat region
	?	NC_000913_3_vch_tpiA	1979279 =	0 (0.000)					noncoding (1979279/4641682 nt)	repeat_region	repeat region

TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_vch_tpiA/1978354‑1978644                                                                                                                                                     |                                                                                                                                               tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTAACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATG                                                                                                                                                >  1:447017‑M1/1‑149 (MQ=255)                  acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                                                               >  2:78685‑M1/1‑132 (MQ=255)                      aTCTAAGATTTTTCCTAAATCGACGCAACCGTACTCGTTACTACACGCACATACATCGGAGGGGGGCTGCTATTTTCAATAATGCTTTATGTATATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                           <  1:47037‑M1/149‑22 (MQ=255)                        cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCATACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  1:248299‑M1/149‑24 (MQ=255)                        cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  2:128167‑M1/149‑24 (MQ=255)                         tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACATAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                        <  1:137613‑M1/149‑25 (MQ=255)                         tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                        >  2:215209‑M1/1‑125 (MQ=255)                                              gCAACTTTACTCGGCACTACACGCACAAACAACGGAGAGGAGCTGCGATTTTCAATAATGCGTAATGCAGATCACACAAAACACTCAATTACATAACACAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:229893‑M1/149‑46 (MQ=255)                                                aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatgtcagtctgttaaaaagt                                                                                                 >  1:537344‑M1/1‑102 (MQ=255)                                                   tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                              >  2:195959‑M1/1‑99 (MQ=255)                                                     tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTGACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                            <  1:553401‑M1/149‑53 (MQ=255)                                                        tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                         >  1:272424‑M1/1‑94 (MQ=255)                                                             aCTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                                                                    >  2:202659‑M1/1‑89 (MQ=255)                                                               tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataattcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                  >  1:458812‑M1/1‑87 (MQ=255)                                                               tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAAATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttatgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                  >  1:8201‑M1/1‑87 (MQ=255)                                                                  acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                               <  2:447017‑M1/149‑66 (MQ=255)                                                                        aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatcttttaaaaagtttttcgctacgctagcaagcaaaa                                                                         >  2:350751‑M1/1‑78 (MQ=255)                                                                             aCGGAGGGGGGCTGCGATTTTCAATAATGCGTAATGCATATCACAAAAAAAACTCCATTAATTAACATAAATGtataagtcatacttttgtttttggtgtatttgaatctgttaaaaagttattagctactctagcatacaataatgaa                                                                    >  1:191969‑M1/1‑73 (MQ=255)                                                                              cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                   <  1:46507‑M1/149‑78 (MQ=255)                                                                              cGGAGGGGGGCTGCGATTGTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                   >  2:341582‑M1/1‑72 (MQ=255)                                                                                 aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgttaaaaagattttcgctacgctagcaagcaaaaatgaaacag                                                                >  2:57799‑M1/1‑69 (MQ=255)                                                                                         tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtaatccatgtaggtaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                        >  1:404315‑M1/1‑61 (MQ=255)                                                                                          gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcg                                                       >  2:188073‑M1/1‑60 (MQ=255)                                                                                             aTTTTCAATAATGCGTGATGGAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttatagagtttttcgctacgctagcaagcaaaaatgaaacaggaagaagcgaaa                                                    <  2:327115‑M1/149‑93 (MQ=255)                                                                                                   aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggactaatcgaaatgggat                                              >  2:45588‑M1/1‑51 (MQ=255)                                                                                                   aataatGCGCGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttattactacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                              >  1:42389‑M1/1‑51 (MQ=255)                                                                                                    ataatTCGATATGGAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                             <  2:102844‑M1/149‑100 (MQ=255)                                                                                                    ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgatatgggtgtatttcaatctgttaaaaagtttttcactaccctggcgagaaaaaatgaaagaagaataatcgaaatgcgctg                                             >  2:377660‑M1/1‑50 (MQ=255)                                                                                                          cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       <  1:327328‑M1/149‑106 (MQ=255)                                                                                                                  aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggatgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacactcat                               >  1:122589‑M1/1‑36 (MQ=255)                                                                                                                   tataacagaaaacACTTAGTTATGTAACATAAATGtataagtcataattgtgttttgggtgtatttcaatcttttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                              <  2:25203‑M1/145‑115 (MQ=255)                                                                                                                     tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtaaagat                            >  2:299952‑M1/1‑33 (MQ=255)                                                                                                                       acacaAAACACTCAATGACTTAACATAAATGtataagtcatacttttgttttgggtgtagttctatctgttaaaaaggttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaacatac                          <  2:510842‑M1/149‑119 (MQ=255)                                                                                                                           aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                      >  2:159489‑M1/1‑27 (MQ=255)                                                                                                                             aaCCCTCAATTACTAAACATAAATTtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                    <  2:42389‑M1/149‑125 (MQ=255)                                                                                                                             aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                    <  1:202659‑M1/149‑125 (MQ=255)                                                                                                                                               aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  <  1:159489‑M1/149‑143 (MQ=255)                                                                                                                                                     |                                                                                                                                               TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_vch_tpiA/1978354‑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: ‑