Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,398,054	Δ7 bp	coding (643‑649/1848 nt)	mutL →	methyl‑directed mismatch repair protein

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4398054	4398060	7	48 [0]	[0] 48	mutL	methyl‑directed mismatch repair protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4398053	0 (0.000)	46 (0.860)	44/292	0.4	100%	coding (642/1848 nt)	mutL	methyl‑directed mismatch repair protein
	?	NC_000913	4398061 =	0 (0.000)					coding (650/1848 nt)	mutL	methyl‑directed mismatch repair protein

AATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGC  >  NC_000913/4397911‑4398209                                                                                                                                                       |                                                                                                                                                     aaaaccgaatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgg                                                                                                                                                      >  2:462917‑M2/151‑151 (MQ=255)      cgaatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCAcc                                                                                                                                                  >  2:640739‑M2/146‑150 (MQ=255)               ccaaattgatgagttcatccgccgcattgcgctggcgcgtttcgacgccacgatcaaactgtcgcataacggtaaaattgtgcgacagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTg                                                                                                                                        >  1:434798‑M2/137‑151 (MQ=255)                  cattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAAc                                                                                                                                     <  2:352765‑M2/18‑1 (MQ=255)                      gatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgagcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAgc                                                                                                                                 >  2:249476‑M2/130‑151 (MQ=255)                         gagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAgcgc                                                                                                                               <  2:458872‑M2/24‑1 (MQ=255)                           gatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTgg                                                                                                                            <  1:331825‑M2/27‑1 (MQ=255)                               atccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGAt                                                                                                                        >  2:353926‑M2/121‑151 (MQ=255)                                 ccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGAtt                                                                                                                       <  1:635369‑M2/32‑1 (MQ=255)                                 ccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTg                                                                                                                      >  1:320382‑M2/119‑151 (MQ=255)                                  cgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTg                                                                                                                      >  2:89214‑M2/118‑150 (MQ=255)                                    ccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGaa                                                                                                                    >  1:643024‑M2/116‑150 (MQ=255)                                     cgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATg                                                                                                                  >  2:104182‑M2/115‑151 (MQ=255)                                        attgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGc                                                                                                                >  2:518397‑M2/112‑150 (MQ=255)                                               tggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGc                                                                                                        >  1:587937‑M2/105‑151 (MQ=255)                                                         cgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGc                                                                                              >  1:186060‑M2/95‑151 (MQ=255)                                                            cgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAc                                                                                           >  2:352080‑M2/92‑151 (MQ=255)                                                            cgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAc                                                                                           <  1:658140‑M2/60‑1 (MQ=255)                                                               cacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAcgcg                                                                                        <  1:89191‑M2/63‑1 (MQ=255)                                                                 cgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGc                                                                                      >  1:547840‑M2/87‑151 (MQ=255)                                                                         ctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGc                                                                              >  2:570368‑M2/79‑151 (MQ=255)                                                                             cgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGa                                                                           <  2:821384‑M2/76‑1 (MQ=255)                                                                              gcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATc                                                                         >  2:319174‑M2/74‑151 (MQ=255)                                                                                  aacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCaa                                                                      >  1:702756‑M2/70‑150 (MQ=255)                                                                                  aacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCaa                                                                      >  1:568395‑M2/70‑150 (MQ=255)                                                                                    cggtaaaattgtgcgtcagtaccgcgcagtgccggaaggggggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAATTc                                                                   <  1:786583‑M2/84‑1 (MQ=255)                                                                                                gcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGccc                                                        <  1:551659‑M2/95‑1 (MQ=255)                                                                                                gcgtcagtaccgcgcagtgccggaaggcaggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGccc                                                        <  2:737482‑M2/95‑1 (MQ=255)                                                                                                    cagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCAc                                                    <  2:527561‑M2/99‑1 (MQ=255)                                                                                                         ccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGCACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAg                                              <  1:138364‑M2/105‑1 (MQ=255)                                                                                                         ccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGc                                                >  1:477488‑M2/47‑149 (MQ=255)                                                                                                          cgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGc                                                >  1:134641‑M2/46‑148 (MQ=255)                                                                                                          cgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCCTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGa                                             <  1:693617‑M2/106‑1 (MQ=255)                                                                                                             gcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAAt                                          <  1:322751‑M2/109‑1 (MQ=255)                                                                                                                    cggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTAt                                   >  2:522475‑M2/36‑151 (MQ=255)                                                                                                                       aaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGc                                >  2:136087‑M2/33‑151 (MQ=255)                                                                                                                        aggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCt                               >  2:500011‑M2/32‑151 (MQ=255)                                                                                                                           cgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACg                            >  1:119886‑M2/29‑151 (MQ=255)                                                                                                                            gggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGt                           >  1:582878‑M2/28‑151 (MQ=255)                                                                                                                               caaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGa                         >  1:809826‑M2/25‑150 (MQ=255)                                                                                                                                 aaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACg                      <  2:320405‑M2/129‑1 (MQ=255)                                                                                                                                   aagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTTAACGGt                    <  1:727840‑M2/131‑1 (MQ=255)                                                                                                                                    agaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTc                   <  1:605381‑M2/132‑1 (MQ=255)                                                                                                                                               taggcgcgGCTCCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATgcgc        >  2:151188‑M2/9‑151 (MQ=255)                                                                                                                                                aggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGcgcg       <  2:643066‑M2/144‑1 (MQ=255)                                                                                                                                                  gcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGAt     <  2:507570‑M2/146‑1 (MQ=255)                                                                                                                                                     cgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGc  >  1:212708‑M2/3‑151 (MQ=255)                                                                                                                                                     cgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGc  >  2:571659‑M2/3‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                     AATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATCGC  >  NC_000913/4397911‑4398209

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