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	57 [0]	[0] 59	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)	57 (1.440)	42/278	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

CATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCG  >  NC_000913/2856964‑2857233                                                                                                                                    |                                                                                                                                           cATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaa                                                                                                                                    <  2:590324‑M1/140‑9 (MQ=255)      aGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcc                                                                                                                               <  2:55941‑M1/140‑14 (MQ=255)      aGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcc                                                                                                                               <  2:42302‑M1/140‑14 (MQ=255)             gCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgta                                                                                                                        >  2:534930‑M1/1‑120 (MQ=255)               cGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaat                                                                                                                      >  1:618332‑M1/1‑118 (MQ=255)                 ttCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataa                                                                                                                    <  1:610338‑M1/140‑25 (MQ=255)                             ccTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgata                                                                                                        <  1:374283‑M1/140‑37 (MQ=255)                                     gATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttg                                                                                                <  2:336874‑M1/140‑45 (MQ=255)                                           cTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcg                                                                                          >  2:67785‑M1/1‑90 (MQ=255)                                              cGGTATCATGTGCGCCTCATGTCATTACAACTAAAATAACAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcggga                                                                                       <  2:322742‑M1/140‑54 (MQ=255)                                                 tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctgtcgggaata                                                                                    >  2:441253‑M1/1‑84 (MQ=255)                                                 tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaata                                                                                    <  2:522375‑M1/140‑57 (MQ=255)                                                 tATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaata                                                                                    >  1:522375‑M1/1‑84 (MQ=255)                                                   tCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACGCCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatatt                                                                                  <  2:618332‑M1/140‑59 (MQ=255)                                                        tgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacga                                                                             >  1:67212‑M1/1‑77 (MQ=255)                                                          cgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaagtcccgatcgtctttggctatcgggaatattaacgaca                                                                           >  1:171404‑M1/1‑75 (MQ=255)                                                           gcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaa                                                                          >  2:554100‑M1/1‑74 (MQ=255)                                                             cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaact                                                                        >  2:602018‑M1/1‑72 (MQ=255)                                                                aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgac                                                                     >  2:472882‑M1/1‑69 (MQ=255)                                                                    gATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaat                                                                 <  2:572537‑M1/140‑76 (MQ=255)                                                                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaata                                                                >  2:548902‑M1/1‑64 (MQ=255)                                                                     aTTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaata                                                                >  2:517485‑M1/1‑64 (MQ=255)                                                                      ttACAACGAAAATAACAACCATGACGCCCGATTTAATTTCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataa                                                               <  1:314442‑M1/140‑78 (MQ=255)                                                                       tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattgccgatagtttttggctgtcgggaatattaaggacaactgacgaataaa                                                              >  2:176506‑M1/1‑62 (MQ=255)                                                                       tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaa                                                              >  2:21183‑M1/1‑62 (MQ=255)                                                                       tACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtctttggctttcgggaatcttaagtacaactgacgaataaa                                                              >  2:308788‑M1/1‑62 (MQ=255)                                                                        gcaACGAATATAAAAACCATCACAACCAATTTAATATCAGAGAACGGGACATAACCCCATGcgcttgaaaagcctggtgtaataacaatccccgatagtcttttgctatcgggaatattaacgacaactgacgaatcaaa                                                             <  2:402835‑M1/139‑80 (MQ=255)                                                                           aCGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaataataacgacaaatgacgaataaaataa                                                          >  1:209190‑M1/1‑58 (MQ=255)                                                                              aaaaTAAAAACCATCGCACCCCGTTCAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaa                                                       <  1:308788‑M1/140‑86 (MQ=255)                                                                                aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                                     <  1:621597‑M1/140‑88 (MQ=255)                                                                                aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                                     <  1:21183‑M1/140‑88 (MQ=255)                                                                                aaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                                     <  2:590750‑M1/140‑88 (MQ=255)                                                                                  tAAAAACCATCACACTCCATTTAATATCATGGAACCGTACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                                   <  1:225089‑M1/140‑90 (MQ=255)                                                                                  tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                                   <  2:336698‑M1/140‑90 (MQ=255)                                                                                  tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                                   <  2:157096‑M1/140‑90 (MQ=255)                                                                                  tAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                                   <  1:67785‑M1/140‑90 (MQ=255)                                                                                   aaaaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccc                                                  >  2:435505‑M1/1‑50 (MQ=255)                                                                                       aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaacttacgaataaaataaaaacaccctgta                                              <  1:504637‑M1/140‑95 (MQ=255)                                                                                                 cccATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaa                                    >  2:110903‑M1/1‑36 (MQ=255)                                                                                                 cccATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtattttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaa                                    >  1:563791‑M1/1‑36 (MQ=255)                                                                                                    aTTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagct                                 <  1:548902‑M1/140‑108 (MQ=255)                                                                                                         atatCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattt                            >  2:648538‑M1/1‑28 (MQ=255)                                                                                                           atCAGGGAACCGGACATCACCCTATGcgcttcaagtgccaggtgtaatcacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                          <  2:192765‑M1/140‑115 (MQ=255)                                                                                                           atCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                          <  1:648538‑M1/140‑115 (MQ=255)                                                                                                              aGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaacaaaaacaccctgtagaataggtaagcttattttacag                       >  1:318273‑M1/1‑23 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagg                      >  2:410391‑M1/1‑22 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagg                      >  1:27214‑M1/1‑22 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagg                      >  2:612239‑M1/1‑22 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatagga                                      >  2:91157‑M1/1‑22 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataatagga                                      <  1:91157‑M1/124‑103 (MQ=255)                                                                                                                ggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaatacaaacaccgtgtgtaataggaacgcttattttacaggg                     >  1:214489‑M1/1‑21 (MQ=255)                                                                                                                      ggACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaac               >  1:475904‑M1/1‑15 (MQ=255)                                                                                                                        aCATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccttgtataaaaggcaaacttattttacagggtaaaacca             >  1:238794‑M1/1‑13 (MQ=255)                                                                                                                          aTAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatg           <  1:612239‑M1/140‑130 (MQ=255)                                                                                                                            aaCCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgcc         >  1:84188‑M1/1‑9 (MQ=255)                                                                                                                                 cATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgccatcta    >  1:50591‑M1/1‑4 (MQ=255)                                                                                                                                   tGcgcttgaagagcccggtgtactatcaactcccgatcgtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgccatctaca  <  1:176506‑M1/140‑139 (MQ=255)                                                                                                                                    |                                                                                                                                           CATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAAACGCGCG  >  NC_000913/2856964‑2857233

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