BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,977,198	A→G	16.5%	S308G (AGT→GGT)	wzxE →	O‑antigen translocase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,977,198	0	A	G	16.5%	292.8 / 4.6	142	S308G (AGT→GGT)	wzxE	O‑antigen translocase
Reads supporting (aligned to +/- strand): ref base A (35/71); new base G (27/3); total (66/76)
Fisher's exact test for biased strand distribution p-value = 1.63e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

TGACCTTGCCTGTTGCTTACATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATG  >  NC_000913/3976962‑3977432

tGACCTTGCCTGTTGCTTACATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGt                                                                                                                                                                                                                                                >  2:976482/1‑239 (MQ=255)

      tGCCTGTTGCTTACATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGC‑CGGGCGC‑GCGGTTTAcc                                                                                                                                                                                                                                          >  2:466568/1‑239 (MQ=255)

            ttGCTTACATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGTCTgg                                                                                                                                                                                                                                    >  2:495828/1‑239 (MQ=255)

            ttGCTTACATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGGCTgg                                                                                                                                                                                                                                    >  2:258981/1‑239 (MQ=255)

                 tACATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGAATATTTGCTGCCCACTTTGTCGCGGCTAACGGAAAAGCGCGATATA‑CCCCGGGCAGGGGTTAAATCGCTGAAATTCGGCTTCCCGGCGGTGGCCGGCCGG‑GAGttgaacagcatgctgg                                                                                                                                                                                                                                 >  1:149301/1‑224 (MQ=255)

                  aCATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACCTTGTCGCGGCTACCGGACACGCCCGAAATC‑CCCCGGGAACTGGTTAAATCCCTGAAATTCTCCTTACCGGCAGTGG‑GCGCCGC‑GGGGTTTACCGCCTGgcgcctt                                                                                                                                                                                                                              >  1:626899/1‑234 (MQ=255)

                      cgtgtagggTAAACTGGTGGGGGCGGAGTTACGCTGGGATGAGGTGGGGAACTGGCAAGGGGTGAGCAGGATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCTTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTg                                                                                                                                                                                                                          <  1:949330/231‑1 (MQ=255)

                        tgatgCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTgtt                                                                                                                                                                                                                        >  2:1190574/1‑237 (MQ=255)

                                      ctgGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTCACGGAAAAGCGCGATATC‑ACCCGGGAACTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGCGG‑CGGCGGC‑GGCTTTTACCGTCTGGCTGCTGAGTGGTTTTGcgttccgtgc                                                                                                                                                                                                          >  1:215965/1‑230 (MQ=255)

                                      ctgGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGAGGTTAAATAGCTGAAATTCGGCTTACCGGCAGCGG‑GGGAGGC‑GGAATGTACCGTCTGggtgtcgaggggtttttggattcgtgg                                                                                                                                                                                                          >  1:962931/1‑216 (MQ=255)

                                       tgGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTACCGGCAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTGGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTGTTACCGTCTGGCGGCTGCGGGGTTTTGTTAttgtgcgg                                                                                                                                                                                                         >  1:822450/1‑232 (MQ=255)

                                       tgGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGTCTGGCTGCTGCGGGATTTTGCTATCTGgccg                                                                                                                                                                                                         >  2:1061314/1‑237 (MQ=255)

                                                  tataGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCAGCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCg                                                                                                                                                                                                                            <  2:131486/209‑1 (MQ=255)

                                                   ataGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTTTTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAAt                                                                                                                                                                                             <  2:162744/239‑1 (MQ=255)

                                                        tggtgtGAGGTGGGGTTGTGGGAAGGGGTGGGCGGTAATTCCGATGGCTACCTGGAATTTATTCCGGCATCGTGCAGCGTATATTTGCTGCCCACGTTTTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑tt                                                                                                                                                                                       <  1:466568/234‑1 (MQ=255)

                                                          ggATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGCAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGG‑CCGTTTTAACCTCTGGTTGTTGTGTGGTTTTTCTATCTTGGGGGTGGTTTCGAATAAA‑att                                                                                                                                                                                      >  1:319336/1‑238 (MQ=255)

                                                          ggATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGGGGTTAAATCGCTGCATTTCGTCCTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGGCTTTCTGTTTCATGGTTTTGTATTCTGGACGCAGTTTTCGataaaaatt                                                                                                                                                                                       >  1:524406/1‑230 (MQ=255)

                                                            aTGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGG‑GAGTTTTACCGTCTGGCTGTTTCGTGATTTTGCTATCTGGCTGCTGTTGTGGGATAAA‑TTTAcc                                                                                                                                                                                   >  2:773749/1‑239 (MQ=255)

                                                            aTGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GCGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑tttta                                                                                                                                                                                    >  2:78436/1‑236 (MQ=255)

                                                                  tGGGGATCTGGCAAGGGGTGAGTAGTATTTCCGATTCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTa                                                                                                                                                                               <  2:206643/237‑1 (MQ=255)

                                                                    gggATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGCAATTCGTCTTACCGGCAGTAG‑CGGCGGC‑GAGGTTTACCGTCTGGCCGCTGCGGGGTTTTGATATCTGGCTGCTGGTTGCGGATAAA‑TTTACCCCTA‑Tgc‑g                                                                                                                                                                         >  2:518582/1‑239 (MQ=255)

                                                                    gggATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGTCTGGCTGCTGCGGGATTTTTCTATCTGGCGGCTTTTTTGGAATAAA‑TTTACCGCTA‑Tgc‑g                                                                                                                                                                         >  2:538652/1‑239 (MQ=255)

                                                                    gggATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑Tgc‑g                                                                                                                                                                         <  2:104556/239‑1 (MQ=255)

                                                                    gggATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATA‑ACCCGGGAAGTGGTTAAATCGATGAAATTCGTCTTACCGGCAGTGG‑CGGCGGG‑GGGGTTTACCGTCCGTGTGCGGCGGGGCTTGTCGATATGGCGGGGGTTTTCGGATAAA‑TTAAACGCTA‑Tgtg                                                                                                                                                                          >  1:680871/1‑237 (MQ=255)

                                                                      gATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGACTTTGCCATATGGCTGCTGTTTTCGAATACA‑TTTACCCCTA‑TGt‑gcg                                                                                                                                                                       >  2:867103/1‑239 (MQ=255)

                                                                       ttttGGCAAGGGGTGGGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGc‑gcg                                                                                                                                                                       <  1:449632/235‑1 (MQ=255)

                                                                       aTCTGGCAAGGGGTGAGCAGTATTTCCAATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGc‑gcg                                                                                                                                                                       <  2:59387/238‑1 (MQ=255)

                                                                          tGGCGAGGGGTGAGGGGGATGTCCGAGGCCGACCTGAAAGGTATGACGTCATCGGTAAGCGTATATTAGGTGCCCAGTGTGGAGGGGTTAACGGAAAAGCGAGATATC‑ACCCGGGAAGTGGTTAAAGCGCTGAAATTGGGCTTACCGGCATGGG‑CGGGGGC‑GAGTGTGACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGGTGGCGAATAAA‑TTTACCGCTA‑TGC‑GCGAgct                                                                                                                                                                   <  1:736661/239‑1 (MQ=255)

                                                                                  gttggtcAGTTTTTCCGTTGCTGTCCTGGAATTTTTTACGGCAGCGTTCGGCGTTTATTTGCTGCCCTCGTTGACGCGGCTAACGGAAAAGCGCGTGCTC‑CCCCGGGAAGTGGTTAAATCGCTGAGTTTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGcc                                                                                                                                                            <  1:167497/232‑1 (MQ=255)

                                                                                   gggAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGGTTTACCGTCTGGCTGCTGCGGGATTTTGCTATCTGGCCTGTGTTTTCGGATAAA‑TTTTCCGCCA‑TGG‑GCGTTCTCTTTTCCT‑g                                                                                                                                                         >  2:1164734/3‑239 (MQ=255)

                                                                                     tagcatgttttCCGATGCCTTCCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑gg                                                                                                                                                        <  2:626899/230‑1 (MQ=255)

                                                                                     taGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCCTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑gg                                                                                                                                                        <  2:696858/237‑1 (MQ=255)

                                                                                      aGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGGTTTTGCTTTCTTGCTGCTGTTGTCGATTAAA‑TTTTCCGCTATTGC‑G‑GGTTTCTTTTCCT‑GGcc                                                                                                                                                      >  2:428070/1‑238 (MQ=255)

                                                                                      aGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGcc                                                                                                                                                      >  2:1091368/1‑238 (MQ=255)

                                                                                           tATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGTCTGGCTGGTGCGTGGTTTTGGTTTCTGGTGGCTGTTTTCGGATAAA‑TTTACCCCCA‑TGC‑GGGATCCTTTTGCCT‑TGGAGtatg                                                                                                                                                 >  1:155081/1‑236 (MQ=255)

                                                                                            aTTTCCGATGCCTACCTGAAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAg                                                                                                                                                 <  2:1320670/238‑1 (MQ=255)

                                                                                             tttCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGCGGGGTTTTCCCTCCGGGTTGCCGGGTG‑TTTTGTTTTCTGGGCTCGGGTGTCAGATAAA‑TTTATCGCTA‑TGC‑GGGGTCTCTTTGCCg‑gggcggttggg                                                                                                                                               >  1:34570/1‑230 (MQ=255)

                                                                                                  ggtgccctccTGCATTTTTTTTCGGCATCGTTCGGCGTTTTGTTGCTGCCGACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑CCCCGGGAAGTGGTTAAATCGCTGAAATTCGTATTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGGTGTCGAATAAA‑TTTGCCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGa                                                                                                                                         <  1:1061314/231‑1 (MQ=255)

                                                                                                  gaggCCTCCCTGGAATTTATTCCGGCCTGGTTCAGTGTATTTTTTGTGGCCACGTGGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGa                                                                                                                                         <  1:105140/236‑1 (MQ=255)

                                                                                                  gATGCCTACCTTCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGa                                                                                                                                         <  2:858309/239‑1 (MQ=255)

                                                                                                  gATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGa                                                                                                                                         <  2:950814/239‑1 (MQ=255)

                                                                                                  gATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGa                                                                                                                                         <  2:1181998/239‑1 (MQ=255)

                                                                                                  gATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGa                                                                                                                                         <  2:597059/239‑1 (MQ=255)

                                                                                                  gATGCCCACCTGCAATTTATTGCGGCACCGTTCAGGGTATATTTGCGGCCCACGTTGTCGCTGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGa                                                                                                                                         <  1:1164734/239‑1 (MQ=255)

                                                                                                     gCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGGAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGTGTTTGCTTTCTGGGTGCGGTTGTCGAATAAA‑TTTTCCCCTC‑TGCGGCGAT‑TCTTTGCCT‑GGcgggtggggggggcggg                                                                                                                                       >  1:615481/1‑223 (MQ=255)

                                                                                                     gCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGGGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑TGCAGTTTTG‑GGGGGAtgt                                                                                                                                      >  2:461663/1‑239 (MQ=255)

                                                                                                        tACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTAACGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCGGTTATT‑GGGTGATGG‑Gtt                                                                                                                                  >  2:138337/1‑239 (MQ=255)

                                                                                                         aCCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTa                                                                                                                                 <  2:1243752/239‑1 (MQ=255)

                                                                                                              gaaTTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGt                                                                                                                            <  2:1233572/238‑1 (MQ=255)

                                                                                                                 ttttttccggcattttttcacGTATTTTTGCTGCCCCCGGTGTCGCGGGTAACGGAAAAGGGCGATATC‑ACCCGGGAAGGGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTGACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTggg                                                                                                                         <  1:1190574/219‑1 (MQ=255)

                                                                                                                 tttATTACGGCATCGTTCAGCGTTTATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGGGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTggg                                                                                                                         <  2:189807/239‑1 (MQ=255)

                                                                                                                 tttATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTggg                                                                                                                         <  2:319336/239‑1 (MQ=255)

                                                                                                                      tttataagagacaggCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCtt                                                                                                                    <  2:416632/225‑1 (MQ=255)

                                                                                                                      taCGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTTTCGAATAAA‑TTTACCGCTA‑TGC‑TCGATCCCTTTGCCT‑GGCAGTTAAT‑GGGTGAAGT‑GTTAAAATGGGGGGCtt                                                                                                                    >  2:755913/1‑239 (MQ=255)

                                                                                                                      taCGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCtt                                                                                                                    <  2:266442/239‑1 (MQ=255)

                                                                                                                           cATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGt                                                                                                                <  2:1094539/238‑1 (MQ=255)

                                                                                                                               gTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTg                                                                                                           <  2:106027/239‑1 (MQ=255)

                                                                                                                                ttCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTgg                                                                                                          <  2:22714/239‑1 (MQ=255)

                                                                                                                                     cGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGGGGTGT‑GTTAAAAGTGGGGGCTTATTTCTTTGGTTAtt                                                                                                     >  2:5770/1‑238 (MQ=255)

                                                                                                                                         tcgtgcgtgcCCACGTTGGCGCGGCCAACGGAAAAGCGCGAAATC‑ACCCGGGAAGTGGTTAAATCGCGGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGt                                                                                                 <  2:302830/232‑1 (MQ=255)

                                                                                                                                          atTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGG‑GAGGTTTACCCTCTTGCTTCTGCGGGATTTTTCTATCTGGGTGCTTTTTTCCGATAAA‑TTTTCCCCTT‑TGT‑GGGGACTCTTTTTCT‑GGCGGGTAGG‑GGTGGGTGTGGTTAAAAG‑GGGGGCTTGTGGTTTTGGTTTTTTgggg                                                                                                >  1:996745/1‑237 (MQ=255)

                                                                                                                                          atTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGGGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTTGT‑TGGTGATGG‑GTTAAAAGTGGGCGCTTATGTCTTTTGTTTTCTGGTg                                                                                                >  2:1303452/1‑239 (MQ=255)

                                                                                                                                               ctgcCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGGCTGGCTTCTTCGTGATTTTTCTCTTTGGCTGCCGTTTTCGAATAAATTTTACCCTTT‑TGC‑‑CGTTCTCTTTGCCT‑GGCAGTTTGT‑GGGGGGTGT‑‑TTAAAAATGGGGGGCTATTTCTTTTGttttttgtgagttgcc                                                                                          >  1:1294572/1‑225 (MQ=255)

                                                                                                                                               ctgcCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGCCCGGGC‑CGGTGGG‑CGGTTTTTTCCTCTAGCTGGTGCGTGTTTTTTTTTTCTGGATTGAGTTTTCGATTAAA‑ATTACCCCTC‑TTC‑CTCGTCTCTTTGCCG‑TGCGGTTAAA‑GGGGGTTAT‑GTTAAAAAGGTGCGCTATTTTCTTTGGATaagtggggaaagc                                                                                           >  1:302830/1‑227 (MQ=255)

                                                                                                                                                tgcCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGc                                                                                           <  2:880038/238‑1 (MQ=255)

                                                                                                                                                   ccACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGTAGTGGTTAAATCGCTGAAATTCGTCTTACCTGCAGTGG‑CGGCCGC‑GGGTTTTACCGGCCTGCTGCCGCGTGATTTTGCTATATGGCTGCCTTTTTCGAAGTAA‑TTTACCCCTATTGC‑CCGTTC‑CTTTTCCTGGGCTGTT‑GG‑GGGTGGTTT‑GTTAAAATGGGGCGCTTTTGTCTTTGGTTTTCttggggttgcccaa                                                                                       >  1:928029/1‑226 (MQ=255)

                                                                                                                                                    cACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGACATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGTTTTACCGGCTGGCTGCTGCGTGCTTTTGCGTGCTTGCTTCCGGTGTCGAATATA‑TTTACTGCTT‑TGC‑GCGTTCTCTTTGCTG‑GGCTGTTAGTGGGGTGGCGG‑TTTAAAAATGGGCGTGTGT‑TCTTTGGGTTTCTgggtgtagctaaga                                                                                      >  1:983898/1‑227 (MQ=255)

                                                                                                                                                    cACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGGTTTACCGGCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGGTTCCTTTTCCT‑GGCGGTTAGT‑GGGTGTTGT‑GTTAAAAGTGGGCGCTTATGTCTTTGTGTATCTgggtgtctcttata                                                                                      >  1:205099/1‑227 (MQ=255)

                                                                                                                                                             gcgGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCCGGCTGCCGCGTGATTTTTCTATCTCGCTGCCGTTGTCGGATAAA‑TTTACCGCTA‑GTC‑GCGGTCTCCTTGTCC‑GGCTGTTTGT‑GGGGGCTGG‑GTTTAAAGAGGGGCCTTCTTTCTTTGGTTATTTGGGTATCGCCAAAAAGTCACAGc                                                                             >  1:2652/1‑239 (MQ=255)

                                                                                                                                                             cgggttAACGGAAAAGCGCGATTTC‑ACCCGGGAAGTGGTTTAATCGGTGAAATTCGTCTTACCGGCAGTGG‑GGGCGGC‑GAGTTTTACACTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGTTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTg                                                                              <  1:40124/233‑1 (MQ=255)

                                                                                                                                                              cgGTTAACGGAAAAGCGCGATTTC‑GCCCGGGAGGTGGTTAAATCGTGGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCg                                                                            <  1:461663/239‑1 (MQ=255)

                                                                                                                                                               ggCTAACGGAAAGGCGCGATATC‑CCCCGGGAAGTTGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCg                                                                            <  2:615481/238‑1 (MQ=255)

                                                                                                                                                                          aaacgcgGTATCTACCCGGGAATTTGTTAAATCGCTTAAATTCCGCTTACCGGCATTGG‑GGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATAtt                                                                  <  1:1199620/235‑1 (MQ=255)

                                                                                                                                                                             cgcgATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGGATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTttatctgtct                                                             >  2:1152588/1‑233 (MQ=255)

                                                                                                                                                                              gcggtttC‑ACCCGGGAGGTGGTTAAATCGCTGAAATTCGTCTTACCGGGAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGc                                                             <  1:91732/232‑1 (MQ=255)

                                                                                                                                                                              gcgATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGc                                                             <  2:545995/238‑1 (MQ=255)

                                                                                                                                                                                   tttcccccGGGGAGTGGTTAAATTCGTGAAATTCGTTTTGCCGGCGGTGG‑CGGGGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAg                                                        <  1:959278/234‑1 (MQ=255)

                                                                                                                                                                                   atC‑CCCTGGGGAGTGGTTAAATCGCTGGAATTCTTCTGACCTGCAGTGG‑CGGCGGC‑GAGTTTGACCGTTTGGCTGCTGCGTGATTTGGCTATCTTGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAg                                                        <  1:218053/238‑1 (MQ=255)

                                                                                                                                                                                       aCCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATTTTCTGGCGGAAGTCa                                                     >  2:777440/1‑238 (MQ=255)

                                                                                                                                                                                        cccGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TTC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGGCTTTGGTTATCTGGTGTTTGCCAAAGCGTCCATGCGGTTTTATATTTTGGCGGAAGTCAg                                                    >  2:759967/1‑238 (MQ=255)

                                                                                                                                                                                        cccGGGAAGGTGTTAAATCCCTGAAATGCGTCGTACCGGCAGTGG‑CGGCTGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTTTCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGc                                                   <  1:138337/239‑1 (MQ=255)

                                                                                                                                                                                             gAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGGTTTACCGTCTGGCTGGTGCTGGATCTTGCTCTCTGGCTGCTTTTGACGAATACA‑TTTAGCCCGA‑TGT‑GCGATCTCCTTGCCT‑TGCAGTTTGT‑GGGAGATGT‑GTTAAAAGTGGGGCGTTTTGCCTTTGATGATATGGTGATAGCGAAAGGGCTACTGGGGGTTTTTATCTTGGTGGAAgcagaaccagt                                              >  1:252683/1‑229 (MQ=255)

                                                                                                                                                                                              aaGTGGTTAACTCGTTGCATTTCGGCTTGCCTGCGGTGG‑CGCGGGA‑GAGTTTTGCCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGGTGTCGACTAAA‑TTTACCGCTA‑GGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATTG‑GTTAAAAGTGGGCGCTCTTGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGtt                                              <  1:518582/238‑1 (MQ=255)

                                                                                                                                                                                              aaGTGGTTAAATCTCGGAAATTCGTTGTGCCGGAAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGtt                                              <  2:1184961/238‑1 (MQ=255)

                                                                                                                                                                                              aaGTGGTTAAATCGGTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGtt                                              <  2:215965/238‑1 (MQ=255)

                                                                                                                                                                                                   gTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GGGGTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACAGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTTGT‑GGGGGTTGT‑TTTAAAAGGGGGCGCCTTTGTCTTTGGTTATCTTGGGATCGGCGAAGAGCAACTTGGGTTTATTATTGTGGGGAGGGTCAGcccagcaagttt                                        >  1:471355/1‑229 (MQ=255)

                                                                                                                                                                                                       aaTCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCATGTCACTTTATTg                                    >  2:917947/1‑239 (MQ=255)

                                                                                                                                                                                                         tCGCTGAAATTCGTCTTGCCGGCAGGGG‑TGGCGGC‑GAGTTTTACCGTTTGGCTTGTGCGTGATTTGGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGAt                                  <  1:389429/239‑1 (MQ=255)

                                                                                                                                                                                                         tCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGAt                                  <  2:680871/239‑1 (MQ=255)

                                                                                                                                                                                                         tCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑AGGTTTTACCGTTTGGCTGCTGCGTGATTTTGCTATCTTGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATTTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGAt                                  <  2:518913/239‑1 (MQ=255)

                                                                                                                                                                                                             tGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCtt                                                                                                                    <  1:8301/153‑1 (MQ=255)

                                                                                                                                                                                                             tGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCtt                                                                                                                    >  2:8301/1‑153 (MQ=255)

                                                                                                                                                                                                                 aTTCGTTTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGc                                                   <  2:150015/214‑1 (MQ=255)

                                                                                                                                                                                                                    gtttttACCGGTAGTGG‑CGGCGGC‑GTGTTTTTCTGTTTGGCTTCTGCGGGATTTTGCTATCGTGCTGCTGGTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTTATGT‑GTTAAAAGTGGGCGCTTATGTCTTGGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGccc                       <  1:1179135/235‑1 (MQ=255)

                                                                                                                                                                                                                     ggcTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGGTGTTTGCGAATAAA‑TTTACCGCTA‑TTC‑GCGGTCTCTTTGCCT‑GGCAGTTAGT‑TGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCa                      <  2:89869/237‑1 (MQ=255)

                                                                                                                                                                                                                     gTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCa                      <  2:34570/239‑1 (MQ=255)

                                                                                                                                                                                                                      tCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTTTCTGGCTGCTGTTGTCGAATAAA‑TTTAACGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCCGTTTGT‑GGGTGTTGT‑GTTAAAAGTGGGGGCTTTTGTTTTTGGTTATCTGGGGGTCTCCAAAGGGTCACGGCGGGTTTATATTTCGGCGGAAGTTAGCCCGTTCCCCTTTTTTGTGGGATTTGcccca                     >  1:1019507/1‑237 (MQ=255)

                                                                                                                                                                                                                           ccGGGAGGGG‑GGGCGGC‑GGTTTTTTCCGTTTGGCTGCTGCGTGATTGTGGTATCTGGCTGGTGTTGTCGAATAAA‑GTTACCGCTA‑TGC‑GCGATCGTTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGGGCTTATGTCTTTGGGTATCTGGTGATCGCAAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCt                <  1:448302/239‑1 (MQ=255)

                                                                                                                                                                                                                            ccgtctgttgtgggtgc‑gTGTTTTTCCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTTGTGTCGAATAAA‑TTTACCGGTA‑TGC‑GCGATCTGTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCt                <  1:657455/224‑1 (MQ=255)

                                                                                                                                                                                                                                           gc‑gAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATg  <  2:962931/238‑1 (MQ=255)

TGACCTTGCCTGTTGCTTACATCATGATGCGTAAACTGCTGGCGGCGCAGTATAGCTGGGATGAGGTGGGGATCTGGCAAGGGGTGAGCAGTATTTCCGATGCCTACCTGCAATTTATTACGGCATCGTTCAGCGTATATTTGCTGCCCACGTTGTCGCGGCTAACGGAAAAGCGCGATATC‑ACCCGGGAAGTGGTTAAATCGCTGAAATTCGTCTTACCGGCAGTGG‑CGGCGGC‑GAGTTTTACCGTCTGGCTGCTGCGTGATTTTGCTATCTGGCTGCTGTTGTCGAATAAA‑TTTACCGCTA‑TGC‑GCGATCTCTTTGCCT‑GGCAGTTAGT‑GGGTGATGT‑GTTAAAAGTGGGCGCTTATGTCTTTGGTTATCTGGTGATCGCCAAAGCGTCACTGCGGTTTTATATTCTGGCGGAAGTCAGCCAGTTCACTTTATTGATGGTATTTGCCCACTGGCTAATCCCTGCGCATG  >  NC_000913/3976962‑3977432

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