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	37 [0]	[0] 37	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)	37 (1.210)	30/266	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

GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_vch_tpiA/1978356‑1978641                                                                                                                                                   |                                                                                                                                            gAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                                                           <  1:445815‑M1/149‑3 (MQ=255)    aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTTAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                        <  1:37959‑M1/149‑6 (MQ=255)    aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGTGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                        <  1:499151‑M1/149‑6 (MQ=255)        gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTAGTCGTCGGTGCGTACACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGTAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                    <  1:446600‑M1/149‑10 (MQ=255)           tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatacgtcatac                                                                                                                                 >  2:537538‑M1/1‑137 (MQ=255)           tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                 >  1:467049‑M1/1‑137 (MQ=255)           tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAAATCCCACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                                 >  2:82557‑M1/1‑137 (MQ=255)                  ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                          <  1:361877‑M1/149‑20 (MQ=255)                  ttATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCAATACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                          >  1:379689‑M1/1‑130 (MQ=255)                            tttgtCCTATATAGTAGCAACTGTATTCGTCACTACACGCACATACAACTGAGGGGTGCTGTTATTTTCAATAATGCGTGATGCACATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                <  2:388256‑M1/149‑30 (MQ=255)                                tCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                            >  1:563846‑M1/1‑116 (MQ=255)                                  cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTCCTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaat                                                                                                          >  1:438081‑M1/1‑114 (MQ=255)                                      aTCGAGTCAACTGTACTCGTAACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAATCACTCAATTACTTAACATAAATGtagcagccagacttttgttttgggtgtatttcaatctgt                                                                                                      <  1:265982‑M1/149‑40 (MQ=255)                                      aTCGACGCAACTGTACTCCTCACTACACGCACATACAACGGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAATTTTCAATTACTTAACATAAATGttcttcacaaacttttgttttgggtgtatttcaatctgt                                                                                                      <  1:266064‑M1/149‑40 (MQ=255)                                                   tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGTGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtgttcgggatagctcaaactattaaattgtttttc                                                                                         >  2:135050‑M1/1‑97 (MQ=255)                                                                acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                            >  2:358634‑M1/1‑84 (MQ=255)                                                                  gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                          >  1:341326‑M1/1‑82 (MQ=255)                                                                    acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGAAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                        <  2:260637‑M1/149‑70 (MQ=255)                                                                          aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAGTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatttgttaaaaagtttttcgctatgctagaaagcaaaaaaga                                                                  >  2:496640‑M1/1‑74 (MQ=255)                                                                             ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                                                               >  2:315474‑M1/1‑71 (MQ=255)                                                                                       tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttattttcagttgatttcaatcttttaaaaaatttttccctacgctagcaagcaaaaatgaaacaggaataatc                                                     <  1:101314‑M1/149‑89 (MQ=255)                                                                                           attattaatgatgcgtgaTGCACATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtgtttcgctacgctagcaagtaaaaatggaacaggaataatcgaaa                                                 <  2:60396‑M1/143‑93 (MQ=255)                                                                                            ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                >  2:565826‑M1/1‑56 (MQ=255)                                                                                            ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaaaggaataatcgaaat                                                >  2:26863‑M1/1‑56 (MQ=255)                                                                                                 aataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat                                           <  2:359200‑M1/149‑99 (MQ=255)                                                                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttttctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                 <  1:229231‑M1/149‑109 (MQ=255)                                                                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgttacgctagtaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                 <  1:147004‑M1/149‑109 (MQ=255)                                                                                                                  aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatatgttaaaaaatttttagctacttcaacaagcaaaaatgaaacagtaattattgaaattggaggttgcccacagtcaaaa                          >  2:16821‑M1/1‑34 (MQ=255)                                                                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                       <  1:561935‑M1/149‑119 (MQ=255)                                                                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                       <  1:358634‑M1/149‑119 (MQ=255)                                                                                                                       acaAAACACTCAATTACTTAACATAAATGtataaatcatacttttgttttgggtgtatttcaatctggtaaaaaatttttcgctacgctagcaagcaaaaatgaaacaggcataatcgaaatgggatgttgcgcacagtcaaaataagt                     >  2:320043‑M1/1‑29 (MQ=255)                                                                                                                              acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta              >  1:343856‑M1/1‑22 (MQ=255)                                                                                                                               cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa             <  2:476190‑M1/149‑129 (MQ=255)                                                                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                          >  1:162522‑M1/1‑18 (MQ=255)                                                                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                          <  2:162522‑M1/133‑116 (MQ=255)                                                                                                                                       cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctgcgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat     >  2:347580‑M1/1‑13 (MQ=255)                                                                                                                                          aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  >  2:510926‑M1/1‑10 (MQ=255)                                                                                                                                                   |                                                                                                                                            GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_vch_tpiA/1978356‑1978641

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