BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	66,946	T→C	Q464Q (CAA→CAG)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	66,946	0	T	C	96.6%	423.4 / NA	154	Q464Q (CAA→CAG)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand): ref base T (0/4); major base C (55/89); minor base A (5/0); total (60/94)
Fisher's exact test for biased strand distribution p-value = 9.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

GATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAG  >  NC_000913/66700‑67189

gATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTACTACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCATTCACCGTGATTTCAATGTCGTGCATCTCGGTGAACTGgc                                                                                                                                                                                                                                                  >  1:515350/1‑251 (MQ=25)

 aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                                                                                                                                                                                 <  1:621931/251‑1 (MQ=25)

 aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                                                                                                                                                                                 <  1:47930/251‑1 (MQ=25)

 aTAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGgcg                                                                                                                                                                                                                                                 <  1:459363/251‑1 (MQ=25)

          gcCAGCGTCGCATCAGGCGTTGAATGCCGGCTGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGAGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCCGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTACTTTTCATGGCGTGCATATCGGCGGAATGGCGCATCTCGtt                                                                                                                                                                                                                                        >  1:756626/1‑251 (MQ=18)

          gcCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGtt                                                                                                                                                                                                                                        <  1:693596/251‑1 (MQ=37)

             aGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGGTCATCTCGGCGAACTGGCGCATATccttt                                                                                                                                                                                                                                       >  1:300319/1‑245 (MQ=25)

               cGTCGCATCTGGCGTTGAAGTCCGGATGCGCTTTGCGTTTCCGGTCTACCAAAGCGCAGCGTGTAGGGCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACATCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCTCAGAGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGTGATTTCAATGTCGGGCATCTCGGCGAACTGGCGCATATCGTTGAGGt                                                                                                                                                                                                                                   <  2:208428/251‑1 (MQ=255)

                gTCGCATCAGGCGTTGAATGCCGGATGCGCTTTTCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCAACGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCCAATCGTTGAGGtt                                                                                                                                                                                                                                  >  1:217463/1‑251 (MQ=25)

                    cATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATAAATCACAGTGGTTTCAATGTAGTGCATCTCGGCGAAATGGGGCATATCGGTGAGGttaagg                                                                                                                                                                                                                              >  1:844684/1‑247 (MQ=21)

                       cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCa                                                                                                                                                                                                                           <  1:929328/251‑1 (MQ=37)

                       cAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCCCCGGGATTTCAATGTCGTGGATCTCCGCGAACTGGCGCATATCGTTGAGGTTCAgcgc                                                                                                                                                                                                                            >  1:429794/1‑247 (MQ=25)

                            gTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCCATCACCGTGATTTCAATGTCGTACATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAgcgccttgct                                                                                                                                                                                                                      >  1:599797/1‑242 (MQ=37)

                                 aTGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGCTGa                                                                                                                                                                                                                    >  1:461821/1‑248 (MQ=37)

                                   gCCGGATGTGCTTTGCTTTTCCGGACTACAAAATCGCAGCGTGTAGGCTTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCGGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACg                                                                                                                                                                                                               <  2:809136/251‑1 (MQ=255)

                                     cggatggtcTTTGATTATCCGGCCTACAAAAACGCAGCGGGTAGGCCCGATAAGACGCGGCAGCGTCGCATCAGTCGTTACACATCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAAGCACCGTGATTTTAATGTCGTGCATCTCGGCGAACTGTCGCATATCGTTGAGGTTCAGCGCATGGCTGAATACGGt                                                                                                                                                                                                             <  2:687803/243‑1 (MQ=255)

                                      ggATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTa                                                                                                                                                                                                            <  1:173395/251‑1 (MQ=38)

                                        attagcTTTGGTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTCCTTAGCGACGAAATTCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATg                                                                                                                                                                                                          <  2:327825/247‑1 (MQ=37)

                                          gcgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATgg                                                                                                                                                                                                         <  1:137351/250‑1 (MQ=38)

                                           cgcTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTGCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTg                                                                                                                                                                                                       <  1:218024/251‑1 (MQ=38)

                                             cTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGGTATCAATCACCGTGGTTTCCATGTCGTGCATCTCGGCGAACTGGCGCGTATCGTTGTGGTTCAGCGCATGGCTGAAGACGGTATTGTgcg                                                                                                                                                                                                     >  1:24781/1‑251 (MQ=37)

                                                         ggCCTAAAAAATCGTAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTGACACACCTGATGCGACTACTCAGCGACGAATTCCGTAATACACTTCGTTCCA‑GCGCATCGCGTCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGa                                                                                                                                                                                         <  2:285197/251‑1 (MQ=255)

                                                          gCCTACAAAATCGCAGAGTGTAGGCCTGATAAGACGCGTAAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGCATCGTTTTCAATCACCGTGATTTCAAATTCGTGCATCTCGGCGAACTGGGGCGTAACGTTGAGGTTCCACGCACGGCTTGAGACGGTATGGTGCGCGccccccgaggg                                                                                                                                                                                        >  2:595543/1‑242 (MQ=255)

                                                           ccTACAAAATCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCCCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCTCATATCGTTGAGGGTCCGAGCATGGCTGCAGACGGTATGGTGCGCGCCCCCCGCGAgg                                                                                                                                                                                       >  2:347618/1‑251 (MQ=255)

                                                              aCAAAAGCGGAGCGTGTAGGCCTGATAAGCCGCGTCAGCAGAGCAGCAGGCGTTACACAACGGAGGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTGCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTGATCAATCACCGTGATTTCAATGTCGTGCATCGCGGCGAACTGGCGAATATCGTTGAGGTTCAGCGCAGGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATc                                                                                                                                                                                    <  2:319158/251‑1 (MQ=255)

                                                                aaaaACGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCa                                                                                                                                                                                  <  2:135556/251‑1 (MQ=255)

                                                                 aaaTCGCCGCGTGTAGGCCTGATAAGACGCGACAGCGTCGCAACAGGCGTTACGCACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑TCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGGTGTCAATCACCGTGATTGCAAGGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGATGATCCAc                                                                                                                                                                                 <  2:205031/251‑1 (MQ=255)

                                                                 aaaTCGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGCTGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTTTCAATCACCGAGATTTATATGTCGGGCATCTCTGCGAACTGGCGCCTATCGTTGGGGGTCAGCGCATGGCTTGAGACGGGATGGCGCGCGccccccgcgcgagtacac                                                                                                                                                                                 >  2:758493/1‑235 (MQ=255)

                                                                     cGCAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACCCTTCGTTCCA‑GCGCATCGCGTCTATAAACGCTGGCAGGCTGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGATGTTCAGCGCATGGCTGAAGACGGTTTGGTGCGCGCCACCAGCGAGGATCCACGCtt                                                                                                                                                                             <  1:584342/251‑1 (MQ=255)

                                                                      gCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTc                                                                                                                                                                            <  1:496030/251‑1 (MQ=255)

                                                                       cAGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCg                                                                                                                                                                           <  1:6759/251‑1 (MQ=255)

                                                                        aGCGTGTATGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGACATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTATAAATAACCGTGATTTCAATGTCGGGCAACCCGGCGAACATGCGCCTATTGTTGAGGTTCAGCGACTGGCCGAAGAAGGACTGGTTCGCGCCACCAGCCAGGacacacgccttcg                                                                                                                                                                          >  2:167525/1‑239 (MQ=255)

                                                                        aGCGTGTAGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAGCGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGAGCCTCGCTTCgg                                                                                                                                                                          <  2:708176/251‑1 (MQ=255)

                                                                               aGGCCTGATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACGTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGt                                                                                                                                                                  <  2:632408/251‑1 (MQ=255)

                                                                                  cctgATAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATAGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTgg                                                                                                                                                               <  1:241907/248‑1 (MQ=255)

                                                                                      aTAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGCACTGGCGCATATCGTTGAGGTTAAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCCGCGAGGATCCACGCTTCGGAA‑GCAGTTGGc                                                                                                                                                              >  1:629216/1‑248 (MQ=255)

                                                                                      aTAAGACGCGTCAGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGc                                                                                                                                                              >  1:702692/1‑248 (MQ=255)

                                                                                                  aGCGTCGCATCAGGCGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCt                                                                                                                                               <  1:269609/251‑1 (MQ=255)

                                                                                                               gCGTTACACCCCTGATGCGTCCACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCAt                                                                                                                                  <  2:519981/251‑1 (MQ=255)

                                                                                                               gCGTTACACACCTGATGCGACTTCTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGCCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCCATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGGTCAGCGCATGGCTTAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTTGCAGATCCGGTTGCGCTTTCCACAGCGCAt                                                                                                                                  <  2:104429/251‑1 (MQ=255)

                                                                                                                cGTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGAGCATGGCTGAAGACGGTATGGGGCGCGGCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCAtt                                                                                                                                 >  1:938169/1‑251 (MQ=255)

                                                                                                                 gTTACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAAATGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGGATTGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTc                                                                                                                                >  1:868891/1‑251 (MQ=255)

                                                                                                                  ttACACACCTGATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTTTCGTTATCAATCACCGTGATTTCAATGTCTTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCCACGCCTGGCTGAAGACGGTATGGTGCGCGCCCCCCGCGGGGATCCACGCTTCGGAA‑GCCGTTTGCAGATTCGGTTGCGCCTTCCACAGCGCATTct                                                                                                                               >  1:349276/1‑250 (MQ=255)

                                                                                                                    acacacctgATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGcc                                                                                                                             <  1:921060/243‑1 (MQ=255)

                                                                                                                            gATGCGACTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCACCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGTATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAg                                                                                                                     <  2:406851/251‑1 (MQ=255)

                                                                                                                              ggcgacttactAGCGACGAAATCCCTAATACACTTTGGTCCAAGC‑CAGGGCGGCTTTAAACGCTGGCAGGCGGGTAGCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTGCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGtt                                                                                                                   <  2:954197/241‑1 (MQ=255)

                                                                                                                                cgacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGTGGTTCAGCGCATGGCTGGAGGCGGGATGGTGCGCGCCCCCAGCGAGGATCCACGCCTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCCTTTGCCACCGGCAGtttt                                                                                                                 >  2:255382/4‑250 (MQ=255)

                                                                                                                                cgacTACTGAGCGACTAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTAGCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTc                                                                                                                 <  2:611681/248‑1 (MQ=255)

                                                                                                                                 gacTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTCTCGTTATCCATCACCGGGATTTCAATGTCGAGCATCCCGGCGAACTGGCGCATCTCGGTGAGGGTAAGCGCATGGCTGAAGACGGTATGGTGCGCGCCCCCAACGGGGGTCCACCCTTCGGAACGGAGT‑GGCAGGTCCGGGTGGGCTTTCCACCGCGCCATCGCCACCAGCaacttcg                                                                                                                >  2:486716/3‑245 (MQ=255)

                                                                                                                                  acTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCgg                                                                                                               <  1:408356/250‑1 (MQ=255)

                                                                                                                                  acTACTTAGCGACGAAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGGATCGTTTTCAATCACCGGGGTTTCCATGTCGTGCCTCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGCAGACGGTATGGTGTCCGCCACCGGCGGAGTGCCCCGCCTCGGGA‑GGAGTTGGCAGAGCCCGTTGCGCTTTCCACCGCCCATTCTCCCCCGGccgattccc                                                                                                               >  2:274558/2‑243 (MQ=255)

                                                                                                                                         aGCGACGAAATCCTGAATACACGTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTTGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGa                                                                                                         <  2:296435/250‑1 (MQ=255)

                                                                                                                                               gaAATCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgt                                                                                                  <  1:420463/251‑1 (MQ=255)

                                                                                                                                               gaAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgt                                                                                                  <  1:36424/251‑1 (MQ=255)

                                                                                                                                                aaatccGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCgggg                                                                                                 >  1:172407/5‑249 (MQ=255)

                                                                                                                                                aaatccGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGgtg                                                                                                 <  1:708766/247‑1 (MQ=255)

                                                                                                                                                  atccGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGGTTTTAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACAAGCGAGGATCCACGCTTCGGAC‑GCAGTTGGAAGATCCGGTTGCGCTTTTCACAGCGCATTCGCCACCGGCAGATTAGGCAGGGAggtgggtggt                                                                                               >  1:361458/3‑242 (MQ=255)

                                                                                                                                                          tACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGGA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCCGTTTCGGCCGGGGGTGCGGTGTTTTCACCgg                                                                                       >  1:332450/1‑250 (MQ=255)

                                                                                                                                                            cacTTCTTTCCA‑GCGCAGCGCGATTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGGATGGTGCGCTCCACCAGCGAGGATCCTCGCTTCGGAA‑GCAGTTGGCATATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCgtgt                                                                                     <  2:559808/251‑1 (MQ=255)

                                                                                                                                                             acTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTc                                                                                    <  1:412330/251‑1 (MQ=255)

                                                                                                                                                               ttCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGGATGGGGCGCGCCCACAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTTCGCTTTCCCCAGCGCATTCGCCACCGGCCGTTTTGGCAGGGGGTGGGGGGTTTTCAcccgggcgg                                                                                  >  2:401323/1‑244 (MQ=255)

                                                                                                                                                                  gTTCCA‑GCGCAGCGCGTCTTTAAACGCGGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGAtac                                                                               <  1:662150/251‑3 (MQ=255)

                                                                                                                                                                     ccA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAg                                                                             <  1:981795/250‑1 (MQ=255)

                                                                                                                                                                         gcgcagcgcGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCCCCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGAGCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTaa                                                                         <  2:504026/251‑1 (MQ=255)

                                                                                                                                                                            cagcgcGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCa                                                                      <  1:741446/251‑1 (MQ=255)

                                                                                                                                                                            cagcgcGTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCa                                                                      <  2:454029/251‑1 (MQ=255)

                                                                                                                                                                                gcgTCTTTAAACGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAgcag                                                                  <  1:652328/251‑4 (MQ=255)

                                                                                                                                                                                        aaactcgggCAGGCGGGTATCTTTATCAATCACCTGGGTTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGAGCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTaa                                                           <  2:513033/243‑1 (MQ=255)

                                                                                                                                                                                        aaaCGCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAAc                                                          <  1:381576/251‑1 (MQ=255)

                                                                                                                                                                                        aaaCGCTGGCAGGCGGGGATCGGTATCAATCACCGAGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCAAATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGGGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGGTGCGCTTTCCCCAGCGCATTCGCCACCGGCAGGTTCGGCAGGGGGGGGGGGGGTTTTACCGTGTCGATAAAGGTAACCAGCAGAAGGTaa                                                           >  1:389155/1‑250 (MQ=255)

                                                                                                                                                                                            gCTGGCAGGTGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATc                                                      <  2:233832/251‑1 (MQ=255)

                                                                                                                                                                                            gCTGGCAGGCGGGTATCTGTAGCAATCACCGGGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTGGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATc                                                      <  2:511647/251‑1 (MQ=255)

                                                                                                                                                                                            gCTGGCAGGCGGGTATCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATc                                                      <  1:285187/251‑1 (MQ=255)

                                                                                                                                                                                               ggcaggctggtatCGTTATCAATCTCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGGTGGCAGATCCGGGTGCGCTGTCCACAGCGCATTCGCCACCGGCAGTTGCGGCATGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGcc                                                   <  1:403323/239‑1 (MQ=255)

                                                                                                                                                                                                      gggtatctttATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTGCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCTCCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTTACCAGCAGACGGTAACGATCGCCGAGatca                                            <  1:394827/243‑1 (MQ=255)

                                                                                                                                                                                                          atCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatc                                         >  1:829869/2‑250 (MQ=255)

                                                                                                                                                                                                          atCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTCGAGGTTCAGCGCCTCGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCCGTTGGCAGATCCGGTTGCGATTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGCGCGGTGTTTTCACCGTATCGATACAGTTACACAACAGACGGTAACGATCGCCGAGatcaatc                                         >  1:589125/2‑250 (MQ=255)

                                                                                                                                                                                                            cGTTATCAATCACCGTTATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatcagg                                      <  2:387439/251‑3 (MQ=255)

                                                                                                                                                                                                            cGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGCGTGCGGTGTTTTAACCGTGTCGGACCAGTTAACCAGCAGACGGAACCGGTCGCCATGatcaatcagg                                      >  2:829669/1‑249 (MQ=255)

                                                                                                                                                                                                            cGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGatcaatcagg                                      <  1:489726/251‑3 (MQ=255)

                                                                                                                                                                                                            cGTTATCAATCACCGTGATTTCAATGTCGTGCATCGCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGagcaatcag                                       <  2:391553/250‑2 (MQ=255)

                                                                                                                                                                                                              ttATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAggg                                                                                                          <  2:690638/181‑1 (MQ=255)

                                                                                                                                                                                                              ttATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAggg                                                                                                          >  1:690638/1‑181 (MQ=255)

                                                                                                                                                                                                                      caCCGTGATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACa                             <  2:120580/250‑1 (MQ=255)

                                                                                                                                                                                                                            gATTTCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTg                      <  1:908744/251‑1 (MQ=255)

                                                                                                                                                                                                                                tCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCCGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                  >  1:525423/1‑248 (MQ=255)

                                                                                                                                                                                                                                tCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                  >  1:709733/1‑248 (MQ=255)

                                                                                                                                                                                                                                tCAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggac                   <  1:575496/250‑3 (MQ=255)

                                                                                                                                                                                                                                 cAATGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                  <  1:145590/250‑4 (MQ=255)

                                                                                                                                                                                                                                 cAATGGCGTGCATCGCGGCGAACGGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTggacc                  <  2:342012/250‑4 (MQ=255)

                                                                                                                                                                                                                                  aaTGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCgg                >  1:955330/1‑251 (MQ=255)

                                                                                                                                                                                                                                  aaTGTCGTGCATCTCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGGTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGGTACAGTTAACCAGCAGACGGGAACGGTCGCCGAGATCAATCAGGCTGGCGACAAACCCTGGACCgg                >  2:310691/1‑251 (MQ=255)

                                                                                                                                                                                                                                   aTGTCGTGCATCTCGGCAAACTGGCGAATATCGTTGAGGTTAAGCGCATGGCTGAAGACGGTATGGAGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATCCAGTTAACCAGCAGACGGTACAGATCGCCGAGATCACTCAGGCTGGCGCCACTCGCTGGACCggg               >  2:630943/1‑250 (MQ=255)

                                                                                                                                                                                                                                    ggTCGTGCATCTCGGCGAACTGGTGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGtt              <  1:652843/250‑1 (MQ=255)

                                                                                                                                                                                                                                     ggcGTGCATCGCGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGtt              <  2:709733/248‑1 (MQ=255)

                                                                                                                                                                                                                                     gaggtgccttccGGCGAATTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGTTATGGTGCTCGTCGCCAGCGAGGATCCTCGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCTCCGGCAGTTTCGGCAGGTAGTGCGGTGCTTTCACCTTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAAGCAGGCTGGCGACAATCGCTGGACCgggtt             <  2:82073/240‑4 (MQ=255)

                                                                                                                                                                                                                                         ggCATCGCGGCGAAGTGGCGCATATCGTTGATGGTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGAGCCTCGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTGCGGCAGGGAGGGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGAGCGCCGAGAGCAATCAGGCTGGCGACAATCGCTGGATCGGTGTGGGt         <  2:168354/250‑1 (MQ=255)

                                                                                                                                                                                                                                              ctcGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGCGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCCCCGGCAGTTTCGGCAGGGAGTGCGGGGCTTTCACCGTGTCGCTCCCGTTAACCAGCAGCCGGTTCCGAGCGCCGTGCTCCATCCGGCTGGCGACACTCGCTGGCCCGGtttccgagttgg    >  2:185765/1‑242 (MQ=255)

                                                                                                                                                                                                                                                cGGCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGGGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGaa   <  2:369577/250‑1 (MQ=255)

                                                                                                                                                                                                                                                 ggCGAACTGGCGCATATCGTTGAGGTTCAGCGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATACAGTTAACCAGCAGACGGTAACGATCGCCGAGATCAATCAGGCTGGCGACAATCGCTGGACCGGTTTGGGTATTGa    >  1:891850/1‑248 (MQ=255)

GATAAGACGCGCCAGCGTCGCATCAGGCGTTGAATGCCGGATGCGCTTTGCTTATCCGGCCTACAAAATCGCAGCGTGTAGGCCAGATAAGACGCGTCAGCGTCGCATCAGGCGTTACATACCGGATGCGGCTACTTAGCGACGAAACCCGTAATACACTTCGTTCCA‑GCGCAGCGCGTCTTTAAACGCTGGCAGGCGTGTGTCGTTATCAATCACCGTGATTTCAATGTCGTGCATCTCGGCGAATTGGCGCATATCGTTGAGGTTCAGTGCATGGCTGAAGACGGTATGGTGCGCGCCACCAGCGAGGATCCACGCTTCGGAA‑GCAGTTGGCAGATCCGGTTGCGCTTTCCACAGCGCATTCGCCACCGGCAGTTTCGGCAGGGAGTGCGGTGTTTTCACCGTGTCGATGCAGTTAACCAGTAGACGGTAACGATCGCCGAGATCAATCAAGCTGGCGACAATCGCTGGGCCGGTTTGGGTATTGAAG  >  NC_000913/66700‑67189

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