Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913_3_vch_tpiA	257,908	Δ776 bp		[crl]	[crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_vch_tpiA	257908	258683	776	27 [0]	[0] 27	[crl]	[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_vch_tpiA	= 257907	0 (0.000)	26 (0.970)	20/268	0.4	100%	noncoding (257907/4641682 nt)	repeat_region	repeat region
	?	NC_000913_3_vch_tpiA	258684 =	0 (0.000)					noncoding (258684/4641682 nt)	repeat_region	repeat region

TTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_vch_tpiA/257763‑258048                                                                                                                                                 |                                                                                                                                              ttGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGGTCAAAAAATTTACCGCACTAGGCCCGAATAATCGTGAAGGTAagtg                                                                                                                                           >  1:194980‑M1/1‑145 (MQ=255)           catTTGCTAAATTTTGCCAAGTTGGTAAAAAAGTTGCATGACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                                                                                                                 <  2:321386‑M1/148‑15 (MQ=255)            aTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatc                                                                                                                                <  2:200744‑M1/149‑16 (MQ=255)                  tAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACACGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAactacaaggataatcaagtct                                                                                                                          >  1:371289‑M1/1‑128 (MQ=255)                           ccAATTTGGTAAAACAGTTGCATCACAACAGGAGAGAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                                                                                                                 <  1:190824‑M1/149‑31 (MQ=255)                                    tAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg                                                                                                        >  1:388282‑M1/1‑110 (MQ=255)                                      aaaCAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgta                                                                                                      >  2:34424‑M1/1‑108 (MQ=255)                                              gCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaa                                                                                              <  1:34424‑M1/149‑50 (MQ=255)                                                        aGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                                                                                                              <  2:94815‑M1/123‑34 (MQ=255)                                                        aGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                                                                                                              >  1:94815‑M1/1‑90 (MQ=255)                                                        aGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTAACCTCACTAGTCCCGTATATTCGTGAAGGAAaaaacaaagataatcgattctttttcgattgtctgactgtatgcgtcaacgtgaaacgg                                                                                    >  1:97057‑M1/1‑90 (MQ=255)                                                           agaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggca                                                                                 <  1:468109‑M1/149‑63 (MQ=255)                                                                 cAAGGACGTTACCGAGTGGACACCCGAAGAGCAGATTGAGCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaa                                                                           <  2:97057‑M1/149‑69 (MQ=255)                                                                   aTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgcttgtctggctgtatgcgccaccgtgacacgggcacgcgaata                                                                         >  2:152671‑M1/1‑79 (MQ=255)                                                                   aTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagt                                                                         >  2:220480‑M1/1‑79 (MQ=255)                                                                   aTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCAATAGGCGCGTATATTCGTGAAGGTAagtgcaaggataatcgattctttttcgattggctggctatgtgcgtcaacatgaaaccggcaccggaagt                                                                         >  2:153071‑M1/1‑79 (MQ=255)                                                                              gAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattct                                                              >  1:315555‑M1/1‑68 (MQ=255)                                                                                               gCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcgaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggag                                             >  1:144100‑M1/1‑51 (MQ=255)                                                                                                   aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagcggg                                         >  1:424312‑M1/1‑47 (MQ=255)                                                                                                                tttACCGCACTAGGCACGTATATTCGTGGAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatc                            <  2:243970‑M1/149‑116 (MQ=255)                                                                                                                             ggccgtATATTTGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggttgatggagctggaagcgcaggaatcccgttttaccttc               <  2:282983‑M1/147‑129 (MQ=255)                                                                                                                              cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaca              <  2:315555‑M1/149‑130 (MQ=255)                                                                                                                              cccGTATAGTCGTGAAGCTAtgtgggaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctaca              <  1:342145‑M1/149‑130 (MQ=255)                                                                                                                                   ttattCGTGAAGGTAagtgcaaagataagcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttac         <  2:144100‑M1/148‑135 (MQ=255)                                                                                                                                   atatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttac         <  1:220480‑M1/149‑135 (MQ=255)                                                                                                                                         gTGAAGGTAagtgcaaagataatcgattcttttttttttttttttctgtatgcgacaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggctggagctggaagcgcaggaatcccgttttgcctacagttaccagttt   >  2:389788‑M1/1‑9 (MQ=255)                                                                                                                                          tcaaggTAagtgcaaagataatcgatgtgttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagtttg  <  1:63213‑M1/147‑142 (MQ=255)                                                                                                                                                 |                                                                                                                                              TTGTGTTAAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_vch_tpiA/257763‑258048

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