BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	484,581	A→C	9.8%	W347G (TGG→GGG)	acrA ←	multidrug efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	484,581	0	A	C	9.8%	372.2 / 3.4	135	W347G (TGG→GGG)	acrA	multidrug efflux system
Reads supporting (aligned to +/- strand): ref base A (69/44); new base C (8/10); total (77/58)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

CATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGA  >  NC_000913/484344‑484816

cATGATGATAATGGCAGTCCCCCACGAAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCAc                                                                                                                                                                                                                                            <  2:195856/239‑1 (MQ=255)

 atgatgATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTAGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTTCCTGGACACCCGGCCGCCCTTTCTGCA‑ACCCCCTTATTCCTACTCGCTCGCCTGATTTCAAACCATCTGTCACCCGCCACt                                                                                                                                                                                                                                           >  1:797706/1‑239 (MQ=255)

   gaggatAATGTCGATCCCCCACGGAAAAATGGGGAGATCAATAAAGAAATTGGGCATTTCTTAACGGTTCGTGTTTAAGTTAAGACTTCGCCTGTTCGGGCTGGGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTACCTTCTTGTGCTTTTACCTGTACCCCAGGCAGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTa                                                                                                                                                                                                                                         <  1:1072250/236‑1 (MQ=255)

    aagatAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTAt                                                                                                                                                                                                                                        <  2:811763/237‑1 (MQ=255)

     agatAAGGGCGATCCCCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATTGCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTAt                                                                                                                                                                                                                                        <  2:671170/237‑1 (MQ=255)

          aTGGGGGTCCCCCACGAAAAAATGGGGCGATTGAAAAAAAAATGAGGCATGTTTAAACGGCTCCGGTTTAAGTTAAGACTTGGGTTGTTCAGGCTGAGCACCGTTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCaa                                                                                                                                                                                                                                  <  1:312265/239‑1 (MQ=255)

           tGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCaaa                                                                                                                                                                                                                                 >  2:1252124/1‑238 (MQ=255)

                aTCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTTCTTGTTATTATCAGCGGTAGCTTCCTGTGCTTTTACCTGGTCACCAGGACGAACTTTCTGCA‑GACCACTTTTTACTACGCGAACGCCTGCTTTCAAACCTTCTGTCACCCGCCACTTATCCCCAAAAggcc                                                                                                                                                                                                                            >  1:726896/1‑236 (MQ=255)

                aTCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCCATAg                                                                                                                                                                                                                               >  2:1022445/1‑236 (MQ=255)

                aTCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCACCCGCTTTCAGACCCTCAGTCACCCACCCCTTAAACCCAATAAccc                                                                                                                                                                                                                            >  1:536790/1‑238 (MQ=255)

                 tCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGAACCCAGGACGCACTTTCTGCA‑GCCCACTTTTTCCTACCCGTACGCCTGCTTTCAGACCTTCTGTCACCCGCCCCTTAACCCCAAAAGccgg                                                                                                                                                                                                                           >  1:863136/1‑237 (MQ=255)

                   aCCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCCGGACGCCCTTTCTGCC‑GCCCCCTTATTACTACGCGATCGCCCTCTTTCAACCCTTCTTTCCCCACCCCCTTATcccccaaaccccg                                                                                                                                                                                                                           >  1:1168660/1‑228 (MQ=255)

                       aCGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCCCTTTCTGCA‑GCCCACTTATTTCTAACCGATACCCTGCTTTCCGACCTTCCTTCCCCAGCCCCTTATCGCCCAAAACCTGGCTTGc                                                                                                                                                                                                                     >  1:71269/1‑239 (MQ=255)

                                 tgggggtttgatAAAGAAATTTGGCATGTTTTAACGGCTCCGTTTTAAGTTAAGATTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTCCCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGa                                                                                                                                                                                                           <  1:199259/230‑1 (MQ=255)

                                 tggggGATCTATAAAGAAATTGGGCATGTCTTAACGGCTCCGGTTTAAGTTAAGGCTGGGGCTGTCTAGGCTGAGCCCCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGa                                                                                                                                                                                                           <  1:435408/238‑1 (MQ=255)

                                                          tttttttACCGTCTCCGTTTTAAGTTAAGCTTTGGCCTTGTCATGCTGAGCACTGCTCGCGGTCTGCTGGTTATTATGAGCAGTAACTTCTTGTGCTTTTCGCTGGACACCGGGACGCACTTTCTGCA‑GCCCTCTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCCCTTTGTCATCcgcg                                                                                                                                                                                   <  1:984449/233‑1 (MQ=255)

                                                           ttttttAAGGGCTCTGTTTTTAGTTAAGGTTTGGACTGTTCAGGCTGGGCACCTGTTGCGGCGTGCTGGTTATTAGTCGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCGTCA‑GCCCACTTATTCCTTCGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGcc                                                                                                                                                                                 <  1:930251/235‑1 (MQ=255)

                                                            ttCTTAACGGGTCCTGTTTAAGTTAAGGCTTTGACTGTTCAGGGTTAGCTACGCTTGCGGCTTGCTGGTTATTACCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCa                                                                                                                                                                                <  2:288775/238‑1 (MQ=255)

                                                            ggttttaCGGCTCTTGTTTAATGTAAGCTTGGGGCTGTTTGGGCGGAGCACCGCTTGCGGCTTGGTGGTTATCATCGGCGGTACCTTCTTGTGCTTTTTCCGGGACCCCAGGACGCACTTTCGGCA‑GCCCCTGTATTACTACGCTATCGCCTGATTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCa                                                                                                                                                                                <  1:971814/233‑1 (MQ=255)

                                                                 aaCGGCTCCTTTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTAGCAGCGGTAACTTCTTGTGCTTTTACCTGGACCCCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTAc                                                                                                                                                                           <  2:856712/239‑1 (MQ=255)

                                                                  aCGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCATTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTAcc                                                                                                                                                                          >  2:1198432/1‑239 (MQ=255)

                                                                    ggCTCCTGTTTAAGTTAAGTCTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCCCTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAg                                                                                                                                                                        <  2:688914/239‑1 (MQ=255)

                                                                      cTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCTGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAAACCTTCTTGCACCCGCCCCTTATCGc                                                                                                                                                                                                                                     >  1:263554/1‑176 (MQ=255)

                                                                      cTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCCCTTATCGCCAAAAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTTTCATCCCCGCCAACAACCaaga                                                                                                                                                                      >  1:319970/1‑236 (MQ=255)

                                                                       tCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGCCGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTAc                                                                                                                                                                     >  2:1046673/1‑239 (MQ=255)

                                                                       tCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGGCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTAc                                                                                                                                                                     >  2:783503/1‑239 (MQ=255)

                                                                       tCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACC‑TCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTTTCATCCGCCCCAACTACCAGTACg                                                                                                                                                                    >  2:822356/1‑239 (MQ=255)

                                                                        ccTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATATCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTTTCTACCGCGCCAACTCCCCGTACg                                                                                                                                                                    >  2:1214305/1‑239 (MQ=255)

                                                                        ccTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGCCGGGTTTCCACTTTGTCATCCGGGCCAACTACCCGTAca                                                                                                                                                                    >  2:887131/1‑238 (MQ=255)

                                                                                    aagACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGcc                                                                                                                                                         >  2:1116929/1‑238 (MQ=255)

                                                                                      ggtttGGACTGTTCATGTTGGGCGCCGCTTGGGGTTTGCTGGTTATTATCAGCTGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCAcg                                                                                                                                                      <  1:1266850/236‑1 (MQ=255)

                                                                                        cTTGGACTTTTCAGGCTGAGCACCGCTTGCGGCTTGCGGGGTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTAACGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCAcg                                                                                                                                                      <  2:1132548/237‑1 (MQ=255)

                                                                                           ggACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCa                                                                                                                                                                                                                    >  2:841388/1‑172 (MQ=255)

                                                                                               tGTTCAGGCTGAGCCCCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTAc                                                                                                                                              <  2:536790/238‑1 (MQ=255)

                                                                                               tGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTAc                                                                                                                                              <  2:919494/238‑1 (MQ=255)

                                                                                                 ttCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGa                                                                                                                                                                                                           >  1:905066/1‑175 (MQ=255)

                                                                                                 ttCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGcc                                                                                                                                                                                 >  1:610799/1‑201 (MQ=255)

                                                                                                 ttCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGcc                                                                                                                                                                                 <  2:610799/201‑1 (MQ=255)

                                                                                                 ttCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGGTTCCACTTTGTCATCCCCGCCAAGTACCAGTCCGAGGGCATCGCCACGCGGCGTAgggg                                                                                                                                           >  2:1005474/1‑239 (MQ=255)

                                                                                                   ccgtgtgtgccccGTTTGCGGTTTTCTTGTTATTATCAGCGGTTACTTCTTGTTCTTTTCCCTGGACCCCAGGACCCCCTTTCTGCA‑GCCCCCTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGtt                                                                                                                                         <  1:1252124/228‑1 (MQ=255)

                                                                                                        tGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGc                                                                                                                                     <  2:911451/238‑1 (MQ=255)

                                                                                                           gCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTCCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCCCTTTGTCAATCGCGCCAACTACCCGGACGACGGCGCCACCACCCGGCGTACgggcaaccccccg                                                                                                                                 >  1:128915/1‑229 (MQ=255)

                                                                                                                gCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGGTTCCACTTTGTCATCACCGCCCACTCCCAGAACGGTGGCATCGCCACCCTGCGGACGGGTACCGCCCCGTTGCg                                                                                                                            >  1:478116/1‑239 (MQ=255)

                                                                                                                 cTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACCCGGCGTACGGGTTACGCCCTGTTGCgg                                                                                                                           >  2:920958/1‑239 (MQ=255)

                                                                                                                    gCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGAc                                                                                                                        >  2:1101900/1‑239 (MQ=255)

                                                                                                                         ttcttGGTTTTTATCGGCGTTCTCTATTTGTGCTTTTACGTGGACACCCGTACGCAC‑TTCTGGATGCCCACTTATTCCTAGTGGATCTCCGGCTTTCAGACCTTCTGTGACCAGCCACTTATCGCCAAGAGCCTGGCTTGCAACGATCGGAGGGGTTTCCACTTTGTCATCCGCGCCAACTTGCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCGGTTGCGGGACTaaaa                                                                                                                   <  1:879310/235‑1 (MQ=255)

                                                                                                                         ttGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTaaaa                                                                                                                   <  2:57849/239‑1 (MQ=255)

                                                                                                                                       tcggggTAACTTCTTGTGCTTTTACCTGGACCCCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTa                                                                                                     <  2:726896/235‑1 (MQ=255)

                                                                                                                                       ctgCGGTAACTTCTTTGGTTTTTACCGGGGCCCCAGGACACACTTTCTGCA‑CCCCACTTATTACTACGCTATCCCATTGTTTCAGACCGTCTGTCACCAGCCCCTTATCGACAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTa                                                                                                     <  1:1116929/237‑1 (MQ=255)

                                                                                                                                        ggtggtcgtgtctttgtcTTTTCCCGGGACCCAAGGGCGCACTTTCTGCT‑CGCCATGTGTTACATGGTGATTGCCTGGTTTAAGACCGTCTGGCGCCAGCCAATGTTCGGCAAGGGCCTAGCTTGCAACGATCGGACGGGTTTCCGCTTTGTCATCCGCGCCAACTACCAGTGCGGTGGCATCGCCGCGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTaa                                                                                                    <  1:336330/222‑1 (MQ=255)

                                                                                                                                         gCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTaa                                                                                                    >  2:344407/1‑238 (MQ=255)

                                                                                                                                          cGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCTCCCACTACCAGTACGGTGGAATCGCCACGCGCCGCAAGGGTTACGCCCTGTTGCAGGACTAAAAACGCGTTTGGATTaaaa                                                                                                  >  1:1001475/1‑237 (MQ=255)

                                                                                                                                          cGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GACCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTAAGGGTTAAGCCCTGTTGCGGGACTAAAATGGCGTTTGGAtacag                                                                                                   >  2:184081/1‑234 (MQ=255)

                                                                                                                                           ggttattttttGTGCTTTTCCCGGGACCCAAGGCCCTCCTTTCGGCG‑GCCCACTTTTTACTACGCAGTCGCCGTGTTTCAGGCCTTCGGTCCCCAGCCCCTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCACCTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGcc                                                                                                 <  1:1046673/230‑1 (MQ=255)

                                                                                                                                            gTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGcc                                                                                                 <  2:1168660/238‑1 (MQ=255)

                                                                                                                                             tAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGCCTAAAATAGCGTTTGGATTAAGCCct                                                                                               >  2:243125/1‑239 (MQ=255)

                                                                                                                                             tAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTGCAGTCCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGGTTCCACTTTGTCATACGCGCCAACTACCAGGACGCTGGCCTCGCCACGCGCCGTACGGGTTACGACCTGTTGGGGGACTCAAATAGCGGTGGGCTTACGGCct                                                                                               >  1:546309/1‑239 (MQ=255)

                                                                                                                                              aaCTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGGTTCCACTTTGTCATCCGCGCCAACTACCAGTACGTGGGCATCGCCCCCCGCCGACAGGGTTACGCCCTCTTGCGGGAGTAAAATAGCCTTTGGATGAAGCCct                                                                                               >  1:99452/1‑238 (MQ=255)

                                                                                                                                               aCTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTCCGGTGGCATCGCCACGCGGCGCACGGGTTACGACCTCTTGCGGGACTAAAAAAGCGTTTGGATTAAGcccc                                                                                               >  2:336330/1‑236 (MQ=255)

                                                                                                                                                  tcttGGTCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCCTTTAGTACTAGGCCATCGCCTTCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGGGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCttcttc                                                                                          <  2:561672/239‑1 (MQ=255)

                                                                                                                                                    ttGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATAGCCAAAAGCCTGGCTTGCAACGATCCGACGGGCTTCCACTTTGTCATCCGCGCCAACTACCCGTACGCTGGCATTCCCCCCCCGCGTACGGGGTCCGCCCTGTTGCGGGCCTCAAAAAGCGCTTTGAATAAGCCCCtcctccc                                                                                        >  1:626570/1‑238 (MQ=255)

                                                                                                                                                           tttACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTgc                                                                                 >  2:1110875/1‑239 (MQ=255)

                                                                                                                                                           tttACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTATTGCGGGACTCAAATAGCGTTTGGATTAAGCCCTTCTTCCATACGTgg                                                                                 >  2:973395/1‑238 (MQ=255)

                                                                                                                                                            ttACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCACCAATAGCCTGGCTTGCAACGATCGGACGGGGTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGCCGTACCGGTTACGCCCCGTTTCGGGCCGAAAATAGACTTTGGATTAACCCCTTCTTCCACACggccg                                                                                >  1:1175147/1‑235 (MQ=255)

                                                                                                                                                              tccgtgcccgccgtaCGCACTTTCTGCA‑GCCCATGTAGTCCTGAGTTATTCCGTGTGTTCAGACCGTCGGTCCCCCGCCCCTTATCGCTAAGAGCCTGGCCTTGAACGATCTGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTTCTGTGGCATCGTCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTgcgcg                                                                              <  1:1214305/225‑1 (MQ=255)

                                                                                                                                                              aCCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGTTTGCAACGATCGGAAGGGATTCCACTTTGTCATCTGCACCCACTACCAGTAAGGAGACATCGCCCCGCGGCGGACGGCTTAAGCCCAGTTGCGGGCTTAAAATAGGTGTCGGATAAAGGCCATCTTACacaagaacggg                                                                              >  1:352717/1‑229 (MQ=255)

                                                                                                                                                                      acCCGGACGCCCTTTCTGCA‑GCACACTTATTACTACGGGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCACCCATAGCATCGCTTGCAACGATCGGGCCGGTTTCCACTTTGCCATCCGCGCCAACTACCAGTACCGTGGAATAGCAACGCGCCGTACGCGTTACGCCCTGTTGCGTGACTTAAATAGCACTCGCATTCAGCCGTTCTTCCAGACGTGGCCGAACGacac                                                                      >  1:100513/1‑236 (MQ=255)

                                                                                                                                                                      acCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACCCGATCGCCTGCTTTCAGACCTTCTGTCTCCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGt                                                                                                                                                <  2:852360/165‑1 (MQ=255)

                                                                                                                                                                        cAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGGCTAACATAACGTTTGGATTAAACCCTTCTTCCAAACGCGCGGGCACGAACac                                                                     >  1:39166/1‑237 (MQ=255)

                                                                                                                                                                        cAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCCCCAACTACCAGTACGGTGGCATCGCCACGCCGCGTACGGGTTACGCCCCGGTGCGGGACTAAAAAAGCGATTGGATTAAGCCCTTCTTCCAGACGTgcgccaacaaaaaa                                                                     >  1:123715/1‑228 (MQ=255)

                                                                                                                                                                             cGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGCCTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCg                                                                >  2:875711/1‑238 (MQ=255)

                                                                                                                                                                                aCTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGGGCGCGCCCGAACTAACCCTgca                                                             >  2:1102432/1‑238 (MQ=255)

                                                                                                                                                                                  tttCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGgcagca                                                          <  2:71269/239‑1 (MQ=255)

                                                                                                                                                                                    tCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCCCTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCAcg                                                                                                                                                      <  1:405472/145‑1 (MQ=255)

                                                                                                                                                                                    tCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCAcg                                                                                                                                                      >  2:405472/1‑145 (MQ=255)

                                                                                                                                                                                             ccACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTCAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGACCATACCCGGCAGCAGGgtggcat                                                 >  2:7438/1‑234 (MQ=255)

                                                                                                                                                                                                 ttattcCTACGGCATCGCCTGCTTTCAGACCTTGTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCggg                                            <  2:478116/239‑1 (MQ=255)

                                                                                                                                                                                                 ttattaCTACGCGATCGCCTGCTTTCAGGCCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCggg                                            <  2:1241939/239‑1 (MQ=255)

                                                                                                                                                                                                 ttattaCTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCggg                                            <  2:863136/239‑1 (MQ=255)

                                                                                                                                                                                                      actacGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGACCTTCTTCCAAACGTGAGCGCCCGAACATACCCGGCAGCAGAGTGTGATACGGGTTCgg                                       >  2:266886/1‑239 (MQ=255)

                                                                                                                                                                                                        ctccgcaTCCCCTGCTTTCAGCCCTTCTGTCCCCACCCCCTTATCGCCATTAGCCGTGCTGGCAACGGTCGGACGGGTTTCCCCTTTGTCATCCGCGCCAACTTCCAGTACGGGGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGa                                     <  1:1022445/233‑1 (MQ=255)

                                                                                                                                                                                                          cgcgATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTTCGCGCACGAACATACCCGGCAGCAGCGTGTGATCCGGGTTCGGGAAg                                   >  2:575445/1‑239 (MQ=255)

                                                                                                                                                                                                            cgATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGAt                                 <  2:319970/239‑1 (MQ=255)

                                                                                                                                                                                                               tCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAgc                              >  2:1232591/1‑239 (MQ=255)

                                                                                                                                                                                                                  ccTGCTTTCAGACCTTCTGTCACCAGCCACTTATAGCCAATAGCCTGGCTTGCAACGATCGGACCGGGTTCCACTTTGTCATCCGCGCCAACTACCAGTCCGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCCTCTTCCAACCCAGCGCGCCCGAACCTACCTCGCAACAGCCGGTAGTCCGGCTTCCGGAagagaggcc                            >  1:242914/1‑235 (MQ=255)

                                                                                                                                                                                                                    tgtgttCAGCCCTGCTGGCACCAGCCTCTTAGCGCCAATATCCTGGCTTGCAACGACCGGACGGGTTGCCACTTTGTCATCCGCGCCAACTACCAGTACGTTGGCATCGCCACGCGGCGTACGGGTGACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTa                         <  2:903020/235‑1 (MQ=255)

                                                                                                                                                                                                                    tggtttCAGACCTTCTTGCCCCAGCCACTTATCGCCAATGGCCTGGCTTGCAACGATCGGACGGGTTTCCACTCTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTa                         <  2:951889/236‑1 (MQ=255)

                                                                                                                                                                                                                    tGCTTTCAGTCCTTCTGTCCCCAGCCACTTAGCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTa                         <  2:1095586/239‑1 (MQ=255)

                                                                                                                                                                                                                         tCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAggggg                    >  2:800960/1‑237 (MQ=255)

                                                                                                                                                                                                                           agcccTTCTGCCACCAGCCCATTATCTCCAATTGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGAt                  <  2:647603/236‑1 (MQ=255)

                                                                                                                                                                                                                              ccTTCTGTCCCCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTGCCAGTACGGGGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGAt                  <  2:797706/236‑1 (MQ=255)

                                                                                                                                                                                                                                ttCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGc                                                                                                                             <  1:1173651/127‑1 (MQ=255)

                                                                                                                                                                                                                                ttCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGc                                                                                                                             >  2:1173651/1‑127 (MQ=255)

                                                                                                                                                                                                                                ttCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCGGACGTGCGCGCACGAAAATACCCGGCCGAAGAGTGTGATCCGGGTTGGGGAAGAAAGGGGGGAGGGCGATAGaac             >  1:923762/1‑237 (MQ=255)

                                                                                                                                                                                                                                 tCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCCGCAGTGCGCGCACGACAATACCCGCCAGCACGGAGTGATCCCGGTTCGGGAAGATAGCCCGTCGGGGGCAAGAccc            >  1:576672/1‑239 (MQ=255)

                                                                                                                                                                                                                                  tcgTCCCCCGCCCCTTTTCGCCCATCGCCTGGCTTTGCACGGTCGGGCGGGTTTCCCTTTTGTCATCCGCGCCAACTTCCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCa           <  1:541956/237‑1 (MQ=255)

                                                                                                                                                                                                                                  cTGTCACCAGCCTCTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCa           <  2:128915/239‑1 (MQ=255)

                                                                                                                                                                                                                                  cTGGCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCa           <  2:1175147/239‑1 (MQ=255)

                                                                                                                                                                                                                                       cccACCCCCTTCTCGCCAATCGCCTGGCCTGCCACGATCGGACGGGTTTCCCTTTTTGCACCCGCGCCAACTACCAGGGCGGTGGCATCGCCACGCGGCGTCCGGCTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGAGCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGt      <  1:1267085/238‑1 (MQ=255)

                                                                                                                                                                                                                                            ccACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGGGCGTAGGGTGATAGACCCAGTGGTCTGa  >  2:734869/1‑238 (MQ=255)

CATGATGATAATGGCGATCACCCACGCAAAAATCGGGCGATCGATAAAGAAATTAGGCATGTCTTAACGGCTCCTGTTTAAGTTAAGACTTGGACTGTTCAGGCTGAGCACCGCTTGCGGCTTGCTGGTTATTATCAGCGGTAACTTCTTGTGCTTTTACCTGGACACCAGGACGCACTTTCTGCA‑GCCCACTTATTACTACGCGATCGCCTGCTTTCAGACCTTCTGTCACCAGCCACTTATCGCCAATAGCCTGGCTTGCAACGATCGGACGGGTTTCCACTTTGTCATCCGCGCCAACTACCAGTACGGTGGCATCGCCACGCGGCGTACGGGTTACGCCCTGTTGCGGGACTAAAATAGCGTTTGGATTAAGCCCTTCTTCCAGACGTGCGCGCACGAACATACCCGGCAGCAGAGTGTGATCCGGGTTCGGGAAGATAGCGCGTAGGGTGATAGACCCAGTGGTCTGA  >  NC_000913/484344‑484816

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