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	56 [0]	[0] 56	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)	56 (1.170)	42/276	0.1	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

GCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGT  >  NC_000913/2856963‑2857234                                                                                                                                     |                                                                                                                                            gCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttga                                                                                                                                      <  1:722671‑M1/140‑8 (MQ=255)     cAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagag                                                                                                                                  <  1:610163‑M1/140‑12 (MQ=255)       aGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcc                                                                                                                                <  2:368340‑M1/140‑14 (MQ=255)       aGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcc                                                                                                                                >  1:319024‑M1/1‑127 (MQ=255)         aaCTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctg                                                                                                                              <  2:319024‑M1/140‑16 (MQ=255)            tCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtg                                                                                                                           >  1:280623‑M1/1‑122 (MQ=255)                    cGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaatatca                                                                                                                   <  2:296560‑M1/140‑27 (MQ=255)                            ccccTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccga                                                                                                           <  1:687897‑M1/140‑35 (MQ=255)                            ccccTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccga                                                                                                           <  1:526820‑M1/140‑35 (MQ=255)                                       aTTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgc                                                                                                >  2:583804‑M1/1‑95 (MQ=255)                                                  tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaata                                                                                     <  2:42534‑M1/140‑57 (MQ=255)                                                       tgtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaac                                                                                >  2:37370‑M1/1‑79 (MQ=255)                                                          gcgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcttggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgac                                                                             >  2:201131‑M1/1‑76 (MQ=255)                                                           cgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgaca                                                                            >  2:587282‑M1/1‑75 (MQ=255)                                                                tATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactga                                                                       >  1:340416‑M1/1‑70 (MQ=255)                                                                 aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgac                                                                      >  2:130226‑M1/1‑69 (MQ=255)                                                                 aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgac                                                                      >  1:421030‑M1/1‑69 (MQ=255)                                                                  tgtgATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacg                                                                     >  1:42866‑M1/1‑68 (MQ=255)                                                                    tgATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaa                                                                   <  1:587282‑M1/140‑75 (MQ=255)                                                                        tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaa                                                               >  1:429632‑M1/1‑62 (MQ=255)                                                                        tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaa                                                               >  2:539150‑M1/1‑62 (MQ=255)                                                                            aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataa                                                           >  2:630885‑M1/1‑58 (MQ=255)                                                                            aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataa                                                           >  1:149461‑M1/1‑58 (MQ=255)                                                                              gAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaa                                                         <  1:801683‑M1/140‑85 (MQ=255)                                                                               aaaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaa                                                        <  2:429632‑M1/140‑86 (MQ=255)                                                                                aaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaac                                                       <  1:130226‑M1/140‑87 (MQ=255)                                                                                aaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaac                                                       >  1:713510‑M1/1‑54 (MQ=255)                                                                                 aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                                      <  2:149461‑M1/140‑88 (MQ=255)                                                                                 aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                                      <  1:37370‑M1/140‑88 (MQ=255)                                                                                   tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgatgaataaaataaaaacacc                                                    <  1:621595‑M1/140‑90 (MQ=255)                                                                                   tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                                    <  1:44369‑M1/140‑90 (MQ=255)                                                                                   tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                                    <  2:213863‑M1/140‑90 (MQ=255)                                                                                   tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                                    <  2:340416‑M1/140‑90 (MQ=255)                                                                                      aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctg                                                 <  2:520506‑M1/140‑93 (MQ=255)                                                                                      aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctg                                                 <  1:583804‑M1/140‑93 (MQ=255)                                                                                      aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctg                                                 <  1:356916‑M1/140‑93 (MQ=255)                                                                                      aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctg                                                 <  2:713510‑M1/140‑93 (MQ=255)                                                                                            tCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataat                                           >  1:682769‑M1/1‑42 (MQ=255)                                                                                             cacaCCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataata                                          <  2:125847‑M1/140‑100 (MQ=255)                                                                                              acacCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatag                                         >  1:678446‑M1/1‑40 (MQ=255)                                                                                               cacCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatagg                                        <  1:519830‑M1/140‑102 (MQ=255)                                                                                                acCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatagga                                       <  2:678446‑M1/140‑103 (MQ=255)                                                                                                            atCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                           <  2:421030‑M1/140‑115 (MQ=255)                                                                                                            atCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                           <  1:665622‑M1/140‑115 (MQ=255)                                                                                                            atCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                           <  2:522724‑M1/140‑115 (MQ=255)                                                                                                              cAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttaca                         >  2:660590‑M1/1‑24 (MQ=255)                                                                                                                gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagg                       >  2:234794‑M1/1‑22 (MQ=255)                                                                                                                gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagg                       >  2:185348‑M1/1‑22 (MQ=255)                                                                                                                   aaCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggta                    >  2:719474‑M1/1‑19 (MQ=255)                                                                                                                   aaCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggta                    <  1:126697‑M1/140‑122 (MQ=255)                                                                                                                      cGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaa                 <  1:539150‑M1/140‑125 (MQ=255)                                                                                                                        gACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaacc               <  2:305444‑M1/140‑127 (MQ=255)                                                                                                                          cATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccat             <  2:280623‑M1/140‑129 (MQ=255)                                                                                                                            tAACCCCATGcgcttgaagagcttggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgc           <  1:201131‑M1/140‑131 (MQ=255)                                                                                                                             aaCCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgcc          <  2:42866‑M1/140‑132 (MQ=255)                                                                                                                                     gcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgccatctacac  >  1:706008‑M1/1‑1 (MQ=255)                                                                                                                                     |                                                                                                                                            GCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCGT  >  NC_000913/2856963‑2857234

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