Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,256,301	G→A	S117S (TCC→TCT)	cyuA ←	putative L‑cysteine desulfidase CyuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,256,301	0	G	A	97.4%	216.7 / ‑5.9	76	S117S (TCC→TCT)	cyuA	putative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (36/38);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCG  >  NC_000913/3256057‑3256539                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               cAAACGGGACTTCATTGACGAACGTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACt                                                                                                                                                                                                                                          <  2:69054/251‑1 (MQ=255)     cGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGTGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTcc                                                                                                                                                                                                                                      >  1:359597/1‑251 (MQ=255)      gggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTCACAGAGACTTTccc                                                                                                                                                                                                                                     >  1:514893/1‑251 (MQ=255)       ggACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTcccc                                                                                                                                                                                                                                    <  2:36021/251‑1 (MQ=255)             cATTGACGAACTTCAGGATCTAAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCACTCTGCCACACACGCCTGCTGGGTACACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCgccgcc                                                                                                                                                                                                                              >  2:613697/1‑251 (MQ=255)                  aCGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGAGCTTAACAGAGACTTTCCCCGCCGCcagcag                                                                                                                                                                                                                         <  1:195600/251‑1 (MQ=255)                         tCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCTTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCtt                                                                                                                                                                                                                    <  1:600743/249‑1 (MQ=255)                          cAGGATCTCAGCCAGCGGCGTTCTGGAAAGCCCCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATAGGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCTTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTgg                                                                                                                                                                                                                 <  1:376397/251‑1 (MQ=255)                            ggATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGc                                                                                                                                                                                                                >  1:53628/1‑250 (MQ=255)                             gaactcagcctgcgccgtcagggAAAGCCCCGTCAGCGGAGCCT‑TTGTCCAGCCCTCTGCCACACCCGCCTGCGGGGTAAACACCCCACCACCGTGCGTCTCGTTATGCACAATGTTGGTATCC‑CCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGAC‑TT‑GGGCGCTTGAGAACAGTATTTCATCGCA‑‑GGTTACC‑TGGATCTGAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCAt                                                                                                                                                                                                              <  2:73773/230‑1 (MQ=255)                             gATCTCAGCCAGCGTCGTTCTGGAAAGCACTGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGAAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGA‑CCGCCGACGCTGGTGACACACGCCCACTTCTCACCGTTACAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTcccc                                                                                                                                                                                                                                    >  1:576721/1‑229 (MQ=255)                             gATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑GGGCGCGTGAGAAGAGGATTTCATCGCA‑ATGTT‑CC‑TGGAGCTTAGCAGAGACGTTcccc                                                                                                                                                                                                                                    <  2:576721/229‑1 (MQ=255)                                 tcAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGc                                                                                                                                                                                                          >  1:77176/1‑251 (MQ=255)                                  cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACCCCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCgg                                                                                                                                                                                                           <  1:199396/249‑1 (MQ=255)                                       agcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGc                                                                                                                                                                                                    >  2:1249/1‑251 (MQ=255)                                          gtcgtTCTGGAAAGCACCTTCAGCGGAGACTCTTG‑CTCGCCCCCGGCCACCCACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTg                                                                                                                                                                                                 <  1:347630/251‑1 (MQ=255)                                                  ggAAAGCCCCGTCAGGGGAGACTCTTC‑CTCGCCCTCTGCCACACCCGCCTGCTGGGTAAACCCCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCCCTTCTCCCCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTg                                                                                                                                                                                         <  2:154150/251‑1 (MQ=255)                                                       gCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCt                                                                                                                                                                                    >  1:118908/1‑251 (MQ=255)                                                       acACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACACTGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCt                                                                                                                                                                                    >  2:28007/2‑251 (MQ=255)                                                           cgcaaGCGGAGCCTCTGT‑CTCGCCCTCTGCCACACCCGCCTGCGGGGTAAACCCCACCCCATCGTGCGGCCCGATATGCACAATGTTGGTAGGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCa                                                                                                                                                                                <  2:257465/247‑1 (MQ=255)                                                            gTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                                                                                                                                                                                 >  1:529912/1‑249 (MQ=255)                                                            gTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGATAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGAGACACACGCCAACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                                                                                                                                                                                 >  1:254800/1‑249 (MQ=255)                                                                 cGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACtt                                                                                                                                                                          >  1:318942/1‑251 (MQ=255)                                                                  ggtgacgctTG‑CTCGCCCTCTGCAACCCCCGCCCGCGGGGTAAACACCACCCCATCGCGCGTCCCGTGATGCACAATGTTGGTA‑GCTCCGCCGACGAGGGTGACACACTCCCGCTTCTCACCGTCACAGGC‑TT‑GGGCGCGGGAGAAGAGGATTCCAGCGCA‑AGGGT‑CC‑TGGCGCGTAACAGAGTCTGTTCGCGTCGTCTGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTAGCGTCTTTCATCACTTc                                                                                                                                                                         <  2:289323/244‑1 (MQ=255)                                                                    ggtctctTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCCCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGAT‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCa                                                                                                                                                                       <  1:52353/248‑1 (MQ=255)                                                                    gagaccctggcccGCCCTCCGCCACACACGCCTGGTGGGTAAACCCCTCCCCACCGGGCGGCTCAATATGCCCATGGGTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGCC‑TT‑GGGCGCGTGAGAAGAGGATTTCAGCGCA‑AGGTT‑CC‑TGGATCTAAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                                                                                                                                                                        <  2:128585/239‑1 (MQ=255)                                                                        tctTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc                                                                                                                                                                   >  1:608521/1‑251 (MQ=255)                                                                        tctTG‑CTCGCCCTCTGCCACACACGACTGCTGGGTAAACACCACACCATCGTGCGTCTAGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACAAGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTCACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc                                                                                                                                                                   >  2:563886/1‑251 (MQ=255)                                                                        tctTG‑ATCGCCCTCTGCCACACACGACTGCTGGGTAAACACCACACCATCGTGCGTCTCGCTATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGCGAAGAGGATTTCATCGCA‑CGGTT‑CC‑TGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGccc                                                                                                                                                                   >  2:611476/1‑251 (MQ=255)                                                                                     tctGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAtt                                                                                                                                                       >  1:99362/1‑251 (MQ=255)                                                                                     cctGCCACCCACTCCTTCTGGGTAAACCCCACCCCTACGTGCGTTTCGATATGCACAATGTGGGTAGGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTCCCAGAT‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCAtt                                                                                                                                                       <  2:234680/250‑1 (MQ=255)                                                                                          cacacaCGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCAc                                                                                                                                                  >  2:376340/1‑251 (MQ=255)                                                                                                     gctgGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGcccc                                                                                                                                        <  1:552245/250‑1 (MQ=255)                                                                                                       tgGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTc                                                                                                                                                                         >  2:526765/1‑215 (MQ=255)                                                                                                       tgGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTATGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTc                                                                                                                                                                         <  1:526765/215‑1 (MQ=255)                                                                                                         ggtaacacccaccccatctggCGTCTGGATATGCACAATTTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTGCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAgcg                                                                                                                                   <  2:629534/231‑1 (MQ=255)                                                                                                                   cccccaTCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAT‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTACCAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCAGCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAAt                                                                                                                         <  1:613697/248‑1 (MQ=255)                                                                                                                    caccaTCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATc                                                                                                                        >  1:34946/1‑251 (MQ=255)                                                                                                                                       atGCACAATTTTGGTATGC‑CCTCCGTCGATGGTGACCCACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGTTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCg                                                                                                     <  1:94308/251‑1 (MQ=255)                                                                                                                                         gCACAATGTGGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTg                                                                                                   <  1:28007/251‑1 (MQ=255)                                                                                                                                            ccattgttgttaggtccgcggACGATGGTGACCCACGCCCCCTTCTCACCGTTC‑AGACCTT‑TGGCGCGTGAGAAGAGGTTTTCACCGCA‑AGGTT‑CC‑TGGATCTAAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCGAACGCCCCCCGCGCCGCCGCAATCGGCAGCCCCACCATGCCCGTGcc                                                                                                 <  1:147695/236‑1 (MQ=255)                                                                                                                                             aaTGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCCTTCCCGTGCCgg                                                                                               >  2:132709/1‑251 (MQ=255)                                                                                                                                                  tGGTATGC‑CCGCCGACGATGGTGACAAACCCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTACCAGAGACTTTCCCCGCCGCCAGCAGTGCTTTCGCATCGGCAATTGCCTGCGCTGCTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTCACGCCCCCAGCGCCGCCGCAATCGGCAGACCCACCATTCCCGTCCCGGGAACg                                                                                          >  2:413286/1‑251 (MQ=255)                                                                                                                                                   ggtagcc‑ccGCCGCCGGTGGTGACACACGCCCACTTCCCCCGGTGCCAGTC‑TT‑TGGCGCGTGAAAAGAGGATTTCTTCGCA‑AGGTT‑CC‑TGGTTCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCATCCCGGCGGTGGCATTTCCACCTAACGCCCCGAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACgg                                                                                         <  1:563886/245‑1 (MQ=255)                                                                                                                                                    gTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACggg                                                                                        >  1:495119/1‑250 (MQ=255)                                                                                                                                                       tGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGAc                                                                                     >  1:47296/1‑251 (MQ=255)                                                                                                                                                               cgacgaTGGTGCCACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGcccc                                                                                                                                        <  2:348798/193‑1 (MQ=255)                                                                                                                                                               cgacgaTGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGcccc                                                                                                                                        >  1:348798/1‑193 (MQ=255)                                                                                                                                                                  cgaTGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATGCCCGTGCCGGGAACGGTGACGCCCAGACCg                                                                           <  1:500112/251‑1 (MQ=255)                                                                                                                                                                       gTGACACACGCCCACTTCTCACCTTTCCAGTC‑TT‑GGGCGCGTGAGAATAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGTTGACGCCCAGACCGttctt                                                                      <  2:53628/251‑1 (MQ=255)                                                                                                                                                                        tGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGttcttc                                                                     >  1:613470/1‑251 (MQ=255)                                                                                                                                                                           cacacGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                  >  1:246396/1‑251 (MQ=255)                                                                                                                                                                            acacGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                                  >  1:288272/1‑250 (MQ=255)                                                                                                                                                                               cGCCCACTTCTCACCGTGCCAGTC‑TT‑GGGCGCGTGAGAAGAGGATTTCATCGCA‑AGTTT‑CC‑TGGATCTTAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGCCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTACCGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTTTTCATCAGAt                                                              <  1:166812/251‑1 (MQ=255)                                                                                                                                                                                  ccACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTg                                                           >  1:62440/1‑251 (MQ=255)                                                                                                                                                                                       tctcACCGTCCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCg                                                        <  1:436091/249‑1 (MQ=255)                                                                                                                                                                                        ctcACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGaa                                                      >  2:545509/1‑250 (MQ=255)                                                                                                                                                                                          cACCGTTCCAG‑C‑TTGTGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGaa                                                      <  2:608521/248‑1 (MQ=255)                                                                                                                                                                                            ccccgttcaatacttggtcgcgTGAGAAGAGGTTTTCA‑CCCAGAGTTT‑CC‑GGGATCTAAACAGAGCCTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATGGCCGGCGCTTTTGCGTCTTTCAGCACTGCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCAGCAGATTTGGCgtatcc                                                   <  2:33745/233‑6 (MQ=255)                                                                                                                                                                                                 ccAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCtt                                            >  1:182940/1‑251 (MQ=255)                                                                                                                                                                                                      c‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGTTT‑CC‑TGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTGTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCg                                                                                                                                                                  <  2:597170/128‑1 (MQ=255)                                                                                                                                                                                                      c‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCg                                                                                                                                                                  >  1:597170/1‑128 (MQ=255)                                                                                                                                                                                                        tt‑tGGCGCGTGAGAAGAGGTTTTCATCGCA‑AGTTC‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTTTTGCGTCTTTCAGCTCTTCCAGCCCGGCGTGGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGTCGTCTac                                        <  2:613470/249‑1 (MQ=255)                                                                                                                                                                                                         ttgggCGCGTGAGAAGAGGATTTCACCGCA‑‑GGTT‑CCCTGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGGTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTaca                                       <  2:194664/247‑1 (MQ=255)                                                                                                                                                                                                            ggCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGt                                    <  1:376340/250‑1 (MQ=255)                                                                                                                                                                                                                  tgaaaagaGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACGGAGACTTTACCCGCCTCCAGCAGTGCTTTAGAATCGGCAATTGCCTGCGCTGTTGCGACTTTCCGCACTTCCCGCACGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCCGCCCCACCATTACCGTGCCGGACACGGGGACGCCCAGACCGTTCTTCCTCCGCTTTGGCGAAACCCAGGCTTCTACACGTTCAACCg                             >  1:365260/5‑251 (MQ=255)                                                                                                                                                                                                                   gagaAGAGGTTTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGCCTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCgg                            <  1:551495/251‑1 (MQ=255)                                                                                                                                                                                                                     aaagtggATTTCACCGCA‑AGGTC‑CC‑TGAATTTTAACAGATCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCTTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGAc                          <  2:580403/246‑1 (MQ=255)                                                                                                                                                                                                                          ggATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGCAACCCAGGCTTCTACACGTTCAACCGGACCTTcc                     >  1:144094/1‑251 (MQ=255)                                                                                                                                                                                                                           gATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTcc                     >  2:14265/1‑250 (MQ=255)                                                                                                                                                                                                                                cATCGCA‑AGGTT‑CC‑TGGATCTTACCAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCCCTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTc                <  1:91117/250‑1 (MQ=255)                                                                                                                                                                                                                                        aGGTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCgcag        <  1:1249/251‑1 (MQ=255)                                                                                                                                                                                                                                          gTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCgcagca      <  1:132709/251‑1 (MQ=255)                                                                                                                                                                                                                                          gTT‑CC‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCATCTACACGTTCAACCGGACCTTCCAGGTCTGCCgcagc       >  2:607999/1‑250 (MQ=255)                                                                                                                                                                                                                                            t‑CC‑TGGATCTTAACAGAGACTTTCCCCGACGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAc    >  1:451459/1‑251 (MQ=255)                                                                                                                                                                                                                                               c‑TGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAAccc  >  2:503412/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               CAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTG‑CTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCACAATGTTGGTATGC‑CCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGAC‑TT‑TGGCGCGTGAGAAGAGGATTTCATCGCA‑AGGTT‑CC‑TGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGCCGCAGCAACCG  >  NC_000913/3256057‑3256539

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