BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	349,717	T→C	4.9%	intergenic (+145/‑295)	prpB → / → prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	349,717	0	T	C	4.9%	507.2 / 8.9	176	intergenic (+145/‑295)	prpB/prpC	2‑methylisocitrate lyase/2‑methylcitrate synthase
Reads supporting (aligned to +/- strand): ref base T (69/92); new base C (4/6); total (73/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

Only 100 of 180 total aligned reads displayed.

GTCTACAACGTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAA  >  NC_000913/349429‑349998

gTCTACAACGTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGt                                                                                                                                                                                                                                                                                             >  1:1194575/1‑289 (MQ=255)

  cTACAACGTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGa                                                                                                                                                                                                                                                                                           <  1:1538103/289‑1 (MQ=255)

     cAACGTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACtt                                                                                                                                                                                                                                                                                        >  2:511681/1‑289 (MQ=255)

       aCGCCCGGCGCCAGTAAGGCACGGAGAAAATCACCAACGACCCCATGCAGACCGGAAACGAGCTGAACGAACGAATCAACTACTCCCAGAACGAAGGGAAGCTCGACAACCTGTGTGCACGTAGCCAGGTGAAATAACAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCATACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCTATGCCTGATGCGCCGCTTACGTGACTTAt                                                                                                                                                                                                                                                                                      <  2:1245841/289‑1 (MQ=255)

         gTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCa                                                                                                                                                                                                                                                                                    <  1:562634/289‑1 (MQ=255)

          tCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGCCTTATCAg                                                                                                                                                                                                                                                                                   >  1:331954/1‑289 (MQ=255)

               ggcCAGAGAGGCAGGCAAAAAACCTGCAGGGACACCAGTCAGGCCCGCAACGAGCTGTACGAAAGCATCAACTGCTCCCAGTACGAAGAGAAGCTCGACAACCTTGTTGGCCGTAGCAAAGTTAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCGGATGCGTCGCTGACGCGACTTATCTGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCatggct                                                                                                                                                                                                                                                                               <  2:966183/287‑6 (MQ=21)

                 ccAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACg                                                                                                                                                                                                                                                                            <  1:844806/289‑1 (MQ=255)

                   aGGAAGGCACGCAGAAAAGCGTCATCGAAACCAGGCAGACCCGCACCGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTGGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                                                                                                                                                          <  2:28786/289‑1 (MQ=255)

                    ggAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCCACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCCGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGCCTTATCAGGCCCACgggg                                                                                                                                                                                                                                                                         >  2:698663/1‑289 (MQ=255)

                      aaGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                                                                                                                                                          <  1:622502/286‑1 (MQ=255)

                      aaGCCACGCAGACAAGGCTCATCCACACCATGCAGACCCGCAACGAGCTGTACGACAGCATAAACGACTCCCAGTACGAAGAGAAGCTCGACAACCTGGTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGCTTGTATTCGACAGCCGATGCCGGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACggg                                                                                                                                                                                                                                                                          <  2:826571/286‑1 (MQ=25)

                             gCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGaa                                                                                                                                                                                                                                                                <  1:913338/289‑1 (MQ=255)

                                  aaagcagaagcGAAACCATGCGGACCCGAAACGAGCAGTACGAAAGCATCAACTACTACCAGTACGAAGGGAAGCTCGACAACCTGTTGGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCGGAGGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCGGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTg                                                                                                                                                                                                                                                             <  2:1251285/277‑1 (MQ=25)

                                    aGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAGGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCACCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAgg                                                                                                                                                                                                                                                         <  1:1314148/289‑1 (MQ=255)

                                    aGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg                                                                                                                                                                                                                                                         <  1:171668/289‑1 (MQ=255)

                                      cggcagcGACACCAGGGACAACCGCAACAGGCTGGACGAAGACATCAACTACTACCAGTAAGAAGATAAGGTCGACAACCGTTTTCCCCGTAGCCAGGTGAAAGAAAAAACGCCCGTGGATTTGATTCTACAGCCGATGCCTGTTGCGTCGCGGAGGCGAGATATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTATGt                                                                                                                                                                                                                                                        <  2:940221/282‑1 (MQ=18)

                                      cGTCCTCGACACCCTGCGCACCCGCAACGAGCGGTACTAAAGCATCAACTACAACCAGAACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                                                                                                                                                                                                                                        <  2:17667/288‑1 (MQ=255)

                                          aTCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                                                                                                                                                                                                                   <  1:1359116/289‑1 (MQ=255)

                                          aTCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                                                                                                                                                                                                                   <  1:188803/289‑1 (MQ=255)

                                             gagaccagccggACCCGCAACTGGCTGTACGAAACATTCAACTCCGACCCGTACGAAGAGAAGCTCGACGACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGGATTCGACAGCCGATGCCTGCTGCATCGCTGACGCGACTTTTCACGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCa                                                                                                                                                                                                                                                                     <  2:786034/257‑1 (MQ=21)

                                             aaaacCATGCAGCCCCGCAACGGGCTGTACGAAAGAATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                                                                                                                                                                                                                <  2:250012/286‑1 (MQ=255)

                                                acCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTTAACTGTAGGTCGGATAAGAc                                                                                                                                                                                                                                              >  2:306751/1‑288 (MQ=255)

                                                acCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGAt                                                                                                                                                                                                                                                   >  1:956292/1‑283 (MQ=255)

                                                  cATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCa                                                                                                                                                                                                                                           <  1:1028145/289‑1 (MQ=255)

                                                    tGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGAAGCATa                                                                                                                                                                                                                                         >  1:250012/1‑289 (MQ=255)

                                                       aGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCCCGACAACCTGTTTGCCCGTCGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGGCAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAgctgactgtaatacacatct                                                                                                                                                                                                                                       >  1:786034/1‑269 (MQ=255)

                                                            cGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCa                                                                                                                                                                                                                                 >  1:1262183/1‑289 (MQ=255)

                                                              cAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGt                                                                                                                                                                                                                                                                                             >  2:1550146/1‑227 (MQ=255)

                                                              cAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                                                                                                                                                                                                                >  2:1422848/1‑288 (MQ=255)

                                                              cAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCAt                                                                                                                                                                                                                                >  1:1133909/1‑288 (MQ=255)

                                                               aaCGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                                                                                                                                                                                                              <  1:178338/289‑1 (MQ=255)

                                                                     cTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACaa                                                                                                                                                                                                                         >  1:427511/1‑288 (MQ=255)

                                                                              aGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGCCGTGACTTACCAGGCCTACGGGGTGACCTGAACTGCAGGTCGGATAAGACGCATAGCGTCGCAGCCGACACCCCATGCCt                                                                                                                                                                                                                >  2:596783/1‑288 (MQ=255)

                                                                                  tCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                                                                                                                                                                                                                <  1:1299446/252‑1 (MQ=255)

                                                                                  tCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAg                                                                                                                                                                                                                                                >  2:955290/1‑252 (MQ=255)

                                                                                  tCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTGATGAGGCCTACGGGGTGCACTGAACGGTAGGTCGGATAAGACGCATAGCGTCGCCTCCGACAACCGCTGCCTGATg                                                                                                                                                                                                            >  1:600928/1‑288 (MQ=255)

                                                                                          cacCAGTACGAAGGGCAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAACAAAAAGCTCCCGTTGATTGGATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACGTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCAGCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGa                                                                                                                                                                                                   <  1:1535751/288‑1 (MQ=255)

                                                                                              ggataaaggaacaGCTCGACAACCTGTTGGCCCTTAGCAATGTGAATTAAAAAACGCCCTTTGATTGTATTCGACAGCCGATGCCGGATGCTGCGCTGACGCGACTTAGCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgc                                                                                                                                                                                                <  2:982370/276‑1 (MQ=17)

                                                                                                tACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTCGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg                                                                                                                                                                                               >  1:577645/1‑287 (MQ=255)

                                                                                                tACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGAcgcg                                                                                                                                                                                               >  1:656426/1‑287 (MQ=255)

                                                                                                tACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGa                                                                                                                                                                                              >  1:952975/1‑288 (MQ=255)

                                                                                                    aaGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGt                                                                                                                                                                                                                                                        >  1:778009/1‑226 (MQ=255)

                                                                                                          aGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAgg                                                                                                                                                                                    >  1:94343/1‑287 (MQ=255)

                                                                                                            cTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCt                                                                                                                                                                                                                                                                               <  2:558231/195‑1 (MQ=255)

                                                                                                            cTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCt                                                                                                                                                                                 >  1:931868/1‑288 (MQ=255)

                                                                                                                aatacatGTATGCCAGTGGCCAGGTGAAAAAAAACACGCCCTTGGATTGAATGCGCCAGCCGATGCCTGATGCGTCGCGGACGCGACGTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACggg                                                                                                                                                                            <  2:395843/283‑1 (MQ=14)

                                                                                                                 cAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACgggg                                                                                                                                                                           >  1:1326625/1‑289 (MQ=255)

                                                                                                                 cAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAACCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTTAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACggg                                                                                                                                                                            <  1:814051/288‑1 (MQ=255)

                                                                                                                     cTGTTTGCCCGTAGCCAGGTGAAATAATAAACGCCCGTTGATTGTATTCGACAGCCGATGCAGGATGCGTCGCTGACGCTACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGTGACGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTTAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCa                                                                                                                                                                       <  1:538680/289‑1 (MQ=255)

                                                                                                                        tttGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTg                                                                                                                                                                    <  1:846180/289‑1 (MQ=255)

                                                                                                                         ttGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGGCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCGGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTg                                                                                                                                                                    <  1:934394/288‑1 (MQ=255)

                                                                                                                         tggCCCGTAGCCAGGTGAAAAAAAAAACGCCCGTTGATTGGATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCGACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTg                                                                                                                                                                    <  2:1133909/286‑1 (MQ=255)

                                                                                                                             cagTAGCCAGGTGAAATAAATAACGCCCGTTGATTTTATTCGACATGCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTAGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACt                                                                                                                                                                <  1:156063/286‑1 (MQ=255)

                                                                                                                                tAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCCGGCCTACGGGGTGCACTGAACTGGAGGTCGGATAAGACGCATCGCGTCGCATCCGACAACCGATGCCTGATg                                                                                                                                                                                                            >  1:329459/1‑242 (MQ=255)

                                                                                                                                tAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAACGGTGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTCCGGGGCGCACTTAACttgag                                                                                                                                                            >  1:459869/1‑285 (MQ=255)

                                                                                                                                  gCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAgg                                                                                                                                                           >  1:750031/1‑288 (MQ=255)

                                                                                                                                       acagaatagaacagccccgtcaatggatttCGGCCGCCGAGTCATGATGCGTCGCTGGCGGAACGGAACAGGGCTACGGTGGGAAGAGAACGGTAGGTCGGATAAGTCGGCTGGCTTCGCAGCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TAAGGCCTACGGGGTGCACTGAACTGTAGGTCGGAt                                                                                                                                                     <  2:810077/261‑1 (MQ=11)

                                                                                                                                                 aaaCGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGCCGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGCTAAGACGCATAGCGTCCCCTCCGACAACCGATGCCTGATGCACCGCTGAACCGACTTACCCAGGCC‑ACGGTGTTCACCTAACTGGAGGTAGGCCAAGACGCATa                                                                                                                                           >  2:1586764/1‑289 (MQ=14)

                                                                                                                                                       ccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑t                                                                                                                                                                                        <  2:754656/238‑1 (MQ=255)

                                                                                                                                                       ccGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑t                                                                                                                                                                                        >  1:754656/1‑238 (MQ=255)

                                                                                                                                                                   ttCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTCGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑Ac                                                                                                                         >  1:637657/1‑288 (MQ=255)

                                                                                                                                                                      gACACCCGTTGCCTGCTGCGTCGCTGACGAGACTTATCAGGCCTACGTGTCGCCCAGAACTGTAGGTCGGATAAGACGGATTGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGGGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGAt                                                                                                                     <  2:1100378/289‑1 (MQ=11)

                                                                                                                                                                       aCAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                                                                                                                                                                                                                   >  1:1463289/1‑148 (MQ=255)

                                                                                                                                                                       aCAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACt                                                                                                                                                                                                                                                                   <  2:1463289/148‑1 (MQ=255)

                                                                                                                                                                       aCAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGAt                                                                                                                     >  1:300190/1‑288 (MQ=40)

                                                                                                                                                                           ccGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTg                                                                                                                >  1:385002/1‑289 (MQ=255)

                                                                                                                                                                               tgccatatgctccggtGACGGGACTGATCGGGCCTACGAGTTGCACTAAACTGTAGATCGAATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATg                                                                                                             <  2:777165/273‑1 (MQ=14)

                                                                                                                                                                                            gCTGACGCGACTTATCTGGTGTACGTAGTGCACTGACCTGGTGGACGTAGAATACGGATGGCGTCGCACGCTGCAACCGATGCGTGGAGCGCCGAGGACGTGACTTATCACGCCTACGGGGTGCACTGAACTGAAGGTCGCATAAGAGTCAGAGCGTCGCATCCGACAACCGATGACTGATGCGCCGCTGACGCGACTGA‑TCAGGCCTACGAGGTGCACTGAACTGTAGGTCTGATAAGACGGATGGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGTGa                                                                                               <  2:933705/289‑1 (MQ=17)

                                                                                                                                                                                                 cgcgACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                          <  1:110754/289‑1 (MQ=255)

                                                                                                                                                                                                 cgcgACTTATCAAGCCTATGATGTGCACTGAACTGTAGGTCTGATAAGACTGATGGCGTCGCATCCGACAACCGATGCGTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCGCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTAt                                                                                          <  1:1354784/289‑1 (MQ=255)

                                                                                                                                                                                                    gACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCa                                                                                        <  1:511681/288‑1 (MQ=255)

                                                                                                                                                                                                       ttATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACCGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGCCCCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGACCTGCACTGTAGGTCGGCTAAGACGCATAGCGTCGCATTCGACC‑ACCCATGCCCGACTCGCCCCTGACGCAACTTATTAGGc                                                                                     >  2:1091480/1‑288 (MQ=14)

                                                                                                                                                                                                          acAGGCGTACGAGGTGCAGTGAACGTTGGGGCGAAAAAGACGGATGGAGTCGCATCCGACAAGCGATGCCTGATGCGCCTCTGCCGGGACTTATCAGGCCTAGGGGGTGCACGTGACTGGAGGTCGTATAAGACGCATAGCGTCGCATCCGACAACCAATGCCGGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGAGGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTAc                                                                                 <  2:926688/288‑1 (MQ=17)

                                                                                                                                                                                                              acgtacgcagTGGACGGAACTGTAGGTCGGATAAGACGGATGGGGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGTGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACggg                                                                              <  2:178762/279‑1 (MQ=255)

                                                                                                                                                                                                               cctcggAGGGGCACTGAACGGTGGGTCGGATAAGACGAGTGGGGTCGCAGCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGt                                                                            <  2:1194575/284‑1 (MQ=255)

                                                                                                                                                                                                                cTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtg                                                                           >  1:987492/1‑289 (MQ=255)

                                                                                                                                                                                                                  aCGATGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaa                                                                         <  1:1134222/289‑1 (MQ=255)

                                                                                                                                                                                                                  aCGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaa                                                                         <  1:685147/289‑1 (MQ=255)

                                                                                                                                                                                                                   cGAGGTGCACTGAACTGTATGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCAGATAAGACGCATACCGTCGCATCCGACC‑ACCCATGGCTGATCCGCCTCTGACTCGACTCTGCAGGCCACAGGGGtgaac                                                                        >  2:342065/1‑289 (MQ=14)

                                                                                                                                                                                                                        tGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGtgaactga                                                                     >  1:228869/1‑287 (MQ=255)

                                                                                                                                                                                                                        tGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGCGCACTGAACTGTAGGTCGGATAAGACGCATAGCATCGCATCCGACA‑ACCGATGCCTGATGAGCCGCGTAACCGACTTACCAGGCCTACGGGGtgacatgaag                                                                   >  2:201394/1‑288 (MQ=14)

                                                                                                                                                                                                                        tGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGACCAACGGGTGCACCGAACTGTAGGTCGGATAAGCAGCATCGCCTCGCATCCGACACACCCA‑GCCTGCTGCGCCGCCGACACGACGTATCCGACCAAACGGGtgaactgaa                                                                    >  2:398232/1‑288 (MQ=14)

                                                                                                                                                                                                                               aaCCGTAGGTCGGATCACACGGATGGCGTCGCAGACGACAACCGAGGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAgg                                                             <  2:1075719/288‑1 (MQ=255)

                                                                                                                                                                                                                                aCTGAAGGTCGGATAAACGGGGTGGGCTAGCATCCGGAAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGt                                                            <  2:129878/288‑1 (MQ=14)

                                                                                                                                                                                                                                     atgTCGGATAAGACGGATGGCGGCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGAt                                                       <  1:711111/286‑1 (MQ=255)

                                                                                                                                                                                                                                      ggTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACa‑a                                                                                                                          >  1:200576/1‑220 (MQ=35)

                                                                                                                                                                                                                                        tCGGATAAGCCGGGGGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGGCGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAg                                                    <  2:656426/288‑1 (MQ=255)

                                                                                                                                                                                                                                        gagaaaagatgtgacgccGTCGCATGCGAGAACAGAGCGCTGCTTCGCGGCTGACGCGACTGATGAGGCGTACGGGGTGCACTGAGCTGTAGGTCGGATACGCCCGATGGCTTCGCATGCGACAACGGATGCCTGATGCGCCGCGGACGCGACTTA‑TGAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCTACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAg                                                    <  2:600928/271‑1 (MQ=11)

                                                                                                                                                                                                                                           gATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTAAGGGGTGCACCGAACTGCAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATCCATGATGCGCCGCTAACCCCACTTACCATGCCTACGGCGTGCACTGCACTCGCAGTCGGATAAGACg                                                 >  2:1008568/1‑288 (MQ=14)

                                                                                                                                                                                                                                                          gtggcaccatcagaCCGATGTCGGATGCGCCGCTGACGTGACTGATCAGGCCTACGGGGTGCCCTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGAGGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATcc                                  <  2:800197/275‑1 (MQ=255)

                                                                                                                                                                                                                                                                       aCCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGa                    <  1:1015252/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                           aTGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTg                                                                                                                                                                    >  2:905323/1‑142 (MQ=34)

                                                                                                                                                                                                                                                                             gCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTGGCATCCGACAACCGATGCCGGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTa               <  2:439059/288‑1 (MQ=255)

                                                                                                                                                                                                                                                                               cggaggggcggcTGAGGTGACGTTGCAGGCATACGGGGGACACTGAACTGGAGGTCGGATAAGAGGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCGGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACa             <  2:300190/278‑1 (MQ=14)

                                                                                                                                                                                                                                                                               cTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACa             <  1:1008568/288‑1 (MQ=255)

                                                                                                                                                                                                                                                                                      accgcTGACTTGACTTATAAGGCCTAAGGGGTGCACTGAACTGTAGGTCGGATAAGGCTGATAGCGTCGCATACGACAACCGATGCCTGATGTGCCGCTGACGCGACCTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCGGATGCGCCGCTGACGCGACTTATCAGGCGTACGGGGTGAACTGAACTATAGGTCGGATAAGACGCATAGCGTCGCATCCGACAAAAATCTCGACCCTACAAATGATa      <  2:788119/287‑1 (MQ=14)

                                                                                                                                                                                                                                                                                        cgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGGACTGAGCTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGGCTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACa   <  2:94343/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                                        cgcTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑ACAGGCCTACGGGGt                                                                                                                                                                          >  2:749340/1‑123 (MQ=17)

                                                                                                                                                                                                                                                                                         gcTGACGCGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACaa  <  1:1091480/289‑1 (MQ=255)

GTCTACAACGTCCTGCGCCAGGAAGGCACGCAGAAAAGCGTCATCGACACCATGCAGACCCGCAACGAGCTGTACGAAAGCATCAACTACTACCAGTACGAAGAGAAGCTCGACAACCTGTTTGCCCGTAGCCAGGTGAAATAAAAAACGCCCGTTGATTGTATTCGACAGCCGATGCCTGATGCGTCGCTGACGCGACTTATCAGGCCTACGAGGTGCACTGAACTGTAGGTCGGATAAGACGGATGGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGTGACTTATCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACCGATGCCTGATGCGCCGCTGACGCGACTTA‑TCAGGCCTACGGGGTGCACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACA‑ACCGATGCCTGATGCGCCGCTGACGCGACTTATCAGGCCTACGGGGTGAACTGAACTGTAGGTCGGATAAGACGCATAGCGTCGCATCCGACAACAATCTCGACCCTACAAATGATAACAA  >  NC_000913/349429‑349998

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