Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,836,957	Δ8,771 bp		[setC]–[adeD]	[setC], yicL, yicU, nlpA, yicS, nepI, yicN, adeQ, [adeD]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3836957	3845727	8771	44 [0]	[0] 43	[setC]–[adeD]	[setC],yicL,yicU,nlpA,yicS,nepI,yicN,adeQ,[adeD]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3836956	0 (0.000)	43 (0.440)	38/492	0.9	100%	coding (4/1185 nt)	setC	putative arabinose exporter
	?	NC_000913	3845728 =	0 (0.000)					coding (1765/1767 nt)	adeD	adenine deaminase

ACGTAACGGCCAGGTACAAAGTCATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTG  >  NC_000913/3845478‑3845968                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 atcaaatcaacttgttttatttttactatcaaaaatagctcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgct                                                                                                                                                                                                                                                  >  1:216285‑M2/251‑251 (MQ=255)  tcaaatcaacttgttttatttttactatcaaaaatagctcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgct                                                                                                                                                                                                                                                  >  2:567735‑M2/250‑250 (MQ=255)        caacttggtttctttttactatcaaaaatagctcattgcattactataattaatatctaaattcggaaattaattcttaacatgcttttactaataatctcaattgctggccatataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaacactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATaaa                                                                                                                                                                                                                                            >  2:468073‑M2/244‑250 (MQ=255)            ttgttttatttttactatccaaaatagctcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATaaaaaa                                                                                                                                                                                                                                         >  1:550428‑M2/240‑249 (MQ=255)              gttttatttttactatcaaaaatagctcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGc                                                                                                                                                                                                                                      <  1:700785‑M2/13‑1 (MQ=255)                     ttttactatcaaaaatagctcattgccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGt                                                                                                                                                                                                                               <  1:554762‑M2/20‑1 (MQ=255)                        tactatcaaaaatagctcattgccttactataattaatatataaattcggtaatgaattattaacatgcgtttaataataatctcaattgctggccctataatattgcactagaattgctttctggatgtatcagcgatactcaaaaattcttaatacaatactcattcgactggaacttatttgtaactcagttatattttttcgcccggtgattcagaataatgcTAATAAAAAAAGCCCGGAGTCATg                                                                                                                                                                                                                           >  2:478659‑M2/228‑251 (MQ=255)                                              ccttactataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTAc                                                                                                                                                                                                     >  1:459550‑M2/206‑251 (MQ=255)                                                  actataattaatatataaattcggtaattaattcttaacatgcgtttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGtt                                                                                                                                                                                                  <  2:550428‑M2/49‑1 (MQ=255)                                                   ctataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGttt                                                                                                                                                                                                 >  1:672975‑M2/201‑250 (MQ=255)                                                   ctataattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGttt                                                                                                                                                                                                 >  1:106470‑M2/201‑250 (MQ=255)                                                     aaaattaatatataaattcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCTATAGTCCCCAGTTACGttt                                                                                                                                                                                                 <  2:526318‑M2/50‑1 (MQ=255)                                                            ataaataaattcggtaattaattcttaccattcttttactaattatctcaattgctggccctataatattgcgctagaattgctttctggttgtatcagcgatacttaaaaattcttaatacaatactcattcgacttgtacttattttgaactcatttatatttttgcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTGTATGCCACTGt                                                                                                                                                                                       <  2:203140‑M2/60‑1 (MQ=255)                                                                   aattcggtaattaattcttaacatgcttttactaataatctcaattgctggacctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattagactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgc                                                                                                                                                                                >  1:115516‑M2/185‑251 (MQ=255)                                                                    attcggtaattaattcttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacactactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAACGTNCCCAGTTACGTTTATGCCACTGTCAActgct                                                                                                                                                                               >  2:366980‑M2/184‑251 (MQ=255)                                                                              taattcttaacatgcttttactaataatctcaattgctggccctataatcttgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatcatcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagactgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTACCCAGTTACGTTTATGCCACTGCCAACTGCTGAAttttttt                                                                                                                                                                     >  2:609135‑M2/174‑251 (MQ=255)                                                                                    ttaacatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg                                                                                                                                                               >  2:428331‑M2/168‑251 (MQ=255)                                                                                       acatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggcg                                                                                                                                                            >  2:677760‑M2/165‑251 (MQ=255)                                                                                       acatgcttttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggcg                                                                                                                                                            >  1:507105‑M2/165‑251 (MQ=255)                                                                                              tttactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAttt                                                                                                                                                       <  1:428331‑M2/92‑1 (MQ=255)                                                                                                tactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCg                                                                                                                                                   >  1:168752‑M2/156‑251 (MQ=255)                                                                                                tactaataatctcaattgctggccctataatattgcgctagcattgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCg                                                                                                                                                   >  1:668824‑M2/156‑251 (MQ=255)                                                                                                          ctcaattgctgggcctataatattgcgctagcattgctttctggttgtatcagcgctaatcaaaaattcttaatacaataatcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGAAAAAGTCAACAGTTACGTTTATGCCACTGTCAACTGCTGACTTTTTTTCTCGCGGCGGATTTTGCGGTCTTCCAtt                                                                                                                                         >  2:331576‑M2/146‑251 (MQ=255)                                                                                                                          ataatattgcgctagcatttctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcatttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGTCACTATCGc                                                                                                                           <  2:512955‑M2/120‑1 (MQ=255)                                                                                                                                          attgctttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAAcc                                                                                                         >  1:291465‑M2/114‑251 (MQ=255)                                                                                                                                               tttctggttgtatcagcgatactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTAtcgccatcaggcagggggagcatttcc                                                                                                    >  1:699767‑M2/109‑238 (MQ=255)                                                                                                                                                  ctggttgtatcagcgatactcaaaacttattaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcgg                                                                                                 >  1:587897‑M2/106‑251 (MQ=255)                                                                                                                                                                 atactcaaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg                                                                                  >  1:314149‑M2/91‑251 (MQ=255)                                                                                                                                                                       aaaaattcttaatacaatactcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAgg                                                                            >  1:521061‑M2/85‑251 (MQ=255)                                                                                                                                                                                       atactctttcgactgggacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACCACAAATCGCGGCGATAACCAGCGCGGCGAAGGTGCCTGCGCAGCcggg                                                            >  2:140818‑M2/69‑250 (MQ=255)                                                                                                                                                                                           tcattcgactggtacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAgg                                                        >  1:409452‑M2/65‑251 (MQ=255)                                                                                                                                                                                            cattcgactggcacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGCTTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATCCCACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAACCGGTAAGGc                                                       >  2:430157‑M2/64‑251 (MQ=255)                                                                                                                                                                                                      gtacttatttgtaactcagttatattttttcgcccggtgattcagcagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCgg                                             >  1:231027‑M2/54‑251 (MQ=255)                                                                                                                                                                                                       tacttatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGc                                            >  1:504022‑M2/53‑251 (MQ=255)                                                                                                                                                                                                           tatttgtaactcagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGAAATCCAGCGCGGCGAAAGGGCCTGCCCATCCGGTAAGGCCGAATACCGGCGTcc                                        >  1:20879‑M2/49‑251 (MQ=255)                                                                                                                                                                                                                      cagttatattttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTCCCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGGgcctgcccagccggtaaggccgaaaccctgcgtaccctccgccatc                             >  2:658126‑M2/38‑229 (MQ=255)                                                                                                                                                                                                                             attttttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTccc                       >  1:68325‑M2/31‑250 (MQ=255)                                                                                                                                                                                                                                 tttcgcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGATCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCCTTCccagccc                  >  2:254000‑M2/27‑248 (MQ=255)                                                                                                                                                                                                                                     gcccggtgattcagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCCGCCGGTAAGGCCGAATACCGGCGGCCCATCGGCAATCCTTCCCAGACCTAAc              >  2:177456‑M2/23‑251 (MQ=255)                                                                                                                                                                                                                                        cggtgattcagaggaatgcTAATAAAAACAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTAACGCCACTATCGCCCTCAGGCAGATACAACCAATCGCGGCGAACACCAGCGCGGCGAAGGTGCCTGCCCCGCCGGGAAGGCCGAATTCCGGCGGCCCATCGGCAATCATTTCCCGACCTCACTTg           >  2:142483‑M2/20‑251 (MQ=255)                                                                                                                                                                                                                                        cggtgattcagaagaatgaTAATAAAAAAAGCCCGGCGACATGCCGGGAAAAAGACACCAGTAACGATTATGCCACTGACAACTGCTGAATTTTTTTCTCGCGGCGGCTTTTGCGTTCTTCCATTACCGCCACTATCGCAATCAAGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCCAAGGTGCATGCCCAGCCGGTACAGCCGCATACCGGCGTCCCCTCGGCAATCATTCCCAGAACTAACTTg           >  2:298516‑M2/20‑251 (MQ=255)                                                                                                                                                                                                                                                cagaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCCGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAgcc   >  1:555789‑M2/12‑250 (MQ=255)                                                                                                                                                                                                                                                 agaagaatgcTAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGACTGCCCAGCCGGTCAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCCGGCCTAACTTGGCAAAGCTg  >  1:692175‑M2/11‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 ACGTAACGGCCAGGTACAAAGTCATCTGCCGTTACCCATTGCCGGGCTGATGAGCACCGACACGGCGCAGTCGCTGGCGGAACAAATTGACGCCTTGAAAGCCGCCGCCCGTGAATGCGGTCCGTTACCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCTGCCCAGCCGGTAAGGCCGAATACCGGCGTCCCATCGGCAATCATTCCCAGACCTAACTTGGCAAAGCTG  >  NC_000913/3845478‑3845968

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