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] 57	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.090)	40/276	0.2	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

CGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACG  >  NC_000913/2859501‑2859773                                                                                                                                      |                                                                                                                                            gcatcaagaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGa                                                                                                                                       >  2:274339‑M2/134‑140 (MQ=255)  catcaagaactcgccgttcgcttcttcccctgaaatgattagcccccgtatcatgtccgccttatgtcattacaacgaaaataaaaaccaacacaccccatttaatatcagggatccgcacataacaccatgCGCCTGaa                                                                                                                                      <  2:781147‑M2/8‑1 (MQ=255)       agaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                                 <  1:720186‑M2/13‑1 (MQ=255)       agaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                                 <  1:294630‑M2/13‑1 (MQ=255)       agaactcgccgttcgcttcttcccctgaaatgattaacaccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                                 <  2:487217‑M2/13‑1 (MQ=255)        gaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcgcaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCt                                                                                                                                >  1:609092‑M2/127‑140 (MQ=255)          actcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTgg                                                                                                                              <  2:342263‑M2/16‑1 (MQ=255)              gccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGta                                                                                                                          >  1:378200‑M2/121‑140 (MQ=255)                  ttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGtaataa                                                                                                                      <  2:498410‑M2/24‑1 (MQ=255)                   tcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAAc                                                                                                                     >  1:646105‑M2/116‑140 (MQ=255)                            cccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa                                                                                                            <  2:378200‑M2/34‑1 (MQ=255)                            cccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa                                                                                                            <  2:110101‑M2/34‑1 (MQ=255)                               ctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccacgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAg                                                                                                         <  2:831273‑M2/37‑1 (MQ=255)                                  aaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCt                                                                                                      >  2:807460‑M2/101‑140 (MQ=255)                                        ttaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa                                                                                                            <  2:309280‑M2/34‑1 (MQ=255)                                        ttaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa                                                                                                            >  1:309280‑M2/95‑128 (MQ=255)                                        ttaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCt                                                                                                <  1:648734‑M2/46‑1 (MQ=255)                                           actccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATc                                                                                             <  1:274339‑M2/49‑1 (MQ=255)                                               cggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGa                                                                                         <  1:431695‑M2/53‑1 (MQ=255)                                               cggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGa                                                                                         <  1:618892‑M2/53‑1 (MQ=255)                                               ccggtatcatgtgcgccttatgtgatttcaacgaaaataaaaccaacacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGa                                                                                         <  2:857809‑M2/53‑1 (MQ=255)                                                       tgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAAc                                                                                 >  1:141995‑M2/80‑140 (MQ=255)                                                                     gattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaat                                                                   <  2:478851‑M2/75‑1 (MQ=255)                                                                     gattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaat                                                                   >  1:85078‑M2/66‑140 (MQ=255)                                                                     gattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaat                                                                   <  1:696851‑M2/75‑1 (MQ=255)                                                                     gattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaat                                                                   <  1:403735‑M2/75‑1 (MQ=255)                                                                     gattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaat                                                                   >  2:216047‑M2/66‑140 (MQ=255)                                                                       ttacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataa                                                                 <  2:340480‑M2/77‑1 (MQ=255)                                                                       ttacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataa                                                                 <  2:324853‑M2/77‑1 (MQ=255)                                                                        tacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaa                                                                >  2:94342‑M2/63‑140 (MQ=255)                                                                         acaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaaa                                                               <  1:216047‑M2/79‑1 (MQ=255)                                                                          caacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaaat                                                              >  1:859463‑M2/61‑140 (MQ=255)                                                                            acgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaaataa                                                            >  1:169877‑M2/59‑140 (MQ=255)                                                                               aaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATaaaaa                                                         <  2:839629‑M2/85‑1 (MQ=255)                                                                                 aataaaaaccatcgcaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  2:609092‑M2/87‑1 (MQ=255)                                                                                 aataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                       <  2:141995‑M2/87‑1 (MQ=255)                                                                                  ataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAacac                                                      <  2:169877‑M2/88‑1 (MQ=255)                                                                                  ataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAacac                                                      <  1:94342‑M2/88‑1 (MQ=255)                                                                                   taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  1:140434‑M2/89‑1 (MQ=255)                                                                                   taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  2:84364‑M2/89‑1 (MQ=255)                                                                                   taaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                     <  1:825451‑M2/89‑1 (MQ=255)                                                                                        accatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGta                                                <  1:807460‑M2/94‑1 (MQ=255)                                                                                                    catttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGc                                    <  1:825549‑M2/106‑1 (MQ=255)                                                                                                      tttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCtt                                  <  2:844142‑M2/108‑1 (MQ=255)                                                                                                           tatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCTGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAtttt                             >  2:400113‑M2/28‑140 (MQ=255)                                                                                                           tatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAtttt                             >  1:423427‑M2/28‑140 (MQ=255)                                                                                                            atcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATTGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTa                            <  1:757468‑M2/114‑1 (MQ=255)                                                                                                            atcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGaaa                                      >  1:491391‑M2/27‑130 (MQ=255)                                                                                                            atcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGaaa                                      <  2:491391‑M2/104‑1 (MQ=255)                                                                                                            atcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTa                            <  2:46639‑M2/114‑1 (MQ=255)                                                                                                             tcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTAc                           <  2:423427‑M2/115‑1 (MQ=255)                                                                                                                gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAgg                        >  2:609543‑M2/23‑140 (MQ=255)                                                                                                                 ggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAggg                       >  1:371073‑M2/22‑140 (MQ=255)                                                                                                                  gaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGt                      <  2:85078‑M2/120‑1 (MQ=255)                                                                                                                                   atgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTac     >  2:303969‑M2/4‑140 (MQ=255)                                                                                                                                    tgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTaca    >  2:177163‑M2/3‑140 (MQ=255)                                                                                                                                      cGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACg  <  2:646105/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                            CGCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACG  >  NC_000913/2859501‑2859773

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