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

AAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATG  >  NC_000913/4397905‑4398200                                                                                                                                                     |                                                                                                                                                    aaaCCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgc                                                                                                                                                   >  2:540820‑M1/1‑149 (MQ=255)   aCCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcac                                                                                                                                                 <  2:153532‑M1/151‑5 (MQ=255)    ccGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcacc                                                                                                                                                <  2:320630‑M1/151‑6 (MQ=255)       aaTTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgct                                                                                                                                             >  2:179372‑M1/1‑143 (MQ=255)       aaTTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgct                                                                                                                                             >  2:70748‑M1/1‑143 (MQ=255)       aaTTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgct                                                                                                                                             <  2:131149‑M1/151‑9 (MQ=255)       aaTTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgct                                                                                                                                             >  2:193155‑M1/1‑143 (MQ=255)       aaTTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgct                                                                                                                                             >  2:353925‑M1/1‑143 (MQ=255)         tttAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgctt                                                                                                                                            <  1:465385‑M1/150‑10 (MQ=255)                acaTTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttga                                                                                                                                     >  2:558634‑M1/1‑134 (MQ=255)                 caTTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaac                                                                                                                                   <  1:488811‑M1/151‑19 (MQ=255)                     gatgaGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagc                                                                                                                               <  2:705443‑M1/151‑23 (MQ=255)                     gatgaGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagc                                                                                                                               <  1:235302‑M1/151‑23 (MQ=255)                     gatgaGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaag                                                                                                                                <  2:649425‑M1/150‑22 (MQ=255)                         agaTCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctg                                                                                                                           <  2:66859‑M1/151‑27 (MQ=255)                             catcCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcg                                                                                                                        <  1:509919‑M1/150‑30 (MQ=255)                                ccgccgCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattg                                                                                                                    >  2:241784‑M1/1‑118 (MQ=255)                                ccgccgCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgatt                                                                                                                     <  1:523764‑M1/150‑33 (MQ=255)                                        ttGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggca                                                                                                             >  2:88725‑M1/1‑110 (MQ=255)                                          gcgcTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaaca                                                                                                          <  1:575062‑M1/151‑44 (MQ=255)                                             cTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacgg                                                                                                       >  1:248801‑M1/1‑105 (MQ=255)                                                  gcgTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatc                                                                                                  >  1:236348‑M1/1‑100 (MQ=255)                                                      ttCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctca                                                                                               <  2:674374‑M1/150‑55 (MQ=255)                                                           cGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctac                                                                                         >  1:241810‑M1/1‑91 (MQ=255)                                                              cACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcg                                                                                      >  2:487217‑M1/1‑88 (MQ=255)                                                                  aTCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctg                                                                                  <  2:722314‑M1/151‑68 (MQ=255)                                                                   tCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgg                                                                                 >  2:360845‑M1/1‑83 (MQ=255)                                                                         tGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggcc                                                                           >  2:530104‑M1/1‑77 (MQ=255)                                                                          gTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccg                                                                          >  2:197045‑M1/1‑76 (MQ=255)                                                                             gCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatc                                                                       <  2:161149‑M1/151‑79 (MQ=255)                                                                                  aCGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaat                                                                  >  2:597137‑M1/1‑68 (MQ=255)                                                                                     gTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcac                                                               >  2:707373‑M1/1‑65 (MQ=255)                                                                                       aaaaTTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacac                                                             >  2:551558‑M1/1‑63 (MQ=255)                                                                                        aaaTTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacac                                                             <  2:525450‑M1/150‑89 (MQ=255)                                                                                          aTTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccac                                                          <  2:527329‑M1/151‑92 (MQ=255)                                                                                             gtgCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcc                                                       <  1:38648‑M1/151‑95 (MQ=255)                                                                                               gCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccg                                                     <  1:237046‑M1/151‑97 (MQ=255)                                                                                               gCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgccc                                                      <  1:283875‑M1/150‑96 (MQ=255)                                                                                                cGTCAGTACCGCGCAGTGCCGGAAGGCGGGAAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgc                                                    <  2:671434‑M1/151‑98 (MQ=255)                                                                                                          gcgcAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcag                                            >  1:344212‑M1/1‑44 (MQ=255)                                                                                                           cgcAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaa                                         <  1:457220‑M1/151‑109 (MQ=255)                                                                                                                   cGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtat                                 >  1:565822‑M1/1‑35 (MQ=255)                                                                                                                    ggAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtatt                                >  1:741230‑M1/1‑34 (MQ=255)                                                                                                                      aaGGCGGGCAAAAAGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgc                              >  1:97205‑M1/1‑32 (MQ=255)                                                                                                                                   aGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtc                 >  1:378443‑M1/1‑19 (MQ=255)                                                                                                                                   aGAACGGCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtc                 <  2:763938‑M1/151‑133 (MQ=255)                                                                                                                                        ggCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatg            <  1:450960‑M1/151‑138 (MQ=255)                                                                                                                                        ggCGCTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatg            >  2:418012‑M1/1‑14 (MQ=255)                                                                                                                                          cgcTTAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgat          >  2:314363‑M1/1‑12 (MQ=255)                                                                                                                                             ttAGGCGCGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgatgcg       >  2:733852‑M1/1‑9 (MQ=255)                                                                                                                                                 gcgcGgcaccgcttttcttgaacaagcgctggcgattgaatggcaacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcagaaattcagtattgctacgtgaacggtcgcatgatgcgcgat   <  2:360024‑M1/151‑147 (MQ=255)                                                                                                                                                  cgcggcaccgcttttcttgaacaagcgctggcgattgaatggcgacacggcgatctcacgctacgcggctgggtggccgatccaaatcacaccacgcccgcactggcataaattcagtattgctacgtgaacggtcgcatgatgcgcgatc  >  2:553563‑M1/1‑4 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AAACCGAATTTAACCACATTGATGAGATCATCCGCCGCATTGCGCTGGCGCGTTTCGACGTCACGATCAACCTGTCGCATAACGGTAAAATTGTGCGTCAGTACCGCGCAGTGCCGGAAGGCGGGCAAAAAGAACGGCGCTTAGGCGCGATTTGCGGCACCGCTTTTCTTGAACAAGCGCTGGCGATTGAATGGCAACACGGCGATCTCACGCTACGCGGCTGGGTGGCCGATCCAAATCACACCACGCCCGCACTGGCAGAAATTCAGTATTGCTACGTGAACGGTCGCATGATG  >  NC_000913/4397905‑4398200

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