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	52 [0]	[0] 52	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)	52 (1.070)	48/292	0.1	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

ATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATC  >  NC_000913/4397912‑4398207                                                                                                                                                      |                                                                                                                                                   aaaccgaatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGc                                                                                                                                                   >  2:540820‑M2/150‑151 (MQ=255)   accgaatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCAc                                                                                                                                                 <  2:153532‑M2/4‑1 (MQ=255)    ccgaatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCAcc                                                                                                                                                <  2:320630‑M2/5‑1 (MQ=255)       aatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCt                                                                                                                                             >  2:179372‑M2/144‑151 (MQ=255)       aatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCt                                                                                                                                             >  2:70748‑M2/144‑151 (MQ=255)       aatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCt                                                                                                                                             <  2:131149‑M2/8‑1 (MQ=255)       aatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCt                                                                                                                                             >  2:193155‑M2/144‑151 (MQ=255)       aatttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCt                                                                                                                                             >  2:353925‑M2/144‑151 (MQ=255)         tttaaccacattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCtt                                                                                                                                            <  1:465385‑M2/9‑1 (MQ=255)                acattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGa                                                                                                                                     >  2:558634‑M2/135‑150 (MQ=255)                 cattgatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAAc                                                                                                                                   <  1:488811‑M2/18‑1 (MQ=255)                     gatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAgc                                                                                                                               <  2:705443‑M2/22‑1 (MQ=255)                     gatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAgc                                                                                                                               <  1:235302‑M2/22‑1 (MQ=255)                     gatgagatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAg                                                                                                                                <  2:649425‑M2/21‑1 (MQ=255)                         agatcatccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTg                                                                                                                           <  2:66859‑M2/26‑1 (MQ=255)                             catccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCg                                                                                                                        <  1:509919‑M2/29‑1 (MQ=255)                                ccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGAtt                                                                                                                     <  1:523764‑M2/32‑1 (MQ=255)                                ccgccgcattgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTg                                                                                                                    >  2:241784‑M2/119‑151 (MQ=255)                                        ttgcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCa                                                                                                             >  2:88725‑M2/111‑150 (MQ=255)                                          gcgctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAaca                                                                                                          <  1:575062‑M2/43‑1 (MQ=255)                                             ctggcgcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACgg                                                                                                       >  1:248801‑M2/106‑151 (MQ=255)                                                  gcgtttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGAtc                                                                                                  >  1:236348‑M2/101‑151 (MQ=255)                                                      ttcgacgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCa                                                                                               <  2:674374‑M2/54‑1 (MQ=255)                                                           cgtcacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAc                                                                                         >  1:241810‑M2/92‑151 (MQ=255)                                                              cacgatcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTAcgcg                                                                                      >  2:487217‑M2/89‑151 (MQ=255)                                                                  atcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTg                                                                                  <  2:722314‑M2/67‑1 (MQ=255)                                                                   tcaacctgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTgg                                                                                 >  2:360845‑M2/84‑151 (MQ=255)                                                                         tgtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGcc                                                                           >  2:530104‑M2/78‑151 (MQ=255)                                                                          gtcgcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCg                                                                          >  2:197045‑M2/77‑151 (MQ=255)                                                                             gcataacggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATc                                                                       <  2:161149‑M2/78‑1 (MQ=255)                                                                                  acggtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAAt                                                                  >  2:597137‑M2/69‑151 (MQ=255)                                                                                     gtaaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATcac                                                               >  2:707373‑M2/66‑151 (MQ=255)                                                                                       aaaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCacac                                                             >  2:551558‑M2/64‑151 (MQ=255)                                                                                        aaattgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCacac                                                             <  2:525450‑M2/88‑1 (MQ=255)                                                                                          attgtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCAcaccac                                                          <  2:527329‑M2/91‑1 (MQ=255)                                                                                             gtgcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGcc                                                       <  1:38648‑M2/94‑1 (MQ=255)                                                                                               gcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGccc                                                      <  1:283875‑M2/95‑1 (MQ=255)                                                                                               gcgtcagtaccgcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCg                                                     <  1:237046‑M2/96‑1 (MQ=255)                                                                                                cgtcagtaccgcgcagtgccggaaggcgggaaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGc                                                    <  2:671434‑M2/97‑1 (MQ=255)                                                                                                          gcgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAg                                            >  1:344212‑M2/45‑149 (MQ=255)                                                                                                           cgcagtgccggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGaaa                                         <  1:457220‑M2/108‑1 (MQ=255)                                                                                                                   cggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTAt                                 >  1:565822‑M2/36‑151 (MQ=255)                                                                                                                    ggaaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTAtt                                >  1:741230‑M2/35‑151 (MQ=255)                                                                                                                      aaggcgggcaaaaagaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGc                              >  1:97205‑M2/33‑151 (MQ=255)                                                                                                                                   agaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTc                 >  1:378443‑M2/20‑151 (MQ=255)                                                                                                                                   agaacggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTc                 <  2:763938‑M2/132‑1 (MQ=255)                                                                                                                                        ggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCatg            <  1:450960‑M2/137‑1 (MQ=255)                                                                                                                                        ggcgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCatg            >  2:418012‑M2/15‑151 (MQ=255)                                                                                                                                          cgcttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCatgat          >  2:314363‑M2/13‑151 (MQ=255)                                                                                                                                             ttaggcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATgcg       >  2:733852‑M2/10‑151 (MQ=255)                                                                                                                                                 gcgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGAt   <  2:360024‑M2/146‑1 (MQ=255)                                                                                                                                                  cgcgGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCGACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCATAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATc  >  2:553563‑M2/5‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                   ATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCGCGATC  >  NC_000913/4397912‑4398207

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