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

GCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACG  >  NC_000913/2859502‑2859773                                                                                                                                     |                                                                                                                                            catcaagaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGaa                                                                                                                                      <  2:590324‑M2/8‑1 (MQ=255)      agaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                                 <  2:42302‑M2/13‑1 (MQ=255)      agaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                                 <  2:55941‑M2/13‑1 (MQ=255)             gccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGta                                                                                                                          >  2:534930‑M2/121‑140 (MQ=255)               cgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGtaat                                                                                                                        >  1:618332‑M2/119‑140 (MQ=255)                 ttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGtaataa                                                                                                                      <  1:610338‑M2/24‑1 (MQ=255)                             cctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATa                                                                                                          <  1:374283‑M2/36‑1 (MQ=255)                                     gattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTg                                                                                                  <  2:336874‑M2/44‑1 (MQ=255)                                           ctccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCg                                                                                            >  2:67785‑M2/91‑140 (MQ=255)                                              cggtatcatgtgcgcctcatgtcattacaactaaaataacaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGa                                                                                         <  2:322742‑M2/53‑1 (MQ=255)                                                 tatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTGTCGGGAata                                                                                      >  2:441253‑M2/85‑140 (MQ=255)                                                 tatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAata                                                                                      <  2:522375‑M2/56‑1 (MQ=255)                                                 tatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAata                                                                                      >  1:522375‑M2/85‑140 (MQ=255)                                                   tcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacgccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATAtt                                                                                    <  2:618332‑M2/58‑1 (MQ=255)                                                        tgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGa                                                                               >  1:67212‑M2/78‑140 (MQ=255)                                                          cgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAAGTCCCGATCGTCTTTGGCTATCGGGAATATTAACGACa                                                                             >  1:171404‑M2/76‑140 (MQ=255)                                                           gccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACaa                                                                            >  2:554100‑M2/75‑140 (MQ=255)                                                             cttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACt                                                                          >  2:602018‑M2/73‑140 (MQ=255)                                                                atgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGAc                                                                       >  2:472882‑M2/70‑140 (MQ=255)                                                                    gattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaat                                                                   <  2:572537‑M2/75‑1 (MQ=255)                                                                     attacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaata                                                                  >  2:517485‑M2/65‑140 (MQ=255)                                                                     attacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaata                                                                  >  2:548902‑M2/65‑140 (MQ=255)                                                                      ttacaacgaaaataacaaccatgacgcccgatttaatttcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataa                                                                 <  1:314442‑M2/77‑1 (MQ=255)                                                                       tacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTGCCGATAGTTTTTGGCTGTCGGGAATATTAAGGACAACTGACGaataaa                                                                >  2:176506‑M2/63‑140 (MQ=255)                                                                       tacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaa                                                                >  2:21183‑M2/63‑140 (MQ=255)                                                                       tacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTGGCTTTCGGGAATCTTAAGTACAACTGACGaataaa                                                                >  2:308788‑M2/63‑140 (MQ=255)                                                                        gcaacgaatataaaaaccatcacaaccaatttaatatcagagaacgggacataaccccatgCGCTTGAAAAGCCTGGTGTAATAACAATCCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaatcaaa                                                               <  2:402835‑M2/79‑1 (MQ=255)                                                                           acgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATAATAACGACAAATGACGaataaaataa                                                            >  1:209190‑M2/59‑140 (MQ=255)                                                                              aaaataaaaaccatcgcaccccgttcaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATaaaaa                                                         <  1:308788‑M2/85‑1 (MQ=255)                                                                                aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  1:21183‑M2/87‑1 (MQ=255)                                                                                aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  1:621597‑M2/87‑1 (MQ=255)                                                                                aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  2:590750‑M2/87‑1 (MQ=255)                                                                                  taaaaaccatcacactccatttaatatcatggaaccgtacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  1:225089‑M2/89‑1 (MQ=255)                                                                                  taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  2:157096‑M2/89‑1 (MQ=255)                                                                                  taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  1:67785‑M2/89‑1 (MQ=255)                                                                                  taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  2:336698‑M2/89‑1 (MQ=255)                                                                                   aaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAccc                                                    >  2:435505‑M2/51‑140 (MQ=255)                                                                                       accatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTTACGAATAAAATAAAAACACCCTGta                                                <  1:504637‑M2/94‑1 (MQ=255)                                                                                                 cccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGaaa                                      >  2:110903‑M2/37‑140 (MQ=255)                                                                                                 cccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTATTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGaaa                                      >  1:563791‑M2/37‑140 (MQ=255)                                                                                                    atttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCt                                   <  1:548902‑M2/107‑1 (MQ=255)                                                                                                         atatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAttt                              >  2:648538‑M2/29‑140 (MQ=255)                                                                                                           atcagggaaccggacatcaccctatgCGCTTCAAGTGCCAGGTGTAATCACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTa                            <  2:192765‑M2/114‑1 (MQ=255)                                                                                                           atcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTa                            <  1:648538‑M2/114‑1 (MQ=255)                                                                                                              agggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAACAAAAACACCCTGTAGAATAGGTAAGCTTATTTTACAg                         >  1:318273‑M2/24‑140 (MQ=255)                                                                                                               gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGa                                        >  2:91157‑M2/23‑124 (MQ=255)                                                                                                               gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGa                                        <  1:91157‑M2/102‑1 (MQ=255)                                                                                                               gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAgg                        >  1:27214‑M2/23‑140 (MQ=255)                                                                                                               gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAgg                        >  2:612239‑M2/23‑140 (MQ=255)                                                                                                               gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAgg                        >  2:410391‑M2/23‑140 (MQ=255)                                                                                                                ggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATACAAACACCGTGTGTAATAGGAACGCTTATTTTACAggg                       >  1:214489‑M2/22‑140 (MQ=255)                                                                                                                      ggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAAc                 >  1:475904‑M2/16‑140 (MQ=255)                                                                                                                        acataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCTTGTATAAAAGGCAAACTTATTTTACAGGGTAAAACCa               >  1:238794‑M2/14‑140 (MQ=255)                                                                                                                          ataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATg             <  1:612239‑M2/129‑1 (MQ=255)                                                                                                                            aaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGcc           >  1:84188‑M2/10‑140 (MQ=255)                                                                                                                                 catgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTa      >  1:50591‑M2/5‑140 (MQ=255)                                                                                                                                   tgCGCTTGAAGAGCCCGGTGTACTATCAACTCCCGATCGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTaca    <  1:176506‑M2/138‑1 (MQ=255)                                                                                                                                     cGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACg  <  1:602018/140‑1 (MQ=255)                                                                                                                                     cGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACg  <  1:110903/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                            GCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACG  >  NC_000913/2859502‑2859773

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