BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,853	C→T	A223A (GCG→GCA)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,853	0	C	T	97.1%	413.1 / ‑2.5	139	A223A (GCG→GCA)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand): ref base C (1/3); new base T (50/83); total (52/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 139 total aligned reads displayed.

CGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACAT  >  NC_000913/73607‑74085

 ggTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCACTGGCCTTACTCAATCGCTGACTCCGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGAGGGCAGGATCGAAGTCGTAACTAAGCGCCTGATAGAGTGCCAGCTAGGAAGTGGTTGCCTg                                                                                                                                                                                                                                      >  2:121873/1‑251 (MQ=255)

 ggTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTg                                                                                                                                                                                                                                      >  1:875593/1‑251 (MQ=255)

  gTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGc                                                                                                                                                                                                                                     >  1:96920/1‑251 (MQ=255)

   tGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCg                                                                                                                                                                                                                                    >  1:119386/1‑251 (MQ=255)

     gcagcaACAGCAGCAGCGCCAGCACAATTCACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCgg                                                                                                                                                                                                                                   >  1:764371/1‑250 (MQ=255)

       agcaACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGAcc                                                                                                                                                                                                                                <  1:371721/251‑1 (MQ=255)

        aaaacaggtgcagcGCCAGCAAAATTAACACCGTGTCGCAAATGCGGCAATGCAGACGCTCTTCCGGCTCACGCCAGCCTTGCAGACGCGTGGTGCCTGGCGCAATGGCCTTACACAATCGCTCACTCAGCAGCACCAACCCAAGGCTACACTCCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGTATCGTAGTCGTAACGCAGCGCCTGATAGATTGCCAGCTCGATAGTGTTTGCCTGCTGACCg                                                                                                                                                                                                                               <  2:443089/242‑1 (MQ=255)

          aaCAGCAGCAGCACCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGcc                                                                                                                                                                                                                             <  1:943347/251‑1 (MQ=255)

           aCAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGccg                                                                                                                                                                                                                            >  1:923249/1‑250 (MQ=255)

             agcagcagcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGccgcc                                                                                                                                                                                                                          >  1:498626/1‑248 (MQ=255)

               cagcagcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGAGCAGCGACAGCATCGCCGCGCGGGAAGGAGCGTAGGCGTAACGCAGCGCCGGATAGACTGCCAGCACGACAGTGGAGGCCTGCGGACCgacgcacc                                                                                                                                                                                                                        >  2:398633/1‑244 (MQ=255)

                  cagcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCg                                                                                                                                                                                                                     <  1:376599/251‑1 (MQ=255)

                   agcGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGCCGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCCCCAACCCGAGGCAACACACCATCTGGATCAGCGCCCGCATCGCCGCGCGGGCAGGATCGTCGTCGTAACCCAGCGCCTGCTAGATTGCCAGCTCGTTAGTGGTTGCCTGCGGACACccccccgaca                                                                                                                                                                                                                     >  2:358905/1‑246 (MQ=255)

                          gCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGTCGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACCCCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACg                                                                                                                                                                                                             <  2:263189/251‑1 (MQ=255)

                            acaATTAACACCGTGTCGAAAATGCTGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGt                                                                                                                                                                                                           <  2:821099/251‑1 (MQ=255)

                                ttAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTgacg                                                                                                                                                                                                       >  1:482847/1‑248 (MQ=255)

                                 tAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCg                                                                                                                                                                                                       >  1:139514/1‑250 (MQ=255)

                                 tAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCg                                                                                                                                                                                                       <  1:92582/250‑1 (MQ=255)

                                 tAACACCGTGTCGCAAAGGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCACCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCg                                                                                                                                                                                                       <  2:577473/250‑1 (MQ=255)

                                  aaCACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGaa                                                                                                                                                                                                     <  1:76892/251‑1 (MQ=255)

                                   acacCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAg                                                                                                                                                                                                    <  2:452417/251‑1 (MQ=255)

                                   acacCGTGTCGCAAATCCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCATGGCGCAATGGACTTACTCAATCGCTGACTCAGCAGCACCAACCCGCGGCAACACACCACCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGCAGCGTAGTCGTAACTCAGCGCCTGATAGATGGCCAGCTGGATAGTGATTGCCTGAGGACCACCGCCCACCGAGAGCACGGGGGCgggg                                                                                                                                                                                                    >  2:801013/1‑248 (MQ=255)

                                     cccGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCt                                                                                                                                                                                                  <  1:483107/250‑1 (MQ=255)

                                      ccGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTg                                                                                                                                                                                                 <  2:273302/251‑1 (MQ=255)

                                          gtCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGa                                                                                                                                                                                             <  1:26516/251‑1 (MQ=255)

                                             gCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGCTAGATTGCCAGCTCGATAGGGGTTGCCTGCGGACCGCCGCCCAGCGGTAGCACGGTGGCGAAGCTGGCGaaa                                                                                                                                                                                           >  2:406638/1‑250 (MQ=255)

                                             gCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGGTTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGaaa                                                                                                                                                                                           >  1:178774/1‑250 (MQ=255)

                                                 aTGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCAAGGTGGCGAAGCTGGCGAAACagag                                                                                                                                                                                      <  1:832727/251‑1 (MQ=255)

                                                     ggCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGCAGCTGGCGAAACAGAGCAt                                                                                                                                                                                   >  1:150563/1‑250 (MQ=255)

                                                      gCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCGGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATaa                                                                                                                                                                                 <  2:838403/251‑1 (MQ=255)

                                                      gCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGGGCCTGGCGCAATGGCCTTACTCAATCCCTGCCTCAGCAGCACCAACCCGAGGCAACACACCATCTGGGTCAGCGCCAGCATCACAGCGCCGGCAGGGTCGTAGTCGTAACTCCGCGCCTGACAGATTGCCAGCTCGCTAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCAAGGCGGCGGAACTGGCGAAACAGAGcg                                                                                                                                                                                    >  2:946775/1‑247 (MQ=255)

                                                        tATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCAACTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGCCCGCCGCCCAGCGATAGCACGGTGGCGAAGATGGAGAAACAGAGCCTaa                                                                                                                                                                                 >  1:475455/1‑249 (MQ=255)

                                                         aTGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAg                                                                                                                                                                               <  1:474229/250‑1 (MQ=255)

                                                           gCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCATCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATa                                                                                                                                                                            <  1:761361/251‑1 (MQ=255)

                                                           gCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGAt                                                                                                                                                                             >  1:67083/1‑250 (MQ=255)

                                                           gCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATa                                                                                                                                                                            <  1:416954/251‑1 (MQ=255)

                                                           gCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCAGCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGAt                                                                                                                                                                             <  2:731936/250‑1 (MQ=255)

                                                              gACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑Agt                                                                                                                                                                        <  1:985110/251‑2 (MQ=255)

                                                              gACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑Agt                                                                                                                                                                        <  2:49464/251‑2 (MQ=255)

                                                              gACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑Agt                                                                                                                                                                        <  1:313888/251‑2 (MQ=255)

                                                                    cGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGAGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTgcagc                                                                                                                                                                   >  1:621961/1‑250 (MQ=255)

                                                                       ccGGGTCGCGCCAGCCTTGCAGCCGCTTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACTCCAGCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAAc                                                                                                                                                                <  2:575363/250‑1 (MQ=255)

                                                                       ccGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAAc                                                                                                                                                                <  1:45197/250‑1 (MQ=255)

                                                                        cGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCTGCGCGGGCAGGATCGTAGTCGTAACTCAGCGACTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGCTAGCACGGTGGCGAAGCTGGCGACACAGAGCATAAAGATA‑AGTGCAGCAACcg                                                                                                                                                              >  1:636149/1‑251 (MQ=255)

                                                                        cGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTAGATAGTGGTTGCATGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATGAAGATA‑AGTGCAGCAACcg                                                                                                                                                              >  1:861808/1‑251 (MQ=255)

                                                                           gTCGCGCCAGCCTTGCAGCAGCGGGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACcggcg                                                                                                                                                           <  1:406638/251‑1 (MQ=255)

                                                                            ccgcgCCAGCCTTGAAGCAGCGTGGTGCCTGGCGCTATGGCCTTACTCAATCGCTGACTGATCAGCACCAACCCTAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGAGTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGGGGCGAAGCTGGTGAAACAGAGCAGAAAGATA‑AGTGCGGCACCcggcgg                                                                                                                                                          <  2:489982/250‑1 (MQ=255)

                                                                                gcCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCggggtt                                                                                                                                                      >  1:205853/1‑248 (MQ=255)

                                                                                gcCAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGAtt                                                                                                                                                      >  1:196406/1‑251 (MQ=255)

                                                                                 ccAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGTGGGCGAAGCTGGCGAAACAGAGCATAGAGATacagggagccaccggcggggttt                                                                                                                                                      >  2:538834/1‑229 (MQ=255)

                                                                                 ccAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAAGCGCTGACTCAGCAGCACCAACCCGAGGCAACACTCCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGAttt                                                                                                                                                     <  1:707355/251‑1 (MQ=255)

                                                                                  cAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTg                                                                                                                                                    <  1:803040/251‑1 (MQ=255)

                                                                                  cAGCCTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTg                                                                                                                                                    <  1:77489/251‑1 (MQ=255)

                                                                                      cTTGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTAATCAATCGCTGACTCAGCAGCACCAACCCGACGCAACACACCATCCGGATCAGCGCCAGCATAGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGCTAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAGGCTGGCGAAACAGAGCATAAAGATACAG‑GCAGCAACCGGCGGGAATTGtacc                                                                                                                                                >  2:941766/1‑248 (MQ=255)

                                                                                       ttGCAGCAGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCATGCGGACCGCCGCCCAGCGATAGCCAGGTGGCGAGGCTGGCGAACCAGAGCATAAAGATA‑AGTGCAGCGACAGGCggggtgttggcgc                                                                                                                                               >  2:351334/1‑241 (MQ=255)

                                                                                         gcagcaGCGTGGTGCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGAGTTGTCGCCg                                                                                                                                              <  1:161854/250‑1 (MQ=255)

                                                                                               gcGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCAc                                                                                                                                       <  1:336811/251‑1 (MQ=255)

                                                                                                  tggtgCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATC‑AGAGCAGCAACCGGCGGGAATTGGCGCCGTAACCAcagc                                                                                                                                    >  2:310878/1‑248 (MQ=255)

                                                                                                   ggtgCCTGGCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGcc                                                                                                                                   >  1:265988/1‑251 (MQ=255)

                                                                                                      gcctggCGCAATGGCCTTACTCAAACGCTGACCCAGCAGCACCAACCCGAGGCAACACGCCATCTGGACCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACCCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCGGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGcaat                                                                                                                                 <  2:302197/246‑4 (MQ=255)

                                                                                                       cctgtcccAATGGCCTTACTCAATCCCTGACTCAGCAGCACCAACCCGAGGCAACACCCCATCTGGACCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTc                                                                                                                               <  2:346974/244‑1 (MQ=255)

                                                                                                          ggCGCAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTcgac                                                                                                                            <  1:793782/251‑1 (MQ=255)

                                                                                                           gcgcAATGGCCTTCCTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTcgacg                                                                                                                           <  2:929817/251‑1 (MQ=255)

                                                                                                             gcAATGGCCTTACTCAATCGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGaa                                                                                                                         <  1:290168/251‑1 (MQ=255)

                                                                                                                    ccTTACTCAAGCGCGGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGaaa                                                                                                                  <  2:902546/251‑1 (MQ=255)

                                                                                                                           caaccGCTGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATGGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCTCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATg                                                                                                             <  2:94273/245‑1 (MQ=255)

                                                                                                                               cGCTGACTCAGCAGCCCCAACCCGAGGCAACACACCATCTGGAGCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCTTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGcc                                                                                                       <  2:577083/251‑2 (MQ=255)

                                                                                                                                  tGACTCAGCAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGccgcg                                                                                                    >  1:594362/1‑247 (MQ=255)

                                                                                                                                       cagcAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGccgcgcat                                                                                                 <  1:607672/246‑8 (MQ=255)

                                                                                                                                        agcAGCACCAACCCGAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAACTGCCAGccgcgcatacc                                                                                              >  1:642671/3‑241 (MQ=255)

                                                                                                                                                caacccGAGGCAATACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTATTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGGTGCCGGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTggg                                                                                      <  2:127998/248‑1 (MQ=255)

                                                                                                                                                      gAGGCAACACACCATCTGGATCAGCGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCaa                                                                                <  1:851094/251‑1 (MQ=255)

                                                                                                                                                          catcacaCCATCTGGAGCAGCGCCAGCTTCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAgctg                                                                            <  2:525244/248‑4 (MQ=255)

                                                                                                                                                               cacCATCTGGATCAGCGCCAGCACCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACgtt                                                                       <  2:862751/251‑3 (MQ=255)

                                                                                                                                                                cccATGTGGACCAGCGCCAGCAGCTCCGCGCGGTCAGGGACGTAGTTGTAACTCAGCGCCTGATAGAGTGCCAGCGCTATAGTGGTTGCCTGCGGACCGCCGCCCAGCGAGAGCACGGTGGCGAAGCTGTCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACgttg                                                                      <  2:482847/250‑4 (MQ=255)

                                                                                                                                                                       gggtcagcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGcc                                                                <  2:464379/246‑1 (MQ=255)

                                                                                                                                                                       ggatcagcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCg                                                               <  1:466743/247‑1 (MQ=255)

                                                                                                                                                                       ggatcagcGCCAGCACCGCCGCGCGGGCAGGATCGCAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGTCCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTc                                                                   <  1:344529/243‑1 (MQ=255)

                                                                                                                                                                          tcagcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCgggg                                                            >  1:767680/2‑251 (MQ=255)

                                                                                                                                                                            agcGCCAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCgggaat                                                          >  1:811311/1‑248 (MQ=255)

                                                                                                                                                                                 cAGCATCGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGtttt                                                     <  1:358905/251‑1 (MQ=255)

                                                                                                                                                                                       cgCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTGGCCTGCGGACCGCCGCCCAGCGATAGCTCGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTg                                               <  1:420721/250‑1 (MQ=255)

                                                                                                                                                                                              cgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAg                                        <  1:193537/251‑1 (MQ=255)

                                                                                                                                                                                              cgGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTAACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTgg                                          >  1:439899/1‑249 (MQ=255)

                                                                                                                                                                                                 cttccgatctAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGta                                      >  1:773032/10‑250 (MQ=255)

                                                                                                                                                                                                   aGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGtaata                                   <  1:402159/251‑1 (MQ=255)

                                                                                                                                                                                                      aTCGTAGTCGTAACTCATCGCCTGATAGATTGCCAGCACGATAGTGGTTGCCTGCGGACCGCAGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAACCAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGc                                <  1:958834/251‑1 (MQ=255)

                                                                                                                                                                                                      aTCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGCCGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGc                                >  1:901332/1‑251 (MQ=255)

                                                                                                                                                                                                             tCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGc                         >  1:766292/1‑251 (MQ=255)

                                                                                                                                                                                                               gTAACTCGTCGCCTGTTTGAGTGCCAGGTAGATAGGGGTTGCCGGCGGACCGCCGCCCAGCGATAGCACGGGGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCa                       <  2:983258/251‑1 (MQ=255)

                                                                                                                                                                                                                   cccAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGAACGGGGGCGAAGCTGGCTAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCgg                   <  2:811289/249‑1 (MQ=255)

                                                                                                                                                                                                                     cAGCGCCTGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCggt                  <  1:395699/250‑2 (MQ=255)

                                                                                                                                                                                                                       gcgcATGATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGa               <  1:11328/251‑1 (MQ=255)

                                                                                                                                                                                                                         gcctggaagtatGCCCGCTCGATAGTGGTCTCCTGCGGACCGCCGCCCAGCGACAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCACGCGGCGGGATTGGTCGCCGTAACCACGGCCATTCGACGAAGCGGATAAAATGCCGGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCTCCATCGGTAGAtt             <  2:144775/240‑1 (MQ=255)

                                                                                                                                                                                                                             gATAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGGGCAACCCAAGCTGGGCGGCAAGCTAACGTTGGTCGCCGGGGATGTTTTCCAGTGCCTGGGGCAATAAGCGGGTCCCCATCggtgggttaaa          >  2:160580/1‑241 (MQ=255)

                                                                                                                                                                                                                             gATAGATGGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaa         <  1:14051/251‑1 (MQ=255)

                                                                                                                                                                                                                               tAGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGACGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGGGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCAGGAGTGATAAGCGGCTCGcccacggttggttaaaaa        >  1:69374/1‑234 (MQ=255)

                                                                                                                                                                                                                               aaGATTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaa       <  2:265988/250‑1 (MQ=255)

                                                                                                                                                                                                                                  aTTGCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCCGGGAATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatac    <  1:121873/251‑4 (MQ=255)

                                                                                                                                                                                                                                     gCCAGCTCGATAGTGGTTGCCTGCGGACCGCCGCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGATCATAAAGATA‑AGTGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCCGCGCATCCCAAGCTGGGCGGCAAGCTGACGTTGTTCGCTGGGAATGTTTTCCAGTGTCAGGAGTAATAAGCGGCTCGCCATCGGTAGATTaaaaaatacgt  <  1:310966/250‑6 (MQ=255)

CGGTGGCAGCAACAGCAGCAGCGCCAGCACAATTAACACCGTGTCGCAAATGCGGCTATGCAGACGATCGTCCGGGTCGCGCCAGCCTTGCAGCAGCGTGGTGCCGGGCGCAATGGCCTTACTCAATCGCTGACTCAACAGCACCAGCCCGAGGCAGCACACCATCTGGAGCAGCGCCAGCATTGCCGCGCGGGCAGGATCGTAGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATA‑AGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCGGGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCAGATTAAAAAACACAT  >  NC_000913/73607‑74085

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