Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,728,118	(C)8→7	coding (1254/1455 nt)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,728,111	0	C	.	87.7%	223.5 / 30.3	66	coding (1261/1455 nt)	xylB	xylulokinase
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base . (29/28);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TCAGAATAAATTCATACTAAATATTAATTAATTACTGAGATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCGGAAAGATAAGGCAGAAACCA  >  NC_000913/3727830‑3728390                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                         tCAGAATAAATTCATACTAAATATTAATTAATTACTGAGATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccg                                                                                                                                                                                                                                                                                   >  2:630021/1‑288 (MQ=255)           tacATACTAAATATTAATTAATTAGTGAGATAGATAGGTGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGGCGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTGGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATc                                                                                                                                                                                                                                                                        <  1:357358/287‑1 (MQ=255)                    aaTATTAATTAATTACTGAGATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGctgctg                                                                                                                                                                                                                                                               >  2:324181/1‑289 (MQ=255)                                                gtgtgAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                                                                                                                                                   <  1:576148/288‑1 (MQ=255)                                                aTGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                                                                                                                                                   <  2:132812/289‑1 (MQ=255)                                                 tgtgAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACTTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                                                                                                                                                   <  1:323361/289‑1 (MQ=255)                                                  gtgAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATCTCCGCCAGCATCTCAccc                                                                                                                                                                                                                                   >  2:365446/1‑285 (MQ=255)                                                  gtgAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGc                                                                                                                                                                                                                                   >  2:202512/1‑287 (MQ=255)                                                     aaTTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGACCTGCTTACCCGTGAACTCCGCCAGCATCTGACGCTAg                                                                                                                                                                                                                                >  2:305102/1‑288 (MQ=255)                                                     aaTTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCC‑TACGGTAATCGATCTGCTGACCGCCGACATGGGCCAGAAtagtaaccaagtc                                                                                                                                                                                                                              >  2:568273/1‑277 (MQ=255)                                                     aaTTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGCGCTGACGACCGCTGATAACCGCCAGCCCCTGACGCCCGTa                                                                                                                                                                                                                              >  2:588719/1‑289 (MQ=255)                                                     aaTTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGTTGATATCCGACAGCATCTGACGCGAGt                                                                                                                                                                                                                               >  2:668566/1‑288 (MQ=255)                                                     aaTTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCTCCAGCCTCTGACGCCAGtc                                                                                                                                                                                                                              >  1:312918/1‑288 (MQ=255)                                                     aaTTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTCCCGCCGGt                                                                                                                                                                                                                               >  2:265271/1‑288 (MQ=255)                                                     aaTTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCAAGCTGCCGACCGCTGATATCCGCCAGCATCTGACGCCAGTa                                                                                                                                                                                                                              >  2:446416/1‑289 (MQ=255)                                                          tCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGGTGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTATCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCGCACAT‑CCCCCCCGTACGGTAAACGAGCTGCTGACCGGTGATTGCCGAGAGCATCTGACGCCAGTACTCAc                                                                                                                                                                                                                         >  1:92917/1‑289 (MQ=255)                                                             cccACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                                                                                                                                                 <  2:501868/278‑1 (MQ=255)                                                             cccACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                                                                                                                                                 >  1:501868/1‑278 (MQ=255)                                                                           caggGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCa                                                                                                                                                                                                        <  1:202512/289‑1 (MQ=255)                                                                             ggggATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                                                                                                                                                                      <  1:227468/289‑1 (MQ=255)                                                                             ggggATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                                                                                                                                                                      <  1:283445/289‑1 (MQ=255)                                                                             ggggATAACGTTTACGCCATTAAGTGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAAt                                                                                                                                                                                                      <  1:378574/289‑1 (MQ=255)                                                                                        ttACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTAACTACGCGCCCCGCCCCCACTCAACGTAacac                                                                                                                                                                                           >  1:333870/1‑289 (MQ=255)                                                                                         tACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTCCTGCGCCCAGTGCTGGCCCCACA‑GCCCCCCCGTACGCTAATCGAGCTGCTGACCGCTCATATCCGCCAGCATCTGACCCCCGTACACACTCGGCCCCCCGCTCTCaacaacacccacaat                                                                                                                                                                                          >  2:636750/1‑276 (MQ=255)                                                                                          aCGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCCGTACTCCCTACGCTCGCCCACCACAATCAAAGTAaccatt                                                                                                                                                                                         >  2:652686/1‑285 (MQ=255)                                                                                          aCGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACtt                                                                                                                                                                                         >  1:381442/1‑289 (MQ=255)                                                                                                ttAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCCCTACGCGCCCCGCCCCCAATCAACGTAACCATTTTCGGt                                                                                                                                                                                   >  1:460435/1‑289 (MQ=255)                                                                                                         gAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGGGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCg                                                                                                                                                                          <  1:568273/289‑1 (MQ=255)                                                                                                          aaGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCGCAGTGCTGGCCGCACAT‑CCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCAAGTACTCAATACGCGCCGCGCCCCCAATCAAAGTCACACTTTTCCGTTTAATACCGc                                                                                                                                                                         >  1:552776/1‑289 (MQ=255)                                                                                                                cTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGCTTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGGCCCGCCCCGAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCac                                                                                                                                                                   >  2:30526/1‑288 (MQ=255)                                                                                                                    tAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGGATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCAc                                                                                                                                                               <  1:355828/289‑1 (MQ=255)                                                                                                                       agGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccgta                                                                                                                                                                                                                                                                                 <  1:617427/173‑4 (MQ=255)                                                                                                                       agGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATcccccccgta                                                                                                                                                                                                                                                                                 >  2:617427/1‑170 (MQ=255)                                                                                                                                                     tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCTCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACGGATGATATCCGCCAGCATCTGACGCCAGGACTCACTCAGCGCCCCGCGCCCAATCAAACTAACATTTTGGGGTTTAATAACGCAGTGCTGCACGACCTC‑AATGCCATCTGCGAGCGGCTAACCCa                                                                                                                              >  1:647004/1‑289 (MQ=255)                                                                                                                                                     tAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGCTATCCGCCAGCATCTGACGCCCGTACTCCCTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGCCATGCACGACACC‑CATGCCACCTGCCCGCGCCTCACCCAc                                                                                                                             >  2:180428/1‑289 (MQ=255)                                                                                                                                                      aaGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTGACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACg                                                                                                                            <  2:312758/289‑1 (MQ=255)                                                                                                                                                       aGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACg                                                                                                                            <  1:345547/288‑1 (MQ=255)                                                                                                                                                           gcctaccGCTGCGCATCTGGTGGATGCGACTGTTGTAACGGTAGTTGCGGCAACAATTCAATGGGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTACCACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTc                                                                                                                       <  1:238126/283‑1 (MQ=255)                                                                                                                                                           gCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCACGCATGCACGAAATC‑CATCCCATCTCCCAGCGGATAACACCCGCCCTc                                                                                                                       >  2:726962/1‑289 (MQ=255)                                                                                                                                                            cATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACCCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCCCGCCATC‑CATTCCATCTGCCAGCGCATAACCCAGCCCttac                                                                                                                      >  2:31437/1‑287 (MQ=255)                                                                                                                                                                     tcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                                                                                                                                                 <  2:165973/173‑1 (MQ=255)                                                                                                                                                                     tcgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCa                                                                                                                                                                                                                                 >  1:165973/2‑174 (MQ=255)                                                                                                                                                                      cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGCGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCCCCCCCAATCAACGTAACACTTTGCGGTTTAATACCCCAGGCATGCACGACATCACAT‑CCATCTGCCAGCGCATAACCACCCCCTTCCACCACTGccc                                                                                                            >  2:64693/1‑287 (MQ=255)                                                                                                                                                                      cgcATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGc                                                                                                              >  2:13767/1‑287 (MQ=255)                                                                                                                                                                         aTCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTcgcg                                                                                                         <  1:324181/289‑1 (MQ=255)                                                                                                                                                                                  aTGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAt                                                                                                <  1:630021/289‑1 (MQ=255)                                                                                                                                                                                   tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGGGCGATTTCTCTGGATTCGCCTCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCAtt                                                                                               <  1:365446/289‑1 (MQ=255)                                                                                                                                                                                   tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGc                                                                                                                                                                         >  2:343964/1‑217 (MQ=255)                                                                                                                                                                                   tGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGc                                                                                                                                                                         <  1:343964/217‑1 (MQ=255)                                                                                                                                                                                      gACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGAAATCCGCCATCATCTGACGCCAGTCCTCACTACGCTCCCCGCCCCCAATCAACGTAACACTTTTCGCTTTAATACCGCCGCCCTGCACGACATC‑CATGCCATCTGCCAGCCCATAACCCACGCCTTCCAGCACTGCTCCCGCCAGTTCATCgg                                                                                             >  1:62373/1‑289 (MQ=255)                                                                                                                                                                                        cTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCCGCCAGCGCCTAACCCACGCCCTCCAGCACTGCCCGCGCCAGTTCATTGGGGc                                                                                          >  2:295112/1‑289 (MQ=255)                                                                                                                                                                                                     aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑GCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGCGt                                                                             >  1:650738/1‑289 (MQ=255)                                                                                                                                                                                                     aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGt                                                                             >  1:190377/1‑289 (MQ=255)                                                                                                                                                                                                     aGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCCTCCAGCACTGCTCGCCCCCGTTCATTGGGGCCATGTTGACGAGt                                                                             >  1:195443/1‑289 (MQ=255)                                                                                                                                                                                                                    gcAATGTGCGATTTCGCGTGATTCGCCCCTACCTCGGCCAGCCTTCCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGaaaa                                                               <  2:552776/287‑1 (MQ=255)                                                                                                                                                                                                                             gtttcgtcgggATTCGCCGCGTTCTGCCCCAGCCTTGCTGCGCCCATTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGc                                                     <  2:381442/280‑1 (MQ=255)                                                                                                                                                                                                                             gATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAGGAAAACCCCCTTCGc                                                     <  2:190377/289‑1 (MQ=255)                                                                                                                                                                                                                                  ctctttaTTCACCGCGATCTGCCCCAGCCTTGCTGCGCCCAGTGCTGTCCCCACAT‑CCCCCCCGTCCCGTAATCCAGCTGCTGACCGCTGTTATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCGATCACCGTAACACTTTGCAGTTTAAGACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGCTGAGTCAAACCAAAGAAAACCCCCTTCGCCt                                                   <  2:92917/280‑1 (MQ=255)                                                                                                                                                                                                                                  ctctGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTggg                                                <  1:722171/289‑1 (MQ=255)                                                                                                                                                                                                                                   gctggtgccgccgcGATCTGCGCCAGCCTTGCTGCGCCCAGGGCTGGCCCAACA‑CCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTggg                                                <  2:312918/280‑1 (MQ=255)                                                                                                                                                                                                                                                             ccTTGCTGCCCCCACTACCTGCCCCCCAT‑CCCCCCCGTACGGTAATCCAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCgg                     <  1:680801/289‑1 (MQ=255)                                                                                                                                                                                                                                                                 gctgcGCCCAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTTGGGCGGACGCTCGCCGGa                    >  2:508838/1‑286 (MQ=255)                                                                                                                                                                                                                                                                     cgcCCAGTCCTGCCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCGGAAAGATa              <  1:668566/288‑1 (MQ=255)                                                                                                                                                                                                                                                                       cccAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTAAtgtg                                      <  1:152549/262‑1 (MQ=255)                                                                                                                                                                                                                                                                       cccAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCACTTCGCCTGGGGATTAAtgtg                                      >  2:152549/1‑262 (MQ=255)                                                                                                                                                                                                                                                                        ccAGTTCTCGCCCCACAT‑CCCCCCCGTACGGTACTCGAGCTGCTGACCGCTGATATCCGCCAGCACCTTACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATTCCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCGGAAAGATAAGGc          <  1:588719/289‑1 (MQ=255)                                                                                                                                                                                                                                                                        ccAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCGGAAAGATAAGGc          <  1:265271/289‑1 (MQ=255)                                                                                                                                                                                                                                                                         cAGTGCTGGCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCGGAAAGATAAGGCa         <  1:726962/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                 gCCCCACAT‑CCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCCACCCCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGCAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCGGACAGATAAGGGAGAAACCa  >  1:188644/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                         TCAGAATAAATTCATACTAAATATTAATTAATTACTGAGATATATAGATGTGAATTATCCCCCACCCGGTCAGGCAGGGGATAACGTTTACGCCATTAATGGCAGAAGTTGCTGATAGAGGCGACGGAACGTTTCTCGTCGTGGCTGATAAGCGGCATAACGCTGCGCATCTGGTAGATGCGACTGTTCTAACGGTAGTTGCGGCAACAATTCAATGAGCGATTTCTCTGGATTCGCCGCGATCTGCGCCAGCCTTGCTGCGCCCAGTGCTGGCCCCACATCCCCCCCCGTACGGTAATCGAGCTGCTGACCGCTGATATCCGCCAGCATCTGACGCCAGTACTCACTACGCGCCCCGCCCCCAATCAACGTAACACTTTGCGGTTTAATACCGCAGGCATGCACGACATC‑CATGCCATCTGCCAGCGCATAACCCACGCCTTCCAGCACTGCTCGCGCCAGTTCATTGGGGCCATGTTGATGAGTCAAACCAAAGAAAACCCCCTTCGCCTGGGGATTATTGTGTGGCGTACGCTCGCCGGAAAGATAAGGCAGAAACCA  >  NC_000913/3727830‑3728390

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