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

AAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCG  >  NC_000913/4397904‑4398202                                                                                                                                                      |                                                                                                                                                      aaaaCCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGg                                                                                                                                                      >  2:462917‑M1/1‑150 (MQ=255)      cGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcacc                                                                                                                                                  >  2:640739‑M1/1‑145 (MQ=255)               ccAAATTGATGAGTTCATCCGCCGCATTGCGCTGGCGCGTTTCGACGCCACGATCAAACTGTCGCATAACGGTAAAATTGTGCGACAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttg                                                                                                                                        >  1:434798‑M1/1‑136 (MQ=255)                  caTTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaac                                                                                                                                     <  2:352765‑M1/151‑19 (MQ=255)                      gatgaGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGAGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagc                                                                                                                                 >  2:249476‑M1/1‑129 (MQ=255)                         gagaTCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgc                                                                                                                               <  2:458872‑M1/150‑25 (MQ=255)                           gaTCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctgg                                                                                                                            <  1:331825‑M1/151‑28 (MQ=255)                               atcCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgat                                                                                                                        >  2:353926‑M1/1‑120 (MQ=255)                                 ccgccgCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattg                                                                                                                      >  1:320382‑M1/1‑118 (MQ=255)                                 ccgccgCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgatt                                                                                                                       <  1:635369‑M1/150‑33 (MQ=255)                                  cgccgcATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattg                                                                                                                      >  2:89214‑M1/1‑117 (MQ=255)                                    ccgcATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaa                                                                                                                    >  1:643024‑M1/1‑115 (MQ=255)                                     cgcATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatg                                                                                                                  >  2:104182‑M1/1‑114 (MQ=255)                                        aTTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggc                                                                                                                >  2:518397‑M1/1‑111 (MQ=255)                                               tGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggc                                                                                                        >  1:587937‑M1/1‑104 (MQ=255)                                                         cGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgc                                                                                              >  1:186060‑M1/1‑94 (MQ=255)                                                            cGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctac                                                                                           >  2:352080‑M1/1‑91 (MQ=255)                                                            cGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctac                                                                                           <  1:658140‑M1/151‑61 (MQ=255)                                                               cACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcg                                                                                        <  1:89191‑M1/151‑64 (MQ=255)                                                                 cGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggc                                                                                      >  1:547840‑M1/1‑86 (MQ=255)                                                                         cTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggc                                                                              >  2:570368‑M1/1‑78 (MQ=255)                                                                             cGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccga                                                                           <  2:821384‑M1/150‑77 (MQ=255)                                                                              gCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatc                                                                         >  2:319174‑M1/1‑73 (MQ=255)                                                                                  aaCGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaa                                                                      >  1:702756‑M1/1‑69 (MQ=255)                                                                                  aaCGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaa                                                                      >  1:568395‑M1/1‑69 (MQ=255)                                                                                    cGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGGGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaattc                                                                   <  1:786583‑M1/151‑85 (MQ=255)                                                                                                gCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgccc                                                        <  1:551659‑M1/150‑96 (MQ=255)                                                                                                gCGTCAGTACCGCGCAGTGCCGGAAGGCAGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgccc                                                        <  2:737482‑M1/150‑96 (MQ=255)                                                                                                    cAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcac                                                    <  2:527561‑M1/150‑100 (MQ=255)                                                                                                         ccGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgcacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcag                                              <  1:138364‑M1/151‑106 (MQ=255)                                                                                                         ccGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggc                                                >  1:477488‑M1/1‑46 (MQ=255)                                                                                                          cgcgCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggc                                                >  1:134641‑M1/1‑45 (MQ=255)                                                                                                          cgcgCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcctttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcaga                                             <  1:693617‑M1/151‑107 (MQ=255)                                                                                                             gcAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaat                                          <  1:322751‑M1/151‑110 (MQ=255)                                                                                                                    cGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtat                                   >  2:522475‑M1/1‑35 (MQ=255)                                                                                                                       aaGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgc                                >  2:136087‑M1/1‑32 (MQ=255)                                                                                                                        aGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgct                               >  2:500011‑M1/1‑31 (MQ=255)                                                                                                                           cGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacg                            >  1:119886‑M1/1‑28 (MQ=255)                                                                                                                            gggCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgt                           >  1:582878‑M1/1‑27 (MQ=255)                                                                                                                               cAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtga                         >  1:809826‑M1/1‑24 (MQ=255)                                                                                                                                 aaaaGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacg                      <  2:320405‑M1/151‑130 (MQ=255)                                                                                                                                   aaGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgttaacggt                    <  1:727840‑M1/151‑132 (MQ=255)                                                                                                                                    aGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtc                   <  1:605381‑M1/151‑133 (MQ=255)                                                                                                                                               tAGGCGCGgctccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgatgcgc        >  2:151188‑M1/1‑8 (MQ=255)                                                                                                                                                aGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgatgcgcg       <  2:643066‑M1/151‑145 (MQ=255)                                                                                                                                                  gcgcGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgatgcgcgat     <  2:507570‑M1/151‑147 (MQ=255)                                                                                                                                                     cggcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgatgcgcgatcgc  >  1:212708‑M1/1‑2 (MQ=255)                                                                                                                                                     cggcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgatgcgcgatcgc  >  2:571659‑M1/1‑2 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AAAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATGCG  >  NC_000913/4397904‑4398202

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