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

GGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCAGTAAAGGCATGGGCTGGCTGAGCGGGATATAGCGCGGTTGTTGGTGGAAATTGTCGTGATATGGTGCGATATCGGCGTCATCCAGGCGTAGCGTCAGGTTGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTC  >  NC_000913/4093744‑4094238                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGGTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCGTCGTTGGTTAACCCCCAAAGAATGTGGGAGCCCGCGCCGTCGCTGTCGACgcat                                                                                                                                                                                                                                                      <  2:342299‑M1/251‑5 (MQ=255) ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTGTTCAAGTACATATGTGAGGGCGATCAGGTTGGTGGCGGGGCAGTGCATGATCTCCCGATCTGTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCGTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCTCCGTCGCTGTCGACgcat                                                                                                                                                                                                                                                      <  2:169563‑M1/251‑5 (MQ=255)   aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgc                                                                                                                                                                                                                                                    >  2:85875‑M1/1‑245 (MQ=255)           caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctc                                                                                                                                                                                                                                            >  1:97704‑M1/1‑237 (MQ=255)            aGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACCTAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctct                                                                                                                                                                                                                                           >  2:274475‑M1/1‑236 (MQ=255)              acacTCTGTGCGGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTTCATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAATTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctag                                                                                                                                                                                                                                         <  1:349523‑M1/251‑18 (MQ=255)                     gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTGTCAGGCACATAGGCGAGGGCGATCAGGTTGGCGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCGCGCAGTGGGAAAGGAAGAGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgac                                                                                                                                                                                                                                  >  2:257973‑M1/1‑227 (MQ=255)                               cccACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtga                                                                                                                                                                                                                         <  2:119961‑M1/250‑34 (MQ=255)                                     gTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgct                                                                                                                                                                                                                  >  2:90005‑M1/1‑211 (MQ=255)                                          cGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTAGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatt                                                                                                                                                                                                             >  1:49454‑M1/1‑206 (MQ=255)                                                aaGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCGCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacagatggtgatcctgctgaatttcatta                                                                                                                                                                                                       <  2:334829‑M1/251‑52 (MQ=255)                                                 aGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattac                                                                                                                                                                                                      <  2:88689‑M1/251‑53 (MQ=255)                                                       cgGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccag                                                                                                                                                                                                <  1:400017‑M1/251‑59 (MQ=255)                                                          tgtCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGACTGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagacgacatatggtgatcctgctgaatttcattacgaccagtct                                                                                                                                                                                             >  1:346086‑M1/1‑190 (MQ=255)                                                                ttgcaagTACATAGGTGAGGGCGATCAGGTTTGTGGCGTGTCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccattctaaaaa                                                                                                                                                                                        <  2:333446‑M1/247‑67 (MQ=255)                                                                          aTAGGTGACGGCGATCAGGTTGGTGGCGTGGCAGTGCCTGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgtatt                                                                                                                                                                             >  2:352010‑M1/1‑174 (MQ=255)                                                                          aTAGGCGAGGGCGCTAAGGTTGGTGGCGTGGCAGTGCCTGATCACACGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGCCGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgcccagtctaaaaagcgcctgaatt                                                                                                                                                                             >  2:265314‑M1/1‑174 (MQ=255)                                                                                           ggTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttac                                                                                                                                                            <  1:122807‑M1/251‑95 (MQ=255)                                                                                                   gcgtggcAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacagga                                                                                                                                                    <  2:132283‑M1/251‑103 (MQ=255)                                                                                                     gtggcAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCGCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgtccttcttgttactgacaggaaa                                                                                                                                                  <  1:319297‑M1/251‑105 (MQ=255)                                                                                                       ggcAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagattcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaa                                                                                                                                                 <  1:291362‑M1/250‑106 (MQ=255)                                                                                                                         gATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagg                                                                                                                              >  2:10844‑M1/1‑127 (MQ=255)                                                                                                                              ttGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggca                                                                                                                           >  2:366022‑M1/1‑122 (MQ=255)                                                                                                                                gCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTCAGCCGGCAGTACGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcatgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaa                                                                                                                           >  1:239243‑M1/1‑120 (MQ=255)                                                                                                                                  cGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatga                                                                                                                     <  1:274475‑M1/251‑134 (MQ=255)                                                                                                                                  cGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatg                                                                                                                      >  1:51383‑M1/1‑118 (MQ=255)                                                                                                                                    ttGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgt                                                                                                                                                               >  2:196545‑M1/1‑116 (MQ=255)                                                                                                                                    ttGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgt                                                                                                                                                               <  1:196545‑M1/207‑92 (MQ=255)                                                                                                                                     tggtggCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaac                                                                                                                   <  2:273123‑M1/250‑136 (MQ=255)                                                                                                                                      ggtggCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgt                                                                                                                 <  1:293008‑M1/251‑138 (MQ=255)                                                                                                                                         ggCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgat                                                                                                              >  1:94884‑M1/1‑111 (MQ=255)                                                                                                                                               aaTGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgt                                                                                                         <  2:372067‑M1/250‑146 (MQ=255)                                                                                                                                                  gcgcTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttca                                                                                                      <  1:208402‑M1/250‑149 (MQ=255)                                                                                                                                                       cGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaattt                                                                                                <  1:173082‑M1/251‑155 (MQ=255)                                                                                                                                                                      aGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtg                                                                                  >  1:296784‑M1/1‑82 (MQ=255)                                                                                                                                                                            ccGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgat                                                                           <  2:176804‑M1/251‑176 (MQ=255)                                                                                                                                                                                         aGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcattt                                                              <  1:265314‑M1/251‑189 (MQ=255)                                                                                                                                                                                         aGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcattt                                                              >  1:179946‑M1/1‑63 (MQ=255)                                                                                                                                                                                                    cTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgaggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaacttctcgttactgacaggcaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggagatgtgatgctcaccgcatttcctgaaaattc                                                   >  1:174942‑M1/1‑52 (MQ=255)                                                                                                                                                                                                        gTTGGATAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctccccgcatttcctgaaaattcacgc                                               >  1:282010‑M1/1‑48 (MQ=255)                                                                                                                                                                                                                   cccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtattttgaa                                    <  1:327217‑M1/251‑215 (MQ=255)                                                                                                                                                                                                                   cccAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgaa                                    <  1:85875‑M1/251‑215 (MQ=255)                                                                                                                                                                                                                          aaTGTGGTAGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgaaaaatcg                              >  1:66066‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                             gtgGTAGCCCGCGCCGTCGCTGTCGACtcatgcggatcctctagagtcgacagatggtgatcctgctgaatttcattacgaccagtctaaaaagcgccggaattcgcgacctgctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgattatgttcacaatttgctgaattgttgtgatgtgatgctctccgcagttcctgaaaagtcacgctgtagcttgaaaaaacgacgt                          <  2:368184‑M1/251‑225 (MQ=255)                                                                                                                                                                                                                                  aGCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgaaaaatcgacgttttt                      <  1:309838‑M1/250‑229 (MQ=255)                                                                                                                                                                                                                                   gCCCGCGCCGTCGCTGTCGACgcatgcggatcctctagagacgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgaaaaatcgacgtttttt                     >  1:256306‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                   gCCCGCGCCGTCGCTGTCGACgcatgcggatccactagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattagcgacctactcgttaatgacaggaaaaagggccattggcaaccagggaaagatgaacgtgatgatcttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgacaaatcgccgtttttt                     >  1:345904‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                      cgcgCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttc                                                                                                       >  2:116611‑M1/1‑18 (MQ=255)                                                                                                                                                                                                                                      cgcgCCGTCGCTGTCGACgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgataatgttc                                                                                                       <  1:116611‑M1/165‑148 (MQ=255)                                                                                                                                                                                                                                               gCTGTCGACgcatgcggatccttaagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccaaggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgaaaaatcgacgttttttacgtggttttccg        >  1:331820‑M1/1‑9 (MQ=255)                                                                                                                                                                                                                                                  ggcgaCgcatgcggatcctctagcgtcgacatatggtgatccttctgaaattcattacgaccagtctaaaaagcgcctgaagtcgcgaccttctcgttactgacaggaaaatgggctattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcacctcatttcctgaaaattcacgctgtatcgtgaaaaatcgacggtttttacttggttgtccgt       <  2:31626‑M1/247‑244 (MQ=255)                                                                                                                                                                                                                                                  gTCGACgcatgcggatcctctagagtcgacatatggtgagcctgctgaagttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgataggaaaatgggccattggcatccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgaaaaatcgacgttttttacgtggttttccgt       <  2:331820‑M1/249‑244 (MQ=255)                                                                                                                                                                                                                                                      aCgcatgcggatcctctagagtcgacatatggtgatcctgctgaatttcattacgaccagtctaaaaagcgcctgaattcgcgaccttctcgttactgacaggaaaatgggccattggcaaccagggaaagatgaacgtgatgatgttcacaatttgctgaattgtggtgatgtgatgctcaccgcatttcctgaaaattcacgctgtatcttgaaaaatcgacgttttttacgtggttttccgtcgaaa  >  1:102954‑M1/1‑2 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         GGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATCTTTGCCGTTGGTGGCTTTAATGCGCTCGCAGTGGGAAAGGAAGTGAGCCGGAAGTTCGGAAGTGGGGACGGCTTCGTTGGTTAACCCCCAAAGAATGTGGTAGCCCGCGCCGTCGCTGTCGACTTTTACGATGCCTAAGTTAGCCGCAGGATCAAGCTGGACGTTACGGAAGAATTTGCCCGAGCCGGTGACAATAAACGGTGTATTTGCCAGTAAAGGCATGGGCTGGCTGAGCGGGATATAGCGCGGTTGTTGGTGGAAATTGTCGTGATATGGTGCGATATCGGCGTCATCCAGGCGTAGCGTCAGGTTGCCGCCGTTGCGCTCATCCCAGCCTTTCAGCCAGGCGTCGGTGGTGGCTTTGATCATTC  >  NC_000913/4093744‑4094238

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