BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,135,441	T→G	11.5%	T311P (ACC→CCC)	pitB ←	phosphate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,135,441	0	T	G	11.5%	288.1 / 2.1	113	T311P (ACC→CCC)	pitB	phosphate transporter
Reads supporting (aligned to +/- strand): ref base T (26/69); new base G (14/3); total (41/72)
Fisher's exact test for biased strand distribution p-value = 2.78e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 113 total aligned reads displayed.

GCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATC  >  NC_000913/3135206‑3135673

cccaggagtatCCCCACCGGGGCGTACTCAATGGTGCTTAACAAATCGCTGGGAAGTTTTGTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGTTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTAtt                                                                                                                                                                                                                                        <  2:548122/230‑1 (MQ=255)

        tttCCACCCCTGGGCGGACCAAAGGGTGCTTACCATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCGTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGGTTTAACGCGCGTAATAGCATCAAAGGTATTTGCCGGAt                                                                                                                                                                                                                                <  2:856186/237‑1 (MQ=255)

        gatCCACCCCGGGGTGTTCTTAATGGTGCTTAACCTATCGCGGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGGTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGAt                                                                                                                                                                                                                                <  2:474394/239‑1 (MQ=255)

               cccGGAGCGTACTCAATGGTGCTTAACATTTCGCTGCGAAGTTTTTTCAGGCGGTTTTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTg                                                                                                                                                                                                                         <  1:700158/238‑1 (MQ=255)

                 cGGTGCGTACTCAATGGTGCTTAACATTTCGCTACGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGaa                                                                                                                                                                                                                       <  1:469186/239‑1 (MQ=255)

                   gAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGGTACATTCCATATTGCCTGGGCGCAT‑CCGTTTTACGCGCCCAAAAACATACAAGGGATTTTCCGGAGGAAAggaaaa                                                                                                                                                                                                                      >  2:68431/1‑233 (MQ=255)

                    aGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGGAACTTTCCATATTGCCTGGCGGCAT‑CCCTTTAAAGCGCGCAATAGCATCAAAAGTATTTTCCGGCTGCCCGTGaaaa                                                                                                                                                                                                                     >  1:1035969/1‑237 (MQ=255)

                    aGCGTACTCAATGCTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACCTTCCCAATTTCCTGGCAGCATCCGTTTTAATG‑ACGCAATTGCCACCAAAGTACTTGCCCGAGGCCAGTTaaaa                                                                                                                                                                                                                     >  2:883289/1‑237 (MQ=255)

                            cAATGGTTCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGGCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACt                                                                                                                                                                                                            <  1:1249841/239‑1 (MQ=255)

                               tGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTTGCAGTTTCGCTAGCCTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGc                                                                                                                                                                                                         <  1:99659/239‑1 (MQ=255)

                                ggtacctAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACGCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCg                                                                                                                                                                                                        <  1:1195351/233‑1 (MQ=255)

                                  tGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCAAATTGCCAGTCAGAGT‑AGCTTTTACGCGCGCAAAAGAATCAAAGGGATTTGCCGGGTGACAGTGAAAATCTGTTACATgatg                                                                                                                                                                                                       >  1:87948/1‑235 (MQ=255)

                                  ggCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGt                                                                                                                                                                                                       <  1:966831/237‑1 (MQ=255)

                                   gCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCCGGCAGCAT‑CGTTTTTACCCGCCCCAAAGGATCAAAGGGATTTTCCGGATGGCAGTGAAACCCTGTTACtttgggg                                                                                                                                                                                                      >  1:534693/1‑233 (MQ=255)

                                    cTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCCAAGGGATTTTCCGGATGACAGTGAAACTCTGTTACTTTCGttc                                                                                                                                                                                                     >  1:752449/1‑236 (MQ=255)

                                     ttAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTTTTTACGGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                                                                                                                                                                                                   <  2:1012485/239‑1 (MQ=255)

                                     ttAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                                                                                                                                                                                                   <  1:1102946/239‑1 (MQ=255)

                                       aaCATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTTCCCGGCAGCAT‑CGTTTTTAAGCCCGCAATAGCATCAAAGGGATTTGCCGGATGTAAGTGGAACCTCGTTAATTCCGTCCCATc                                                                                                                                                                                                 >  1:101199/1‑239 (MQ=255)

                                         cATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACGTTCCCTATTTCCTGGCAGCAA‑CGCTTTAACGGCCGCAAAAGCATCCAAGGAATTTTACGGATAGCAGGGAAAATCTGTTACTTGCGTGccttacg                                                                                                                                                                                               >  1:1161994/1‑234 (MQ=255)

                                            atCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCCGCAT‑ACTTTTTAAGACCGCAATACCCTCAAAAGGATTTGCCGGTTGAACGGGAAAATCTGTTACTTTCGTGCCacacgacc                                                                                                                                                                                            >  2:854218/1‑232 (MQ=255)

                                             taagagaccggTTTTTTCAGGAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGGTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGa                                                                                                                                                                                           <  1:685675/229‑1 (MQ=255)

                                              cGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTTCCTGGCAGCAT‑CGGTTTAACCCGCGCAAAAGCATCCAAGGGATTTTCCGGATGACCGTGAAACTCTGTTACTTGCCTGCCCCCAGCCGAt                                                                                                                                                                                          >  1:813005/1‑239 (MQ=255)

                                              cGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTGCCTGGGAGCAT‑CGTTTTAAAGCGCGCACAAGCATCAAAGGGATTTGCCGGATGACAGTGAAAACCTGTTACTTTCGGGCCATCAGTccgt                                                                                                                                                                                          >  1:504197/1‑236 (MQ=255)

                                                 tgctatgTTTTTTACGCGGGTTTTGGTCTTCTTTCCTGCCGCCTGGCCGTTTCGCTAGCTTCGCGGGGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgc‑t                                                                                                                                                                                      <  2:585543/233‑1 (MQ=255)

                                                 tGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACATCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCCTGACCTTCCCTATTTCCTGGCAGCCC‑ACCTTTTACACGCGCCAAAGCATCTATGCTAGTTGCCGGATGAAAAGGAAACACCGCTACTCCCGTACCAACACTCCAtgca                                                                                                                                                                                       >  2:558780/1‑238 (MQ=255)

                                                  gggAAGTTTTTGCAGCGGGGTCTGGTCTTCTTTACGGACGCCTGGCAGTTTCGGTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgc‑tg                                                                                                                                                                                     <  1:364021/237‑1 (MQ=255)

                                                      aGTTTTTTGAGTTGGTTATGGTCTTTTTTACTGACGCCTGGCAGTTTCACTAGCTTCGGGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCGCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGTCAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAgg                                                                                                                                                                                 <  1:305143/239‑1 (MQ=255)

                                                                             ttttttGCTGGCCCCTGGGAGGTTCGCGCGGTTCGGGGGGGTATCAGAGATGCACTGCATAAGGGGGCCAAGCTGGCTGCGTTGAGTCTGCCTTAACGGCTCGTAACTTGCCATATTGCCTGGCAGCAT‑CGGTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCa                                                                                                                                                          <  2:417822/236‑1 (MQ=255)

                                                                               cttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGTCCGGCCGCAT‑CGTTTTTAAACGCGCCACAGCATCAAAGGGATTTTCCGGATGACAGTTGAAATCTTTTACTTGCGTGCCCTCAGCCGAGGC‑TGCCGGCAAAGGGGGGTcccaccacataaa                                                                                                                                                         >  2:517902/1‑227 (MQ=255)

                                                                                            tGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTTCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTCCGGCGGt                                                                                                                                            >  1:54550/1‑238 (MQ=255)

                                                                                              gCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTTCCTGGCAGCAT‑CGTTTTAACGCGCGCAAAAGCATCAAAGGGATTTTCCGGATGACAGGTAAACACTTTTTATTGCGTGCCATCCAACGAGTC‑TGCCGGCAAAGGGGGGTCCATACCAATCAACCTCTTCGGCGGTTCa                                                                                                                                         >  2:1251427/1‑239 (MQ=255)

                                                                                                 ttttCCCTAGTTTCCGGGGGGTATCCGGAATTCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCTCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGa                                                                                                                                      <  2:478918/238‑1 (MQ=255)

                                                                                                   ttccctggtttCGCGGTGTTTTCAGAGATGCACAGCTTAAGGCGGCGCGGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATg                                                                                                                                    <  2:24537/230‑1 (MQ=255)

                                                                                                   ttCGCTAGCTTCGCGGAGGTTTCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATg                                                                                                                                    <  1:656882/239‑1 (MQ=255)

                                                                                                   ttCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAAAGCGCGCAATAGCATCAAAGGTATTTTCCGGGTGACAGTGAAACTCTGTTATTTTCGTGCCATCAATCGATGC‑TGCAGGCAATAGGGGTTCCATCGCAAACAACTTTTGCGGcgagtacgggag                                                                                                                                    >  2:1325071/1‑228 (MQ=255)

                                                                                                     cGCTGGCTTCGGGGAGGTATCAGAGATGCACAGAATAATGCGGCGCAGCTGGCTGCGCTGACTCACGTTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCt                                                                                                                                  <  2:54550/239‑1 (MQ=255)

                                                                                                        tAGCGTCGCGGAGGTAACAAGGATGCCCAGGATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGtt                                                                                                                               <  1:1122062/239‑1 (MQ=255)

                                                                                                         aGCTTCGGGGGGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTg                                                                                                                              <  1:93061/239‑1 (MQ=255)

                                                                                                              cgcgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGt                                                                                                                         <  1:287419/239‑1 (MQ=255)

                                                                                                                   aGGTATCAGGAAGGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAATCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGa                                                                                                                                      <  1:780760/221‑1 (MQ=255)

                                                                                                                     gTATCAGATATGCAACGCATAATGCGGCCCAGCTTGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCAGCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGtt                                                                                                                   <  2:981128/238‑1 (MQ=255)

                                                                                                                                 cccaGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACGTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTaa                                                                                                       <  1:944500/236‑1 (MQ=255)

                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                                                                                                                                                                                                   >  1:511355/1‑142 (MQ=255)

                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTTGCCGCAT‑CGTTTTAAAGCGCGCAATAGCATCCAAAGTATTTGCCGGATGACAGGGAAACTCTGTTACTTGCGTGCCACCAGTCGGTGCATGC‑GGGAAAAGAGGGTCCATCGCAATCAACATCTTCCGCAGTTTACGAGGCCGTTGCGGGTAGGGTTCGGAGTTGGTAAAggaca                                                                                                 >  2:1066973/1‑236 (MQ=255)

                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGGATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCCTCGCAATCAACTTCTGCGGCCGTTCAGGATGCTGTTGCAGGTAGTGTTCCCAGTTGGGAAAGGCa                                                                                                  >  1:1080076/1‑238 (MQ=255)

                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAC‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                                                                                                                                                                                                   <  2:511355/142‑1 (MQ=255)

                                                                                                                                       aTAATGGGGCCCAGGTGGCTGGTGGTACTAACGCTTAGCGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATc                                                                                                <  2:361870/239‑1 (MQ=255)

                                                                                                                                        tAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATcg                                                                                               <  1:517902/239‑1 (MQ=255)

                                                                                                                                            gtctcccAGCTTGTTGCGCTTCCTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGGAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGGCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcggg                                                                                           <  2:375253/233‑1 (MQ=255)

                                                                                                                                             cggcgCAGCCGGCTTCGCTGACTCCCCTTTACCGGCTGCTAAATTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGGGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggt                                                                                          <  2:1261051/239‑1 (MQ=255)

                                                                                                                                              ggcgCAGCTGGCTGCGCTTCCTTCCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                         <  2:1080076/239‑1 (MQ=255)

                                                                                                                                               gcgcAGCTTGCCGCGCGTACTCCCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                         <  1:607718/238‑1 (MQ=255)

                                                                                                                                               gcgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                         <  1:68431/238‑1 (MQ=255)

                                                                                                                                                cgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAATCTGTTACTTGCGTGCCATCCTTCGCTGC‑TGCAGGCCATGGAGGTTCCATCGCAATCAACTTCTGCGGCCGTTCAGGATGCTGTTGCAAGTGTTGTTCGTAGTTGGTAACGGCATCCcgggtacg                                                                                       >  1:285144/1‑239 (MQ=255)

                                                                                                                                                  cAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAgg                                                                                                                          <  1:597624/202‑1 (MQ=255)

                                                                                                                                                    gctggctgCGCTGGCTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggtacgggta                                                                                   <  1:1325071/239‑1 (MQ=255)

                                                                                                                                                                  cTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATGGCCGGa                                                                     <  1:1293425/239‑1 (MQ=255)

                                                                                                                                                                  cTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                     <  1:1251427/239‑1 (MQ=255)

                                                                                                                                                                   tCACGCTTAACGGTTCTTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGAc                                                                    <  1:1138426/239‑1 (MQ=255)

                                                                                                                                                                   ctaCGCTTAACGGCTCGTAACTTGCCATTGTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCACCAAAGGGAGTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                     <  1:70117/236‑1 (MQ=255)

                                                                                                                                                                    ctactttAACGGCTCGTAACTTTCCATATTGCCTGGCCGCAT‑CGTTTTAACGCGCGCAATAGCAACAAAGGTAGTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTTCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACg                                                                   <  1:540978/234‑1 (MQ=255)

                                                                                                                                                                            aCGGCTCGTAACTTTTCATATTTCCTTGCCGCAT‑CGTTTTAACGCGCCCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAATCTGTTACTTTCGTGCCATCAATCGATTC‑TGCCGGCAATGGGGGTTCCATCGCAATCCACTTTTGCGGCAGTTCAGGAATCTGTTGCAGGTAGAGTTCGACGTTGGTAACGGCACCGCGGGGTCGGGGAATTTCATACCCGGACGCAtttcaa                                                           >  1:370145/1‑235 (MQ=255)

                                                                                                                                                                              ggCTCGGAACTTTCCATATTGCCTGGCAGCAT‑CGTTTGAACGCGCGAAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACGTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATAtt                                                         <  1:883289/239‑1 (MQ=255)

                                                                                                                                                                                   gTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGGGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTgacga                                                    <  1:220688/239‑1 (MQ=255)

                                                                                                                                                                                     aaCCTTCCATTTTGCCTGGCCGCCT‑CGGTTTAACGGCGGCAATAGGATCCAAGGGAATTGCCGGATGGCAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGG‑TGCAGGCAAGGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGacgacg                                                  <  1:678718/239‑1 (MQ=255)

                                                                                                                                                                                           ccATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                     <  1:835550/214‑1 (MQ=255)

                                                                                                                                                                                             atatTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagc                                          <  1:219483/239‑1 (MQ=255)

                                                                                                                                                                                             atatTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagc                                          <  1:991944/239‑1 (MQ=255)

                                                                                                                                                                                              atttGCCTGGCAGGAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagca                                         <  1:383010/237‑1 (MQ=255)

                                                                                                                                                                                                tttccgtgcagcaT‑CGTTTTAACCCGCGCAATTGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagc                                          <  1:404848/229‑1 (MQ=255)

                                                                                                                                                                                                  gCCTGGGCGCAT‑GGGTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGGTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAgggg                                     <  2:1090623/239‑1 (MQ=255)

                                                                                                                                                                                                  gCCTGGCAGTAT‑CGTTTTAACGCGCGAAATAGCCACAAAGGGATTGGCCGGGTGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAgggg                                     <  2:50911/239‑1 (MQ=255)

                                                                                                                                                                                                  gCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGCTGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCCGGTAGGGTTCGAAGTTGGCAACGGCATTCCGGGTACGGGTACTTTCATAGCCGGACGCATTATTATTGAGCACGAACCCCGCAgggg                                     >  1:611019/1‑239 (MQ=255)

                                                                                                                                                                                                      ggCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCaa                                  >  1:1104051/1‑238 (MQ=255)

                                                                                                                                                                                                        cagcaT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATc                                >  1:65154/1‑238 (MQ=255)

                                                                                                                                                                                                           caT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAc                                                                                                                                                                                                                     <  1:23399/55‑1 (MQ=255)

                                                                                                                                                                                                           caT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAc                                                                                                                                                                                                                     >  2:23399/1‑55 (MQ=255)

                                                                                                                                                                                                            ct‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATccc                              <  2:752449/235‑1 (MQ=255)

                                                                                                                                                                                                             t‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGCGTTCGAAGTTGGTAACGGCATCGCGGGGACGGGTAATTTCATAGCCGGACGCCTTCATATTTCCGACGAAGCCAGCAAGGGCAATCCccacc                          >  1:343141/1‑239 (MQ=255)

                                                                                                                                                                                                             ctgtttttAACGCCCGCAATAGCATCCAAGGTTGTGGCCGGTGTGCAGGGAAACTCTGTTCCTGGCGGCCCTTCCGTCGGTGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCccac                           <  2:269785/235‑1 (MQ=255)

                                                                                                                                                                                                             actgTTTAAACCCGCGAAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGGGCCAGCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTGCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCccac                           <  2:285144/236‑1 (MQ=255)

                                                                                                                                                                                                                 ttttAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACATGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGt                       <  1:473656/239‑1 (MQ=255)

                                                                                                                                                                                                                  tttAACGCGCGCAATTGCATCAAAGGTTTTTGGCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTa                      <  2:534693/239‑1 (MQ=255)

                                                                                                                                                                                                                  tttAACGCGAGCAATAGCCTCAAAGGTATTTGCCGGATGACAGTGAAATGTTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTa                      <  1:311124/239‑1 (MQ=255)

                                                                                                                                                                                                                      ccccgcgcAATAGCATCAAAGGTATTTGCCGGAGGACAGTGAAACTTTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGGAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAg                  <  2:1144226/236‑1 (MQ=255)

                                                                                                                                                                                                                      aCGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCGGCAGGGGCAATCCCCACCAGTACCa                   >  1:1144226/1‑238 (MQ=255)

                                                                                                                                                                                                                           cgcAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTTTTCGAAGTTGGTAACGGCATCCCGGGTACGGGGAATTTCATAGCCCGACGCATTCATATTTACCACGAAGCCAGCAAGGGCAAACCCCACCAATACCaaccata             >  2:412959/1‑233 (MQ=255)

                                                                                                                                                                                                                           cgcAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCCTACTGACGACGAAGCCAGCAGGGGAAATCCCCACCAGTACCACCAtt              >  1:202767/1‑238 (MQ=255)

                                                                                                                                                                                                                           cgcAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCACCAtt              >  1:781300/1‑238 (MQ=255)

                                                                                                                                                                                                                           cgcAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGGGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCCCCAGTCCCCGCCtt              >  1:52309/1‑238 (MQ=255)

                                                                                                                                                                                                                                     tCAAAGGTTTTTTCTGGATGACATTGAAACTCTGTTATCTGCGTGCCATCAGGGGATGC‑TGAGGGAAAATGAGGTTTCATCGGAATCAACGTACGCGGCAGTTCAGGATGTTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt   <  2:370145/239‑1 (MQ=255)

                                                                                                                                                                                                                                     tCAAAGGTATTTGGCGGATGACAGTGAAACTCTGTTACTTGCGTTCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt   <  1:505810/239‑1 (MQ=255)

                                                                                                                                                                                                                                     tCAAAGGTATTTGGCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt   <  1:1229772/239‑1 (MQ=255)

                                                                                                                                                                                                                                      cAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTTACGACGAAGCCAGACGGGGAGATCCCCACCAGGACCAGCATTCCCAggacact   >  1:693344/1‑233 (MQ=255)

                                                                                                                                                                                                                                        aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                     >  1:457005/1‑170 (MQ=255)

                                                                                                                                                                                                                                        aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                     <  2:457005/170‑1 (MQ=255)

                                                                                                                                                                                                                                        aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATc  >  1:765069/1‑237 (MQ=255)

                                                                                                                                                                                                                                        aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGAAATTTCATAGCCGGACGCATTCATATTGACGAAGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCCGGCCGATc  >  1:888180/1‑237 (MQ=255)

                                                                                                                                                                                                                                        aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGGACGGGTAATTTCATAGCCGGACGCATTCATATTTACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCCGCATTACCAGGCCGATc  >  1:838691/1‑237 (MQ=255)

                                                                                                                                                                                                                                        aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGGGTTCGAAGTTGGTAACGGCACCGCGGGTACGGGTAATTACATAGCCGGACGCATTAATATTGACGCCGAAGACAGCAAGGTCAAACCCCCCCCGAACCGGCATTAcacgccagacc  >  2:965445/1‑227 (MQ=255)

GCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCAT‑CGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGC‑TGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATC  >  NC_000913/3135206‑3135673

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