Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,093,991	3461 bp→GCATGCGGATCCTCTAGAGTCGACATATG		[rhaD]–rhaB	[rhaD], rhaA, rhaB

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4093991	4097451	3461	53 [0]	[0] 53	[rhaD]–rhaB	[rhaD],rhaA,rhaB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4093990	0 (0.000)	46 (1.300) +GCATGCGGATCCTCTAGAGTCGACATATG	42/422	0.1	100%	coding (283/825 nt)	rhaD	rhamnulose‑1‑phosphate aldolase
	?	NC_000913	4097452 =	0 (0.000)					intergenic (‑4/‑284)	rhaB/rhaS	rhamnulokinase/DNA‑binding transcriptional activator RhaS

GGTATCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTC  >  NC_000913/4097206‑4097693                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  cccacagttggcgagtgaagaccgcggtgtcgttttcaagtacataggtgagggcgatcaggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGa                                                                                                                                                                                                                                                <  2:119961‑M2/4‑1 (MQ=255)       gttggcgagtgaagaccgcggtgtcgttttcaagtacataggtgagggcgatcaggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCctgct                                                                                                                                                                                                                                         >  2:90005‑M2/241‑251 (MQ=255)            cgagtgaagaccgcggtgtcgttttcaagtacataggtgagggcgatcaggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttaggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAAtt                                                                                                                                                                                                                                    >  1:49454‑M2/236‑251 (MQ=255)                  aagaccgcggtgtcgttttcaagtacataggtgagggcgatcaggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgcgcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacagatgGTGATCCTGCTGAATTTCATTa                                                                                                                                                                                                                              <  2:334829‑M2/22‑1 (MQ=255)                   agaccgcggtgtcgttttcaagtacataggtgagggcgatcaggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTAc                                                                                                                                                                                                                             <  2:88689‑M2/23‑1 (MQ=255)                         cggtgtcgttttcaagtacataggtgagggcgatcaggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAg                                                                                                                                                                                                                       <  1:400017‑M2/29‑1 (MQ=255)                            tgtcgttttcaagtacataggtgagggcgatcaggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagactgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagacgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCt                                                                                                                                                                                                                    >  1:346086‑M2/220‑251 (MQ=255)                                  ttgcaagtacataggtgagggcgatcaggtttgtggcgtgtcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCATTCTaaaaa                                                                                                                                                                                                               <  2:333446‑M2/37‑1 (MQ=255)                                            ataggtgacggcgatcaggttggtggcgtggcagtgcctgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGTAtt                                                                                                                                                                                                    >  2:352010‑M2/204‑251 (MQ=255)                                            ataggcgagggcgctaaggttggtggcgtggcagtgcctgatcacacgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggccggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGCCCAGTCTAAAAAGCGCCTGAAtt                                                                                                                                                                                                    >  2:265314‑M2/204‑251 (MQ=255)                                                             ggttggtggcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTAc                                                                                                                                                                                   <  1:122807‑M2/65‑1 (MQ=255)                                                                     gcgtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGa                                                                                                                                                                           <  2:132283‑M2/73‑1 (MQ=255)                                                                       gtggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgcgcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGTCCTTCTTGTTACTGACAGGaaa                                                                                                                                                                         <  1:319297‑M2/75‑1 (MQ=255)                                                                         ggcagtgcatgatcacccgatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagattcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGaaaa                                                                                                                                                                        <  1:291362‑M2/76‑1 (MQ=255)                                                                                           gatctttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAgg                                                                                                                                                     >  2:10844‑M2/157‑251 (MQ=255)                                                                                                ttgccgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACcagggca                                                                                                                                                  >  2:366022‑M2/152‑247 (MQ=255)                                                                                                  gccgttggtggctttaatgcgctcgcagtgggaaaggaagtcagccggcagtacggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCATGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGaa                                                                                                                                                  >  1:239243‑M2/150‑247 (MQ=255)                                                                                                    cgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATg                                                                                                                                             >  1:51383‑M2/148‑250 (MQ=255)                                                                                                    cgttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGa                                                                                                                                            <  1:274475‑M2/104‑1 (MQ=255)                                                                                                      ttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGt                                                                                                                                                                                      <  1:196545‑M2/62‑1 (MQ=255)                                                                                                      ttggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGt                                                                                                                                                                                      >  2:196545‑M2/146‑207 (MQ=255)                                                                                                       tggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAAc                                                                                                                                          <  2:273123‑M2/106‑1 (MQ=255)                                                                                                        ggtggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGt                                                                                                                                        <  1:293008‑M2/108‑1 (MQ=255)                                                                                                           ggctttaatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGtgat                                                                                                                                     >  1:94884‑M2/141‑251 (MQ=255)                                                                                                                 aatgcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGt                                                                                                                                <  2:372067‑M2/116‑1 (MQ=255)                                                                                                                    gcgctcgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTca                                                                                                                             <  1:208402‑M2/119‑1 (MQ=255)                                                                                                                         cgcagtgggaaaggaagtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAAttt                                                                                                                       <  1:173082‑M2/125‑1 (MQ=255)                                                                                                                                        agtgagccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTgtggtg                                                                                                         >  1:296784‑M2/112‑250 (MQ=255)                                                                                                                                              ccggaagttcggaagtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGgtgatgtgat                                                                                                  <  2:176804‑M2/146‑1 (MQ=255)                                                                                                                                                           agtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttt                                                                                     >  1:179946‑M2/93‑251 (MQ=255)                                                                                                                                                           agtggggacggcttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCAttt                                                                                     <  1:265314‑M2/159‑1 (MQ=255)                                                                                                                                                                      cttcgttggttaacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgaggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGAACTTCTCGTTACTGACAGGCAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGAGATGTGATGCTCACCGCATTTCCTGAAAATTc                                                                          >  1:174942‑M2/82‑251 (MQ=255)                                                                                                                                                                          gttggataacccccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCCCCGCATTTCCTGAAAATTCACGc                                                                      >  1:282010‑M2/78‑251 (MQ=255)                                                                                                                                                                                     cccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATTTTGaa                                                           <  1:327217‑M2/185‑1 (MQ=255)                                                                                                                                                                                     cccaaagaatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGaa                                                           <  1:85875‑M2/185‑1 (MQ=255)                                                                                                                                                                                            aatgtggtagcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCg                                                     >  1:66066‑M2/60‑250 (MQ=255)                                                                                                                                                                                               gtggtagcccgcgccgtcgctgtcgactcatgcggatcctctagagtcgacagatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCGGAATTCGCGACCTGCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATTATGTTCACAATTTGCTGAATTGTTGTGATGTGATGCTCTCCGCAGTTCCTGAAAAGTCACGCTGTAGCTTGAAAAAACGACGt                                                 <  2:368184‑M2/195‑1 (MQ=255)                                                                                                                                                                                                    agcccgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGttttt                                             <  1:309838‑M2/199‑1 (MQ=255)                                                                                                                                                                                                     gcccgcgccgtcgctgtcgacgcatgcggatcctctagagacgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGtttttt                                            >  1:256306‑M2/51‑250 (MQ=255)                                                                                                                                                                                                     gcccgcgccgtcgctgtcgacgcatgcggatccactagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTAGCGACCTACTCGTTAATGACAGGAAAAAGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATCTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGACAAATCGCCGtttttt                                            >  1:345904‑M2/51‑250 (MQ=255)                                                                                                                                                                                                        cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTc                                                                                                                              >  2:116611‑M2/48‑165 (MQ=255)                                                                                                                                                                                                        cgcgccgtcgctgtcgacgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATAATGTTc                                                                                                                              <  1:116611‑M2/118‑1 (MQ=255)                                                                                                                                                                                                                 gctgtcgacgcatgcggatccttaagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAAGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCg                               >  1:331820‑M2/39‑251 (MQ=255)                                                                                                                                                                                                                    gtcgacgcatgcggatcctctagagtcgacatatgGTGAGCCTGCTGAAGTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGATAGGAAAATGGGCCATTGGCATCCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGt                              <  2:331820‑M2/214‑1 (MQ=255)                                                                                                                                                                                                                    ggcgacgcatgcggatcctctagcgtcgacatatgGTGATCCTTCTGAAATTCATTACGACCAGTCTAAAAAGCGCCTGAAGTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCTATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCTCATTTCCTGAAAATTCACGCTGTATCGTGAAAAATCGACGGTTTTTACTTGGTTGTCCGt                              <  2:31626‑M2/214‑1 (MQ=255)                                                                                                                                                                                                                        acgcatgcggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGaaa                         >  1:102954‑M2/32‑250 (MQ=255)                                                                                                                                                                                                                                ggatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGt                >  1:51607‑M2/24‑251 (MQ=255)                                                                                                                                                                                                                                 gatcctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTa               <  2:97704‑M2/229‑1 (MQ=255)                                                                                                                                                                                                                                    ccgctagagtcgacatatgGGGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCGTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGGTCACAATTTGCTGAATTGTGGTGATGTGAGGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTGGCCGTCGAAAAGTTAAGGTAAg             <  2:168627‑M2/231‑1 (MQ=255)                                                                                                                                                                                                                                     ctctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGa            >  2:285595‑M2/19‑250 (MQ=255)                                                                                                                                                                                                                                       ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGaa           <  1:90005‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                       ctagagtcgacatatgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATGGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAAGCGACGTTGTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGaa           <  2:248709‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                              cgacatttgGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTc  >  1:118761‑M2/10‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  GGTATCAATCCCAATGCTATCGATACGAATCCCTTCCTCGCACACCTTGTTTAATCCAAGGCGAATGGCACTTTCAAGGCTATCCACATCCCAGGTGACATAGCCGTTCTGACTATGCAGCCCATTGTTAAAACGATGGATTTCGCGCAGCGTCAGGCTGCGGCATTCACGCTCGTAACGCGCCAGCATCACGCGCCCACTGGATGCGCCGAGATCGACGGCGACACAATTGCGAAAGGTCATAATGTGATCCTGCTGAATTTCATTACGACCAGTCTAAAAAGCGCCTGAATTCGCGACCTTCTCGTTACTGACAGGAAAATGGGCCATTGGCAACCAGGGAAAGATGAACGTGATGATGTTCACAATTTGCTGAATTGTGGTGATGTGATGCTCACCGCATTTCCTGAAAATTCACGCTGTATCTTGAAAAATCGACGTTTTTTACGTGGTTTTCCGTCGAAAATTTAAGGTAAGAACCTGACCTC  >  NC_000913/4097206‑4097693

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