BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,083,186	C→T	V204V (GTG→GTA)	speB ←	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,083,186	0	C	T	99.1%	403.2 / ‑4.2	115	V204V (GTG→GTA)	speB	agmatinase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (52/62); total (52/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 115 total aligned reads displayed.

GCAGTGATTTCCGACTGATCGTATGCCGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCAC  >  NC_000913/3082941‑3083430

      aTTTCCGACTGATCGTATGCCGGGGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGc                                                                                                                                                                                                                                           >  2:54728/1‑251 (MQ=255)

          ccGACTGATCGTATGCCGGAGCCACTTCCTCTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAATCAGGATCCATGCAGTCGATATCAAAAGTCAGGTAAACCGTCATATCACCCACAATCGGTGTCACTGGGGCGATAACGTCATCCTCGCTGCGATCGTTTACCTGGCAGGCGTc                                                                                                                                                                                                                                        <  2:393299/250‑1 (MQ=255)

                gATCGTATGCCGTAGCCACTTCCACTACGTCCATACCAACAATGTTGAGATCTTACAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCGCTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGTCATATCACCCTCAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCAc                                                                                                                                                                                                                                  <  1:34700/250‑1 (MQ=255)

                gATCGTATGCCGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCAc                                                                                                                                                                                                                                  <  1:323702/250‑1 (MQ=255)

                gATCGTATGACGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACg                                                                                                                                                                                                                                 >  2:80103/1‑251 (MQ=255)

                 aTCGTATGCCGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGTCGTCCAGCACgg                                                                                                                                                                                                                                >  2:371033/1‑251 (MQ=255)

                   cGTATGCCGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGt                                                                                                                                                                                                                               <  1:514307/250‑1 (MQ=255)

                   cGTATGCCGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGt                                                                                                                                                                                                                               <  1:317744/250‑1 (MQ=255)

                        gCCGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGcc                                                                                                                                                                                                                         <  2:507319/251‑1 (MQ=255)

                         ccGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCg                                                                                                                                                                                                                        <  1:316374/251‑1 (MQ=255)

                          gggAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTGTAATAGCGCGATCGGAGGTCATGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGAGATTAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTGGTGCGATAACGTCATCCACGCTGCGATCGTTTACCTGTCAGGCGTCCATCACGGTagatacgt                                                                                                                                                                                                                       <  2:414020/250‑8 (MQ=255)

                           ggAGCAACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGtt                                                                                                                                                                                                                      >  2:648422/1‑251 (MQ=255)

                            gAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTg                                                                                                                                                                                                                     <  1:8435/251‑1 (MQ=255)

                             ggCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGt                                                                                                                                                                                                                    <  2:205057/250‑1 (MQ=255)

                                   ttCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTAt                                                                                                                                                                                                              >  1:322877/1‑251 (MQ=255)

                                      cACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCa                                                                                                                                                                                                            >  1:429693/1‑250 (MQ=255)

                                         tACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACt                                                                                                                                                                                                        >  1:438160/1‑251 (MQ=255)

                                          aCGTCCATCCCAACAATGTTTAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGTCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCATGTAAACCGGCATATCTCCCACAATCTGTTTCACTTGGGCGATAACGTCCTCCACGCTGCAATCGTTTACCTGTCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTc                                                                                                                                                                                                       <  2:275089/251‑1 (MQ=255)

                                          aCGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTc                                                                                                                                                                                                       <  1:390069/251‑1 (MQ=255)

                                           cGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCg                                                                                                                                                                                                      <  2:322877/251‑1 (MQ=255)

                                              ccATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTa                                                                                                                                                                                                   <  1:355688/251‑1 (MQ=255)

                                              ccANCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCGCCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCGGTTGTCTTTATCAAACTCGGTa                                                                                                                                                                                                   <  2:584731/251‑1 (MQ=255)

                                               cATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTAc                                                                                                                                                                                                  >  1:143786/1‑251 (MQ=255)

                                               cATCCCAACAATGTTGAGATCTGTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTAc                                                                                                                                                                                                  <  2:48221/251‑1 (MQ=255)

                                               cATCCAAACAATGTTGAGATCTTTCAGGCCACGTAACAGGTTAATAGAGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGAGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGACAGGTAAACCGGCATATCACCCAAAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCAAGCACGGTAAAGCAGTTGTATTTCTCAAACTCgnaa                                                                                                                                                                                                   >  2:609930/1‑247 (MQ=255)

                                                     aacaaTGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTTCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCTGCATATCACCCACAAGCTGTTTCACGTGGGCGATAACGTCATCCACGCTGCGATCTTGTACCTGGCAGGCCTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATAc                                                                                                                                                                                            <  1:80103/251‑1 (MQ=255)

                                                     aacaaTGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATAc                                                                                                                                                                                            <  1:537161/251‑1 (MQ=255)

                                                             ggagaTCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTg                                                                                                                                                                                     <  2:310020/249‑1 (MQ=255)

                                                              gagaTCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGc                                                                                                                                                                                    >  2:448563/1‑250 (MQ=255)

                                                               agaTCGTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCAc                                                                                                                                                                                  <  1:394874/251‑1 (MQ=255)

                                                                gaTCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACg                                                                                                                                                                                 >  1:536324/1‑251 (MQ=255)

                                                                 aTCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACgg                                                                                                                                                                                >  2:268506/1‑251 (MQ=255)

                                                                   cTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCGGCACGGa                                                                                                                                                                               >  2:532279/1‑250 (MQ=255)

                                                                     ttCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGaa                                                                                                                                                                              <  1:268506/249‑1 (MQ=255)

                                                                      tCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATACAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATg                                                                                                                                                                            >  1:630983/1‑250 (MQ=255)

                                                                       cAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGAt                                                                                                                                                                          >  2:254404/1‑251 (MQ=255)

                                                                           ccGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCgg                                                                                                                                                                      <  2:490983/251‑1 (MQ=255)

                                                                               gTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCg                                                                                                                                                                  >  2:430362/1‑251 (MQ=255)

                                                                               gTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCg                                                                                                                                                                  >  1:184760/1‑251 (MQ=255)

                                                                                 aCCAGGTTAATAGCGCGATCGGAGGTCAGGCCTCCAATCACTGGCGTACAGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACATTCTGTTTCACTTGGGCGATAACGTCATCCACTCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGAt                                                                                                                                                                <  1:124578/251‑1 (MQ=255)

                                                                                        tAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACct                                                                                                                                                         >  1:490766/1‑251 (MQ=255)

                                                                                         aaTAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACctt                                                                                                                                                        >  1:266020/1‑251 (MQ=255)

                                                                                         aaTAGCGCGATCCGAGGTCAGGTCGCCAATCACTGGCGTACCGGTGCCTTGTGCAAAAGCCGGAACCAGGCAGTCCATATCAAAAGTCAGGTAAACCCGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACTTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAGTACCAATCTGCCCGGAATGATTCGGGTCGATCAGACct                                                                                                                                                         <  2:410449/250‑1 (MQ=255)

                                                                                          aTAGCGCGAGCGGAGGTCAGTCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                                                                                                                                       <  1:597717/251‑1 (MQ=255)

                                                                                            aGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttctt                                                                                                                                                     <  1:595732/251‑1 (MQ=255)

                                                                                               gcgATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGTATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCg                                                                                                                                                  <  2:191520/251‑1 (MQ=255)

                                                                                                cgATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCgg                                                                                                                                                 >  2:231123/1‑251 (MQ=255)

                                                                                                 gATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgc                                                                                                                                                <  1:648422/251‑1 (MQ=255)

                                                                                                   tCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGCATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgcgc                                                                                                                                              >  2:116614/1‑251 (MQ=255)

                                                                                                   tCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgcgc                                                                                                                                              >  2:531788/1‑251 (MQ=255)

                                                                                                      gAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGt                                                                                                                                           <  2:135878/251‑1 (MQ=255)

                                                                                                           cAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCGTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGa                                                                                                                                      <  1:129914/251‑1 (MQ=255)

                                                                                                              gccgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACa                                                                                                                                   <  1:527303/251‑1 (MQ=255)

                                                                                                                  ccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGt                                                                                                                               >  2:5576/1‑251 (MQ=255)

                                                                                                                     aTCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGcc                                                                                                                            <  1:333246/251‑1 (MQ=255)

                                                                                                                       cACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGt                                                                                                                          <  1:437026/251‑1 (MQ=255)

                                                                                                                          tGGCGCACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGACTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGAGCGGTATAGAACATAGTGCCGTGGt                                                                                                                       >  2:108531/1‑251 (MQ=255)

                                                                                                                                          ggTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGtt                                                                                                       <  2:266020/251‑1 (MQ=255)

                                                                                                                                             gCAAAAGCATGCTCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTGTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCTCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGc                                                                                                    <  2:630983/251‑1 (MQ=255)

                                                                                                                                             gCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCg                                                                                                     <  2:184760/250‑1 (MQ=255)

                                                                                                                                               aaaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCAt                                                                                                  <  1:288154/251‑1 (MQ=255)

                                                                                                                                                  aGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAgg                                                                                               >  1:478286/1‑251 (MQ=255)

                                                                                                                                                      ggATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATc                                                                                           >  1:363100/1‑251 (MQ=255)

                                                                                                                                                      ggATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATc                                                                                           >  1:512714/1‑251 (MQ=255)

                                                                                                                                                      ggATCCAGGCAGTCGATATCAAAAGGCAGGTAAACCGGCATATCACCCTCAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGGCCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATc                                                                                           <  2:572799/251‑1 (MQ=255)

                                                                                                                                                          ccAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGgtg                                                                                       <  1:430362/251‑1 (MQ=255)

                                                                                                                                                                aGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCg                                                                                 <  2:213718/251‑1 (MQ=255)

                                                                                                                                                                        tCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAgtg                                                                         >  2:216737/1‑251 (MQ=255)

                                                                                                                                                                          aaaaGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAgtgt                                                                        >  1:184502/1‑250 (MQ=255)

                                                                                                                                                                                 aGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGgtg                                                                                       <  1:89121/228‑1 (MQ=255)

                                                                                                                                                                                      aaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCAttt                                                           <  1:639880/251‑1 (MQ=255)

                                                                                                                                                                                            cATATCACCAACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTc                                                                                                                                                                                                                   <  1:103158/93‑1 (MQ=255)

                                                                                                                                                                                            cATATCACCAACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTc                                                                                                                                                                                                                   >  2:103158/1‑93 (MQ=255)

                                                                                                                                                                                             atatCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGa                                                     >  2:521571/1‑250 (MQ=255)

                                                                                                                                                                                             atatCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGa                                                     >  1:236956/1‑250 (MQ=255)

                                                                                                                                                                                               atCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAAt                                                  <  1:531788/251‑1 (MQ=255)

                                                                                                                                                                                               atCACCCACAATCTGTTTCACTTGGGCGATAACGTCATACACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAAt                                                  >  1:124091/1‑251 (MQ=255)

                                                                                                                                                                                                  aCCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCt                                               <  2:552646/251‑1 (MQ=255)

                                                                                                                                                                                                  aCCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTAGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGAAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCt                                               >  1:658315/1‑251 (MQ=255)

                                                                                                                                                                                                       caATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCa                                          <  2:157509/251‑1 (MQ=255)

                                                                                                                                                                                                                tCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACg                                  <  1:254404/250‑1 (MQ=255)

                                                                                                                                                                                                                 cACTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgc                                 <  1:231123/250‑1 (MQ=255)

                                                                                                                                                                                                                  ccTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgcag                               <  1:549580/250‑1 (MQ=255)

                                                                                                                                                                                                                   cTTGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgcagc                              >  2:465625/1‑251 (MQ=255)

                                                                                                                                                                                                                    ttGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgcagc                              >  1:542220/1‑250 (MQ=255)

                                                                                                                                                                                                                     tGGGCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCGGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgcagca                             <  1:532279/250‑1 (MQ=255)

                                                                                                                                                                                                                      gggCGATAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGAACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCagcagc                           >  2:658431/1‑251 (MQ=255)

                                                                                                                                                                                                                            tAACGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcagc                     >  1:366225/1‑251 (MQ=255)

                                                                                                                                                                                                                             aaCGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcagcg                    <  1:5576/251‑1 (MQ=255)

                                                                                                                                                                                                                             aaCGTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGACGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcagcg                    >  1:643259/1‑251 (MQ=255)

                                                                                                                                                                                                                                gTCATCCACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGAATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTa                  >  2:359767/1‑250 (MQ=255)

                                                                                                                                                                                                                                 tcatcaACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAAc                >  2:268896/1‑251 (MQ=255)

                                                                                                                                                                                                                                      cACGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGt           >  2:136287/1‑251 (MQ=255)

                                                                                                                                                                                                                                       aCGCTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTg          <  2:143786/251‑1 (MQ=255)

                                                                                                                                                                                                                                          cTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgcagc                              >  1:482118/1‑228 (MQ=255)

                                                                                                                                                                                                                                          cTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgcagc                              <  2:482118/228‑1 (MQ=255)

                                                                                                                                                                                                                                          cTGCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTc       >  1:146388/1‑251 (MQ=255)

                                                                                                                                                                                                                                            gCGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTc       <  2:425622/249‑1 (MQ=255)

                                                                                                                                                                                                                                             cGATCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcacc    <  1:130583/251‑1 (MQ=255)

                                                                                                                                                                                                                                               aTCGTTTACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcaccac  >  1:382340/1‑251 (MQ=255)

GCAGTGATTTCCGACTGATCGTATGCCGGAGCCACTTCCACTACGTCCATCCCAACAATGTTGAGATCTTTCAGGCCGCGTACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCAC  >  NC_000913/3082941‑3083430

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