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

CATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGAGTGCAATAGAAAATTTCGACGCCCATACGCCCATGATGCAGCAGTATCTCAGGCTGAAAGCCCAGCATCCCGAGATCCTGCTGTTTTACCGGATGGGTGATTTTTATGAACTGTTTTATGACGACGCAAA  >  NC_000913/2856964‑2857226                                                                                                                                    |                                                                                                                                    cATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGc                                                                                                                                    >  2:23168‑M1/1‑132 (MQ=255)        aaCTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATCTAATATCAGGGAACCGGACATAACCCCATGcgcttgaa                                                                                                                             <  2:73578‑M1/133‑9 (MQ=255)           tCGCCGTTCGCTTCTTCCCCTGCAATGATTCACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagag                                                                                                                          <  1:214063‑M1/133‑12 (MQ=255)           tCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagag                                                                                                                          >  2:475982‑M1/1‑122 (MQ=255)            cGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagc                                                                                                                         <  1:632091‑M1/133‑13 (MQ=255)             gCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcc                                                                                                                        >  1:12368‑M1/1‑120 (MQ=255)             gCCGTTCGATTCTTCGCCTGAAGTGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcc                                                                                                                        <  1:294083‑M1/133‑14 (MQ=255)                gTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctgg                                                                                                                     <  1:463399‑M1/133‑17 (MQ=255)                   cGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCCGGGAACCGGACAAAACGCCATGcgcgtgaagagcctggtgt                                                                                                                  >  2:195042‑M1/1‑114 (MQ=255)                        cttcCCCTCAAATGACTAGCTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataa                                                                                                             <  1:297344‑M1/133‑25 (MQ=255)                        cttcCCCTCAAATGACTACCTCGGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataa                                                                                                             <  1:245970‑M1/133‑25 (MQ=255)                                 aaaTGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccg                                                                                                    >  1:512329‑M1/1‑100 (MQ=255)                                  aaTGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccga                                                                                                   <  1:23168‑M1/133‑35 (MQ=255)                                    tGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgata                                                                                                 <  2:12368‑M1/133‑37 (MQ=255)                                                  aTCATGTGCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcg                                                                                   <  1:26465‑M1/133‑51 (MQ=255)                                                      tgtgCGCCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaa                                                                               >  1:334231‑M1/1‑79 (MQ=255)                                                         gcgcCTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatat                                                                            >  1:155281‑M1/1‑76 (MQ=255)                                                             cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaac                                                                        <  1:488798‑M1/133‑62 (MQ=255)                                                             cTTATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaac                                                                        >  1:35453‑M1/1‑72 (MQ=255)                                                               tATGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacga                                                                      >  1:425175‑M1/1‑70 (MQ=255)                                                                aTGTGATTACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgac                                                                     <  1:50487‑M1/133‑65 (MQ=255)                                                                      ttACAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactga                                                               >  1:613391‑M1/1‑63 (MQ=255)                                                                         cAACGAAAATAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacga                                                            >  2:642903‑M1/1‑60 (MQ=255)                                                                             gAAAATAAAAACCATAACCCCCCATTTAATGTCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataa                                                        <  1:263089‑M1/133‑78 (MQ=255)                                                                              aaaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaattttaacgtcaactgccgaataaa                                                       >  1:229805‑M1/1‑55 (MQ=255)                                                                               aaaTAAAAACCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaa                                                      <  2:155281‑M1/133‑80 (MQ=255)                                                                                     aaaCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaa                                                <  2:613391‑M1/133‑86 (MQ=255)                                                                                       aCCATCACCCCCCATCTACTATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  1:195042‑M1/133‑88 (MQ=255)                                                                                       aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  1:641386‑M1/133‑88 (MQ=255)                                                                                       aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  1:358558‑M1/133‑88 (MQ=255)                                                                                       aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  1:398871‑M1/133‑88 (MQ=255)                                                                                       aCCATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  1:437924‑M1/133‑88 (MQ=255)                                                                                       aCCATCACACCCCATTTAATAACAAGGACCCGCACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaaca                                              <  2:201727‑M1/133‑88 (MQ=255)                                                                                         cATCACACCCCATTTAATTTCACGGAACCGTACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                            <  1:227595‑M1/133‑90 (MQ=255)                                                                                         cATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                            <  2:334231‑M1/133‑90 (MQ=255)                                                                                         cATCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacacc                                            <  2:425175‑M1/133‑90 (MQ=255)                                                                                          aTCACACCCCATTTAATATCAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccc                                           <  2:512329‑M1/133‑91 (MQ=255)                                                                                                             cAGGGAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagctta                        >  2:254751‑M1/1‑24 (MQ=255)                                                                                                               gggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccaatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatt                      >  1:347291‑M1/1‑22 (MQ=255)                                                                                                                ggAACCGGACATAACCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattt                     >  1:526342‑M1/1‑21 (MQ=255)                                                                                                                  aaCCGGACATTACCCCATGcgctttaagagccgggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttatttta                   <  2:229805‑M1/133‑115 (MQ=255)                                                                                                                            aaCCCCATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaa         <  1:124444‑M1/133‑125 (MQ=255)                                                                                                                              ccccATGcgcttgaagagcctggtgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaacc       >  1:350564‑M1/1‑7 (MQ=255)                                                                                                                                 cATGcgcttgaagagcctggtgtaataacaattcccaatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatg    <  2:347291‑M1/133‑130 (MQ=255)                                                                                                                                   tGcgcttgaagagcctgtaataacaattcccgatagtcttttgctatcgggaatattaacgacaactgacgaataaaataaaaacaccctgtataataggaaagcttattttacagggtaaaaccatgccatc  >  2:405063‑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 < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑