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	45 [0]	[0] 46	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)	45 (1.240)	37/264	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

GCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTA  >  NC_000913/2859502‑2859769                                                                                                                                     |                                                                                                                                        catcaagaactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgc                                                                                                                                         >  2:23168‑M2/133‑133 (MQ=255)        aactcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatctaatatcagggaaccggacataaccccatgCGCTTGaa                                                                                                                                  <  2:73578‑M2/8‑1 (MQ=255)           tcgccgttcgcttcttcccctgcaatgattcactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAagag                                                                                                                               <  1:214063‑M2/11‑1 (MQ=255)           tcgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAagag                                                                                                                               >  2:475982‑M2/123‑133 (MQ=255)            cgccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGc                                                                                                                              <  1:632091‑M2/12‑1 (MQ=255)             gccgttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                             >  1:12368‑M2/121‑133 (MQ=255)             gccgttcgattcttcgcctgaagtgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGcc                                                                                                                             <  1:294083‑M2/13‑1 (MQ=255)                gttcgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTgg                                                                                                                          <  1:463399‑M2/16‑1 (MQ=255)                   cgcttcttcccctgaaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatccgggaaccggacaaaacgccatgCGCGTGAAGAGCCTGgtgt                                                                                                                       >  2:195042‑M2/115‑133 (MQ=255)                        cttcccctcaaatgactagctccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGtaataa                                                                                                                  <  1:297344‑M2/24‑1 (MQ=255)                        cttcccctcaaatgactacctcgggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGtaataa                                                                                                                  <  1:245970‑M2/24‑1 (MQ=255)                                 aaatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCg                                                                                                         >  1:512329‑M2/101‑133 (MQ=255)                                  aatgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGa                                                                                                        <  1:23168‑M2/34‑1 (MQ=255)                                    tgattaactccggtatcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATa                                                                                                      <  2:12368‑M2/36‑1 (MQ=255)                                                  atcatgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCg                                                                                        <  1:26465‑M2/50‑1 (MQ=255)                                                      tgtgcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGaa                                                                                    >  1:334231‑M2/80‑133 (MQ=255)                                                         gcgccttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAatat                                                                                 >  1:155281‑M2/77‑133 (MQ=255)                                                             cttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAAc                                                                             <  1:488798‑M2/61‑1 (MQ=255)                                                             cttatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAAc                                                                             >  1:35453‑M2/73‑133 (MQ=255)                                                               tatgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGa                                                                           >  1:425175‑M2/71‑133 (MQ=255)                                                                atgtgattacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGAc                                                                          <  1:50487‑M2/64‑1 (MQ=255)                                                                      ttacaacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGa                                                                    >  1:613391‑M2/64‑133 (MQ=255)                                                                         caacgaaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGa                                                                 >  2:642903‑M2/61‑133 (MQ=255)                                                                             gaaaataaaaaccataaccccccatttaatgtcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataa                                                             <  1:263089‑M2/77‑1 (MQ=255)                                                                              aaaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATTTTAACGTCAACTGCCGaataaa                                                            >  1:229805‑M2/56‑133 (MQ=255)                                                                               aaataaaaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGaataaaa                                                           <  2:155281‑M2/79‑1 (MQ=255)                                                                                     aaaccatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATaaaaa                                                     <  2:613391‑M2/85‑1 (MQ=255)                                                                                       accatcaccccccatctactatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                   <  1:195042‑M2/87‑1 (MQ=255)                                                                                       accatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                   <  1:641386‑M2/87‑1 (MQ=255)                                                                                       accatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                   <  1:358558‑M2/87‑1 (MQ=255)                                                                                       accatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                   <  1:437924‑M2/87‑1 (MQ=255)                                                                                       accatcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                   <  1:398871‑M2/87‑1 (MQ=255)                                                                                       accatcacaccccatttaataacaaggacccgcacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAaca                                                   <  2:201727‑M2/87‑1 (MQ=255)                                                                                         catcacaccccatttaatttcacggaaccgtacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                 <  1:227595‑M2/89‑1 (MQ=255)                                                                                         catcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                 <  2:334231‑M2/89‑1 (MQ=255)                                                                                         catcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAcc                                                 <  2:425175‑M2/89‑1 (MQ=255)                                                                                          atcacaccccatttaatatcagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACAccc                                                <  2:512329‑M2/90‑1 (MQ=255)                                                                                                             cagggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTa                             >  2:254751‑M2/25‑133 (MQ=255)                                                                                                               gggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCAATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAtt                           >  1:347291‑M2/23‑133 (MQ=255)                                                                                                                ggaaccggacataaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTAttt                          >  1:526342‑M2/22‑133 (MQ=255)                                                                                                                  aaccggacattaccccatgCGCTTTAAGAGCCGGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTa                        <  2:229805‑M2/114‑1 (MQ=255)                                                                                                                            aaccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTaaaa              <  1:124444‑M2/124‑1 (MQ=255)                                                                                                                              ccccatgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAAcc            >  1:350564‑M2/8‑133 (MQ=255)                                                                                                                                 catgCGCTTGAAGAGCCTGGTGTAATAACAATTCCCAATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATg         <  2:347291‑M2/129‑1 (MQ=255)                                                                                                                                   tgCGCTTGAAGAGCC‑‑‑TGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCatc    >  2:405063‑M2/3‑132 (MQ=255)                                                                                                                                     cGCTTGAAGAGCC‑‑‑TGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTa  <  1:405063/133‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        GCCGCTACGCAAGTGGATGGTACGCAAATGTCTTTGCTGTCAGTACCAGAAGAAACTTCGCCTGCGGTCGAAGCTCTGGAAAATCTTGATCCGGATTCACTCACCCCGCGTCAGGCGCTGGAGTGGATTTATCGCTTGAAGAGCCTGGTGTAATAACAATTCCCGATAGTCTTTTGCTATCGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTA  >  NC_000913/2859502‑2859769

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