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	58 [0]	[0] 58	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)	58 (1.640)	45/264	0.0	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

CATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAA  >  NC_000913/2856964‑2857226                                                                                                                                    |                                                                                                                                    cATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGc                                                                                                                                    <  1:600353‑M1/133‑2 (MQ=255)   tCAAGAACGCGCCGTCCGGTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgc                                                                                                                                  <  1:230287‑M1/133‑4 (MQ=255)       gAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttga                                                                                                                              <  2:99246‑M1/133‑8 (MQ=255)       gAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttga                                                                                                                              >  1:536251‑M1/1‑126 (MQ=255)        aaCTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaa                                                                                                                             <  2:609635‑M1/133‑9 (MQ=255)           tCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagag                                                                                                                          >  1:504907‑M1/1‑122 (MQ=255)            cGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagc                                                                                                                         >  1:501480‑M1/1‑121 (MQ=255)             gCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcc                                                                                                                        >  1:559398‑M1/1‑120 (MQ=255)              ccGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcct                                                                                                                       <  2:499147‑M1/133‑15 (MQ=255)                 ttCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggt                                                                                                                    >  1:540633‑M1/1‑116 (MQ=255)                   cGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaa                                                                                                                             <  1:482820‑M1/122‑9 (MQ=255)                   cGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaa                                                                                                                             >  2:482820‑M1/1‑114 (MQ=255)                       tcttcCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaata                                                                                                              >  1:517721‑M1/1‑110 (MQ=255)                           ccccTGAAATGATTAACTCCGGTACCTTGTGCCCCTTATGTTATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaa                                                                                                          <  2:293192‑M1/133‑28 (MQ=255)                                  aaTGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccga                                                                                                   <  2:504907‑M1/133‑35 (MQ=255)                                  aaTGATTAACTCCGGCATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGCACCGGACATAACCCCATGctcttgaagagcctggtgaaataacaattcccga                                                                                                   >  2:162024‑M1/1‑99 (MQ=255)                                              cGGTATCATGTGCGCTTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgct                                                                                       <  2:236634‑M1/133‑47 (MQ=255)                                                     aTGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcggga                                                                                <  2:501480‑M1/133‑54 (MQ=255)                                                      tgtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaa                                                                               >  2:327922‑M1/1‑79 (MQ=255)                                                        tgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaata                                                                             >  1:92323‑M1/1‑77 (MQ=255)                                                           gcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctattgggaatatta                                                                          >  2:270119‑M1/1‑74 (MQ=255)                                                           gcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatatta                                                                          >  1:69429‑M1/1‑74 (MQ=255)                                                             cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaac                                                                        <  2:536251‑M1/133‑62 (MQ=255)                                                             cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaac                                                                        >  2:529814‑M1/1‑72 (MQ=255)                                                               tATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacga                                                                      >  2:604640‑M1/1‑70 (MQ=255)                                                                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactg                                                                >  1:600020‑M1/1‑64 (MQ=255)                                                                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactg                                                                >  2:391802‑M1/1‑64 (MQ=255)                                                                       tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgac                                                              >  2:94828‑M1/1‑62 (MQ=255)                                                                       tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtctcttgctatcggggatattaacgccaactgac                                                              >  1:306602‑M1/1‑62 (MQ=255)                                                                         cAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacga                                                            >  2:444525‑M1/1‑60 (MQ=255)                                                                           aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaat                                                          <  2:69429‑M1/133‑76 (MQ=255)                                                                           aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgatatcgggcatagtaacgacaactgacgaat                                                          >  1:302638‑M1/1‑58 (MQ=255)                                                                                    aaaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaa                                                 <  2:565218‑M1/133‑85 (MQ=255)                                                                                       aCCATCACTCCCCATTTACTATTAGGGAACTGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcggtaatattaacgacaactgacgaataaaataaaaaca                                              <  1:162024‑M1/133‑88 (MQ=255)                                                                                       aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  2:600020‑M1/133‑88 (MQ=255)                                                                                       aCCATCACACCCCATATAATATCAAGGACCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  2:230316‑M1/133‑88 (MQ=255)                                                                                         cttcaCCCCACATTTAATATCAGGGAGCGGGAAATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                            <  2:306602‑M1/131‑90 (MQ=255)                                                                                         cATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                            <  2:592425‑M1/133‑90 (MQ=255)                                                                                               acCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtat                                      <  1:529814‑M1/133‑96 (MQ=255)                                                                                                   cATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataata                                  >  2:405325‑M1/1‑34 (MQ=255)                                                                                                     tttAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatagg                                <  1:46533‑M1/133‑102 (MQ=255)                                                                                                     tttAATATCAGGGAAACGGACATATCCCCATGcgcttgaagtgcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatagg                                <  2:191469‑M1/133‑102 (MQ=255)                                                                                                      ttAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatagga                               <  1:7727‑M1/133‑103 (MQ=255)                                                                                                        aaTATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaa                             <  1:391802‑M1/133‑105 (MQ=255)                                                                                                            tCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagctt                         <  1:270119‑M1/133‑109 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaccccctgtatgatagaaaagcttatt                      >  2:159975‑M1/1‑22 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatt                      >  1:10577‑M1/1‑22 (MQ=255)                                                                                                                  aaCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                   <  1:444525‑M1/133‑115 (MQ=255)                                                                                                                  aaCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                   <  2:540633‑M1/133‑115 (MQ=255)                                                                                                                  aaCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                   <  2:559398‑M1/133‑115 (MQ=255)                                                                                                                  aaCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                   <  2:92323‑M1/133‑115 (MQ=255)                                                                                                                  aaCCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                   <  2:610075‑M1/133‑115 (MQ=255)                                                                                                                    ccGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaataccctgtctaataggaaagcttattttaca                 >  1:275850‑M1/1‑17 (MQ=255)                                                                                                                           ttaccCCATTcgcttgaagagcctggtgaaatcacaattccagatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaa          <  2:302638‑M1/131‑124 (MQ=255)                                                                                                                              ccccATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaacc       >  1:549406‑M1/1‑7 (MQ=255)                                                                                                                                ccATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccat     <  1:94828‑M1/133‑129 (MQ=255)                                                                                                                                  aTGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgc   <  2:549406‑M1/133‑131 (MQ=255)                                                                                                                                   tGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgcc  >  2:622392‑M1/1‑2 (MQ=255)                                                                                                                                    |                                                                                                                                    CATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAA  >  NC_000913/2856964‑2857226

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