BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	693,319	T→C	10.7%	E174G (GAG→GGG)	Y75_RS03425 ←	PhoH family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	693,319	0	T	C	10.7%	342.6 / 4.1	134	E174G (GAG→GGG)	Y75_RS03425	PhoH family protein
Reads supporting (aligned to +/- strand): ref base T (83/32); new base C (3/15); total (87/47)
Fisher's exact test for biased strand distribution p-value = 1.04e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

ATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATC  >  NC_007779/693085‑693549

aTAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTcc                                                                                                                                                                                                                                         <  2:1677643/237‑1 (MQ=255)

        gtcgtTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCaa                                                                                                                                                                                                                               <  1:734981/239‑1 (MQ=255)

        gtcgtTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCa                                                                                                                                                                                                                                >  1:1013302/1‑238 (MQ=255)

         tcgtTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCa                                                                                                                                                                                                                                >  1:34708/1‑237 (MQ=255)

                gCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCaa                                                                                                                                                                                                                               >  2:69184/1‑231 (MQ=255)

                gCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCaa                                                                                                                                                                                                                               <  1:69184/231‑1 (MQ=255)

                 cGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCa                                                                                                                                                                                                                       >  2:515365/1‑238 (MQ=255)

                  gTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGc                                                                                                                                                                                                                     <  1:1663690/239‑1 (MQ=255)

                  gTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGc                                                                                                                                                                                                                     <  1:96942/239‑1 (MQ=255)

                  gTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGc                                                                                                                                                                                                                     <  1:1608507/239‑1 (MQ=255)

                    aCGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCa                                                                                                                                                                                                                   <  2:1604636/239‑1 (MQ=255)

                      gACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGcc                                                                                                                                                                                                                    >  1:312371/1‑236 (MQ=255)

                       aCCACGAACATAGGACAGCGCAGCGACTACAAAACCGGTGCGCAAATTCAGTTCCTCATCTTTCTCAAAGCCAAGCCTTTCAAACAGCGCGTCGTACAGCGGGCTGA‑GATATGGGTCTACTTTCTGGCTTAAATCGCCAGCAAGGAAGCCCAGTTTCTCCCC‑GGC‑TTCTTCTGCAGAACGAGTCAGCGGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACt                                                                                                                                                                                                                 <  2:119205/238‑1 (MQ=255)

                           cGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTg                                                                                                                                                                                                                <  2:938436/235‑1 (MQ=255)

                           cGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCag                                                                                                                                                                                                            >  1:1120461/1‑239 (MQ=255)

                            gCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCagg                                                                                                                                                                                                           <  1:1639804/239‑1 (MQ=255)

                            gCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCagg                                                                                                                                                                                                           <  1:515365/239‑1 (MQ=255)

                                  aGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGtt                                                                                                                                                                                                     >  1:843447/1‑239 (MQ=255)

                                      cAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCCTCAACTGCCGCAGCCACTGCCAGGTGGGGTTta                                                                                                                                                                                                  >  2:1652642/1‑238 (MQ=255)

                                      cAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTta                                                                                                                                                                                                  >  1:712027/1‑238 (MQ=255)

                                      cAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTta                                                                                                                                                                                                  >  1:153074/1‑238 (MQ=255)

                                       aGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtac                                                                                                                                                                                                 >  1:1138470/1‑238 (MQ=255)

                                       aGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtac                                                                                                                                                                                                 >  1:1739672/1‑238 (MQ=255)

                                        gCGGTGCGACTTCACTAACGGTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCTTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTAAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAACCACTGCCAGGTAGGTTTtaccc                                                                                                                                                                                               >  2:828524/1‑239 (MQ=255)

                                        gCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtaccc                                                                                                                                                                                               >  1:1334170/1‑239 (MQ=255)

                                         cGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTtacccg                                                                                                                                                                                              >  1:531350/1‑239 (MQ=255)

                                              cGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGAAGGTTTtacccgtaccc                                                                                                                                                                                         >  1:461749/1‑239 (MQ=255)

                                              cGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACCGCCAGGTAGGTTTtacccgtacc                                                                                                                                                                                          >  2:1115570/1‑238 (MQ=255)

                                                           gTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAAACGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACg                                                                                                                                                                             >  2:33753/1‑238 (MQ=255)

                                                            ttGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGcc                                                                                                                                                                           >  1:1785516/1‑239 (MQ=255)

                                                             tGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGGGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGcc                                                                                                                                                                           <  1:1674702/238‑1 (MQ=255)

                                                              gcgcTCAATCAGTTTTCCGACTTTCTCAAACCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGGGCA‑GATACGGGTCTACTTTCGGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGa                                                                                                                                                                         <  1:1841032/239‑1 (MQ=255)

                                                              gcgcTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGGGGGCGCA‑GATACGGGTCTCCTTTCGGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAAAATAC‑GGCGAATTTCCGGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGa                                                                                                                                                                         <  1:1260884/239‑1 (MQ=255)

                                                              gcgcTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTATTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCATTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCGGACGCCCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGa                                                                                                                                                                         <  1:1707786/239‑1 (MQ=255)

                                                               cgcTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGaa                                                                                                                                                                        >  2:686125/1‑239 (MQ=255)

                                                                          tttCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCAt                                                                                                                                                             >  1:1207111/1‑239 (MQ=255)

                                                                              tcGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCagg                                                                                                                                                                                                           >  2:1397169/1‑189 (MQ=255)

                                                                              tcGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTc                                                                                                                                                         >  1:1589756/1‑239 (MQ=255)

                                                                              tcGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTc                                                                                                                                                         >  1:752620/1‑239 (MQ=255)

                                                                              tcGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTc                                                                                                                                                         >  1:154081/1‑239 (MQ=255)

                                                                              tcGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTCAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCagg                                                                                                                                                                                                           <  1:1397169/189‑1 (MQ=255)

                                                                                    ttCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCCACGCCGAAGGTAATGTCATGGTCGAGaac                                                                                                                                                   >  2:1864052/1‑238 (MQ=255)

                                                                                         aatgCCCATCTTTACAAACGGCGGGTCGTCAAGCGGG‑GCAGGATTCGGTTCTCTGTTCTGGCTAAAATGGCCAGGCAGTTACCCCTTTTTCCCACC‑GTC‑TCCTCCGGCCGGACGGGTCATCATACCCC‑GTCTAATTCCCTGACCCCCCGTGCCTTCATCTCCCGCAGCCACAGACATTTAGTTTTTACCGGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTgg                                                                                                                                              <  2:1207111/236‑1 (MQ=255)

                                                                                           aGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCg                                                                                                                                            <  1:820393/239‑1 (MQ=255)

                                                                                                aGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGt                                                                                                                                        >  1:603209/1‑238 (MQ=255)

                                                                                                              cgcgTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTACATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACGGCCAGGTAGGTTTTACCCGTACCCGCCCGGCCAACGCCGACGGTAATGGCATGGTCGAGAATATTGGCGATGTACCGCGCCTGGCTTg                                                                                                                         >  2:114883/1‑239 (MQ=255)

                                                                                                              cgcgTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCaggtaggt                                                                                                                                                                                                      <  1:933588/162‑1 (MQ=255)

                                                                                                              cgcgTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTg                                                                                                                         >  1:1069380/1‑239 (MQ=255)

                                                                                                              cgcgTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTg                                                                                                                         >  1:1520239/1‑239 (MQ=255)

                                                                                                              cgcgTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGAAGTACTGGGCCTGGTTTg                                                                                                                         >  2:878344/1‑239 (MQ=255)

                                                                                                                    gtACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGTCGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGACTGGTTTGGCGTgc                                                                                                                   >  1:1390863/1‑239 (MQ=255)

                                                                                                                      aCAGCGGGCGCA‑GATGCGGGTCTACTTTCTGGCTAAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTgcg                                                                                                                  <  1:960404/238‑1 (MQ=255)

                                                                                                                            ggCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTa                                                                                                           >  1:1883365/1‑239 (MQ=255)

                                                                                                                             gcgcA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTaa                                                                                                          >  1:1637870/1‑239 (MQ=255)

                                                                                                                             gcgcA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTaa                                                                                                          >  1:344656/1‑239 (MQ=255)

                                                                                                                                        gggTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCCCCAGGGCATCAACTGTCGCAACCACTGCCCGGTAGGTTTTACCCGCACCCGCCGGGCCAACGCCGAAAGCAATGTCCTGGTACAGAATATTGGCGATGTCCTGCGCCTCGTTTTGCGAGTCCGGCTTAAttctgtcgcg                                                                                                >  2:1123378/1‑231 (MQ=255)

                                                                                                                                              aCTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTggg                                                                                        >  1:1651245/1‑238 (MQ=255)

                                                                                                                                              aCTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGt                                                                                        >  1:476021/1‑239 (MQ=255)

                                                                                                                                              aCTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGt                                                                                        >  1:364260/1‑239 (MQ=255)

                                                                                                                                                ttttcggCTTAAATCCCCCGGACAGACGCCCGATTTTTCCCC‑GGCTTTCTACTGCGGGCGGAGAC‑GCAAAAAACGGGCAAATTCCCTGCGGCCCCGGGGCAACACCGGCAGCGGCCACGGCCAGGTAGGTTTTACCCTTCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGtt                                                                                       <  2:296113/234‑1 (MQ=255)

                                                                                                                                                  tCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGAAATGTCATGGTCGAGACTATTGGCGATGTACTGCGCCTGGTTGGGCGTGCGCGGCTTACTTCCCCCGCGg‑tttggtttg                                                                                     >  2:1565077/1‑238 (MQ=255)

                                                                                                                                                    tGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGATTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGAACCGCGCCTGGTTTGGCGTGCGGGGCGTAATTAATCCGC‑GTTTTGGttttat                                                                                   >  2:1233235/1‑237 (MQ=255)

                                                                                                                                                          aagcccccgggcgggAAGCCATGTTTC‑CACCGGGC‑TTCTCCGGCGGAAGGAGCCACCAGAATCC‑GGGGATTTTCCGGCCGCCCCGGGGCATCACCGCCCCCGCCCCCTGCAAGGTAGGTTTTACCCGTCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGAc                                                                            <  2:1883365/230‑1 (MQ=255)

                                                                                                                                                               gcgaggagggACGCCCATTTTCCCCCC‑GGC‑GTCTACGGCGGGACGGCCCAGCAGAACAC‑GGCGATTTTCCGGCCCCCCCGGGGCATAAACGGCCGCAGCCGCTGCCAGGTAGGTTTTACCAGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCtt                                                                       <  2:1391402/231‑1 (MQ=255)

                                                                                                                                                                ccAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCttt                                                                      >  2:626082/1‑239 (MQ=255)

                                                                                                                                                                ccAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGGTTTGATATTGACCGCttt                                                                      >  1:1034510/1‑239 (MQ=255)

                                                                                                                                                                ccAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGGTTTGATATTGACCGCttt                                                                      >  1:538191/1‑239 (MQ=255)

                                                                                                                                                                     caggAGGCCCATTTTCTCCCG‑GGT‑TTCTACTGCCGCAGCAGTCAGCACAACAG‑GGCAAATTCCGTACCGCCCCGGGGCATCAACGCCCCCGCCCACTGCCGGGTAGGTTTTACCCGTCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCg                                                                  <  2:116510/238‑1 (MQ=255)

                                                                                                                                                                           gCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCg                                                           >  1:598206/1‑239 (MQ=255)

                                                                                                                                                                               aGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTATGCCGC‑G‑GTTGGTGTTGATATTTATCGCCTTGGCGCAGTCCGG‑cg                                                       >  2:935809/1‑237 (MQ=255)

                                                                                                                                                                                    ctcACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCCTAATTACGCCCC‑G‑GTTGGGTTTTGTATTGACCGCTTTGCCGGACTCCGGCCACGct                                                   >  1:861451/1‑238 (MQ=255)

                                                                                                                                                                                     tcACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAAAATAC‑GGCGAATTCCCTGCCGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCtctc                                                <  1:192100/239‑1 (MQ=255)

                                                                                                                                                                                      cccc‑GGT‑TTCTACTGCCGGACGAGTCAGCAGAATCC‑GGCGAATTTCCTGCCGCTCCAGGGCACCAACGGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTcc                                               <  1:1204324/237‑1 (MQ=255)

                                                                                                                                                                                      cACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAACCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGGAATGTCATGGTCGGGGATATTGGCGATGGACTGCGCCTGGTTTGGCGTGCCCGGCCTTATTACGCCGC‑G‑TTTGGTTTTGTTATTGACCCCTTTGCCGTATCCCGG‑CACCCtctc                                                >  2:146860/1‑238 (MQ=255)

                                                                                                                                                                                         c‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGCAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGGTCAATTACGCCGC‑G‑TTTGGATTTGATGTTGACCGCTTTGCCGTCCTACGG‑CCCCGTCTCcgc                                             >  2:756890/1‑238 (MQ=255)

                                                                                                                                                                                         c‑GGC‑TTCTACGGCCGGACGAGTCAGCAGAATCC‑GGCGAATTCCCTGACGCCCCAGGGCATCAACGGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCcgcg                                            <  1:166064/239‑1 (MQ=255)

                                                                                                                                                                                               ttCTCCTGCCGGACGAGCCGCAAGAAAAG‑GGGGAATTTCCTGACCCCCCGGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCttt                                        <  1:935809/239‑1 (MQ=255)

                                                                                                                                                                                                 cTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCttt                                        >  1:859081/1‑237 (MQ=255)

                                                                                                                                                                                                   aCTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTgcg                                                                                                                  <  1:89781/164‑1 (MQ=255)

                                                                                                                                                                                                   aCTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTgcg                                                                                                                  >  2:89781/1‑164 (MQ=255)

                                                                                                                                                                                                   aCTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCt                                     >  1:634437/1‑238 (MQ=255)

                                                                                                                                                                                                       ccGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTTCAGt                                >  1:848059/1‑239 (MQ=255)

                                                                                                                                                                                                       ccGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGt                                >  1:789846/1‑239 (MQ=255)

                                                                                                                                                                                                                 cagcagAATCC‑GGCGAATTCCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGct                       <  1:1782150/238‑1 (MQ=255)

                                                                                                                                                                                                                        ttgc‑GGCGAATTCCCTGCCGCCCACGGGCATCACCGGCCGCAGCCACGGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTa                 <  1:1038822/234‑1 (MQ=255)

                                                                                                                                                                                                                         tAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGGTTTGGTATTGACCGCTTTGCAGTACTCCCG‑CACGCCCCCCGCGCTTTGCTCCAGGACCCCCGCCTCTTTTAt               >  2:1655069/1‑238 (MQ=255)

                                                                                                                                                                                                                         tAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGCAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTTACCGCTTTGCCGTCCTCCGG‑CACGCTCCCCGCGCTTTGCTCCAGTACCCGCGCCTCTTTAAt               >  1:308940/1‑238 (MQ=255)

                                                                                                                                                                                                                             ggCGAACTGCCGGACGCCACGGGGCATCAACTGCCGCGGCCGCTGCCAGGTCGTTTTTACCCGGCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCa           <  2:1872548/239‑1 (MQ=255)

                                                                                                                                                                                                                              gCGATTTCCCGGCCGCCCCGGGGCATCACCTCCCGCACCCACTGCCAGGTGGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCACGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCa           <  2:344656/238‑1 (MQ=255)

                                                                                                                                                                                                                               cGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCCC‑G‑TTTTGGTTTTATATTTACCGCTTTGTCGTACTCCGG‑CACGCCCTCCGCGCCTTTCTTCAGGACCCCCGCTTCTTTAATCGCAAg         >  2:498564/1‑239 (MQ=255)

                                                                                                                                                                                                                               cGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCCCGGCTTAATTACGCCCC‑G‑GTTTGTTTTGATATTGACCGCTTTTCCGTAACCCGG‑CACGCCCCCCGCGCTTTGCTTCAGTACCCCCCGCTCTTTAATAGCaat         >  2:1722811/1‑238 (MQ=255)

                                                                                                                                                                                                                                ggaTTTCCGGGCGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAgg        <  1:1766468/237‑1 (MQ=255)

                                                                                                                                                                                                                                  ttttCCTGCCCCCCCCGGGCATCCACGGCCGCGCCCACGGCCGGGTGGTTTTTGCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAgg        <  2:1069380/236‑1 (MQ=255)

                                                                                                                                                                                                                                  aTTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGGCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCCAgggg      >  1:508531/1‑237 (MQ=255)

                                                                                                                                                                                                                                  aTTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTg      >  1:1603295/1‑239 (MQ=255)

                                                                                                                                                                                                                                   tttCCTGACGCTCCAGGGCATCACCTGCCGCAGCCACTGCCAGGTAGTTTTTCCCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTgg     <  1:1787644/239‑1 (MQ=255)

                                                                                                                                                                                                                                     tCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGAt   <  1:686125/239‑1 (MQ=255)

                                                                                                                                                                                                                                     tCCTCACGCCCCAGGGCACCACGTGCCACACCCAGTGCCAGTTAGGTTTTCCCCTTCCCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGAt   <  2:1722711/239‑1 (MQ=255)

                                                                                                                                                                                                                                      ccTGACGCCCCAGGCCATCAGCGGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATc  <  1:245688/239‑1 (MQ=255)

                                                                                                                                                                                                                                       cTGACGCGCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTTGGTTTGATATTGACCCAGTTGCCGGACTCCGG‑CACGCCCTGCGCTCTTGTCTACAGTACTCGCGCTGGTTTTATCGCacgaggga    >  2:1781573/1‑229 (MQ=255)

ATAAACGCGTCGTTCAGCGTACGACCACGCATATAGGCCAGCGGTGCGACTTCAATAACGTTGCGCTCAATCAGTTTCTCGACTTTCTCAAAGCCCAGCATTTCAAACAGCGCGTCGTACAGCGGGCGCA‑GATACGGGTCTACTTTCTGGCTTAAATCGCCAGGCAGGAAGCCCAGTTTCTCACC‑GGC‑TTCTACTGCCGGACGAGTCAGCAGAATAC‑GGCGAATTTCCTGACGCTCCAGGGCATCAACTGCCGCAGCCACTGCCAGGTAGGTTTTACCCGTACCCGCCGGGCCAACGCCGAAGGTAATGTCATGGTCGAGAATATTGGCGATGTACTGCGCCTGGTTTGGCGTGCGCGGCTTAATTACGCCGC‑G‑TTTGGTTTTGATATTGACCGCTTTGCCGTACTCCGG‑CACGCTCTCCGCGCTTTGCTCCAGTACCCGCGCTTCTTTAATCGCAAGGTGGATC  >  NC_007779/693085‑693549

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