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	1979278	776	48 [0]	[0] 48	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)	47 (1.040)	40/270	0.2	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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_vch_tpiA/1978357‑1978643                                                                                                                                                   |                                                                                                                                              aaaTGGGTCTTTACAATTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtat                                                                                                                                            >  2:692744‑M1/1‑146 (MQ=255)        tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcat                                                                                                                                     >  2:678546‑M1/1‑139 (MQ=255)        tCTTTACACTTATCTAAGATTTATTCTAGATCGCACCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGAAGATCACACAAAACAATCAATTACTTAAACTA‑AATGtatacgtact                                                                                                                                     >  2:356723‑M1/1‑139 (MQ=255)              cacTTATCTAAGAATTTTCCTAAATCGACGCAACTGTACTCTTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTCACATA‑AATGtataagtcatactttt                                                                                                                               <  2:35873‑M1/149‑17 (MQ=255)                 ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgtt                                                                                                                            >  2:656146‑M1/1‑130 (MQ=255)                    tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttg                                                                                                                         >  2:474693‑M1/1‑127 (MQ=255)                     cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGATGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgg                                                                                                                        <  2:734527‑M1/149‑24 (MQ=255)                     cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgg                                                                                                                        <  2:438071‑M1/149‑24 (MQ=255)                             tttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtattt                                                                                                                <  2:365075‑M1/149‑32 (MQ=255)                             tttCCTAAATCGACGCAACTGTACGCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtattt                                                                                                                <  2:365074‑M1/149‑32 (MQ=255)                                     aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                        <  1:678546‑M1/149‑40 (MQ=255)                                             aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagt                                                                                                >  1:659988‑M1/1‑102 (MQ=255)                                                     tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggcgtatttcaatctgttaaaaaatttttcgct                                                                                        >  2:45704‑M1/1‑94 (MQ=255)                                                       gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                      >  2:116304‑M1/1‑92 (MQ=255)                                                               acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                              >  2:738111‑M1/1‑84 (MQ=255)                                                                   acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                          >  1:739065‑M1/1‑80 (MQ=255)                                                                       aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                      >  1:179614‑M1/1‑76 (MQ=255)                                                                          aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                   >  1:715983‑M1/1‑73 (MQ=255)                                                                           cgtgagggggctgcgATTTTCAATAATGCGTTATGAAGATCACATAAAACACTTAATTACTTAACATA‑AATGtataagtcataattttgttttgggtgtattgccatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                  <  2:391245‑M1/144‑78 (MQ=255)                                                                               ggggggCTGCGATTTTCAATAATGCGTGATGCCGATCACACAAAACACTCCATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                                                              <  1:490611‑M1/149‑82 (MQ=255)                                                                                    gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                         >  1:800981‑M1/1‑63 (MQ=255)                                                                                     cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                        <  1:238912‑M1/149‑88 (MQ=255)                                                                                     cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                        >  1:362718‑M1/1‑62 (MQ=255)                                                                                     cTGCGATTTTCAATAATGCGTACTGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgataaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                        >  2:395583‑M1/1‑62 (MQ=255)                                                                                            tttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                 >  1:45669‑M1/1‑55 (MQ=255)                                                                                                   aaagcgtGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AAGGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                          <  1:70767‑M1/146‑102 (MQ=255)                                                                                                     tGCTTGATGCAGGTCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                        <  2:316926‑M1/149‑104 (MQ=255)                                                                                                         tGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgctttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgagatgggatgtgtcgccc                                    >  2:269185‑M1/1‑42 (MQ=255)                                                                                                             gCAGATCACACAAAACACTCAATTACTTAACATACAATGtttggccatacttttgttttggatgtacttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgagatgggctgtttggaacagtc                                 >  2:181061‑M1/1‑39 (MQ=255)                                                                                                                gATCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaccaggaataagcgaaatgggatgttgcgcacagtcaaa                             >  1:442227‑M1/1‑35 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                            <  1:303241‑M1/149‑116 (MQ=255)                                                                                                                 aTCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcaagcaagcaaaaatgaaacagaaataatcgaaatgggatgttgcgcacagtcaaaa                            >  1:750536‑M1/1‑34 (MQ=255)                                                                                                                  tCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                       <  2:260370‑M1/89‑57 (MQ=255)                                                                                                                  tCACACAAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                       >  1:260370‑M1/1‑33 (MQ=255)                                                                                                                    acacaAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                         <  1:570317‑M1/149‑119 (MQ=255)                                                                                                                     cacaAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcgcagtcaaaataac                        >  1:250621‑M1/1‑30 (MQ=255)                                                                                                                      acaAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                       >  1:557892‑M1/1‑29 (MQ=255)                                                                                                                      acaAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                       >  1:545478‑M1/1‑29 (MQ=255)                                                                                                                      acaAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                       >  2:355140‑M1/1‑29 (MQ=255)                                                                                                                      acaAAACACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtattacaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                       >  1:705698‑M1/1‑29 (MQ=255)                                                                                                                        aaaaCACTCAATTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                     <  2:63964‑M1/149‑123 (MQ=255)                                                                                                                                 aaTTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccttaaata            >  2:473364‑M1/1‑18 (MQ=255)                                                                                                                                 aaTTACTTAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata            <  1:738111‑M1/149‑132 (MQ=255)                                                                                                                                  atgacttAACATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa           <  2:634143‑M1/146‑133 (MQ=255)                                                                                                                                    taattaacATA‑AATGaataagtcatacttttgtttagagtgtatttcaatctgttaaaaagtttttcgctacgataaaaaaaaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc         <  2:779937‑M1/146‑135 (MQ=255)                                                                                                                                         aaCATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg    <  1:470210‑M1/149‑140 (MQ=255)                                                                                                                                         aaCATA‑AATGtataagtcatacttttgttctgggtgtatttcaatctgtcaaaaagtttttcgctacgctagcaagcaaaaatgagacaggaataatcgacatgggttgttgcgcacagtcaaaataatccaacgtaaatactcaactt    >  2:763927‑M1/1‑10 (MQ=255)                                                                                                                                           cATA‑AATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct  >  1:416632‑M1/1‑8 (MQ=255)                                                                                                                                                   |                                                                                                                                              AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATA‑AATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_vch_tpiA/1978357‑1978643

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