BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,394	A→G	S522G (AGC→GGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,394	0	A	G	99.1%	346.6 / NA	115	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand): ref base A (0/1); major base G (42/71); minor base T (0/1); total (42/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 116 total aligned reads displayed.

AACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACA  >  NC_000913/1337148‑1337639

aaCCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCg                                                                                                                                                                                                                                                       >  1:838048/1‑250 (MQ=255)

         gtgCTGATGGCCTCAGGCTTGCTGGCGAAAAAAGCCGTAACTCGGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCg                                                                                                                                                                                                                                              <  2:455613/250‑1 (MQ=255)

          tgCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                                                                            >  1:107409/1‑251 (MQ=255)

          agCTGATGGCCGCAGGCTTGCTGGCGAAAAACGACGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATATCGACGAACTGGGGTTTCACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCCGGGCCGATGCCCGTTCCTATCGAACCGGCAAACAAAAAAGGCGCATTAACCGTc                                                                                                                                                                                                                                            >  2:318749/2‑251 (MQ=255)

           gCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTc                                                                                                                                                                                                                                            <  2:897731/250‑1 (MQ=255)

                    cGCAGGCTTGCTGTCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTACGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCggtgcggtgc                                                                                                                                                                                                                                  <  2:972784/251‑1 (MQ=255)

                                  cGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGt                                                                                                                                                                                                                    >  1:514949/1‑251 (MQ=255)

                                  cGAAAAAAGCCGTAACGCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAATTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCGTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTCCGGCAACCGt                                                                                                                                                                                                                    <  2:640169/251‑1 (MQ=255)

                                   gAAAAAAGCCGTAACTCTGGGCCTCACGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGATACGCTGCCCGATCCTATCAAAACGGCCATCAAAACAGGTGACTTAACCGTCGGTTAGGTGCttgcaggacaccgac                                                                                                                                                                                                                   >  2:148226/1‑237 (MQ=255)

                                    aaaaaaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGACCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCACTCAAACAAGGCGATTTAACCGTCGGTGCGGGGCTGTCCGGCAACCGTa                                                                                                                                                                                                                   >  2:324468/1‑250 (MQ=255)

                                    aaaaaaGCCGTAACGCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTTGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACTCCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTaa                                                                                                                                                                                                                  <  1:951148/251‑1 (MQ=255)

                                       aaaGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACtt                                                                                                                                                                                                               <  1:667711/251‑1 (MQ=255)

                                         ggcggTAACTCTGGGCCTTTAGCAGCAACCATAGGTCAAAGCCCCGCTGGCCCCGGGTTCGAAAGACGTTTCTGATTATCTTGCAAAAGTGAAACTGACACCGTAGTTCGACGAACTGGGGTTTAACCTTGTGGGAGACGGTTGTACCCCC‑TGT‑ATTGGTAACTCTGGGCCGGTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATCTAACCGTCGGTGCGGTGCTGTCCGACAAACGTAACTTTg                                                                                                                                                                                                             <  2:823169/247‑1 (MQ=255)

                                            cGTATCTCTGGTCCTCAAGCGGCAACCATGGGTAAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTTATTATCTGGCAAAAGCGAAACTGACACCTGATCTCGACGAACTGTGGTTTAACCTTGTGGGATACGGTTGTACCACC‑GGT‑ATTGGTAACTCTGGGCCGCTGCCCGAGCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGCCCGGCACCCTTAACTTTGaat                                                                                                                                                                                                          <  2:909546/251‑2 (MQ=255)

                                            cGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAg                                                                                                                                                                                                          <  1:53776/251‑1 (MQ=255)

                                             gTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAgg                                                                                                                                                                                                         >  1:432875/1‑251 (MQ=255)

                                               aaCTCGGGGCCCCAAGCTGCAACCATGGGTCACAGCGTCGTTGGCACCGGGTTCGAAAGTCGTTTTTGATTATCGGGCAAAAGCGAAACTGACATCGTATCTCGACGAACTGGGTTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGcc                                                                                                                                                                                                       <  2:248996/251‑1 (MQ=255)

                                                   ctGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTTGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTat                                                                                                                                                                                                   >  1:523483/1‑251 (MQ=255)

                                                     gggCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatc                                                                                                                                                                                                  <  1:968637/250‑1 (MQ=255)

                                                      ggCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGTTTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTTGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGTCAACCGTAACTTTGAAGGCCGTatcca                                                                                                                                                                                                <  2:566740/251‑1 (MQ=255)

                                                      ggCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatcca                                                                                                                                                                                                <  1:113993/251‑1 (MQ=255)

                                                       gcatCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACAGGGTTCGAAAGTCGTTTCTGATTATCTGGCAACAGCGAACCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGCAACGGCAATCAAAAAAGGAGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGGatccc                                                                                                                                                                                                >  2:626519/4‑249 (MQ=255)

                                                       cccTCATGCGGCCTCCATGGGTCAAAGCTTCGCTGGCCCCTTGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAATCTGACACCGTATCTCGACGAACTGGTGTTTACCCTTGTGGGTTACTGTTGTACCTCC‑TGT‑ATTTGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAATGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccat                                                                                                                                                                                               <  1:52264/250‑1 (MQ=255)

                                                       accTCAAGCTGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGCATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGGCCGGCAACCGAAACTTTGACGGCCGTatccat                                                                                                                                                                                               >  2:53217/2‑251 (MQ=255)

                                                        ccTCAAGAGGCAACCATGGGTCAAAGCGTCGCTGGCCCCGGGTTCGAAAGTCGTTTCTGCTTATCTGGCAAAAGCGACACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATAAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTACATTTGAAGGCCGTatccatc                                                                                                                                                                                              >  2:110885/1‑251 (MQ=255)

                                                               cGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGACGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGCTTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGt                                                                                                                                                                                       >  1:858657/1‑251 (MQ=255)

                                                                ggCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGtt                                                                                                                                                                                      <  2:710415/251‑1 (MQ=255)

                                                                       aggggTCAAAGCGATGCTGGCACCGGGTTCGAAAGTCGTTTCTGAGTAGCTGGCAAAAGCGAAACGGACACCGTATCTCGACGAACTGGGGTGTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGCCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAaacta                                                                                                                                                                               <  2:293962/249‑1 (MQ=255)

                                                                         gggTCAAAGCGTCGGTGGCACAGGGTGCGATAGTCGGTTCTGAGTATTTGGCACAAGCGAAACGGTCTCCGTAGCGTGACGAACTGGGTGTTAACCTTGTGGGACACGGTTGTACCACC‑GGT‑TTTGGGAACTCTGGGCTTCTGCCCGAACCTATCGAAACGGCTATCAAAAAGGGCGATCTAACCGTCGGGGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGGATCCATCCGCTGGTTAAaactaac                                                                                                                                                                             <  2:781961/251‑1 (MQ=255)

                                                                             cAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTGGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggc                                                                                                                                                                         <  2:971463/251‑1 (MQ=255)

                                                                               aaGCGTCGCTGGCACCGGTTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCTAAACTGACCCCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATCCGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAActggctg                                                                                                                                                                       <  2:792766/251‑1 (MQ=255)

                                                                                aGCGTCGCTGGCACCGGGTTCGAAAGTCGTTACTGATTCTCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGA‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCACCCGAAACTTTGAAGGCCGTATCCTTCCGCTGGTTAAAACTAAatggctgg                                                                                                                                                                      >  2:750100/1‑251 (MQ=255)

                                                                                 gCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACtggctggc                                                                                                                                                                     <  1:274470/251‑1 (MQ=255)

                                                                                  cgccGCTGTCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCGATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGcc                                                                                                                                                                    <  2:712284/248‑1 (MQ=255)

                                                                                   gtctcTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCCCC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCt                                                                                                                                                                   <  1:559128/247‑1 (MQ=255)

                                                                                     cGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcg                                                                                                                                                                 <  1:112855/251‑1 (MQ=255)

                                                                                        tGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTCTCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGaa                                                                                                                                                                                                           >  1:793951/1‑206 (MQ=255)

                                                                                         ggccccGGGTTCTAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTAGCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTACCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTcgcc                                                                                                                                                               <  2:771698/245‑1 (MQ=255)

                                                                                             ccGGGTTCGAAATCCGTTTCTGATTATCTGGCAAAAGCAAAACTGACACCGTATCTCGACGAACGGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAAACAAAAAAGGCGATTTACCCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCt                                                                                                                                                         <  2:424909/251‑1 (MQ=255)

                                                                                             ccGGGTTCGAAAGCCGTTGCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTGTAAACTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCGGCCCGATCCTATCGAAACGGCAATCAAAAAAGTCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCACCCGCTGGTTAAAACTAACGGGCTGGCCTcgccgccg                                                                                                                                                           <  2:263720/249‑1 (MQ=255)

                                                                                                 gTTCGAAAGTCGTTTCTGCTTCTCTGGCAAAAGCGAAACTGACACCGTATCTCGCCGAACTGGGGTTTAACCTTGTGGGCTACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTCTCGACACGGCACTCAAAAACGGCGCTTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTTACCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                     >  2:958088/1‑251 (MQ=255)

                                                                                                 gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                     <  2:493495/251‑1 (MQ=255)

                                                                                                 gTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCtggtg                                                                                                                                                     >  1:441009/1‑251 (MQ=255)

                                                                                                     gaaatgtggttatggttaacttggcaaaGCGGAACTGAAACCGTATCCCCGCGAACTGGGGTTTAACCTTGTGGGGTTCGGTTGGTCCACCCTG‑‑AATGGGAACTCTGGGCCGCGTCCCGGGCCTAACCAAAACGCAAATAAAAAAAGCGATTTAACCGGCGGTGGGGTTCTTGCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTg                                                                                                                                                 <  2:841162/226‑1 (MQ=255)

                                                                                                      aaaGTCGTTTCAGCTTATCTGGAAAAAGCGAAACGGACACCGTATCTTGACGAACTGGGGTTTAACCTTGGGGGATACGGTTGTTCCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTACCGAACCGGCAATCAAAAAAGGCGATTTAACAGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGc                                                                                                                                                <  2:310457/251‑1 (MQ=255)

                                                                                                       aaGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGGGGTTGcc                                                                                                                                               <  1:784306/251‑1 (MQ=255)

                                                                                                         gTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTa                                                                                                                                             <  1:675689/251‑1 (MQ=255)

                                                                                                          acGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGATCTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCGCGCCGCCGCTGGTGGTTGCCTAt                                                                                                                                            >  2:474653/2‑251 (MQ=255)

                                                                                                                  gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTACCGCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTgg                                                                                                                                     <  1:540901/250‑1 (MQ=255)

                                                                                                                      aTCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCggg                                                                                                                                 >  1:583761/1‑250 (MQ=255)

                                                                                                                               aaGCGAAACTGACACCGTATCTCGTCGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAat                                                                                                                      <  1:291648/251‑1 (MQ=255)

                                                                                                                                   gAAACTGACATCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGGTGTACCACC‑TGT‑ATTGGTAACTCGGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCTCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑Ca                                                                                                                 <  2:90420/251‑1 (MQ=255)

                                                                                                                                   gAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑c                                                                                                                  >  1:68869/1‑250 (MQ=255)

                                                                                                                                    aaaCTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACA‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCACCCGTAACTTTAAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCAGCCGCTGGTGGTTGCCTATGAGCTGGCgggtactctcaatttca                                                                                                                   >  2:231596/1‑236 (MQ=255)

                                                                                                                                        tGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑GAATATGAATAT‑CAACCTg                                                                                                            >  1:595264/1‑251 (MQ=255)

                                                                                                                                         gagacCGTATCGCGACGAAGGGGGGTTTAACCTTGGGGGATACGGTTGGACCACG‑GGT‑ATTGGTAACTCGGGGCCGCGGCCGGAGCCTATCGAAAGGGCAATCAAAAAAGGGGATTGAACCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTAACGGGCTGGCCGCGCCGCCGCTGGTGGGTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTg                                                                                                            <  2:913095/247‑1 (MQ=255)

                                                                                                                                            cccgtattccgacgaACTGGGGTTTAACCTTTTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGTCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCtt                                                                                                        <  1:888528/242‑1 (MQ=255)

                                                                                                                                                   ctcGACGAACTGGGGTTTCACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTCAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCCTCTGaacc                                                                                                 >  2:362584/1‑248 (MQ=255)

                                                                                                                                                   ctcGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGcc                                                                                                 >  1:379042/1‑251 (MQ=255)

                                                                                                                                                    ccgtcgaACCGGGGTTTAACCTTGGGGGCTAAGGTTGTTCCCCC‑TGT‑ATTGGTAACTCTGTGCCGCTGCCCGAGCCTACCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTTCTGACCGGCAACCTTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGT‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCt                                                                                                <  2:70498/247‑1 (MQ=255)

                                                                                                                                                     cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAg                                                                                                   <  1:120579/247‑1 (MQ=255)

                                                                                                                                                     cgacgaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGc                                                                                                  <  1:777635/248‑1 (MQ=255)

                                                                                                                                                     cgaagaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGGGCTGTCCGGCAACCGTAACGTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGc                                                                                                  <  2:665354/244‑1 (MQ=255)

                                                                                                                                                       acgaACGGGGGTTTAAACGTGTGGGATACGGTTGTACCCAC‑TGT‑ATTGGTAACACCGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTTCGGCACCCGTAACTTTGAAGGCCGTAGCCATCCGCTGGTTAAAACGAACGTGATGGCCTGGCCGCTGCTGGGGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATc                                                                                             <  2:499721/251‑1 (MQ=255)

                                                                                                                                                         gaACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCgg                                                                                           >  1:600465/1‑251 (MQ=255)

                                                                                                                                                               gttgttAATCTTGTGGGGTCCGGTTGTTCCACC‑‑GTGATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAAGCAAAAAAGGCGATTTAACCGTCGGTGCGTTGCTGTCCGGCAACCGTAACTTTTAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATTCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGa                                                                                     <  2:360219/247‑1 (MQ=255)

                                                                                                                                                               gggTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGa                                                                                     <  1:953774/251‑1 (MQ=255)

                                                                                                                                                                gtgttAACCTTGTGGGATACGGTTGTTCCACC‑TGT‑AATGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAAGCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGAt                                                                                    <  2:878399/248‑1 (MQ=255)

                                                                                                                                                                ggTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGTTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAATAGGCGATTTAACCGCCGGGGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGAt                                                                                    <  2:428454/251‑1 (MQ=255)

                                                                                                                                                                ggTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGAt                                                                                    <  1:860662/251‑1 (MQ=255)

                                                                                                                                                                      aCCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTTAGCCTATCGGCCATGTTCGCaa                                                                               >  2:585691/1‑250 (MQ=255)

                                                                                                                                                                       ccTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAg                                                                             <  2:838048/251‑1 (MQ=255)

                                                                                                                                                                           gtgGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCg                                                                          >  1:90185/1‑250 (MQ=255)

                                                                                                                                                                            gggggtcagttgtttctcac‑TGT‑TTTGGTAATTCTGGGCCGCAGCCTGATCCTAACGATACGGCAATTAAAAAAGGCGATTCAACCTTCTGTTCGTTGCTGTCCGGCAACCGTAACTTTGAAGTCCGTATCCATCCGCTGGTGAAAACTAACTGGCTGGCCTCGCCGCCGCTGTTGGTTGCCTGTGCGCTGGCGTG‑AAATATGAAGAT‑CAACCTGTCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGAt                                                                        <  2:804876/232‑1 (MQ=255)

                                                                                                                                                                               gATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCg                                                                     <  1:270358/251‑1 (MQ=255)

                                                                                                                                                                                     gTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAACCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGGAAAT‑TGAATAT‑CAACCTGGCTTCTGAGCCTAACGGCCATGATCGCAAAAGCGGTCCGGTTTAt                                                               >  2:376617/1‑251 (MQ=255)

                                                                                                                                                                                     gTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGGGGTTGCCTATGCGCTGGCGGG‑AAAAATGAATAT‑CAACCTGGCTTCTGAGCCCACCGGCCaagctagcaaacgcagtccggttttt                                                               >  2:738285/1‑226 (MQ=255)

                                                                                                                                                                                       tGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCACTGGCGGG‑AAACATTAATAT‑CAACCTTGCTTCTGAGCCTATCGGCCATGATCGCAAAAACGATTCGGTTTATCt                                                             >  2:593855/1‑251 (MQ=255)

                                                                                                                                                                                            cacc‑TGT‑ATTGGTACCTCTGGGCCGTGGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAg                                                        <  2:758333/251‑1 (MQ=255)

                                                                                                                                                                                             ccc‑TGT‑AGTGGTACCTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCACCCGTAATGTTGAAGGACGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATAGGAATAT‑CAACCTGGCTTCTGAGCCGATCGGCCATGAGCGCAAAGGCGTTCCGGTTTATCTGAAAGa                                                       <  2:432875/250‑1 (MQ=255)

                                                                                                                                                                                                 tGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGatat                                                    >  1:292826/1‑251 (MQ=255)

                                                                                                                                                                                                 cgggaTTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTTTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGata                                                     <  2:799348/247‑1 (MQ=255)

                                                                                                                                                                                                      ttGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAAACTGAATATACAACC‑GGCTTCTGAGCCTATCGGCCATGATCGCAAAAGCGATCCAGTTTATCTGAAAGATAT‑ATgg                                               >  2:594779/1‑251 (MQ=255)

                                                                                                                                                                                                          tAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCAt                                           >  1:605849/1‑251 (MQ=255)

                                                                                                                                                                                                           acctctGGGCCGGTGCCCGTGCCTCTCTAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTATCTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACGGGCTGGCCGCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTGCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAATATAT‑CTGGCCATc                                          <  2:673688/249‑1 (MQ=255)

                                                                                                                                                                                                                 gggCCGCTGCCCGATCCTACCGAACCGTCAACCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGTCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTGGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGca                                      <  2:514949/249‑1 (MQ=255)

                                                                                                                                                                                                                     cGCTGCCCGATCCGAGCGAAACGGCAATCAAAAAAGGCGATTTAACCGCCGGTGCGGTGCTGTCCGGCAACCGTAAATTGGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCGGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGaa                                <  2:416486/251‑1 (MQ=255)

                                                                                                                                                                                                                      gCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCCCAAGaa                                >  1:891128/1‑250 (MQ=255)

                                                                                                                                                                                                                       cTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAAt                              >  1:128807/1‑251 (MQ=255)

                                                                                                                                                                                                                         gCCCGATCCTATCGCAACGGCAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTTCTGTCCGGCAACCGTATCTGTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGTCTCGCCGCCGCTGGTGGTTGCCTATGCGCGGGCGTG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCTGTTTAGCTGAAAGATAT‑CTGGCCATCGTCACAACAAAt                              <  1:248781/249‑1 (MQ=255)

                                                                                                                                                                                                                         gCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGTTCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAATAAAtt                             <  1:275462/250‑1 (MQ=255)

                                                                                                                                                                                                                         acccGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAAAATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTTTCTGAAAGA‑ATACTGGCCCTCGGACCAAGAAAtgt                            >  2:665989/2‑249 (MQ=255)

                                                                                                                                                                                                                             gATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCg                        <  1:376617/251‑1 (MQ=255)

                                                                                                                                                                                                                                ccTATCGAAACGGCAAACAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCGTg                      <  1:952027/250‑1 (MQ=255)

                                                                                                                                                                                                                                        aaCGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAAc             <  1:377394/251‑1 (MQ=255)

                                                                                                                                                                                                                                           ggCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAg          <  1:466311/251‑1 (MQ=255)

                                                                                                                                                                                                                                                tCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATTAATAT‑CAACCTTGCTTCTGAGCCTATCGGCCATTATCGCAAAGGCCATCCGGTTTTTCTTATAAATAT‑CTTGCCCTCCGCCCCAGTAATTGCCCCTGCGGTAGAACACGTCTcc     >  2:392388/1‑251 (MQ=255)

                                                                                                                                                                                                                                                 cAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGACCAAGTCTcc     <  1:204573/250‑1 (MQ=255)

                                                                                                                                                                                                                                                 cAAAAAAGGCGATTTAACCGTCGGGGCGGGGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTcc     <  1:174818/250‑1 (MQ=255)

                                                                                                                                                                                                                                                   aaaaagGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTCTGCGCTGGCGGG‑CAATATGAATAT‑CAACATGGCTTCTGAGCCTATCGGACATGACACCAAAGGCGATCCGGTTTATCTGGAAATTAT‑CTGTCCATCGGCCCAAGAAATTGCCCGTGGGGGAGAAAAAGTCTccccc  >  2:404397/1‑248 (MQ=255)

AACCCAAGTGTGCTGATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACC‑TGT‑ATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGG‑AAATATGAATAT‑CAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATAT‑CTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCTCCACA  >  NC_000913/1337148‑1337639

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