BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,812,243	T→C	4.8%	T273T (ACT→ACC)	emrA →	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,812,243	0	T	C	4.8%	328.6 / 3.0	107	T273T (ACT→ACC)	emrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base T (49/52); new base C (4/2); total (53/54)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 107 total aligned reads displayed.

CCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCAATACCACTA  >  NC_000913/2812013‑2812469

gtcgggtgtttttaagaggacgccGAGGTGCGTACCGCCTGGTTGGCTGTGGAGCGTACTCG‑TTTTATCAGCCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑AGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATcacca                                                                                                                                                                                                                                 <  2:916090/216‑1 (MQ=255)

 cagccGTGCACCAGGCTGCCACCGAAGTGCGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TAGTAGCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATcac                                                                                                                                                                                                                                   <  1:273760/235‑1 (MQ=255)

 aagccGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATcacca                                                                                                                                                                                                                                 <  1:468160/236‑1 (MQ=255)

        gCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCCCCAATATTTGGGGGGATGCCAACTTTAAAGAGACGCAGATTTCCAATAACCGTATCGGTCAGCCGGGCACCACCACCCCGGAAAttt                                                                                                                                                                                                                        >  2:646020/1‑239 (MQ=255)

              ggcttgcACCGAAGTACGTAACGCCTGGCTGGCTCTGGACCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCCCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCg                                                                                                                                                                                                                  <  1:593626/233‑1 (MQ=255)

                cTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACCATCACCACGGATATTTACGGCga                                                                                                                                                                                                                 <  1:646020/238‑1 (MQ=255)

                 cccccccgaatgccGTACCGCTTGGGTGGCGGTGG‑GGGTTCCCGTTATTACCCGTCCGATGACCGGTTATGTCACCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCgatg                                                                                                                                                                                                               <  2:319111/226‑1 (MQ=255)

                  gCCACCGAAGTGCGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCgatg                                                                                                                                                                                                               <  1:158087/238‑1 (MQ=255)

                    ccccGAAGTACGTAACGCCTGGCTGGCGTTGGAGCGTACTCG‑TATTATCAGCCCGATGCCCGGTTATGTCCCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCgatgat                                                                                                                                                                                                             <  1:180340/236‑1 (MQ=255)

                    cACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAAATTGCCAATATGGGTATCGGGCAGCCGGTCCCTAACACCCCGGAAATTTACGGGgaggatt                                                                                                                                                                                                            >  2:124645/1‑238 (MQ=255)

                    cACCGAAGTACGTAACGCCGGGCTGGCGCTGGAGCATACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCgatgat                                                                                                                                                                                                             <  1:662629/238‑1 (MQ=255)

                           gTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAAATTCCCAATATGCGTAACGGTCAGCCGGTCACTACCACCACGGATATTTACCGCCGTGAGGTGgaaaa                                                                                                                                                                                                     >  2:869331/1‑237 (MQ=255)

                             aCGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGAATCGGTCAGCCGGTCACTATCACCCCGGATATTTACGGCAATGATGTTGAATac                                                                                                                                                                                                    >  1:809393/1‑238 (MQ=255)

                                 aaCGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGGTGCCCACTTTAAAGAGACGCCGATTGCCCATATGCGTATCGGTCAGCCGGTCCCTATCACCACGGAAATTTCAGGCGAGGATGTGAAAAAAACCgg                                                                                                                                                                                               >  2:22870/1‑239 (MQ=255)

                                 aaCGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAGAAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCCCTATCACCACGGGTATTTAcctggtctttttacacatccccg                                                                                                                                                                                                >  1:916090/1‑216 (MQ=255)

                                           tggcGCTGGAGCGTACTCG‑TATTACCAGTCCGATGACCGGTGATGGCGCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGtggt                                                                                                                                                                                      <  1:596194/238‑1 (MQ=255)

                                               gcTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGtggttggt                                                                                                                                                                                  <  1:124645/238‑1 (MQ=255)

                                               gcTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGtggttggt                                                                                                                                                                                  <  1:206335/238‑1 (MQ=255)

                                               gcTGGAGCGTACGCG‑TATTAACAGTCCGATGCCCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGtggttggt                                                                                                                                                                                  <  1:698587/238‑1 (MQ=255)

                                                gtGGAGCGTTCTCG‑GATTATCCGTCCGATGAACGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGG‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTc                                                                                                                                                                                 <  1:535915/237‑1 (MQ=255)

                                                gctgggggcttcccgtttttgcggTCCGAATACCGGTTTGGCCCCCCGGCGCCGGGGAAAGGCTGGGGCGCAAATTAGGCCAACGACCGCCGTTG‑TGGCGGTCGTTCCAGGCACCAATATGTGGGGTGATGCCAACTTTAAAGAGACCCAGATTTCCAATATGGGGATTGGGCAACCGGTCACTAACACCCCGGATATTTACGGCGATGATGTGAAATACACCGGTTAAGtggttggt                                                                                                                                                                                  <  2:708184/215‑1 (MQ=255)

                                                cTGGGGCGTACTGG‑TATTAGCAGTCCGATGACCGGTTATGGCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTc                                                                                                                                                                                 <  1:22870/238‑1 (MQ=255)

                                                   gAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCAACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTAACTATCACCACGGCTATTTACGGCGATGATTTGTAAAACACCGGTAAAGTGGTTGGCCTg                                                                                                                                                                               >  1:596451/1‑237 (MQ=255)

                                                      ggTACTTG‑TTTTATCAGTCCGATGACCGGTTATTTCCCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGata                                                                                                                                                                           <  2:28741/237‑1 (MQ=255)

                                                      cggaCTCG‑TATTACCAGTCCGATGACCGTTTATGTCCACCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGata                                                                                                                                                                           <  1:574995/235‑1 (MQ=255)

                                                            tcggttgcgcacgcctgtggcgggTTAGGGCTCCCGCCGCGGGGTACAGCCTTGGGAGAAAATTAGCCCAACGA‑CGCCGCTGATGGCGGGCGATCCAGCCACCAATATGGGGGTGGATGCCAACGTGAAAGAGACGCAGAGTGCCAATATGCGGATCGGTCAGCCGGTCACTATCACCACGGATAGTTACGGCGATGATGTGAAATACTCCGGTAAAGTGGTTGGTCTGGATATGGGca                                                                                                                                                                     <  2:596451/217‑1 (MQ=255)

                                                             g‑TATTATCAGTCCGATGACCGGTTATGTCCCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGtggt                                                                                                                                                                                      <  1:153144/220‑1 (MQ=255)

                                                             g‑TATTATCAAGCCGAGCACCGGTTATGCCGCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCTCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATa                                                                                                                                                                                                     <  1:237326/205‑1 (MQ=255)

                                                                    acAGTCCGATGACCGGTTAAGGCCCCCGCCGCGGGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑gc                                                                                                                                                           <  1:80433/238‑1 (MQ=255)

                                                                      aGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGCTGCCACATTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGCAATTTACGGCGATGATGTGAAAAACACCGGGAAAAGGGTTGGTCTGGATATGGGccaagggaaggc                                                                                                                                                           >  2:554071/1‑228 (MQ=255)

                                                                        tCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGCCTGGAAATGGGctgtcccttttaa                                                                                                                                                          >  1:32153/1‑226 (MQ=255)

                                                                        tCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGAATTTTCCGGCGGTGATGTGAAATCCACCGGAAAAGGGGTTGTGCTGGAAATGGGGCAA‑GGTA‑ACGCGt                                                                                                                                                       >  2:80184/1‑239 (MQ=255)

                                                                        tCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATAACCACGGATATTTACGGCGATGATATGAAAAACACCGGCAAAAAGGGTGGTCTGGAGATGGGCAAA‑gagaaggagc                                                                                                                                                        >  2:399041/1‑230 (MQ=255)

                                                                        tCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGGGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGGATCGGTCAGCCGGTCACTATCACCACGGCTATTTACGGCGATGATGTGAAATACACCGGTAAAAGGGGTGGTTTGGatactggttcttaaaaacaac                                                                                                                                                         >  2:153144/1‑221 (MQ=255)

                                                                        tCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTAGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGCCGCCGATTGCCAATATGCGTATCGGGCAGCCGGTCACTATCAACACGGAATTTTCCGAGGAAGATGGTAAATAAAACGTAACAGTGGTGGGTTTGGAGAAGGGcacaagcaaagata                                                                                                                                                         >  2:569115/1‑229 (MQ=255)

                                                                            aTGACCGGTTATGTCTCCCGCCGCGCGGGACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGGGGATGCCAACTTTAAAGAGACGCAGATTGCCAAAATCGGTATCGGTCACCCGGCTACTATCGCCACGGATTATTCCGGGCGTGATTTGGAATAAACCGGGAAAGAGTGTGGTCTGGATATGGGCACA‑AGGA‑Gcgcgcttc                                                                                                                                                    >  2:499959/1‑234 (MQ=255)

                                                                             tGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGAGGATGTGAAATCAACCGGTAAAGTGGTTGGCTTGGATATGGGCCCA‑GGTA‑ggcggttctt                                                                                                                                                   >  1:610395/1‑229 (MQ=255)

                                                                                           tcCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑CCGCGTTTT‑CACTGCTTCCA‑gcgcc                                                                                                                                  >  1:30778/1‑238 (MQ=255)

                                                                                              cgccgcGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGtggt                                                                                                                                                                                      >  2:323981/1‑188 (MQ=255)

                                                                                              cgccgcGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGtggt                                                                                                                                                                                      <  1:323981/188‑1 (MQ=255)

                                                                                              cgccgcGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCaaaa                                                                                                                               >  1:52394/1‑239 (MQ=255)

                                                                                              cgccgcGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGCGGTTGGTCTGGATATggg                                                                                                                                                                       >  1:627424/1‑203 (MQ=255)

                                                                                              cgccgcGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTAAGGCGATGATGTGAAATACACCGGTCAATGGGTTGGTCTGGACATGGGCACC‑GGGA‑Gcgcggttccccccgtccccacgcaaaa                                                                                                                                 >  2:600673/1‑216 (MQ=255)

                                                                                                          cAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTCTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGGTGATGTGAAATACACCGGTCAACGGGTTGGTCTGGGTATGGGCACA‑GGTA‑AGGCGGTCC‑CACTGCTTCCA‑GCGGAAAAAGCGACCGGGaa                                                                                                                   >  2:22370/1‑239 (MQ=255)

                                                                                                                ggggCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCAAA‑GGTA‑GCGCGTTTT‑CACTGCTTCCA‑CCGAAAAAtagaaccgggaaagggac                                                                                                             >  2:847251/1‑222 (MQ=255)

                                                                                                                  ggCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCCCGGCTATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGGTGGTCTGGGAATGGGCACAGGGCA‑ACG‑GTTCTACAC‑ACTTCCA‑GGGCAAAATGCGCCCGGGAAcgtgataa                                                                                                           >  2:407557/1‑232 (MQ=255)

                                                                                                                  ggCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCCCGGATATTTACGGCGATGATGGGAAATACACCGGTAAAGTGGTTGGACTGGATATGGGCACA‑GGAA‑GCCCGTTCT‑CCCTGCTTCCA‑GGCCAAAAAGGCACCGGAAAcggaataa                                                                                                           >  2:154416/1‑232 (MQ=255)

                                                                                                                  ggCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATAAACCGGTAAAGGGGGTGGTCTGGAGATGGGCACA‑GGGA‑GCGCGTTCT‑CACTGATTCCA‑GCGCAAAATgacgacgggaaagggaaaa                                                                                                           >  2:54265/1‑221 (MQ=255)

                                                                                                                  ggCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGGGGAATACACCGGTAAAGGGGTTTGTGTGGTAATTGGCACA‑GGCAGGCGCTTTC‑‑CCCGGCTTCCACGGGAAAAATG‑GACCGGGAAATGGacaa                                                                                                           >  2:224532/1‑236 (MQ=255)

                                                                                                                  ggCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACCATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTTT‑CCCTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCa                                                                                                           >  1:681645/1‑239 (MQ=255)

                                                                                                                       aaaTTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGGGACCGGTAACTGGATCAAAGTc                                                                                                      >  1:82206/1‑239 (MQ=255)

                                                                                                                              cccAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACg                                                                                                                                                                                                                               >  1:325457/1‑115 (MQ=255)

                                                                                                                              cccAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACg                                                                                                                                                                                                                               <  2:325457/115‑1 (MQ=255)

                                                                                                                                cAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGt                                                                                             <  1:332072/239‑1 (MQ=255)

                                                                                                                                         ccgcTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATTTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CCCTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCttt                                                                                     >  1:293124/1‑237 (MQ=255)

                                                                                                                                                    ggTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGaa                                                                           <  1:473989/238‑1 (MQ=255)

                                                                                                                                                      tCGTTCCAGCCACCAATATGGGGGTGGATGCCAACTTTAAAGGGCCGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACCATCCCCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACt                                                                         <  2:681645/238‑1 (MQ=255)

                                                                                                                                                                caccaATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGa                                                                 <  1:154416/236‑1 (MQ=255)

                                                                                                                                                                    aaTATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCCTCGTCTGCCGTTGCGTATCGAACTGGCCCA‑AAACaa                                                            >  2:506361/1‑235 (MQ=255)

                                                                                                                                                                     atatGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTg                                                         <  1:80184/239‑1 (MQ=255)

                                                                                                                                                                          tgggtggATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCa                                                    <  1:224532/239‑1 (MQ=255)

                                                                                                                                                                             gtggATGCCAACTTTAAAGGGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAat                                                  <  1:207034/238‑1 (MQ=255)

                                                                                                                                                                             gtggATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCGGCGTCTGCCTGCCCGTATCGAACTGGa                                                                      >  1:407075/1‑218 (MQ=255)

                                                                                                                                                                             gtggATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATAGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGGAACTGGATCAAGGTCGTTACGCGTTTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGGGCAat                                                  >  1:66037/1‑238 (MQ=255)

                                                                                                                                                                               ggttGCCAACTTTAAAGGGGCGCAGATTGCCAATTTGGGGTACGGTTAGCCGGTCACTACCCCCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATc                                               <  2:52394/236‑1 (MQ=255)

                                                                                                                                                                                gATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATggg                                                                                                                                                                       >  1:672687/1‑122 (MQ=255)

                                                                                                                                                                                gATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATggg                                                                                                                                                                       <  2:672687/122‑1 (MQ=255)

                                                                                                                                                                                gATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATACGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCCAAATTCGACCGGTAACTGGATCAAAGTCGTTCCGCTTCCTCCTGTCCGTATCCACCGGGACCCACAACCCCTGGGCCAATATc                                               >  2:875402/1‑238 (MQ=255)

                                                                                                                                                                                 aTGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGACGTTCAGCGTCTGCCTGTGCGTATCGAAACGGACAAGAAAAAGCTGGAGCAAAAACCg                                             >  2:335663/1‑239 (MQ=255)

                                                                                                                                                                                  tGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAATTGGACCAGAAACACCTGGGGCAATATcccc                                            >  1:673877/1‑237 (MQ=255)

                                                                                                                                                                                   gCCCACTTTAAAGAGTCCCTGATTGCCAATTGTCGTATAGGGCAGCCGGTACCTATCACCACGGATATTTACGGCGGTGATGTGAAATACACCGGTCAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCTGTTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCt                                           <  2:231824/239‑1 (MQ=255)

                                                                                                                                                                                   gCCAACTTTAAAGAGACGAAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATTCACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGc                                                     <  1:87648/229‑1 (MQ=255)

                                                                                                                                                                                        cTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGGCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTa                                      <  1:847251/239‑1 (MQ=255)

                                                                                                                                                                                           tAAAGAGACGCAGATTGCCAATATGCGTCTCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGCTGTGAAATCCACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGGAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGGACGGAACTGGACCCAAAACAGCAGGAGCacaatacgcgtcggttcg                                   >  2:193875/1‑222 (MQ=255)

                                                                                                                                                                                            aaaGAGGCGCAGATTGCCAATTTGCGTATCGGTCAGCCGGTCACTTTCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCgg                                  <  1:577237/239‑1 (MQ=255)

                                                                                                                                                                                             aaGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGACAAGTGGTTGGTCTGGATATGGACACA‑GGTA‑GACCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCATGGTCTGCCTTTGCGTAGCGAACCGAACCAAAAAAAGCCGGAGCAAAAACAGCTgaagttcata                                 >  2:87648/1‑230 (MQ=255)

                                                                                                                                                                                                  aCGCAGATTGCCAATATGCGTATCGGTCAGCCGTTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGAAACTGGATAAAAGTCGTCCAGCGTCGGCCTGTGCGAATCGAACGGGACCAGAAACAGCTGGAGCAATATCCGCTGctaaccggatggg                             >  2:679684/1‑226 (MQ=255)

                                                                                                                                                                                                  aCGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCCAAGTCGTTCAGCGTCTTCCTGTGCGTATCGAACTGGACCCGAAACAGCTGGGGCAATAACCGCTGCGTATCGGttttt                             >  1:447138/1‑236 (MQ=255)

                                                                                                                                                                                                   cGCAGATTGCCAATTGGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGt                             <  1:335663/237‑1 (MQ=255)

                                                                                                                                                                                                            ccAATATGCGTTTCGGTTAGCCGGTCAGTATCCCCACGGATATTTACGGCGATGATGTGAAATACACCGGGAAAGTTGTTGGTCTGGATAGGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTg                  <  2:447138/239‑1 (MQ=255)

                                                                                                                                                                                                            ccAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGGGCAATATCCGCTGCGAATCGGTTTTTCCACCCTgtga                  >  1:289376/1‑236 (MQ=255)

                                                                                                                                                                                                            ccAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACg                        <  1:221956/233‑1 (MQ=255)

                                                                                                                                                                                                            ccAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTg                  <  1:554071/239‑1 (MQ=255)

                                                                                                                                                                                                            ccAATATGCGGAGCGGTCCGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGGGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTg                  <  1:91693/239‑1 (MQ=255)

                                                                                                                                                                                                             caagttgcctttCGGTCAGCCGGTCACTATCGCCACGGATATTGACGGCGAGGATGTGAAATCCAACGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTACCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGAGTATCGAACTGGACCAGAAACAGCTGGAGCAATAGCCGCTGCGTATCGGTTTGTCCACGCTGGTg                  <  1:169493/227‑1 (MQ=255)

                                                                                                                                                                                                             cAATATGCGTATCGGTCAGCCGGTCACTATCCCCACGGATATTTACGGCGATGATATGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTg                  <  1:399041/238‑1 (MQ=255)

                                                                                                                                                                                                               attggcttATCGGTCAGCCGGTCCCTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CTCTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGc               <  1:556900/232‑1 (MQ=255)

                                                                                                                                                                                                                 tatgcagtCGGTTAGCCGGTCACTTTCTCCACGGATATTTACGGCGATGATGTGAAATTCACCGGTAAAGTGGTTGGGCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATAGCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGt             <  2:23210/232‑1 (MQ=255)

                                                                                                                                                                                                                   tggttttggTCAGCCGGGTCCTATCACCCAGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GGGGGTTCT‑CACTGCTTCCA‑GCGGAAAATGCGACCGGTAACTGGATCAAAGGCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCa           <  2:809393/232‑1 (MQ=255)

                                                                                                                                                                                                                   gCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCa           <  1:600673/239‑1 (MQ=255)

                                                                                                                                                                                                                    ggtttCGGTCCGCCGGTACCTGTCACCCGGGATAGTTACTGCGATGAGGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGTCCA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCaa          <  1:276122/235‑1 (MQ=255)

                                                                                                                                                                                                                      tATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGGGAGCGTCAAt         >  1:371649/1‑238 (MQ=255)

                                                                                                                                                                                                                      tATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGAATCGAACTGGACCAGAAACAGCTTGAGCAGTATCCGCAGCGTATCGGTTTGTCCACGCTGGTGAGCgaccat         >  1:717131/1‑233 (MQ=255)

                                                                                                                                                                                                                      tATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGGGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCCCTCGGTCTCGGGTTGTCCACGCTGGGAAGCggcttg         >  2:221956/1‑233 (MQ=255)

                                                                                                                                                                                                                      tATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGCGCGTATCGAACTGGACCAGCAACCGCTGGAGCCATATCCGCTGCCTAACGGttttgccccgcgggggcggggaaaaa        >  2:468091/1‑216 (MQ=255)

                                                                                                                                                                                                                      tATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGAA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGGCCAAAATCGTTCCGGGGCGGCCTGTGCGCATCGAACTGGACCATAAAAAACTTGAGCAAGTGCCGCTGCGGACCGGTTTGTTCACACTGGTGaaggtgaaa         >  2:169493/1‑230 (MQ=255)

                                                                                                                                                                                                                       aTCGGGCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATTCACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCAATAc       <  1:773944/239‑1 (MQ=255)

                                                                                                                                                                                                                           gTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCAATACCACt   >  1:546590/1‑239 (MQ=255)

                                                                                                                                                                                                                            tCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCAATACCACTa  >  1:503231/1‑239 (MQ=255)

                                                                                                                                                                                                                                   gTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATAACCGCTCCGTATCGGCTTGTCCACGCTGGGGGGCGTCAATACCACt   >  1:23210/1‑231 (MQ=255)

CCAGCCGTGCAACAGGCTGCCACCGAAGTACGTAACGCCTGGCTGGCGCTGGAGCGTACTCG‑TATTATCAGTCCGATGACCGGTTATGTCTCCCGCCGCGCGGTACAGCCTGGGGCGCAAATTAGCCCAACGA‑CGCCGCTGATGGCGGTCGTTCCAGCCACCAATATGTGGGTGGATGCCAACTTTAAAGAGACGCAGATTGCCAATATGCGTATCGGTCAGCCGGTCACTATCACCACGGATATTTACGGCGATGATGTGAAATACACCGGTAAAGTGGTTGGTCTGGATATGGGCACA‑GGTA‑GCGCGTTCT‑CACTGCTTCCA‑GCGCAAAATGCGACCGGTAACTGGATCAAAGTCGTTCAGCGTCTGCCTGTGCGTATCGAACTGGACCAGAAACAGCTGGAGCAATATCCGCTGCGTATCGGTTTGTCCACGCTGGTGAGCGTCAATACCACTA  >  NC_000913/2812013‑2812469

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