BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	1,568,023	G→T	14.3%	Q278K (CAA→AAA)	dosC ←	diguanylate cyclase DosC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	1,568,023	0	G	T	14.3%	344.8 / 8.3	120	Q278K (CAA→AAA)	dosC	diguanylate cyclase DosC
Reads supporting (aligned to +/- strand): ref base G (32/69); new base T (16/3); total (48/72)
Fisher's exact test for biased strand distribution p-value = 4.52e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

CTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTA  >  NC_007779/1567789‑1568261

cTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑Ga                                                                                                                                                                                                                                               <  1:1547895/238‑1 (MQ=255)

cTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑Ga                                                                                                                                                                                                                                               <  1:893668/238‑1 (MQ=255)

 tttACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTTC‑GA‑Ta                                                                                                                                                                                                                                            <  2:674531/239‑1 (MQ=255)

 tttACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑Ta                                                                                                                                                                                                                                            <  1:837844/239‑1 (MQ=255)

  ttACGCATAATTTCATCACCAGTGTTATGGCGCCACGTATCGTTAATCTCTTTGAATTTACCGACGTCAATAATCAGCGCTGACAGCTGTGTACCGTTCCGGTTGGCATGTGAAATTTCGCGTTTGCAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATCTGC‑GA‑tcc                                                                                                                                                                                                                                           <  2:1839824/239‑3 (MQ=255)

  ttACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TAc                                                                                                                                                                                                                                           <  2:881783/239‑1 (MQ=255)

    aCGTAGAATTTCATCGCCAGTGTTAGGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CCT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACgg                                                                                                                                                                                                                                         <  2:1084647/239‑1 (MQ=255)

     cGCAGAATTTCATCACCAGTGTTATGG‑CCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTa                                                                                                                                                                                                                                       <  1:1873433/239‑1 (MQ=255)

                        tgtTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATT‑CCTACCCACTTCGTGGCGCGATA‑CTT‑CTTTAAACAATTCACGCAGCAAGGTAATAATTTTC‑GA‑TACGGTATTTCTTTTCTTTAATaaa                                                                                                                                                                                                                     >  1:1538408/1‑239 (MQ=255)

                                   cgcttATCGTTTATCTCTTTGAGTTTATCAGGTTCAATAATCAGCACTGAGAGCGGTGTGCCGGTCCGGTTGGCATTCACAGTTTCACGTTTGAAGATATTCGGTAGGAAGCAGCGGTTAAGTTATTTCGTCACTACATC‑CATGCCGTCTTCGTGGCGCGATG‑CTT‑CTTCAAACTATTCACGCAGCAGGATAATAATTTTC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCCCCCg                                                                                                                                                                                                          <  2:1141296/235‑1 (MQ=255)

                                   cACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCg                                                                                                                                                                                                          <  1:1302889/239‑1 (MQ=255)

                                   cACGTATCATTGATATCTTTGAATATATAAACGTCAATAATCAGCACTGACAGCTGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGCTAGACGGTATGAAACGGCGTTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGTCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTTTAATAGAAATTTCACCCg                                                                                                                                                                                                          <  2:1746840/239‑1 (MQ=255)

                                    aCGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGGACATC‑CCTACCCACTTCGTGGCGCGATA‑CTT‑TTTCCAACAATTCACCCCGCAAGGGAATAATTTTC‑GA‑TACGGTATTTCTTATCTGTAATAAAAAATTTAcccc                                                                                                                                                                                                          >  1:501334/1‑237 (MQ=255)

                                     cGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGGACATC‑CATAACGACGTCTTGGGGCGATA‑CTT‑TTTCAAACAAATCACGCAGCAAGGTAATAATTTTC‑GA‑TACGGTATTTTTTATCTGTAATAAAAAATTCACCCCCa                                                                                                                                                                                                        >  1:587532/1‑239 (MQ=255)

                                     cGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCa                                                                                                                                                                                                        <  1:372472/239‑1 (MQ=255)

                                     cGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCa                                                                                                                                                                                                        <  2:515535/239‑1 (MQ=255)

                                          cGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGTAAGTAATTTTGTTAGTAAATC‑CATACCCACCTCGTGGCGCGATA‑CTT‑CTTTCAAAAATTTACCCAGCCAAGGAATAATTTTCTGA‑T‑CGGGATTTTTTTTCTTTAATAAAAAATTTACCCCCCGAcc                                                                                                                                                                                                    >  1:1002275/1‑237 (MQ=255)

                                              gATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTTCGTACATA‑CATTCCCAATTCGTGGCGCGGGA‑CTT‑TCTTAAAAAATTCACACAGCCAGGTAATAATTTGT‑CA‑TACGGTCTTTCTTTTTTTTAATTAAAAATTTACCCCCCGACTTCttt                                                                                                                                                                                               >  2:856177/1‑237 (MQ=255)

                                                 ctctTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAATAATTTTGTTAGTACATC‑CATAACCACTTCGGGGCGGGATA‑CTT‑CTTCAAACAAATCAAGCAGCAAGGGAATAATTTTC‑GA‑TACGGGATTTTTTTTCTTTAATAAAAAATTCACCCGCCGACTTCTGtttt                                                                                                                                                                                            >  1:926270/1‑237 (MQ=255)

                                                  tctTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTAGCGGTCCGGGTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTAtt                                                                                                                                                                                           <  2:871714/239‑1 (MQ=255)

                                                            tATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTGTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCGCGCAGCAAGGTAATAATTTGC‑GA‑TCCGGTATTTCTTATCTGTAATAACAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTg                                                                                                                                                                                 <  2:1682742/239‑1 (MQ=255)

                                                                  cGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGCTAAGTAAATTTCTCACTAACTT‑CATACCCACTTCTGGGCGCGGTACCTT‑TTTTCAAAAATTCAC‑CACCAAGGGAATAAATTTC‑TG‑TACCGGCTTTTTTTTTTCTAATAAAAATTTTACCCGCACACCTCTTTTATTTTAATTTCTTGTGttt                                                                                                                                                                           >  2:739557/1‑238 (MQ=255)

                                                                           tCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGtt                                                                                                                                                                 <  1:935756/239‑1 (MQ=255)

                                                                             aGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGGAATTTTGTCAGTACATC‑CATACCGACTTCGTGGGGCGGTA‑ATT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTttg                                                                                                                                                               >  1:1450714/1‑239 (MQ=255)

                                                                             aGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATC‑CATACCGACTTCGTGGCGCGGTA‑CTT‑CTTTAAACAATTTACGCAGCAAGGGAATAATTTGC‑GA‑TACGGTATTTTTTATCTGTTATAAAAAATTTACCCGCAGACCTCTGTTTTTTAAATTTTTTGTGTTACGCC‑TGGtttt                                                                                                                                                               >  1:610948/1‑238 (MQ=255)

                                                                             aGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGGTAGTAATTTTTTCAAGAACTT‑CAAACCGACCTTTCGGCGCGAAA‑CTT‑TTTTAAAAAATTTACTCCCAACGGGAATAATTTTC‑GAGAACGGTATTT‑TTTTCTTTAATAAAAAATTTAACCACAGAATTCTGTTATTTTAATTTCTGGTGTTTAGCATTGGttt                                                                                                                                                                >  2:1871638/1‑239 (MQ=255)

                                                                              gCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGGAATAATTTTC‑GA‑TACGGTATTTTTTATTTGTAATAAAAAATTTAACCCCAGACTTCTGTTTTTCAAATTTCTGGTGTTACGCA‑TGGtttta                                                                                                                                                              >  1:1928453/1‑237 (MQ=255)

                                                                               cACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTGCCCGTACATT‑CATAACCACCTTGCGGGGCGGTA‑ATT‑CTTCCAAAAAATAACCCAGCAAGGGAATAATTTGT‑GA‑TACGGTATTTCTTATTTGTAATAAAAAATTCACCCCCCGACTTTTCTTATTTAAATTTTTTTGGTTAAGCA‑TGGtttttt                                                                                                                                                             >  2:1795583/1‑236 (MQ=255)

                                                                               cACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAAAAATTTTGTCAGTACAAC‑CATACCCACTTCGTGGCGGGGTA‑CTT‑TTTTAAAAAAATTACCCAACAAAGGAATAATTTGC‑GA‑TACGGGATTTTTTTTCTTTTATAAAAAATTTACCCGCAGCACTCTGTTTTTTAAATTTTTGGTGTTAATCA‑TGGGttgat                                                                                                                                                             >  1:815991/1‑239 (MQ=255)

                                                                                   gACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTg                                                                                                                                                                                 <  1:914004/216‑1 (MQ=255)

                                                                                    aCAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAAATTCGTCAGTAAATC‑CCTACCGACTTCCTGGGGCGGTA‑GTT‑CTTTAAAAAAATCACGCAGCAAGGTAATAATTTTC‑GA‑TACGGTATTTCTTATCTGTAATAAAAAATTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCT‑TGGTTTGATTTaaa                                                                                                                                                        >  1:1398744/1‑239 (MQ=255)

                                                                                    aCAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAAATTTGTCAGGACATC‑CCTACCGACTTCCTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGAAAGGGAATAATTTTC‑GA‑TACGGTATTTTTTTTTTGTAATAAAAAATTCACCCCCAGACTTCTTTTTTTTAAATTTTTTGTGTTATCCA‑TTGttttttttaaa                                                                                                                                                        >  1:1203803/1‑231 (MQ=255)

                                                                                              aCCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAATTTCTTCCGTACATT‑CCTAACCACCTCTTGGCGGGCGA‑CTTACTTCCATCAATT‑ACGCACCAAGGCAAAAATATGC‑GG‑CAAGGTGTTTTTTTTCTGTAATAAAAAATTTACTCCCCCACCTCTTTTTTTCCAAATTCTTGTTTTATTCACTGGTTT‑ATTGaaaaaaaacgcta                                                                                                                                              >  2:1627616/1‑230 (MQ=255)

                                                                                              aCCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCg                                                                                                                                               <  1:1862481/238‑1 (MQ=255)

                                                                                              aCCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGTAAGTAATTTCGTCCGTAAAAC‑CATACCGACGTCGTGGCGGGGTA‑CTT‑CTTTAAAAAAATTACCCAGCAAGGGAATAATTTTC‑GA‑GACGGGATTTCTTTTTTTTATTAAAAAATTTACCCCCCCAATTCTTTTATTTAAATTTTTGTTGTTATGCA‑TGGGTTGCTTGAAAATACCggcca                                                                                                                                              >  2:1721168/1‑235 (MQ=255)

                                                                                                cggtccggtTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGGTTAGGAATTTTCTCAGTAGAACACAAACC‑CCCTCGGGGGGGGGTG‑ATT‑CTTCCAAAAATTTAACCAGGAAGGGAAATATTTTC‑GG‑TGAGGTATTTGTTATTTTTAATAAAAAATTTAACCGCCCACTCTTGTTATTTTAATATCTGTTGTTATGCA‑TTGTTTTTTTTaaaaaaaaatacaaa                                                                                                                                            >  2:825028/1‑228 (MQ=255)

                                                                                                       gtTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTTAGTAAATTCGTCAGTACATC‑CCTAACGACTTCGTGGCGCGAGA‑CTT‑CTTCAAAAAAATCACGCAGCAAGGGAATAATTTGC‑GG‑TACGGGATTTTTTATCTGTAATAAAAAATTCACCCGCAGACTTCTGTTATTTAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATAACGTCGAAATCCTGaa                                                                                                                                     >  1:795355/1‑239 (MQ=255)

                                                                                                       gtTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACgg                                                                                                                                                                                                                                         <  1:709630/140‑1 (MQ=255)

                                                                                                       gtTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGGAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGaa                                                                                                                                     >  1:524812/1‑239 (MQ=255)

                                                                                                                        cgcgTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATgg                                                                                                                    <  2:1002275/239‑1 (MQ=255)

                                                                                                                         gcgTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGc                                                                                                                   <  1:739557/239‑1 (MQ=255)

                                                                                                                          cgTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTTGTTAGTAAATC‑CATACCGACTTCGTGGGGCGGTA‑CTT‑TTTTAAACAATTCACGCAGCCAGGTAATAATTTTC‑GA‑TACGGTATTTCTTTTCTGTAATAAAAATTTTACCCGCAGACTTCTGTTATTTAAATTTTTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGTATCAGACGGGAGATATgg                                                                                                                    >  1:884469/1‑237 (MQ=255)

                                                                                                                           gTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGcc                                                                                                                  <  1:1675995/238‑1 (MQ=255)

                                                                                                                             ttGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCTACGTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTAc                                                                                                               <  2:67993/239‑1 (MQ=255)

                                                                                                                             ttGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGTTAATAATTTGC‑GA‑TACGGTATTTCTTAGCTGTAATAAGAATTTCGCCCGCAGACTTCGGTTATTCAAAGTTCTGGTGTTACGCG‑TGGTTTGATTGAAAACTGCGTCGAAATCCTCAGTCAGGCGGGAGGTAGGGCCTag                                                                                                               <  2:1001051/239‑2 (MQ=255)

                                                                                                                             ttGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTAc                                                                                                               <  1:560784/239‑1 (MQ=255)

                                                                                                                             ttGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTAc                                                                                                               <  1:177961/239‑1 (MQ=255)

                                                                                                                              tGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACt                                                                                                              <  2:926270/239‑1 (MQ=255)

                                                                                                                               gAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACtt                                                                                                             <  1:1643416/239‑1 (MQ=255)

                                                                                                                                        cGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                                                                                                      <  2:795355/237‑1 (MQ=255)

                                                                                                                                        cGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                                                                                                      <  1:97315/237‑1 (MQ=255)

                                                                                                                                        cGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                                                                                                      <  1:1435477/237‑1 (MQ=255)

                                                                                                                                        cGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAAt                                                                                                      <  1:73061/237‑1 (MQ=255)

                                                                                                                                        cGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATAc                                                                                                    <  1:856177/239‑1 (MQ=255)

                                                                                                                                         ggTATGAAACGGCGGTTAAGTAATTTCGTCAGTACCTC‑CATGCCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGTTAAGGCCTACTTCGGCAATAc                                                                                                    <  2:1584002/238‑1 (MQ=255)

                                                                                                                                            aGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCAc                                                                                                 <  1:862708/238‑1 (MQ=255)

                                                                                                                                                         tAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGa                                                                                   <  2:1398744/239‑1 (MQ=255)

                                                                                                                                                         tAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGa                                                                                   <  2:689938/239‑1 (MQ=255)

                                                                                                                                                              aaTTTAGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGGAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGCGAGATATGGCCTACTTCGGCAATACCACGAAAATAATGTCGACCttt                                                                              <  2:587532/239‑1 (MQ=255)

                                                                                                                                                                tttCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTAt                                                                            <  1:827033/239‑1 (MQ=255)

                                                                                                                                                                tttCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTAt                                                                            <  1:1156490/239‑1 (MQ=255)

                                                                                                                                                                  tCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATgg                                                                          >  2:1088320/1‑239 (MQ=255)

                                                                                                                                                                  tCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATgg                                                                          >  1:693731/1‑239 (MQ=255)

                                                                                                                                                                     tCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTa                                                                       >  1:1632356/1‑239 (MQ=255)

                                                                                                                                                                     tCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTa                                                                       >  1:1753289/1‑239 (MQ=255)

                                                                                                                                                                       aGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTaa                                                                      >  1:695646/1‑238 (MQ=255)

                                                                                                                                                                           cATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACcac                                                                 <  2:1203803/239‑1 (MQ=255)

                                                                                                                                                                                  tACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCa                                                           >  2:982766/1‑239 (MQ=255)

                                                                                                                                                                                    ccGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCaaa                                                         >  1:1850393/1‑239 (MQ=255)

                                                                                                                                                                                    ccGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCaaa                                                         >  1:1255866/1‑239 (MQ=255)

                                                                                                                                                                                    ccGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCaaa                                                         >  1:29368/1‑239 (MQ=255)

                                                                                                                                                                                       aCTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTAt                                                                            <  2:239785/217‑1 (MQ=255)

                                                                                                                                                                                       aCTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTAt                                                                            >  1:239785/1‑217 (MQ=255)

                                                                                                                                                                                         ttCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAg                                                    <  1:1643518/239‑1 (MQ=255)

                                                                                                                                                                                                  cgATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTa                                           >  2:1582694/1‑239 (MQ=255)

                                                                                                                                                                                                   gATA‑CTT‑CTTCAAACAATTCACGCAGCGAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTaa                                          <  1:1755946/239‑1 (MQ=255)

                                                                                                                                                                                                   gATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTaa                                          <  1:878350/239‑1 (MQ=255)

                                                                                                                                                                                                   gATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTaa                                          <  1:1721168/239‑1 (MQ=255)

                                                                                                                                                                                                         tt‑cttcAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCa                                     <  2:884469/239‑1 (MQ=255)

                                                                                                                                                                                                         tt‑cttcAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCa                                     <  2:1467055/239‑1 (MQ=255)

                                                                                                                                                                                                            cttcAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCaa                                    >  1:943469/1‑238 (MQ=255)

                                                                                                                                                                                                               cAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAacta                                <  1:1493499/239‑1 (MQ=255)

                                                                                                                                                                                                               cAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAacta                                <  1:1003711/239‑1 (MQ=255)

                                                                                                                                                                                                                 aaCAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAACCCACAGGCAAAAATCAGCCTGGCTTAAAGCCAAtctact                              >  2:1459719/1‑239 (MQ=255)

                                                                                                                                                                                                                 aaCAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAact                                 >  1:1001534/1‑236 (MQ=255)

                                                                                                                                                                                                                    aaTTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGcc                           <  1:1184978/239‑1 (MQ=255)

                                                                                                                                                                                                                      ttCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATCTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCTTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTa                         >  2:327552/1‑239 (MQ=255)

                                                                                                                                                                                                                           gcagcaAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCa                    <  1:872823/239‑1 (MQ=255)

                                                                                                                                                                                                                              gcaAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGaa                 >  1:1421531/1‑239 (MQ=255)

                                                                                                                                                                                                                                aaGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGaa                 <  2:1450714/237‑1 (MQ=255)

                                                                                                                                                                                                                                    taataaTTTTC‑GA‑TACTGTATTTCTTATCTGTAATAAAAATTTCTCCCCCAGACTTCTGTTATTCAAATTTCTGGTGTGACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGTCCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAACCTACTGCCTAAATCAGAATCcag            <  1:1047818/238‑1 (MQ=255)

                                                                                                                                                                                                                                    taataaTTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAAGCAGAATCcag            <  2:501334/238‑1 (MQ=255)

                                                                                                                                                                                                                                     aataatTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCcagca          <  1:1878912/239‑1 (MQ=255)

                                                                                                                                                                                                                                     aataatTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCcagca          >  2:1125885/1‑239 (MQ=255)

                                                                                                                                                                                                                                         atTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGAt      >  2:11219/1‑239 (MQ=255)

                                                                                                                                                                                                                                            tGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGTTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGAtttt   >  1:794083/1‑239 (MQ=255)

                                                                                                                                                                                                                                             gC‑GA‑TTCGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTta  <  2:1928453/239‑1 (MQ=255)

                                                                                                                                                                                                                                             gC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGgtgt                                                                                                                                                                             <  1:1560630/69‑1 (MQ=255)

                                                                                                                                                                                                                                             gC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGgtgt                                                                                                                                                                             >  2:1560630/1‑69 (MQ=255)

                                                                                                                                                                                                                                             gC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGTCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTta  <  2:524812/239‑1 (MQ=255)

CTTTACGCAGAATTTCATCACCAGTGTTATGGCCCCACGTATCGTTGATCTCTTTGAATTTATCAACGTCAATAATCAGCACTGACAGCGGTGTACCGGTCCGGTTGGCATGGGCAATTTCGCGTTTGAAGATAGTCGGTAGGAAACGGCGGTTAAGTAATTTCGTCAGTACATC‑CATACCGACTTCGTGGCGCGATA‑CTT‑CTTCAAACAATTCACGCAGCAAGGTAATAATTTGC‑GA‑TACGGTATTTCTTATCTGTAATAAAAATTTCACCCGCAGACTTCTGTTATTCAAATTTCTGGTGTTACGCA‑TGGTTTGATTGAAAATACCGTCGAAATCCTGAATCAGACGGGAGATATGGCCTACTTCGGCAATACCACTAAAATAATGTCGACCTTTATGGTTAAACCACAGGCCAAAATCAGCCTGGCTTAAAGGCAAACTACTGCCTAAATCAGAATCCAGCAGGATTTTA  >  NC_007779/1567789‑1568261

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