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	20 [0]	[2] 22	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)	19 (0.810)	17/270	0.4	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

TTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_vch_tpiA/1979157‑1979424                                                                                                                           |                                                                                                                                                  gatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaagggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGatgt                                                                                                                         >  2:84419‑M2/123‑149 (MQ=255)     ttcgctaaatcgacgtacgtttacacgttactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                     <  2:330566‑M2/31‑1 (MQ=255)         ctaaatcgacgcaactgtactcgtcactacacgcacatataacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactgaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                 <  1:123010‑M2/35‑1 (MQ=255)         ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGCTGTATTTCAAt                                                                                                                 >  1:72237‑M2/115‑149 (MQ=255)              tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                                                                            >  2:220966‑M2/110‑149 (MQ=255)              tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgagtttcaataatgcgggatgccgatcacccaaaacactcaattacttaacataaatgTGTAATTCATAATTTTGTTTTGGGTTTATTTCAATCTGtt                                                                                                            >  1:162327‑M2/110‑149 (MQ=255)                 acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                                         >  1:220492‑M2/107‑149 (MQ=255)                 acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTaaa                                                                                                         >  1:220542‑M2/107‑149 (MQ=255)                      actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGGCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                    <  1:261202‑M2/48‑1 (MQ=255)                               gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgatac                                                                                           >  2:294211‑M2/93‑149 (MQ=255)                                    tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacatacaattacttaacataactgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                      >  1:397126‑M2/88‑149 (MQ=255)                                       acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                   <  2:266874‑M2/65‑1 (MQ=255)                                              tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                            <  1:203247‑M2/72‑1 (MQ=255)                                                 aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                         >  1:293290‑M2/75‑149 (MQ=255)                                                  acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                        >  1:106068‑M2/74‑149 (MQ=255)                                                              tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                            >  2:39664‑M2/62‑149 (MQ=255)                                                                                     gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                                     <  1:39664‑M2/111‑1 (MQ=255)                                                                                                            tacttaacATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGTCCACCGGAAACTGGTCTCATAATTATGTGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc              <  1:213271/142‑1 (MQ=255)                                                                                                                 aacataaatgTATAAATCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATGATCGAAATGGGATGTTGGGCACAGTCAAAATAACTCACCGTATATAATCATCTg         >  2:332055‑M2/11‑149 (MQ=255)                                                                                                                  acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCCAACATGAAACAGGAATAATCGAAATGGAATGTTGCGCACAGTCACAATAACTCACCGTAACTAATCAACTGc        >  2:280631‑M2/10‑149 (MQ=255)                                                                                                                     tAAATGCATAAGTCAAACTTTTGTTTTGGGTGTATTTCAAACAGTGAAAAAGTTTATCGCTACGCTAGCGAGCAAAAATGAAAAAGGAATAATCGACATGGGATGTTCCGCACAGTCAACAAAACACACAGTAGATAATCATCTGCtat     >  2:6697/1‑149 (MQ=255)                                                                                                                        atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa  <  1:294211‑M2/146‑1 (MQ=255)                                                                                                                           |                                                                                                                                                  TTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_vch_tpiA/1979157‑1979424

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