Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,619,478	Δ9 bp	coding (1250‑1258/1323 nt)	deoA →	thymidine phosphorylase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4619478	4619486	9	52 [0]	[0] 52	deoA	thymidine phosphorylase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4619477	0 (0.000)	52 (1.040)	38/190	0.2	100%	coding (1249/1323 nt)	deoA	thymidine phosphorylase
	?	NC_000913	4619487 =	0 (0.000)					coding (1259/1323 nt)	deoA	thymidine phosphorylase

GATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAG  >  NC_000913/4619360‑4619597                                                                                                                                |                                                                                                               gacaccatcgattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTg                                                                                                     >  1:307208‑M2/128‑139 (MQ=255)        tcgattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATaa                                                                                              >  1:1076953‑M2/121‑139 (MQ=255)          gattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAg                                                                                            <  2:307208‑M2/21‑1 (MQ=255)                  cgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGaa                                                                                    >  1:652259‑M2/111‑139 (MQ=255)                     cggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGc                                                                                 >  2:1179516‑M2/108‑139 (MQ=255)                       gctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGcac                                                                               <  1:730123‑M2/34‑1 (MQ=255)                          ttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACAcc                                                                            <  2:628937‑M2/37‑1 (MQ=255)                              tgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGTAAAGCACACCAACt                                                                        >  2:182646‑M2/99‑139 (MQ=255)                                  atggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGc                                                                                                    <  2:225535‑M2/13‑1 (MQ=255)                                  atggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGc                                                                                                    >  1:225535‑M2/95‑107 (MQ=255)                                  atggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCt                                                                    <  1:741722‑M2/45‑1 (MQ=255)                                    ggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTAt                                                                  <  2:361436‑M2/47‑1 (MQ=255)                                     gcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATc                                                                 >  1:1116315‑M2/92‑139 (MQ=255)                                        cgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGcc                                                              >  2:1077823‑M2/89‑139 (MQ=255)                                        cgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGcc                                                              >  1:58624‑M2/89‑139 (MQ=255)                                         gtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCg                                                             <  1:1077823‑M2/52‑1 (MQ=255)                                          tctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                          >  2:188330‑M2/87‑125 (MQ=255)                                          tctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                          <  1:188330‑M2/39‑1 (MQ=255)                                          tctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGt                                                            >  1:368623‑M2/87‑139 (MQ=255)                                                 gaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCg                                                     <  2:652259‑M2/60‑1 (MQ=255)                                                   ccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCGGCGaa                                                   <  1:925102‑M2/62‑1 (MQ=255)                                                   ccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   >  1:32033‑M2/78‑139 (MQ=255)                                                   ccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   <  2:58624‑M2/62‑1 (MQ=255)                                                      ggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                <  2:32033‑M2/65‑1 (MQ=255)                                                      ggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                <  1:1311986‑M2/65‑1 (MQ=255)                                                       gtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   <  2:1255461‑M2/62‑1 (MQ=255)                                                       gtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   >  1:1255461‑M2/74‑135 (MQ=255)                                                        tagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                              <  1:1179516‑M2/67‑1 (MQ=255)                                                        tagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                              <  2:134290‑M2/67‑1 (MQ=255)                                                          gacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGt                                            <  2:349523‑M2/69‑1 (MQ=255)                                                             ggtcagcgtccgctggcggttatccacgcgaaagactaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGtata                                         >  1:719349‑M2/68‑139 (MQ=255)                                                               tcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa                                                 <  1:436820‑M2/64‑1 (MQ=255)                                                               tcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa                                                 >  2:436820‑M2/66‑129 (MQ=255)                                                                 agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActga                                     >  1:1203584‑M2/64‑139 (MQ=255)                                                                     tccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   <  2:11905‑M2/62‑1 (MQ=255)                                                                     tccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   >  1:11905‑M2/60‑121 (MQ=255)                                                                       cgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctgat                               <  2:1241378‑M2/82‑1 (MQ=255)                                                                         ctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCa                             >  1:431573‑M2/56‑139 (MQ=255)                                                                                       cgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGc               <  2:1116315‑M2/98‑1 (MQ=255)                                                                                         cgaaagacgaaaaccactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCTCCGTATCAGCGAATAACGGTATACTGTTCTGATCATTTAAATTCGAAGCAc             <  1:823454‑M2/100‑1 (MQ=255)                                                                                          gaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAg                                                                                            <  1:710398‑M2/21‑1 (MQ=255)                                                                                          gaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAg                                                                                            >  2:710398‑M2/39‑59 (MQ=255)                                                                                          gaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACt            >  2:1192280‑M2/39‑139 (MQ=255)                                                                                             agacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAg         <  2:1076953‑M2/104‑1 (MQ=255)                                                                                                  aaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGcac                                                                               <  2:539771‑M2/34‑1 (MQ=255)                                                                                                  aaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGcac                                                                               >  1:539771‑M2/31‑64 (MQ=255)                                                                                                  aaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                              <  1:770636‑M2/67‑1 (MQ=255)                                                                                                  aaaacaactggcaggaaacggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTATCTATCGCCGTATCAGCGAATAACg                                              >  2:770636‑M2/31‑97 (MQ=255)                                                                                                    aacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag  <  2:368623‑M2/111‑1 (MQ=255)                                                                                                    aacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag  <  2:27239‑M2/111‑1 (MQ=255)                                                                                                             caggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   >  1:718755‑M2/20‑81 (MQ=255)                                                                                                             caggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                   <  2:718755‑M2/62‑1 (MQ=255)                                                                                                                                |                                                                                                               GATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAG  >  NC_000913/4619360‑4619597

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