BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	1,674,232	T→C	5.6%	Y187H (TAC→CAC)	Y75_RS08420 →	serine protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	1,674,232	0	T	C	5.6%	459.6 / 10.0	145	Y187H (TAC→CAC)	Y75_RS08420	serine protease
Reads supporting (aligned to +/- strand): ref base T (65/71); new base C (4/4); total (69/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01

Only 100 of 146 total aligned reads displayed.

CTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGAC  >  NC_007779/1673995‑1674468

cTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCt                                                                                                                                                                                                                                              <  1:107931/238‑1 (MQ=255)

 tttGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGCCGTAAAGTGACTCAGGCAGGCTa                                                                                                                                                                                                                                             >  1:89094/1‑238 (MQ=255)

   tGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTAcc                                                                                                                                                                                                                                           >  2:598177/1‑238 (MQ=255)

    ggCGCTATGGGAGCCCCGACATGGAAGGCCGCGTTGATCCGCCCCTGGGAAAGCGGTTAAAAGCAAATGGGGATGGTTGGATTGTACCTCCCGCCGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAACCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCt                                                                                                                                                                                                                                         <  2:77177/239‑1 (MQ=255)

          aTGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGat                                                                                                                                                                                                                                   <  1:61987/239‑1 (MQ=255)

            gagaTCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatca                                                                                                                                                                                                                                 <  1:973012/239‑1 (MQ=255)

              gaTCCACGCCAAGGAGGGCCGGGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGGTTGTACCCCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTGGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGatcat                                                                                                                                                                                                                                <  2:1130425/238‑1 (MQ=255)

                 ccACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGAGGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCg                                                                                                                                                                                                                            <  1:992609/239‑1 (MQ=255)

                 caaagACATGGAAGGCCGCGTGGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCg                                                                                                                                                                                                                            <  2:562683/235‑1 (MQ=255)

                   tcgaggtAGAAGGCCGGGTTGATCCGACTCGGGGAAAGCGGTTAAAGGCAGAGGGGGATGGTTGGGTTGTCCCGCCCGGAGCCGGGCCGGGGGACTTCGGATTGATTGTGTTACGTAATCCCCTTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCGGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCCACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  2:1889975/233‑1 (MQ=255)

                   gggACATGGAGGGCGGGTTTGTTCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGTATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTTTGCTACGTAATCCCCCTTCTGGCATTACTCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  2:881574/237‑1 (MQ=255)

                   aCGACTTGGAAGGCCCCGTGGCTCCTCCACTGGGAAAGCGGGTAAAAGCGGATGGGGTTTGTTGGGTGATCCCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGGTACGTAATCCCCCTTCTGGCATTACCTCGTTGCCGTTATTTGGGGGAGATAAAGCCGCTTTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  2:756024/239‑1 (MQ=255)

                   aCGACATAGAAGGCCGCGTTGGTCCGACCCTGGGAAAGCGGTTAAAAGCAGGTGGGGATGGTTGGATTGTACCTCTCGCGGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTCCTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGAACAGTGCCTCCGGCAGGGTACCCTGac                                                                                                                                                                                                                                      >  2:1605247/1‑226 (MQ=255)

                   aCGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  1:253467/239‑1 (MQ=255)

                   aCGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  1:1110126/239‑1 (MQ=255)

                   aCGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  1:1436560/239‑1 (MQ=255)

                   aCGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  2:907669/239‑1 (MQ=255)

                     gACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAACGGGAGTCAGCCAGGCTACCCTGACGATCACCTCGATAc                                                                                                                                                                                                                        >  2:1845473/1‑239 (MQ=255)

                      aagtAGAAGGGCGGGTTGATCCGACGCTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCCCCCGCAGCCGCGCCGGGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAGGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGAt                                                                                                                                                                                                                          <  2:809891/233‑1 (MQ=255)

                         tAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGGGACCCAGGCAGGCTACCCGGAAGATAATCCCGATACGtag                                                                                                                                                                                                                    >  2:1327858/1‑237 (MQ=255)

                          aGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGt                                                                                                                                                                                                                   >  1:1325866/1‑239 (MQ=255)

                             aGGGCGCGTTTGTCCGCACCTGGAAAGGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCCGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                                                                                                                                                                                                                <  2:1600440/239‑1 (MQ=255)

                               gCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                                                                                                                                                                                                                >  1:881574/1‑237 (MQ=255)

                               gCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCCACCCTGAAGATCATCTCGATACGTTGTACa                                                                                                                                                                                                                <  1:566982/237‑1 (MQ=255)

                                 cgcgTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtca                                                                                                                                                                                                            >  2:1403245/1‑239 (MQ=255)

                                 cgcgTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtca                                                                                                                                                                                                            >  2:972203/1‑239 (MQ=255)

                                      tGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCCGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCCGAAGATCATCTCGATACGTTGTACAGCCATCaaa                                                                                                                                                                                                       >  1:1245406/1‑239 (MQ=255)

                                              cacTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGaa                                                                                                                                                                                               >  1:974859/1‑239 (MQ=255)

                                                   ggAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATa                                                                                                                                                                                                                         >  1:652360/1‑208 (MQ=255)

                                                   ggAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                                                                                                                                                                                                                >  1:461115/1‑217 (MQ=255)

                                                         cGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCCACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGActggctg                                                                                                                                                                                    >  1:1205038/1‑239 (MQ=255)

                                                             tAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGgcg                                                                                                                                                                                <  1:773787/239‑1 (MQ=255)

                                                                     gATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGCAAAGTGACCCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTc                                                                                                                                                                        >  2:1136816/1‑239 (MQ=255)

                                                                      aTGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCg                                                                                                                                                                       >  1:1150872/1‑239 (MQ=255)

                                                                        ggggATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGt                                                                                                                                                                     <  1:1671103/239‑1 (MQ=255)

                                                                        ggggATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGt                                                                                                                                                                     <  1:1836973/239‑1 (MQ=255)

                                                                        ggggATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTATTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGt                                                                                                                                                                     >  2:495634/1‑239 (MQ=255)

                                                                         ggttgtgTTGGATTGTCCCTCCCGACGCCGCGCCTTGGGATGTCGGTTTGATTTTGCTACGTAACCCCCCTTCGGCCATTACCCCTTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACCGGCTGGGCGCAAACGTCGGt                                                                                                                                                                     <  2:32392/232‑1 (MQ=255)

                                                                         cggATGGTTTGATTGTCCCCCCCGCACCCGCGCGGTGGGACTTCGGTTTGTTTGTGATACGTAACCCCCCTTCTGGCATTACGCCGTTGCCTTTATTTGGGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGt                                                                                                                                                                     <  2:101775/237‑1 (MQ=255)

                                                                            aTGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                                                                                                                                                                                                                >  1:331140/1‑192 (MQ=255)

                                                                            aTGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                                                                                                                                                                                                                <  2:331140/192‑1 (MQ=255)

                                                                             tGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATg                                                                                                                                                                 >  1:1590145/1‑238 (MQ=255)

                                                                                  gggttGTCCCTCCCGCCGCCGGGCCGTGGGACTTCGGATTGATTTGGCTACGTACTCCCCTTTCTGGCATTACGCCGTTGCGGTTATTTGGGGGGGAAAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACAt                                                                                                                                                           <  2:1446306/236‑1 (MQ=255)

                                                                                    aTTGTCCCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATc                                                                                                                                                          <  1:1661497/238‑1 (MQ=255)

                                                                                    aTTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCa                                                                                                                                                         >  1:1217785/1‑239 (MQ=255)

                                                                                        tACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACAt                                                                                                                                                           <  1:1049569/233‑1 (MQ=255)

                                                                                        tACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGGGAAGTGCCTGGCTGGGCGCAAACGTCGGAGATGTCCCATCAGTGc                                                                                                                                                     >  1:348334/1‑239 (MQ=255)

                                                                                           ctcccccgattgcgcCTTGGGTCTTCGGAGTGATTCTGCTGGCGAACTCCCCTTTTCGGTCTGCGGCGTTGCCGTTATTTGAGGGAAATAAAGCCGCGCTCACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTGACATCAGTGCGAt                                                                                                                                                  <  2:1159528/228‑1 (MQ=255)

                                                                                           cgcccGCGGCGGCGCCGTGGGACTTCGGTTTGATTGTGCTACGTAATCCCCCTTCTGGCATTCCGCCGTTGCCGTTATTGGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGAt                                                                                                                                                  <  2:187025/237‑1 (MQ=255)

                                                                                             cccGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCg                                                                                                                                                                                                                            >  1:344113/1‑163 (MQ=255)

                                                                                                  ggCCGGGCCGTGGGGCTTCGGATTGATTGTGCTGCGAAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAACGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGc                                                                                                                                           <  2:547081/238‑1 (MQ=255)

                                                                                                  aGCCGCGCGGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACg                                                                                                                                                                                                                       <  2:1547619/163‑1 (MQ=255)

                                                                                                  aGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACg                                                                                                                                                                                                                       >  1:1547619/1‑163 (MQ=255)

                                                                                                     cgcgCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCg                                                                                                                                         <  1:1608311/238‑1 (MQ=255)

                                                                                                       cgcCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGt                                                                                                                                                                                                              <  2:179486/167‑1 (MQ=255)

                                                                                                       cgcCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGt                                                                                                                                                                                                              >  1:179486/1‑167 (MQ=255)

                                                                                                         ccGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTg                                                                                                                                     >  1:1741378/1‑238 (MQ=255)

                                                                                                            tGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACa                                                                                                                                  >  1:1165185/1‑238 (MQ=255)

                                                                                                                 cTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCCACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGtt                                                                                                                            <  1:263476/239‑1 (MQ=255)

                                                                                                                    cGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCg                                                                                                                          >  1:1718285/1‑238 (MQ=255)

                                                                                                                    aggATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCAGTTGCCGTTATTTGAGGGAGATAAACCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGGGAATCAGGCAGGATACCCTGAAGATCATCTCGATACGTTGTACAGTCATTAAAACTGTGAAGTGACTGGCTGGGCGGAAACGGCGGGGATGGCACATCAGTGGTATACCTTGGCGGGTGAAAGCggctag                                                                                                                          >  2:973908/2‑234 (MQ=255)

                                                                                                                         tgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTg                                                                                                                    <  1:1142010/239‑1 (MQ=255)

                                                                                                                         tgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTg                                                                                                                    <  2:1154497/239‑1 (MQ=255)

                                                                                                                         tgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTg                                                                                                                    <  1:936212/239‑1 (MQ=255)

                                                                                                                         tgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTg                                                                                                                    <  1:502169/239‑1 (MQ=255)

                                                                                                                         gtatgtggatccGTAATCCCCCTTCTGGAATGACGCCGTTGCCGTTTTTGGAGGGAGAAAAAGCCCCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGAGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTg                                                                                                                    <  2:355519/228‑1 (MQ=255)

                                                                                                                          gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCCACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGGGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGa                                                                                                                   >  2:517823/1‑239 (MQ=255)

                                                                                                                           attgTGCTCGGTAATCCCCCGTCTGGCATTCCGCCGTTGCCGTTTTTGGAGGGGGAAAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGa                                                                                                                   <  2:1824373/238‑1 (MQ=255)

                                                                                                                                 cTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAGAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTg                                                                                                                                                      >  1:1860432/1‑197 (MQ=255)

                                                                                                                                    cGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTAACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATa                                                                                                          >  1:214905/1‑238 (MQ=255)

                                                                                                                                         tCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATg                                                                                                    >  1:966744/1‑239 (MQ=255)

                                                                                                                                              cttgtagaataaCGCCGTTGCCGTTCTTTGAGGGGGAGAAAGCCGCGCTTACCGCCGCATTAAAAGCGTCAGGTCGTAAAGTGACTCAGGCAGGCTTCTCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACGGGCTGGGCGCAAAATTCGGTGATGTCACATCAGTTCGATACCTTGCCGTGTGACAGCGGTTCGCGTCTGAGGTTGCATACCGATGACGGc                                                                                               <  2:1456521/228‑1 (MQ=255)

                                                                                                                                              cTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACgg                                                                                                >  2:1780516/1‑238 (MQ=255)

                                                                                                                                              cTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACg                                                                                                 >  1:793375/1‑237 (MQ=255)

                                                                                                                                              cTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACg                                                                                                 >  1:1809219/1‑237 (MQ=255)

                                                                                                                                                 cTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTgg                                                                                            <  1:1308516/239‑1 (MQ=255)

                                                                                                                                                   ggaaTTACGCCTTTGCCGTGATTTTGGGGAGTTAACACCCCGTTGACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCa                                                                                          <  2:1106916/236‑1 (MQ=255)

                                                                                                                                                   ggCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGCCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGGGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCa                                                                                          >  1:1273644/1‑239 (MQ=255)

                                                                                                                                                        tgtggcgttgccgttATTAGAGGGAGATAAAGCCGCGCTTTCCGCAGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaa                                                                                     <  2:1614860/234‑1 (MQ=255)

                                                                                                                                                        tACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaa                                                                                     <  1:1403245/239‑1 (MQ=255)

                                                                                                                                                         aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACg                                                                                                                                                                                                                       <  2:1794193/108‑1 (MQ=255)

                                                                                                                                                         aCGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACg                                                                                                                                                                                                                       >  1:1794193/1‑108 (MQ=255)

                                                                                                                                                               ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCACTTAATTGGGGt                                                                              >  1:571897/1‑239 (MQ=255)

                                                                                                                                                               ttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGt                                                                              >  1:1504503/1‑239 (MQ=255)

                                                                                                                                                                   ggtgattgaggggAGATAACGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACGCAGGCAGGCTACCCTGAAGATCATCGCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCa                                                                           <  2:231660/228‑1 (MQ=255)

                                                                                                                                                                   cgtttgttGGGGGGGATAAGGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCa                                                                           <  2:1718285/232‑1 (MQ=255)

                                                                                                                                                                            aGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACAt                                                                                                                                                           <  1:1062241/149‑1 (MQ=255)

                                                                                                                                                                            aGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACAt                                                                                                                                                           >  2:1062241/1‑149 (MQ=255)

                                                                                                                                                                            aGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGCTGACGGACGGCAATTAATTGGGGGGCAAAGtt                                                                      >  2:1218146/1‑234 (MQ=255)

                                                                                                                                                                                agaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTcac                                                                                                                                                             <  1:146005/143‑1 (MQ=255)

                                                                                                                                                                                      aGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATTATCTCGATAGGTTGTACAGTCATCAAAACTGTGAAGTAACTGGCTGGGCGCAAACGTCGGTGATATCACATCCGTGCGATACGTTGCCGGGGGACAGCGGTTCGCCTCTGAAGTGTCATATCGATGACGGCAGGCCATGAACTGGGGGGGAGAGTTCCGCTCCCGcactga                                                       >  2:228731/1‑234 (MQ=255)

                                                                                                                                                                                           cgcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGGCGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGa                                                   >  1:1755162/1‑238 (MQ=255)

                                                                                                                                                                                                   gccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGgcg                                          <  1:1811768/239‑1 (MQ=255)

                                                                                                                                                                                                        aTTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCg                                     >  1:92450/1‑239 (MQ=255)

                                                                                                                                                                                                                 ggCAGGTCGTAAAGTGACTCAGGCAGGCCACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCg                              <  1:1871102/237‑1 (MQ=255)

                                                                                                                                                                                                                       tCGTAAGGTGATCCAGGCGGGCTACCTTGAAGCTCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCCAACGTCGGTGATGTCCCATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCAtt                       <  2:1581691/238‑1 (MQ=255)

                                                                                                                                                                                                                                 aCTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTgg                                                                                            >  2:1430117/1‑159 (MQ=255)

                                                                                                                                                                                                                                 aCTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTgg                                                                                            <  1:1430117/159‑1 (MQ=255)

                                                                                                                                                                                                                                     aggcagtcAACCCTGAAGATCAGCTCGATCCGTTGTACAATCATCAAAACTGGGAAGGGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATCCCGATGACGGCTGTCAATTAATTGGGGTGCAAATTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGtttt        <  2:1165185/239‑1 (MQ=255)

                                                                                                                                                                                                                                           gcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCa                                                                                          <  1:1459213/151‑1 (MQ=255)

CTTTGGCGCTATGAGATCCACGACATAGAAGGCCGCGTTGATCCGACACTGGGAAAGCGGTTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGCCGCGAAAGATCGCTGGCGCGCCGATAACCGGGCCATTTCTGTTACCGGTTTTCGCGAC  >  NC_007779/1673995‑1674468

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: `‑`