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	35 [0]	[0] 35	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)	32 (0.750)	26/274	0.8	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

GCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_vch_tpiA/1979142‑1979424                                                                                                                                          |                                                                                                                                                  tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacccaaaacactcaattacttaacataaattTATAAGTCATAc                                                                                                                                        >  1:62895‑M2/138‑149 (MQ=255)   tacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtt                                                                                                                                      >  2:79254‑M2/136‑149 (MQ=255)     cacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                    <  1:194549‑M2/16‑1 (MQ=255)            ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                             <  2:277408‑M2/23‑1 (MQ=255)             taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                            >  2:441096‑M2/126‑149 (MQ=255)             taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                                                                            >  1:716287‑M2/126‑149 (MQ=255)             taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacaaacaacggagggagacagcgattttcaataatgcacgatacagatcacacaaaacactaaattactcaacataaatgTAAAAATCATACTCTTGCTTCggg                                                                                                                            >  2:51130‑M2/126‑149 (MQ=255)                   ttttcctaaatcgacgcaactgtactcgtcactaaaaacaaatacaacggaggggggctgcgattttcaatactgcgtgatgcagatcacacaaaacactcaattacttaacataactgTATAAGTCATACTTTTCTTTTGCGTGTAtt                                                                                                                      >  1:621341‑M2/120‑149 (MQ=255)                    tttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                     <  1:79254‑M2/31‑1 (MQ=255)                             tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgagattttcaataatgcgtgatgcggatcacacaaaaccctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCGGtt                                                                                                            >  2:44417‑M2/110‑149 (MQ=255)                                     actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                    <  2:252932‑M2/48‑1 (MQ=255)                                       tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                  <  1:77947‑M2/50‑1 (MQ=255)                                            tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaattTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTggct                                                                                             >  1:508766‑M2/95‑149 (MQ=255)                                            tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                             >  2:697363‑M2/95‑149 (MQ=255)                                              gtcactacacgcacatacaacggaagggggctgcgatttacaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                           <  1:282725‑M2/57‑1 (MQ=255)                                                   tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacataacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGATAAAAAGTTTTTCGCTACGCTag                                                                                      >  1:661465‑M2/88‑149 (MQ=255)                                                            atacaacggaggggggctgcgattttcaataatgcataatgcaggtcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                             >  1:464003‑M2/79‑149 (MQ=255)                                                                  cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTGTAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                       <  1:749312‑M2/77‑1 (MQ=255)                                                                             tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                            >  2:425387‑M2/62‑149 (MQ=255)                                                                              gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                           >  1:204288‑M2/61‑149 (MQ=255)                                                                                         taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGCATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                                <  1:697363‑M2/100‑1 (MQ=255)                                                                                          aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                               >  2:216424‑M2/49‑149 (MQ=255)                                                                                                           acagaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGGTTTTCGCCACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                              <  2:148027‑M2/118‑1 (MQ=255)                                                                                                           acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                              <  2:204288‑M2/118‑1 (MQ=255)                                                                                                               aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGCATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                          >  1:516783‑M2/28‑149 (MQ=255)                                                                                                                    actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTa                     >  1:588694‑M2/23‑149 (MQ=255)                                                                                                                        aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                       <  2:61607‑M2/77‑1 (MQ=255)                                                                                                                        aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                       >  1:61607‑M2/19‑95 (MQ=255)                                                                                                                          ttacttaacataaacgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTCaataat               <  2:247327‑M2/133‑1 (MQ=255)                                                                                                                            acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatca             <  1:463234‑M2/135‑1 (MQ=255)                                                                                                                             cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTTAGTATGTTTAACAGTTTTTTGCTACGCTAGCGAGCAAAAATGAATCAGGAATAATCGAAATGGGGTGTTGCGCACAGTCAAAATGACTCAGAGTAAATAatcat            >  2:330088‑M2/14‑149 (MQ=255)                                                                                                                             cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACATGAATAATCGAAATGGGATGTTGCGCACAGTCACAATAACTAACCGTAAATAatcat            >  2:329857‑M2/14‑149 (MQ=255)                                                                                                                                   ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta      >  1:355483‑M2/8‑149 (MQ=255)                                                                                                                                   ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTATTCGCTACGCTAGCAAGCAAAAATTAAACCGGAATAATCGAAATGGCATTTTGCGCACAGTCAAAATAACTCACCGTaaataatcatccgata      >  2:150872‑M2/8‑144 (MQ=255)                                                                                                                                       atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa  >  1:409785‑M2/4‑149 (MQ=255)                                                                                                                                          |                                                                                                                                                  GCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_vch_tpiA/1979142‑1979424

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