BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	72,367	C→T	S187S (TCG→TCA)	thiQ ←	thiamine ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	72,367	0	C	T	94.7%	335.9 / 3.8	117	S187S (TCG→TCA)	thiQ	thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand): ref base C (3/3); new base T (70/38); total (74/43)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 117 total aligned reads displayed.

GTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAG‑ACAACGCGCTAACGCCACTCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCT  >  NC_000913/72122‑72611

gTCTTGTGCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGCGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGacac                                                                                                                                                                                                                                                    >  1:130820/1‑251 (MQ=255)

    tgtgCGGGTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACAcc                                                                                                                                                                                                                                                   <  1:496612/248‑1 (MQ=255)

           gTTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGACGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGCCACCATCAATaacg                                                                                                                                                                                                                                         >  2:422700/1‑248 (MQ=255)

             tGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTc                                                                                                                                                                                                                                       >  1:402160/1‑251 (MQ=255)

                     gCGACGCTCACCGCGCCTTATCCGGCCTACAAAACGCACTGCCCGTAGGTCGGATAAGTTTTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACATCCCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATTCGCGCGGCATCTTCCACACTGTGTGACACCAGCAATAACGTCATTTTTtg                                                                                                                                                                                                                               <  1:207806/251‑1 (MQ=255)

                         cgctgaccGCGTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTtgctgc                                                                                                                                                                                                                           <  1:82102/245‑1 (MQ=255)

                                cgcgTCTTATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGc                                                                                                                                                                                                                    <  1:391518/251‑1 (MQ=255)

                                       tATCAGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGc                                                                                                                                                                                                             <  1:649887/250‑1 (MQ=255)

                                        aTCAGGCCTACAAAACGCACTACCCGTAGGTCGGCTAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCACCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCCTTTTTTGCTGCTGGCAGCTC‑GGGCTc                                                                                                                                                                                                           >  2:325396/1‑251 (MQ=255)

                                           aGGCCTACAACACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTACCCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGACGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGCGTGACACCATAAATAACGTCATTTTTTGCTGCTGGCACCTC‑GGGCTAAcc                                                                                                                                                                                                        >  2:533882/1‑251 (MQ=255)

                                           aGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCAcc                                                                                                                                                                                                        >  1:253871/1‑251 (MQ=255)

                                           aGGCCTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGGCATTTTTTGCTGCTGGCAGCTC‑GTGCTCAcc                                                                                                                                                                                                        >  2:172565/1‑251 (MQ=255)

                                               cTACAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCg                                                                                                                                                                                                    >  1:552869/1‑251 (MQ=255)

                                                  cAAAACGCACTGTCCGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCa                                                                                                                                                                                                 <  1:538629/251‑1 (MQ=255)

                                                   caaaCGCACTGTCCGTAGGCCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGCCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTc                                                                                                                                                                                                  >  2:87633/2‑249 (MQ=255)

                                                        gcactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCAtct                                                                                                                                                                                           >  1:508501/8‑251 (MQ=255)

                                                         cactgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCCTCTTCCACACTGTGTGACACCATCCATAACGTCCTTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCAtctc                                                                                                                                                                                          >  2:44122/7‑251 (MQ=255)

                                                            tgtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCCTCAATAACGCAATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCACCTg                                                                                                                                                                                       >  2:149741/4‑251 (MQ=255)

                                                            tACCCGTAGGTCAGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAACATCTCCTg                                                                                                                                                                                       >  1:55711/1‑251 (MQ=255)

                                                             gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCt                                                                                                                                                                                        <  1:228813/247‑1 (MQ=255)

                                                             gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCt                                                                                                                                                                                        <  1:144894/247‑1 (MQ=255)

                                                             gtccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTg                                                                                                                                                                                       >  1:664746/3‑250 (MQ=255)

                                                               ccGTAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACg                                                                                                                                                                                    >  2:522392/1‑251 (MQ=255)

                                                                  tAGGTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTa                                                                                                                                                                                  >  1:215659/1‑250 (MQ=255)

                                                                    ggTCGGATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACg                                                                                                                                                                               >  1:588556/1‑251 (MQ=255)

                                                                        ggATAAGGTGTTCACGCCGCATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCCGCATCTCCTGACGTAACGcagg                                                                                                                                                                           >  2:485785/1‑248 (MQ=255)

                                                                          aTAAGGCTTTCACGCCGCACCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGAt                                                                                                                                                                         <  2:659108/251‑1 (MQ=255)

                                                                             aggtgTTCACGCCGCATCTGACAGTGCATACTAACCAGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCCCCAGCGTCAGCATCTCCTGACGTACCGCCGGATCg                                                                                                                                                                       >  2:389903/5‑250 (MQ=255)

                                                                                          gcATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACgg                                                                                                                                                         >  2:236088/1‑251 (MQ=255)

                                                                                          gcATCTGACAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACg                                                                                                                                                          <  1:236088/250‑1 (MQ=255)

                                                                                                  cAGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCACAGTCAGAATCCCCTGACGTAACGCCGGATCAAGCGAAGCGAACGTTTCATCg                                                                                                                                                  >  1:491846/1‑250 (MQ=255)

                                                                                                   agagCATACTAACCCGTACTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGCTGCGCCCGTCGGCGACCACCACCGAGCGCGTGTCGATCCGCGCGACATCTTCCCCACAGTGTGACACCCTCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGAAGTAACTCCGGCTCGAGCGCCGAGAACGGTTCATCg                                                                                                                                                  >  2:122261/4‑249 (MQ=255)

                                                                                                   aGTGCATACTAACCCGTACTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGAAGAAAACGGTTCATCgcg                                                                                                                                                >  2:68298/1‑249 (MQ=255)

                                                                                                   aGTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATCACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCgcg                                                                                                                                                >  1:392351/1‑249 (MQ=255)

                                                                                                      gCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCa                                                                                                                                              <  1:146419/250‑1 (MQ=255)

                                                                                                        aTACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATa                                                                                                                                           <  1:207359/251‑1 (MQ=255)

                                                                                                             aaCCCGTAATCCCAAATAGAGCTGAAGCACTCGCCTTACCGCTCACCAACTCCTTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCGCGCGGCGATCCGCGCGGCATCTTCCACACTGTGTGCCACCATCAATACCGTCATTTTTTGCTGCTTGCAG‑TCAGTGCTCACCCGCGTCAGCATCTCCTGACGTCACGCCGCATCGCGCGCAgagacagcttcatctagcactcaaatc                                                                                                                                      >  2:207164/1‑228 (MQ=255)

                                                                                                               cccGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTAATCGAGCAATacattcgg                                                                                                                                    >  2:707109/1‑244 (MQ=255)

                                                                                                                     aTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTc                                                                                                                              >  2:272069/1‑251 (MQ=255)

                                                                                                                        cccAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCACTAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcc                                                                                                                           >  2:166643/1‑250 (MQ=255)

                                                                                                                        cccAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAAACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGAATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTcgcg                                                                                                                           >  1:348263/1‑251 (MQ=255)

                                                                                                                            atagggctgAAGCACTAGCCTTACCGCTCAACAACTCATAGGTCTTACCATGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACAACATCAATACCGTCATTTTTTGCTGATGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCAGACGTAACGCCGGATCGAGCGCAGAGAACGGTACATCGAGCAATAAAATCGGCTGTTcgcgcccc                                                                                                                       >  2:22573/9‑247 (MQ=255)

                                                                                                                            aTAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCAcc                                                                                                                       >  1:446767/1‑251 (MQ=255)

                                                                                                                                     aaGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGCTGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGCCGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAAcgc                                                                                                             >  2:229188/1‑251 (MQ=255)

                                                                                                                                     aaGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAAcgc                                                                                                             >  1:33074/1‑251 (MQ=255)

                                                                                                                                     aaGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAAcgc                                                                                                             >  1:92005/1‑251 (MQ=255)

                                                                                                                                           tCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAAc                                                                                                       >  1:545725/1‑251 (MQ=255)

                                                                                                                                            cGCCTTACCGCTCACCAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAAc                                                                                                       <  1:663763/250‑1 (MQ=255)

                                                                                                                                              ccTTACCGCTCAACAACTCAGTGGTATTACCCTGTCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGCTCCGCGCGGCATCTCCCACACTGTGGGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGcc                                                                                                    <  2:130820/251‑1 (MQ=255)

                                                                                                                                                 tACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCAcg                                                                                                 >  1:203184/1‑250 (MQ=255)

                                                                                                                                                 tACCGCGCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCa                                                                                                   <  2:290148/249‑1 (MQ=255)

                                                                                                                                                    cGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGAGGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGCGACAGGTTCATCGAGCAATAAAATCGGCTGTTAGCGCACCAG‑ACACCGCGCTAACGccccgcac                                                                                              >  1:75190/1‑245 (MQ=255)

                                                                                                                                                     gCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑t                                                                                            >  1:377045/1‑251 (MQ=255)

                                                                                                                                                       tCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCCTTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGCGCACTAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTACCGCCCCgcgc‑tgt                                                                                          >  2:209334/1‑251 (MQ=255)

                                                                                                                                                       tCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑tgt                                                                                          >  1:473237/1‑251 (MQ=255)

                                                                                                                                                       ccaacacatcctggctcttACCCTGGCTGGCGATGCGCCCGTCGGCGACCACCACCGAGCGTGTGGCGAGCCGCGCGGCATCTTCCACACCGTTTGACACCATCTATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGATAACGGTTCATCGAGCAATAAAATCGGCTGGTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑tg                                                                                           <  2:380863/232‑1 (MQ=255)

                                                                                                                                                          acaacTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑tgtcgc                                                                                       >  1:449623/1‑251 (MQ=255)

                                                                                                                                                           caacTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAATTCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑tgtcgct                                                                                      >  1:256999/1‑251 (MQ=255)

                                                                                                                                                            cacTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGACCCCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACTTAACGCCGGATCGAGCGCAGAGAAAGGTTCATCGAGCAATAAAATCGGCTGTTCCCGAACAAG‑ACAACGCGGTAACGCCACgcca‑tgtcgctg                                                                                     >  2:401387/2‑251 (MQ=255)

                                                                                                                                                            aacgcATTGGTTGTACCCGGCCAGGCGAGGCGCCCGTCGGCGACCGCCACCGAGCGCGTGGCGATCCGCGCGGCATCGTCCACACTGTGGGACACAATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCGGTTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑tgtcgctg                                                                                     <  2:126087/247‑1 (MQ=255)

                                                                                                                                                                     gtcttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGa                                                                            <  2:377045/247‑1 (MQ=255)

                                                                                                                                                                      tcttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGCAAAAC‑CGCTAACGCCCCCCCC‑TGTCGCTGCCCCCCGGa                                                                            >  2:433216/4‑250 (MQ=255)

                                                                                                                                                                        ttccccTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑tgtcgctg                                                                                     <  2:320021/236‑1 (MQ=255)

                                                                                                                                                                        ttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACgcgc‑tgtcgctg                                                                                     >  1:320021/2‑239 (MQ=255)

                                                                                                                                                                         ttcccTGCCAGGCGATGTGCCCGTCGGCGACCACCCCCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGATACCATCAATATCGTCATTTTTTGCTGCTGGCAGCTC‑GCGCTCACCAGCGTTAGCACCTCCTGACGTAACGCCGGACCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGTTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAg                                                                         <  1:258520/248‑1 (MQ=255)

                                                                                                                                                                         tACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGaaa                                                                          <  1:533882/249‑1 (MQ=255)

                                                                                                                                                                         tACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCAGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGaaa                                                                          <  1:189778/249‑1 (MQ=255)

                                                                                                                                                                          aCCCTGCCAGGCGATGCGCCCGTCGGCGCCCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGCCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTAATCGAGCAATAAAATCGGCTGTTAGCGCACCGG‑ACAACGCTCTACCGCCACACGCCTG‑CGCTGACCGCCGGAAAggc                                                                       >  2:317465/1‑249 (MQ=255)

                                                                                                                                                                                 cAGGCGATGCGCCAGTCGGCGAACACCACCGAGAGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGACGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACAGCCGGAAAGCTCGccagg                                                                >  2:316903/1‑248 (MQ=255)

                                                                                                                                                                                      gATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑CCAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCCGGTAAcc                                                           >  2:505551/1‑251 (MQ=255)

                                                                                                                                                                                       aTGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCACTAAAATCGGCTGTTCGCGCACCCG‑ACAACCCGCTAACGCCACGCCC‑TGTACCTGACCGCCGGaaatcccgccagggaaccc                                                          >  2:399025/1‑235 (MQ=255)

                                                                                                                                                                                         gcgcCCGCCGGCGACCACCACCGAGCGCGGGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACcgcg                                                        <  1:490892/251‑1 (MQ=255)

                                                                                                                                                                                              cGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCAtt                                                   >  2:100214/1‑251 (MQ=255)

                                                                                                                                                                                              cGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCAtt                                                   <  1:348239/251‑1 (MQ=255)

                                                                                                                                                                                               gTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGccctt                                                   >  2:677608/1‑247 (MQ=255)

                                                                                                                                                                                                  ggCGACCACCACCGCGCGCGTGGCGATCCGCGCGGCATCTGCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGCTCGAGCGCAGAGACCGGTTCATCGAGCAATAACATCGGCTGTTCGCACACCAG‑ACCACGCGCTCACGCCACCCGC‑TGTAGCTGACCTCCGGAAAGCTCGCCAGAGAACCGCGcccttcat                                                >  2:453347/1‑244 (MQ=255)

                                                                                                                                                                                                     gaccaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTAt                                            >  1:311618/8‑251 (MQ=255)

                                                                                                                                                                                                      accaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGa                                                                                    <  1:687268/204‑1 (MQ=255)

                                                                                                                                                                                                         actacCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAG‑ACAACGCGCTAACGCCACTCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAAt                                        <  1:325396/251‑1 (MQ=255)

                                                                                                                                                                                                            acCGAGCGCGTGTCGAGCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCCCCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATccc                                     <  2:449623/251‑1 (MQ=255)

                                                                                                                                                                                                            acCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGCTCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCCG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAAGTAACCGCGCCCTTAAATTATCAATccc                                     >  2:367126/1‑251 (MQ=255)

                                                                                                                                                                                                                aGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATc                                 >  1:433131/1‑251 (MQ=255)

                                                                                                                                                                                                                 gcgcGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATACCGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCt                                <  1:118723/251‑1 (MQ=255)

                                                                                                                                                                                                                     gTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCTGAAAGCGCGACAGGCAACCGCGCCATTAAATTATCAATCCCCATCTggcg                            >  2:15894/1‑251 (MQ=255)

                                                                                                                                                                                                                           aTCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGAt                      <  1:433216/251‑1 (MQ=255)

                                                                                                                                                                                                                            cccccaccgcATCTGCCACGCGGGGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGCACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAG‑ACAACGCGCTATCGCCGCTCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTggcgggcgtta                     <  2:23102/245‑4 (MQ=255)

                                                                                                                                                                                                                             ccGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGGCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTTGCGGGCGAt                      >  2:336952/1‑249 (MQ=255)

                                                                                                                                                                                                                              cgcgcgGCATCTTCCTCACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCCCCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGc                   <  1:270353/251‑1 (MQ=255)

                                                                                                                                                                                                                                cgcggcATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGCTAGCGt                 >  1:369299/5‑251 (MQ=255)

                                                                                                                                                                                                                                   ggcATCTTCCACTCTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGTCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTTACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCa              <  1:49239/250‑1 (MQ=255)

                                                                                                                                                                                                                                    gcATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACAAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAAAGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATAAACTTATCAATCCCCATCTGGCGGGCGATAGCGTGCAt             >  1:377407/1‑251 (MQ=255)

                                                                                                                                                                                                                                     cATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAtt            >  2:530579/1‑251 (MQ=255)

                                                                                                                                                                                                                                      aTCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGCGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAttt           >  1:432114/1‑251 (MQ=255)

                                                                                                                                                                                                                                         ttCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGACTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCACCAGGTAACCGCGCCATTAAATTATCAATCCACATCTGGCGGGCGATAGCGTGCATTTTcc        >  2:513750/1‑251 (MQ=255)

                                                                                                                                                                                                                                          tCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTct        <  1:485996/250‑2 (MQ=255)

                                                                                                                                                                                                                                          tCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTct        <  1:100214/250‑2 (MQ=255)

                                                                                                                                                                                                                                          tCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGGCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTct        <  1:149741/250‑2 (MQ=255)

                                                                                                                                                                                                                                           ccACCCTGTGTGACACCATCAATAACTTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTctcc      <  2:253871/251‑4 (MQ=255)

                                                                                                                                                                                                                                           ccACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTctcc      >  2:607253/1‑248 (MQ=255)

                                                                                                                                                                                                                                               acTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAG‑ACAACGCGCTAACGCCACGCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCctgct  <  1:272069/250‑1 (MQ=255)

GTCTTGTGCGGGTTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACAAAACGCACTACCCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAGTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTC‑GTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAG‑ACAACGCGCTAACGCCACTCGC‑TGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCT  >  NC_000913/72122‑72611

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