Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_010473	4,296,510	671 bp→GAGCGGCCGCCACCGCGGGAAGTTCCTATACTTTCTAGAGAATAGGAACCTTCGCATGCACGCAGCATATGC	coding (2‑672/672 nt)	nfi →	deoxyribonuclease V

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_010473	4296510	4297180	671	38 [0]	[0] 38	nfi	deoxyribonuclease V

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_010473	= 4296509	0 (0.000)	27 (1.100) +GAGCGGCCGCCACCGCGGGAAGTTCCTATACTTTCTAGAGAATAGGAACCTTCGCATGCACGCAGCATATGC	26/254	0.1	100%	coding (1/672 nt)	nfi	deoxyribonuclease V
	?	NC_010473	4297181 =	0 (0.000)					intergenic (+1/‑42)	nfi/yjaG	deoxyribonuclease V/YjaG family protein

GGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAG  >  NC_010473/4296310‑4296691                                                                                                                                                                                                           |                                                                                                                                                                                       ggATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑CAATC‑ACTGT‑TG‑AAAAATATCACCGCTAAGGAGTGATTAg                                                                                                                                                                                       <  1:24681S26‑M1/200‑2 (MQ=255)  gATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAg                                                                                                                                                                                       <  1:4794S41‑M1/200‑2 (MQ=255)     ccGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCGATGCC‑ATCGACTG‑CTG‑AACAGTATCGCCGCTAAGGAGTGATTAgagc                                                                                                                                                                                    >  1:27946S14‑M1/1‑196 (MQ=255)         cGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggcc                                                                                                                                                                                >  1:22427S94‑M1/1‑192 (MQ=255)                      gccgccTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaa                                                                                                                                                                   >  1:29000S20‑M1/1‑179 (MQ=255)                          ccTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgaagccgcggccgccaccgcgggaa                                                                                                                                                               <  1:1479S11‑M1/200‑26 (MQ=255)                               ccGCATTGAAGAAGAAGTAGCGACTAT‑GTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTCGAGGCGCTG‑C‑AT‑GACTGTCTG‑‑ACATTATCACCACTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactt                                                                                                                                                       >  1:7595S11‑M1/1‑167 (MQ=255)                                      gaagaagaaGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTCT‑‑GTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcggggagttctgatacttttgaga                                                                                                                                                 <  1:18514S55‑M1/200‑40 (MQ=255)                                           agaagTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaa                                                                                                                                              <  1:22336S13‑M1/200‑43 (MQ=255)                                             aagTAGCGACGATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGAT‑agcggccgccaccgcgtagttctgatactttctagagaataggaa                                                                                                                                            <  1:477S12‑M1/200‑46 (MQ=255)                                                               cAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTGAC‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggagttcaggttcccagagaataggaaccttcgcatgcacgcag                                                                                                                           <  1:2485S90‑M1/200‑62 (MQ=255)                                                                           aCGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgctt                                                                                                              <  1:14614S4‑M1/200‑75 (MQ=255)                                                                             ggCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttca                                                                                                            >  1:2590S55‑M1/1‑124 (MQ=255)                                                                             ggCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGC‑‑GTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaatagggaaccttcgcatgcacgcagcatatgcttcag                                                                                                          <  1:18753S1‑M1/200‑79 (MQ=255)                                                                                  aGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaa                                                                                                       <  1:15710S7‑M1/200‑82 (MQ=255)                                                                                         gTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagaaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcgga                                                                                                <  1:2793S16‑M1/200‑89 (MQ=255)                                                                                           cTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaattttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcgga                                                                                              >  1:29883S26‑M1/1‑110 (MQ=255)                                                                                               aaCCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgaggataatt                                                                                          >  1:26951S29‑M1/1‑106 (MQ=255)                                                                                                aCCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataattt                                                                                         >  1:4193S46‑M1/1‑105 (MQ=255)                                                                                                       tCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatacttttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtat                                                                                  >  1:31892S58‑M1/1‑98 (MQ=255)                                                                                                         aCGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatacttttagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttga                                                                                >  1:18012S10‑M1/1‑96 (MQ=255)                                                                                                            gCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatctttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaa                                                                             >  1:20808S17‑M1/1‑93 (MQ=255)                                                                                                                tcatcaGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctgagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactca                                                                         >  1:32217S55‑M1/1‑89 (MQ=255)                                                                                                                   tcaGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaac                                                                      <  1:2096S4‑M1/200‑115 (MQ=255)                                                                                                                         tgtgCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagaatagggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttac                                                                <  1:30126S70‑M1/200‑121 (MQ=255)                                                                                                                           tgCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttaca                                                              >  1:13063S53‑M1/1‑78 (MQ=255)                                                                                                                                        tGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgaaaccaacatgttacaaaacccaattcatct                                                 <  1:551S55‑M1/200‑136 (MQ=255)                                                                                                                                         gCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTGAAACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactgacatgttacaaaacccaattcatc                                                 <  1:11998S30‑M1/200‑136 (MQ=255)                                                                                                                                           tGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctg                                              >  1:30468S105‑M1/1‑62 (MQ=255)                                                                                                                                                 gtTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggagttctgactccagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctgcgtctggagcg                                        <  1:30468S34‑M1/200‑145 (MQ=255)                                                                                                                                                    cGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTCaagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctgcgtctggag                                     >  1:16858S7‑M1/1‑53 (MQ=255)                                                                                                                                                      tGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctagtactctcagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctgcgtctggagcgcc                                   <  1:14569S15‑M1/200‑150 (MQ=255)                                                                                                                                                          ggCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctttttatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctgcgtctggagcg                               >  1:26250S6‑M1/1‑47 (MQ=255)                                                                                                                                                                  c‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctgcgtctggagcgcctgaaagctg                        <  1:16421S6‑M1/200‑161 (MQ=255)                                                                                                                                                                    aTCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttagaaaacccaattcatctgcgtctggagcgcctggaaagctg                       <  1:22175S25‑M1/200‑162 (MQ=255)                                                                                                                                                                            tCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaagcctcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttaaaaaactggagcatgtaccaaacccaattcatctgcgtctggagcgcctggaaagctgggccaca               >  1:14568S2‑M1/1‑31 (MQ=255)                                                                                                                                                                            tCTG‑AACAGTATCACCGCTAAGGAGTGATTAgagcggccgccaccgcgggaagttcctatactttctagagaataggaaccttcgcatgcacgcagcatatgcttcaggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctgcgtctggagcgcctggagctgcagcatgtca               <  1:26255S4‑M1/200‑170 (MQ=255)                                                                                                                                                                                          aCCGCTAAGGAGTGATTAgagcggccgccaccgcggggagttctgatactttctcagagaataggaaccttcgcatgcacgcagcatatgcttcagggtaaactgcggataatttccgtatttgagaactcaacatgttacaaaacccaattcatctgcgtctggagcgcctggaaagctggcagcatgtcactttcatg  <  1:1195S32‑M1/200‑183 (MQ=255)                                                                                                                                                                                                           |                                                                                                                                                                                       GGATCCGTCGATGCTGTACGCGCCGCCTGCCCGCATTGAAGAAGAAGTAGCGACTATACTTGCAGGTTTCGGTCACGGCGAAGGTCATGTCTTTAACCTTGGTCACGGCATTCATCAGGATGTGCCGCCAGAACATGCTGGCGTGTTCGTGGAGGCAGTG‑C‑ATCGACTGTCTG‑AACAGTATCACCGCTAAGGAGTGATTATGGATCTCGCGTCATTACGCGCTCAACAAATCGAACTGGCTTCTTCTGTGATCCGCGAGGATCGACTCGATAAAGATCCACCGGATCTGATCGCCGGAGCCGATGTCGGGTTTGAGCAGGGCGGAGAAGTGACGCGAGCGGCGATGGTGCTGCTGAAATATCCCTCGCTTGAGCTGGTCGAG  >  NC_010473/4296310‑4296691

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 47 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑