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	28 [0]	[0] 28	[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)	28 (1.010)	27/270	0.1	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

AAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913_3_vch_tpiA/257770‑258047                                                                                                                                          |                                                                                                                                             aaaaaTTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATCGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                                                                                                                   >  1:459197‑M1/1‑138 (MQ=255)    agtttgcTAACTTGTGCCAAGTTGGTAAAAGAGTTGCAGCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATGCGTGAAGGTAagtgcaaagataat                                                                                                                                <  2:212406‑M1/147‑15 (MQ=255)       ttGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcga                                                                                                                             >  1:395390‑M1/1‑132 (MQ=255)       ttGATAAATTTAGCCAATTTGGTAAAACTGTTGCATCACAAAAGCAGATAGCAATGACGTTACTGAGTGGACACCCGAGAAGTAGATTGATCAAAATATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcga                                                                                                                             <  2:464415‑M1/149‑18 (MQ=255)          cTAAATTTTGCGAAGTCGGTAAAACAGTTACATCACAACAGGAGATAGCGATGACGTTACCGCGTGGACACCCGAAGAGCAGATTCATCAAAAAATTTACTGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                                                                                                                          >  2:464691‑M1/1‑129 (MQ=255)                     cAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattg                                                                                                               >  2:213120‑M1/1‑118 (MQ=255)                      aaTTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGCATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                                                                                              <  2:342658‑M1/149‑33 (MQ=255)                          tGGTAAAACAGTTGCAACACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctgg                                                                                                          <  1:195097‑M1/149‑37 (MQ=255)                                        cATCACAACAGGAGATAGCAATGACGTTACCGAGTTGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaac                                                                                            <  1:399384‑M1/149‑51 (MQ=255)                                          tCACAACAGGGGATAGCAATGACGTTACCGAGTGGACACTCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                          <  2:459197‑M1/149‑53 (MQ=255)                                          tCACAACAGGAGAAAGAAATGACGTTACCGAGTGGACACCCGAATAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATTTTCGGGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgt                                                                                          <  2:214995‑M1/149‑53 (MQ=255)                                                    agaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTCTATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                                                                                                  >  1:236970‑M1/1‑87 (MQ=255)                                                    agaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagatattcgattctttttcgattgtctggctgtatgc                                                                                                  <  2:236970‑M1/131‑45 (MQ=255)                                                          cAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaa                                                                          <  1:105004‑M1/149‑69 (MQ=255)                                                           aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaag                                                                         >  2:206451‑M1/1‑80 (MQ=255)                                                                 gTTACCGATTGGACACTCGAAGAGTAGATTGATCAAAAAATTTACCGCACTAGGTCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtg                                                                   <  2:451153‑M1/149‑76 (MQ=255)                                                                      cGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattc                                                              >  1:27235‑M1/1‑69 (MQ=255)                                                                         gTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgg                                                           <  1:213120‑M1/149‑84 (MQ=255)                                                                                    aaGAGTAAATTGATCAAAAAATTTACCGCACTAGGCTCGTATAGTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggat                                                <  2:156023‑M1/149‑95 (MQ=255)                                                                                            aTTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctgg                                        >  2:408693‑M1/1‑47 (MQ=255)                                                                                             ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcg                                                                     >  1:362692‑M1/1‑46 (MQ=255)                                                                                             ttGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcg                                                                     <  2:362692‑M1/119‑74 (MQ=255)                                                                                              tGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaa                                      <  1:67259‑M1/149‑105 (MQ=255)                                                                                                 tCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaaagtgaaaccggcatcggcagtgcttgaatgctggggctggtggttggagcgggaagcg                                   >  2:287656‑M1/1‑42 (MQ=255)                                                                                                   aaaaaaTTTACCGCACTAGGCCCGTATAATCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgca                                 <  1:126422‑M1/149‑110 (MQ=255)                                                                                                                      gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattattttgcgattgtctgggtgtatgcgtcaacgttaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctac              >  2:494189‑M1/1‑21 (MQ=255)                                                                                                                       gccgTATATTCCTGAAAGTAagtgaaaagataatcgattctttttctattgactggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgtggctggtggatggagctggaagcgcaggaatcccgttttacctaca             <  2:120769‑M1/148‑130 (MQ=255)                                                                                                                                  gTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatggagctggaagcgcaggaatcccgttttacctacagttaccagttt  >  2:249192‑M1/1‑9 (MQ=255)                                                                                                                                          |                                                                                                                                             AAAAATTTGCTAAATTTTGCCAATTTGGTAAAACAGTTGCATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCC  >  NC_000913_3_vch_tpiA/257770‑258047

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