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	33 [0]	[0] 33	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)	32 (0.580)	30/274	0.6	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

GGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCA  >  NC_000913/2856956‑2857224                                                                                                                                            |                                                                                                                                  ggAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATg                                                                                                                                   <  2:921862/140‑1 (MQ=255)         cATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaa                                                                                                                           <  2:622968‑M1/140‑9 (MQ=255)             aaGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagc                                                                                                                       <  1:157682‑M1/140‑13 (MQ=255)                   tCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtg                                                                                                                 >  1:605194‑M1/1‑122 (MQ=255)                      ccGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaa                                                                                                              <  1:685488‑M1/140‑22 (MQ=255)                          tCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataa                                                                                                           <  1:289077‑M1/139‑25 (MQ=255)                          tCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataa                                                                                                           >  2:289077‑M1/1‑115 (MQ=255)                             cttcttCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaat                                                                                                       >  2:190857‑M1/1‑112 (MQ=255)                                     ccTGAAATGATTAACTCCGGTATCATGTGCGTCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGCAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgata                                                                                               <  2:787257‑M1/140‑37 (MQ=255)                                               ttAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgct                                                                                     <  1:190857‑M1/140‑47 (MQ=255)                                                              tgtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaac                                                                      >  2:14982‑M1/1‑79 (MQ=255)                                                                tgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacga                                                                    >  1:390279‑M1/1‑77 (MQ=255)                                                                  cgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgaca                                                                  >  2:150219‑M1/1‑75 (MQ=255)                                                                   gcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaa                                                                 >  2:405917‑M1/1‑74 (MQ=255)                                                                     cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaact                                                               >  1:837873‑M1/1‑72 (MQ=255)                                                                            gATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaat                                                        <  1:405917‑M1/140‑76 (MQ=255)                                                                               tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaa                                                     >  2:354016‑M1/1‑62 (MQ=255)                                                                               tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaa                                                     >  2:272364‑M1/1‑62 (MQ=255)                                                                                     gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaa                                               <  1:150219‑M1/140‑85 (MQ=255)                                                                                     gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaa                                               <  2:837873‑M1/140‑85 (MQ=255)                                                                                        aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                            <  1:272364‑M1/140‑88 (MQ=255)                                                                                         aTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacac                                           <  1:696335‑M1/140‑89 (MQ=255)                                                                                             aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctg                                       <  2:325827‑M1/140‑93 (MQ=255)                                                                                                   tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataat                                 >  1:28990‑M1/1‑42 (MQ=255)                                                                                                    cacaCCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataata                                >  1:688667‑M1/1‑41 (MQ=255)                                                                                                          ccATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgta                                     >  2:348589‑M1/1‑35 (MQ=255)                                                                                                          ccATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgta                                     <  1:348589‑M1/129‑95 (MQ=255)                                                                                                              ttAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagctta                      >  1:1830‑M1/1‑31 (MQ=255)                                                                                                                   atCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaatcttatttta                 >  1:587822‑M1/1‑26 (MQ=255)                                                                                                                    tCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttac                <  1:14982‑M1/140‑116 (MQ=255)                                                                                                                             cGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaa       <  2:605194‑M1/140‑125 (MQ=255)                                                                                                                               gACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaacc     >  1:231003‑M1/1‑14 (MQ=255)                                                                                                                                  aTAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgatgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatg  <  2:390279‑M1/140‑130 (MQ=255)                                                                                                                                            |                                                                                                                                  GGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCA  >  NC_000913/2856956‑2857224

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