BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,621,090	C→T	R633R (CGC→CGT)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,621,090	0	C	T	97.8%	256.1 / NA	143	R633R (CGC→CGT)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand): ref base C (2/0); major base T (75/61); minor base G (1/2); total (80/63)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCG  >  NC_000913/3620843‑3621334

    ggaccctgcccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTCCCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCAGGAAACGCGGTATGAATACAACAGCGCCGTTGACCTGACCGCCGGCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTATTATgg                                                                                                                                                                                                                                               <  2:567537/242‑1 (MQ=11)

       ccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCAGTGGACCGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATAGCCGGTGAACTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGAGGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGGCGACGCGAAGTATGGAAt                                                                                                                                                                                                                                            >  1:345309/1‑251 (MQ=11)

           gaccGCTTCGGGCAGATGACGGCGGTGCACCGCGCGGAAGGGCTGAGTCAGTAGCGCGCATACGACACCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAACCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACg                                                                                                                                                                                                                                         >  1:459886/1‑250 (MQ=9)

           gaccGATTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCCGTACCGCGCATACGCCAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGCCGGCAGCAGAAACGGGACAAAGTACGATGCGAGGGGACAGGCCGTCCGGACCACGCAGGGCGGGCCGACGCGTAGTATGGACTACg                                                                                                                                                                                                                                         >  2:585359/1‑250 (MQ=11)

            gccGCTGCGGGCAGAGGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGGACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCTGGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGCCCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAt                                                                                                                                                                                                                                       <  2:159272/250‑1 (MQ=11)

             ccGCTTCGGGCAGATGCCGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtg                                                                                                                                                                                                                                      <  1:200175/251‑1 (MQ=14)

             ccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtg                                                                                                                                                                                                                                      >  1:110119/1‑251 (MQ=14)

             ccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAAACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtg                                                                                                                                                                                                                                      >  2:516089/1‑251 (MQ=14)

               ccgtCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCGGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGGATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtgct                                                                                                                                                                                                                                    <  2:34063/248‑1 (MQ=9)

                cTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtgctg                                                                                                                                                                                                                                   >  1:138652/1‑251 (MQ=14)

                            gACGGCGGTGCACCCAGACGAAGGGCTGACTCAGTACCGCGCATACGACAGCCGTGGAACGTTCAT‑TCCCGCGCAAGAAACGCAGGGCCATGACACGCGGTATGAATACCACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACCGAAAGATGCCTGGGGAACGGCCGCCCGTACCCCGCAGGGCGGGCTGCCGCGTAGTATGGAATACGATGCTGCCGGACggggcc                                                                                                                                                                                                                       >  2:241894/1‑248 (MQ=11)

                               ggcggTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGGAGTATGGAATACGATGCTGCCGGCCGGGTCATcc                                                                                                                                                                                                                    >  2:7755/1‑251 (MQ=14)

                                  ggTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAGACGGGACACAGTACGATGCGTGGGGAAAGTCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGGATGGAATACGATGCGGCCGGACGGGTCATCCGcc                                                                                                                                                                                                                 >  2:427166/1‑251 (MQ=9)

                                   gTGCACCGCGAGGAAGGGCTGAGTCAGTCCCGCGCATACTACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGTCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGGCCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCt                                                                                                                                                                                                                <  1:131118/251‑1 (MQ=9)

                                   gTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAGTACGATGCTGCCGGACGGGTCATCCGCCt                                                                                                                                                                                                                >  1:207680/1‑251 (MQ=14)

                                     gCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                                                                                                                                                                                                               >  1:120695/1‑250 (MQ=14)

                                      cACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAc                                                                                                                                                                                                             >  2:270150/1‑251 (MQ=14)

                                         cgcgAGTAAGGGCGTAGGCATTACCGGGCATACGACAGCCGTGGAAAGT‑AATGTGCCGTGAAAGACCCGCAGGGCCGTGAAACGCGGTACGAAGACAACATCGCCGGTAGCCTGACCGCCGTCAGTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGGAGTATGGAATACGATGCTGCCGTACGAGTCATCCGCCTGACCAg                                                                                                                                                                                                          <  2:391377/251‑1 (MQ=11)

                                         cgcgAGGAAGGGCTGGGTCAGTCCCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCAGGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAg                                                                                                                                                                                                          <  2:499579/251‑1 (MQ=9)

                                          gcgAGGAAGGGCTGAGTAAGGACCGCGCATACGACAGCCGGGGACAGGTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAATCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGTACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCAGTACCAGGCAGGGCGGGCTGACGCTTAGTATGGAATACGATGCTGTCGGACGGGTCATCCGAATGACCAGt                                                                                                                                                                                                         <  1:396188/251‑1 (MQ=11)

                                                aaGGGCGGAGTAAGTACCGCGCATACGACAGCCGGGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTTACCAGTGAAAAc                                                                                                                                                                                                   <  1:152611/251‑1 (MQ=9)

                                                 aGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg                                                                                                                                                                                                  >  1:293167/1‑251 (MQ=14)

                                                   ggCTGAGTCAGTCCCGCGCATCCGCCAGCCGTGGTCAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGc                                                                                                                                                                                                <  1:96734/251‑1 (MQ=9)

                                                   ggCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAACCGGc                                                                                                                                                                                                >  2:100571/1‑251 (MQ=14)

                                                    gCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCa                                                                                                                                                                                               <  1:223179/251‑1 (MQ=14)

                                                           caggaCCGCGCCAACGACCGCCGTGGACAGTAAAT‑TGCCGGGAAAGACACGCAGGGCCAAGAAACGCGGTATGAATACAACCACGCCGGTGACCGGCCCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCacac                                                                                                                                                                                        <  1:51850/247‑1 (MQ=11)

                                                              ccccGCGCATGCGACAGCCGTGGACAGTTAAT‑GCCCGTGAAAGACACGCGGGGCCATGAAACGCGGTATGAACACAACAACGCCGGTGACCGGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAcaccac                                                                                                                                                                                     <  2:229326/249‑1 (MQ=11)

                                                                ccGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCt                                                                                                                                                                                   >  1:395138/1‑251 (MQ=14)

                                                                   cgcAGACGACAGCCGTGTACAGTTAAT‑GGCCGTGAAAGACACGCAGGGCCAGGAAACGCGTTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTTCCCCGGACGGCAGCAGAAACGGGACACAGTACGAGGCGTGGGGAAAGGCCGGCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTcc                                                                                                                                                                                <  1:241894/251‑1 (MQ=11)

                                                                     cATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGt                                                                                                                                                                              >  1:299172/1‑251 (MQ=14)

                                                                      aTACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATCCAACATCGCCGTTGACCTGACCGCCGTCATTGCCCCGGACGGCATCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGtt                                                                                                                                                                             <  1:567469/251‑1 (MQ=9)

                                                                      aTACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACCCGGTATGAATACAACATCGCCGGTGCCCTGACCGCCGTCATTGCCCCGGACGGCCGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCCCGCAGGGCGGGCTGACGCGTAGACTGGCATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGATAACGGCAGCCCCACCCCCCTCCGtt                                                                                                                                                                             >  2:541133/1‑251 (MQ=11)

                                                                       tACGAAAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACAAGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACCCCACCTTCCGtt                                                                                                                                                                             >  2:386900/1‑250 (MQ=9)

                                                                        aCGACAGCCGTGGACTGTTAAT‑TGCCGGGAAAGACACGCAGGGCCATGAACCGCGGTATGAATACACCATCGCCGGTGACCTGACCGCCGTCAGTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTa                                                                                                                                                                            <  2:15967/250‑1 (MQ=255)

                                                                          gACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGTCAGCCACACCACCTTCCGTTACg                                                                                                                                                                          >  1:322075/1‑250 (MQ=14)

                                                                           aCAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGAt                                                                                                                                                                        <  1:19466/251‑1 (MQ=14)

                                                                                 gTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGGACCACGCAGGGCGGGCTGACGCGTAGGATGGAATACGATGCTGCCGGACGGGTCATCCGCCTAACCAGTGAAAACGGCAGCCACACCACCTTCCGGTACGATGTACtt                                                                                                                                                                  >  2:260241/1‑251 (MQ=9)

                                                                                 gTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGGGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACtt                                                                                                                                                                  <  1:48739/251‑1 (MQ=14)

                                                                                  tggccAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGCATGAATACAACATCGCCGGTGCCCTGACCTCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGCCGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCCGAGAAAACGGCCGCCACACCACCTTCCGTTAAGAGGTACTTg                                                                                                                                                                 >  2:153256/5‑251 (MQ=11)

                                                                                    gACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGAc                                                                                                                                                               >  1:547639/1‑251 (MQ=14)

                                                                                    gACAGTTAAT‑TGCCGGGAAAGACACGCAGGGCAATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGAc                                                                                                                                                               <  2:421758/251‑1 (MQ=14)

                                                                                       aGTTAAT‑TGCCGAGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGCGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACCCAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCCCACCACCTTACGACACGATTAACTTGACCgg                                                                                                                                                            >  2:474760/1‑251 (MQ=11)

                                                                                                 ccGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCCGCCGAAACGGGACCCAGTACGATGCGCGGGGAAAGCCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGCTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGCTGTACTTGACCGGCTGATAc                                                                                                                                                     >  1:60774/1‑249 (MQ=11)

                                                                                                   gTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAgg                                                                                                                                                  >  1:398114/1‑250 (MQ=14)

                                                                                                     gaaatacgcGCAGTGCCATGAAACGCGGTATGATTACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAGCGGCAGCCACACCCCCTTCCGTTACGATGTACTTGACCGGCTGATACAGGaaa                                                                                                                                               <  1:44594/243‑1 (MQ=11)

                                                                                                     gAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCCTCCGTTACGAAGCACTTGCCCGGCTGATAAAGGaa                                                                                                                                                >  2:44797/1‑250 (MQ=255)

                                                                                                      aaaGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGAGACAGGAAAc                                                                                                                                              >  2:180984/1‑251 (MQ=14)

                                                                                                       aaGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCCCCTTCCGCTACGATGTACTTGCCCGGCTGATACAGGAAAcc                                                                                                                                             >  2:472506/1‑251 (MQ=14)

                                                                                                           cccGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCt                                                                                                                                         <  2:538097/249‑1 (MQ=14)

                                                                                                              gCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTg                                                                                                                                      <  1:179240/251‑1 (MQ=14)

                                                                                                                aGGGCCATGAAACGCGGTATGAATACACCATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGAc                                                                                                                                    <  1:4742/251‑1 (MQ=14)

                                                                                                                     caggAAACGCGTGATGAATACAACAGGGCCGGTGACCGGACCGCCGTCAGTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCAGCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCg                                                                                                                               <  2:385285/248‑1 (MQ=11)

                                                                                                                     cATGAAACGCGGTATGGATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGGCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGGACTTGACCGGCTGATACAGGAAACCGGCTTTgccggccg                                                                                                                               >  2:153135/1‑251 (MQ=9)

                                                                                                                      aTGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCTTTGCCCCGGACGGCAGCAGACACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGACGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGCCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTACGTTACGATGACCTTGACCGGCTGATAAAGGAAACCGGCTTTGACGGCCGc                                                                                                                              >  2:301956/1‑251 (MQ=11)

                                                                                                                          aaCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCCGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGCTGTACTTGACCGGCTGATACCGGAAACCGGCTTTGACGGCCgaccac                                                                                                                          >  2:564756/1‑246 (MQ=9)

                                                                                                                           aCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcacaca                                                                                                                         <  1:58200/251‑1 (MQ=14)

                                                                                                                            cgcggaaggAATACAACACCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAg                                                                                                                        <  2:288945/243‑1 (MQ=9)

                                                                                                                               ggTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGCAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCATCCACACCACCTTCAGTTCCGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACCCAGCGt                                                                                                                     >  2:155324/1‑251 (MQ=255)

                                                                                                                                  aTGAAGACACCATCGCCGGTGACCGGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCCCCTTCCGTTACGATGTACTTGACCGGCTGATACAGTAAACCGGCTTTGACGGCCGCACACAGCGTTAt                                                                                                                  <  1:424303/251‑1 (MQ=255)

                                                                                                                                       tACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATcaccc                                                                                                             >  2:514107/1‑250 (MQ=14)

                                                                                                                                        acaacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTTATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATcaccac                                                                                                            <  1:240072/251‑1 (MQ=14)

                                                                                                                                          aacatagacggagaactgaccGCCGTCATTCCCCCGGTCGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGACCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGCAAACGGCAGCCACACCACCTTCCGTTACGATGGACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACCCAGCGTTATcaccccgg                                                                                                          >  2:390589/16‑247 (MQ=11)

                                                                                                                                           acaTCGCCGGTGCCCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGAc                                                                                                         <  2:110119/251‑1 (MQ=14)

                                                                                                                                            caTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGAcc                                                                                                        >  1:52706/1‑251 (MQ=14)

                                                                                                                                            caTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGAcc                                                                                                        >  1:238557/1‑251 (MQ=14)

                                                                                                                                                 ccGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACTATGCTGCCGGACGGGTCATCCGCCTGACCAGTGACAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGATTTGAAGGCCGCACCCAGCGTTATCCCCACGAACTGAcc                                                                                                   >  1:586873/1‑251 (MQ=11)

                                                                                                                                                   ggTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGAc                                                                                                                                    >  1:344262/1‑216 (MQ=14)

                                                                                                                                                   ggTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCGTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGAc                                                                                                                                    <  2:344262/216‑1 (MQ=14)

                                                                                                                                                     tgacctgaccGCCGTCATTGCCCCGGCCGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAAGAAACCGGCGTTGACGGCCGCCCACAGCGTTATCACCCCGACCTGACCGGc                                                                                                >  2:390729/1‑250 (MQ=255)

                                                                                                                                                         ctgaccGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCCCGACCTGACCGGCAAACt                                                                                           >  2:156095/1‑251 (MQ=14)

                                                                                                                                                          tgaccGCCGTCATTGCCCAGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACtt                                                                                          >  1:465152/1‑251 (MQ=14)

                                                                                                                                                            accGCCGGCATGGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTAt                                                                                        <  1:7755/251‑1 (MQ=14)

                                                                                                                                                             ccgccgTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATc                                                                                       <  1:140732/251‑1 (MQ=14)

                                                                                                                                                                  gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGc                                                                                    >  1:221666/1‑249 (MQ=14)

                                                                                                                                                                  gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGc                                                                                    >  1:424568/1‑249 (MQ=14)

                                                                                                                                                                       tGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTGCGACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACGTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAg                                                                              <  2:209150/250‑1 (MQ=9)

                                                                                                                                                                              gACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGCCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCCTCCGCCTGACCAGTGAAAACGGCAGCCACACCACCGTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCGGACCGGCAAACTTTTCCGCAGCGAGGATGAggg                                                                      <  2:537520/251‑1 (MQ=255)

                                                                                                                                                                               ccgtcaccagAAACGGGACACAGTACTATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAGGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGGTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAACCTTATCCGCAGCGAGGATGAGGGt                                                                     <  2:340122/244‑1 (MQ=11)

                                                                                                                                                                                 ggCAGCAGAAACGGGACACAGTACGAGGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCt                                                                   <  1:189052/251‑1 (MQ=14)

                                                                                                                                                                                  gcagcaGAAACGGGACACAGTACGATGCGTGGGGAAAGGACGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCTAGGATGAGGGTCTg                                                                  <  1:445791/251‑1 (MQ=14)

                                                                                                                                                                                          aaCGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCa                                                          <  2:299172/251‑1 (MQ=14)

                                                                                                                                                                                            cGGGACCCAGTACGATGCGTGAGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCCCACAGCGTTATCACCACGCCCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGGCACCCAcc                                                        >  2:554826/1‑250 (MQ=255)

                                                                                                                                                                                                acacAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGc                                                     >  1:36338/1‑250 (MQ=14)

                                                                                                                                                                                                    attACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTAt                                                >  2:251319/3‑251 (MQ=14)

                                                                                                                                                                                                      tACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATg                                               >  1:493574/1‑250 (MQ=14)

                                                                                                                                                                                                              gTGGGGCAAGGCCGGCCGTACCACGCAGGGCGGACTGACGCGCAGTATGGCATACGATGCTGCCGGACGGGTCCTCCGCCTGCCCAATGAAACCGGCAGCCACACCCCCTTCAGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGATTCGACGGCCGCACAAAGCGTTCTCACCAACACCTGACCGGCAAACTTATCCGCAGCGATGATGAGGATATGGTCCCCCCCTGGCCCTATGCCGACGCa                                       >  2:269429/1‑250 (MQ=11)

                                                                                                                                                                                                              gTGGGGACAGGCCGTCCGGACCACGCAGGGCGGGCTGACGCGTAGTATGGGATACGATGCTGCCGGAAGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGGACAACGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGACGCTGAGGGTCTGGTCACCCACTGGCACTATGACGAAga                                        >  2:356686/1‑248 (MQ=11)

                                                                                                                                                                                                                   gAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGc                                 <  1:270150/251‑1 (MQ=14)

                                                                                                                                                                                                                    aaaGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCAGCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGcc                                <  2:276801/251‑1 (MQ=14)

                                                                                                                                                                                                                          cgtccgtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgc                          >  1:359628/1‑251 (MQ=14)

                                                                                                                                                                                                                              cgtACCACACAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTACCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgcaccg                      >  1:277458/1‑251 (MQ=14)

                                                                                                                                                                                                                                tACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgcaccgc                     >  1:56481/1‑250 (MQ=14)

                                                                                                                                                                                                                                 aCCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGTACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgcaccgca                    <  2:522844/250‑1 (MQ=14)

                                                                                                                                                                                                                                    aCGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGcaccgcaccgg                >  1:295544/1‑250 (MQ=14)

                                                                                                                                                                                                                                     cGCAGTGCGGGCTGACGCGTAGTATGGAATACGAGGCTGCCGGACGGGTCCGCCGCCTGATCAGTGAAAACGGCAGCCACACCACCTTCCGTGACGATGGACTTGACCGGCTGATACAGGAAGCCGGCGTTGACGGCCGCACACAGCGTTAGCACCACGACCTGACCGGCAAACGTATCCGCAGCGAGGATGAGGGTCGGGTCACCCTCTTGCACTATGACGAATCAGACGGCCTCAcgagcctcacagtg               <  2:67247/251‑13 (MQ=11)

                                                                                                                                                                                                                                        aGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAg            <  1:156095/251‑1 (MQ=14)

                                                                                                                                                                                                                                          ggcgggcTGACGCGGAGTATGGAATGCGATGCGGCCTGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAgtg          <  2:465152/251‑3 (MQ=9)

                                                                                                                                                                                                                                                 tGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGCAGCAGACCGCCTCACGCACCGCACCGGGAAGGGTGAAAcc   >  2:432148/1‑251 (MQ=11)

                                                                                                                                                                                                                                                  gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACCCCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATAACCACGAACTGACCGGCAAACTTATCCGCAGCGAGGGTGCGGGTCTGGTCCCCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGACGGGTGCAACCg  >  2:373654/1‑251 (MQ=11)

                                                                                                                                                                                                                                                  gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTAACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAAcc   <  1:260241/250‑1 (MQ=11)

GTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAAT‑TGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCG  >  NC_000913/3620843‑3621334

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