Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	2,857,096	Δ2,538 bp	coding (4‑2541/2562 nt)	mutS →	DNA mismatch repair protein MutS

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	2857096	2859633	2538	22 [0]	[0] 22	mutS	DNA mismatch repair protein MutS

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 2857095	0 (0.000)	22 (0.890)	19/274	0.3	100%	coding (3/2562 nt)	mutS	DNA mismatch repair protein MutS
	?	NC_000913	2859634 =	0 (0.000)					coding (2542/2562 nt)	mutS	DNA mismatch repair protein MutS

TACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACA  >  NC_000913/2859507‑2859749                                                                                                                                |                                                                                                                    gaactcgccgttcgcttcttcccctgaaatgattaactccggatatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                         <  2:326738‑M2/13‑1 (MQ=255)       cgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGgtgt                                                                                                   <  1:310210‑M2/19‑1 (MQ=255)        gccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGta                                                                                                  >  2:209438‑M2/121‑140 (MQ=255)              cgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACa                                                                                            <  2:285720‑M2/26‑1 (MQ=255)                                       tccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCgg                                                                   <  1:369877‑M2/51‑1 (MQ=255)                                           gtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCtt                                                                             <  2:372902‑M2/41‑1 (MQ=255)                                           gtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCtt                                                                             >  1:372902‑M2/86‑126 (MQ=255)                                                        cttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACt                                                  >  2:221635‑M2/73‑140 (MQ=255)                                                        cttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACt                                                  >  2:205577‑M2/73‑140 (MQ=255)                                                                 ttacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataa                                         >  2:168761‑M2/64‑140 (MQ=255)                                                                    caacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaaat                                      >  2:369465‑M2/61‑140 (MQ=255)                                                                     aacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataa                                         >  2:98018‑M2/60‑136 (MQ=255)                                                                     aacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataa                                         <  1:98018‑M2/77‑1 (MQ=255)                                                                       cgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAataaaataaa                                   >  2:159102‑M2/58‑140 (MQ=255)                                                                        gaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAAtaaaataaaa                                  <  1:175484‑M2/84‑1 (MQ=255)                                                                             taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                             <  2:160513‑M2/89‑1 (MQ=255)                                                                             taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                             <  1:369465‑M2/89‑1 (MQ=255)                                                                                aaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTg                          <  1:55935‑M2/92‑1 (MQ=255)                                                                                aaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTg                          <  1:221635‑M2/92‑1 (MQ=255)                                                                                          accccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGa                <  2:65052‑M2/102‑1 (MQ=255)                                                                                                    atatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAttt      >  2:138512‑M2/29‑140 (MQ=255)                                                                                                        cagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACa  >  2:306845‑M2/25‑140 (MQ=255)                                                                                                                                |                                                                                                                    TACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACA  >  NC_000913/2859507‑2859749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑