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	79 [0]	[0] 79	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)	77 (1.400)	47/194	0.1	100%	coding (1249/1323 nt)	deoA	thymidine phosphorylase
	?	NC_000913	4619487 =	0 (0.000)					coding (1259/1323 nt)	deoA	thymidine phosphorylase

CTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATAT  >  NC_000913/4619349‑4619604                                                                                                                                 |                                                                                                                                cTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaact                                                                                                                       >  2:1037820‑M1/1‑129 (MQ=255)  tGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaactt                                                                                                                      <  1:1027869‑M1/139‑12 (MQ=255)       ccATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccga                                                                                                                 >  1:827946‑M1/1‑123 (MQ=255)          tCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataa                                                                                                              >  2:1269960‑M1/1‑120 (MQ=255)          tCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataa                                                                                                              >  1:420910‑M1/1‑120 (MQ=255)           cGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaa                                                                                                             >  2:1155065‑M1/1‑119 (MQ=255)           cGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaa                                                                                                             >  2:73502‑M1/1‑119 (MQ=255)                    cGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaa                                                                                                    <  2:827946‑M1/139‑30 (MQ=255)                    cGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaag                                                                                                            >  1:1366652‑M1/1‑110 (MQ=255)                    cGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaag                                                                                                            <  2:1366652‑M1/131‑22 (MQ=255)                      tCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaag                                                                                                  <  1:151423‑M1/139‑32 (MQ=255)                        ggCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagca                                                                                                <  1:145438‑M1/139‑34 (MQ=255)                                    atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtct                                                                                    >  1:1101784‑M1/1‑94 (MQ=255)                                    atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtct                                                                                    >  1:1287622‑M1/1‑94 (MQ=255)                                      ggCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagca                                                                                                <  2:297134‑M1/125‑34 (MQ=255)                                      ggCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagca                                                                                                >  1:297134‑M1/1‑92 (MQ=255)                                         gcgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgc                                                                               >  2:447262‑M1/1‑89 (MQ=255)                                            tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgt                                                                            >  1:411478‑M1/1‑86 (MQ=255)                                            tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgata                                                                                                               <  1:101586‑M1/104‑19 (MQ=255)                                            tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgata                                                                                                               >  2:101586‑M1/1‑86 (MQ=255)                                            tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgc                                                                                                                    >  1:374327‑M1/1‑86 (MQ=255)                                            tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgc                                                                                                                    <  2:374327‑M1/99‑14 (MQ=255)                                                ggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaa                                                                                                    <  2:189185‑M1/111‑30 (MQ=255)                                                ggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaa                                                                                                    >  1:189185‑M1/1‑82 (MQ=255)                                                     ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   <  2:1287622‑M1/139‑63 (MQ=255)                                                     ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   <  2:411478‑M1/139‑63 (MQ=255)                                                     ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   <  1:1037820‑M1/139‑63 (MQ=255)                                                      cAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaat                                                                  <  2:420910‑M1/139‑64 (MQ=255)                                                       aGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaata                                                                 <  1:442054‑M1/139‑65 (MQ=255)                                                          tAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                                              >  2:198156‑M1/1‑72 (MQ=255)                                                               ggTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtata                                                         >  2:1259969‑M1/1‑67 (MQ=255)                                                                gTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   <  2:783775‑M1/128‑63 (MQ=255)                                                                gTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   >  1:783775‑M1/1‑66 (MQ=255)                                                                 tCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaat                                                                  <  2:98195‑M1/128‑64 (MQ=255)                                                                 tCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaat                                                                  >  1:98195‑M1/1‑65 (MQ=255)                                                                   aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactga                                                     >  2:1104094‑M1/1‑63 (MQ=255)                                                                   aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   >  2:1018283‑M1/1‑63 (MQ=255)                                                                   aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   <  1:1018283‑M1/125‑63 (MQ=255)                                                                   aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaag                                                                                                  >  2:1372892‑M1/1‑63 (MQ=255)                                                                   aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaag                                                                                                  <  1:1372892‑M1/94‑32 (MQ=255)                                                                   aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgc                                                                                                                    <  2:1233525‑M1/76‑14 (MQ=255)                                                                   aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgc                                                                                                                    >  1:1233525‑M1/1‑63 (MQ=255)                                                                           cTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                                              >  1:158390‑M1/1‑55 (MQ=255)                                                                           cTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                                              <  2:158390‑M1/122‑68 (MQ=255)                                                                           cTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaat                                                                  >  1:616568‑M1/1‑55 (MQ=255)                                                                           cTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaat                                                                  <  2:616568‑M1/118‑64 (MQ=255)                                                                             ggCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtgtcagcgaata                                                                 >  2:1248959‑M1/1‑53 (MQ=255)                                                                             ggCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtgtcagcgaata                                                                 <  1:1248959‑M1/117‑65 (MQ=255)                                                                                    aTCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                                >  1:1322375‑M1/1‑46 (MQ=255)                                                                                    aTCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                                <  2:1322375‑M1/111‑66 (MQ=255)                                                                                     tCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttg                                   <  1:1104094‑M1/139‑95 (MQ=255)                                                                                        aCGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtata                                                         <  2:780796‑M1/114‑73 (MQ=255)                                                                                        aCGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtata                                                         >  1:780796‑M1/1‑42 (MQ=255)                                                                                        aCGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   <  2:1159619‑M1/104‑63 (MQ=255)                                                                                        aCGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                                   >  1:1159619‑M1/1‑42 (MQ=255)                                                                                         cgcgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagc                               <  2:1101784‑M1/139‑99 (MQ=255)                                                                                         cgcgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagc                               <  1:73502‑M1/139‑99 (MQ=255)                                                                                         cgcgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaata                                                                 >  1:502094‑M1/1‑41 (MQ=255)                                                                                         cgcgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaata                                                                 <  2:502094‑M1/105‑65 (MQ=255)                                                                                           cgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcac                             <  1:1259969‑M1/139‑101 (MQ=255)                                                                                            gAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcact                            <  1:447262‑M1/139‑102 (MQ=255)                                                                                            gAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcact                            <  1:1269960‑M1/139‑102 (MQ=255)                                                                                                    aaaaCAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacgg                    <  1:198156‑M1/139‑110 (MQ=255)                                                                                                      aacaacTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagca                              <  2:547094‑M1/127‑100 (MQ=255)                                                                                                      aacaacTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagca                              >  1:547094‑M1/1‑28 (MQ=255)                                                                                                           gtgGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctcatcatttaaatttgaaccac                             <  1:972090‑M1/122‑101 (MQ=255)                                                                                                           cTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcac                             >  2:972090‑M1/1‑23 (MQ=255)                                                                                                              gCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacggag                  >  1:853084‑M1/1‑20 (MQ=255)                                                                                                              gCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacggag                  <  2:853084‑M1/131‑112 (MQ=255)                                                                                                              gCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgaacatttaaatttgaagcactgagtacggag                  <  1:151686‑M1/131‑112 (MQ=255)                                                                                                              gCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgaacatttaaatttgaagcactgagtacggag                  >  2:151686‑M1/1‑20 (MQ=255)                                                                                                                   aaGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcac                             <  2:1088262‑M1/115‑101 (MQ=255)                                                                                                                   aaGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcac                             >  1:1088262‑M1/1‑15 (MQ=255)                                                                                                                    aGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacggagaacatatgaaacgt    >  1:68172‑M1/1‑14 (MQ=255)                                                                                                                      cggcgAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacggagaacatatgaaacgtgc  <  2:812191‑M1/139‑128 (MQ=255)                                                                                                                        gcgAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                                <  2:372765‑M1/75‑66 (MQ=255)                                                                                                                        gcgAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                                >  1:372765‑M1/1‑10 (MQ=255)                                                                                                                             aGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                                              <  2:606994‑M1/72‑68 (MQ=255)                                                                                                                             aGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                                              >  1:606994‑M1/1‑5 (MQ=255)                                                                                                                                 |                                                                                                                                CTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATAT  >  NC_000913/4619349‑4619604

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