BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,638,874	A→T	10.9%	L232F (TTA→TTT)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,874	0	A	T	10.9%	487.7 / 7.8	171	L232F (TTA→TTT)	creD	inner membrane protein
Reads supporting (aligned to +/- strand): ref base A (71/73); new base T (18/2); total (91/80)
Fisher's exact test for biased strand distribution p-value = 5.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTT  >  NC_000913/4638638‑4639111

ggtgtggtgAAAGCCGCTGAAGTTAACGGAACGGCGCTGACCAGTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAgg                                                                                                                                                                                                                                                <  2:321133/239‑1 (MQ=255)

  tgtggtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATTGCCCTGAATTTAAACGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGGGTAAATGCGAAATGACCTTAACCA‑GCCAC‑CGGCCCCATCCCAGTTTTTTaaggg                                                                                                                                                                                                                                              >  1:608365/1‑235 (MQ=255)

      gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTTAATATTGCCCTGAATTTTAGCGGGACCGGCGATCTTTTTTGGGTGCCTGGCCGGGGGAATAGCGGAAATACCCTAACCC‑GCAACTTGGCCGCCCCCCGGTTTTTGGGGGGAtt                                                                                                                                                                                                                                           >  2:1307159/1‑239 (MQ=255)

      gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTTTGGTGCCTGGCGGGCGCAATAGCGAAATGGCCCTCACCA‑GCAAC‑CGGCCCCCCCCCCGTGTTTTTGGGGAttt                                                                                                                                                                                                                                          >  2:162730/1‑239 (MQ=255)

      gtgAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAACTGAATATGGCCCTGAATTTAAGCGGTACCGGCGACTTTTCATGGGTGCCTTGCGGGCCGAAAAGCGAAAAGACCTTTACCC‑ACACA‑TGGCCGCCTCCACGTGTTTTTGGGGAttt                                                                                                                                                                                                                                          >  2:827594/1‑239 (MQ=255)

           aGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGGAGCAGAACCTGAAGCTGAAAATGGCCCTGAATTTAAGCGGTAACGGCGAACCTTTTGTGGTGCCCGGCCGGCTGAAAAGGGAAAAGACCTTTACCA‑CCCAC‑GGGCCGCCTGCCCGGTTTTTTGGGGGTTTTTTA‑c                                                                                                                                                                                                                                    >  2:716466/1‑239 (MQ=255)

           aGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTTTGGTGCCTGGCGGGCGGAATAAGGAAATGACCTTAACCA‑GCAAC‑ATGACGCATCCCGGGTTTTTAGGGGATTTTTTa                                                                                                                                                                                                                                      >  1:474937/1‑238 (MQ=255)

               ccTTAAGTTAACGGGACGGGGGCGACCATTGAACCCGGCCCCGGGTTAGAGCAAGGGGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGc                                                                                                                                                                                                                                <  2:685323/239‑1 (MQ=255)

                         acggaacggCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGGAGCTGAATATGGCCCTGAATTTAACCGGTACCGGGGATCTTTCTTTGGGGCCCTGCCGGCGTGAAAGCGCAATGACCTTTACCC‑GCAAA‑CGGCCGCCCCCCCGGTTTTTTGGGGGTTTTTTACCCAGCAAAACgggg                                                                                                                                                                                                                       >  2:1245914/1‑238 (MQ=255)

                                       aCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGGATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAACGAAATGAACTTTACCA‑ACAAC‑TGGCCGCCTCCCAGTTTTTTTGGGGATTTTTTA‑CCCACCCAAAGGG‑GA‑GTTAACGGGTGAgg                                                                                                                                                                                                      >  2:1277048/1‑239 (MQ=255)

                                       aCCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGGACCGGCGATCTTTCTGTGGTGCCTGGCGGGGGTAAAAACGAAAAGACCTTAACCA‑GCAAC‑TGGCCCCATCCCAGTTTTTTAGGTGATTTTTTC‑CCAGCCAAAAGGG‑GA‑GATAAGGAGACAgg                                                                                                                                                                                                      >  2:246305/1‑239 (MQ=255)

                                         cATTGACCCCGGCGCCGGGTTAGGGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGtt                                                                                                                                                                                                    <  2:1021051/239‑1 (MQ=255)

                                         cATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAg                                                                                                                                                                                                           <  1:1102213/232‑1 (MQ=255)

                                                  cggacctgggTTAGAGGAGGGTGGGCAGGGTGTGCATATCCCGTTTCCTGAAAGGGACTGGGGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGGGTGGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGgcgc                                                                                                                                                                                           <  2:339302/232‑1 (MQ=255)

                                                   ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGGAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCCGTTTTTTAGGGGATTTTCTA‑CCCGCCAAACGGG‑AA‑GGTAGCGAGTCAGGTTTTCAGgcgc                                                                                                                                                                                           >  1:924053/1‑238 (MQ=255)

                                                   ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTAACGGCGATCTTTCTGTGGGGCCTGGCGGGCGTAATAACGAAATGACCTTAACCC‑GCAAA‑CGGCCGCCTCCCAGTTTTTTTGGGGGTTTTTTA‑CCAGCCCAACGGG‑AA‑GTTAGCGAGTCAGGGTTTTAGgcgg                                                                                                                                                                                           >  2:1296522/1‑237 (MQ=255)

                                                   ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAAAATGGCCCTGAATTTAAACGGTACCGGCGATCTTTTTGTGGGGCCTGGCGGGGGGAATAACGCAACGACCTTAACCC‑GCACC‑TGGCCGCCTCCCCGCTTTTTTGGGGATTTTCTA‑CCCGCCAAACAGG‑AA‑GGTAGCggggaggagtttccatgtc                                                                                                                                                                                           >  2:1102213/1‑220 (MQ=255)

                                                   ggCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATACCTTTACCTGAAGGGGACTGGCGGGAGCCGAACCTGAAGCTGCATATGGCCCTGAATTTAAACGGTACCGGGGATCTTTCTTTGGTGCCCTGCGGGGCTAATAGCGAAATGGCCTTTACCAGGCAAC‑‑GGGCGCATCCCCGGTTTTTTGGTGAGTTTTTC‑CCCACCAAACGGGCAG‑GTTACGGAGGCAGGTTTTAGGgcg                                                                                                                                                                                            >  2:892659/1‑238 (MQ=255)

                                                   cgccccGGGTTGGAGGAAGGCGGGCAGGGCGGGCATAGCCATTTACCCGAAGGGGGCTGGGGGAAGCAAAACCTGAAGCTGAATATGGCCCTGAATTTAAGCTGTACCGGCGATCTTTCTGTGGTGGCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAA‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCa                                                                                                                                                                                          <  2:1189699/235‑1 (MQ=255)

                                                              agggCAAGGCGGGCAGGGCGGGAATATCCCTTTGCCGGAAGGGGACTGGCGGAAGCAGAAACTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAagc                                                                                                                                                                                <  1:1289309/235‑1 (MQ=255)

                                                                     ggCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGGAGCCGAACCTGAAGCTGAATATGGCCCTTGATTTTAGCGGTACCGGCGATCTTTCTGTGGTGCCCGGCGGGAGTAATAGCGAAATGAACTTAACCA‑GCAAC‑TGGCCGCCACCCCATTTTTTCTGTGATTTTCTC‑AACTCCAAACGGG‑AA‑GCAAAAGAGGCAAGTTTTTCGGCGCACCGGGAAaacaacagggg                                                                                                                                                                        >  2:621672/1‑232 (MQ=255)

                                                                     ggCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTTTGTGGTGCCTGGCGGGCGGAATAACGAAAGAGCCTTAACCA‑GCAAC‑TTGCCCCCTCCCAGGTTTTTTGGTGATTTTTTA‑CCAACCAAAAGGG‑GA‑GGTAGCGAGTCGGGGTTTTTTGGGCACTGGCAAAACAACTgggt                                                                                                                                                                        >  2:1094443/1‑237 (MQ=255)

                                                                     ggCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTTGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCCAC‑TGGCCGCATCCCCGTTTTTTTGGCGGTTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGGTTTCCGGGGCACTGGCCAAACAGCGGGtt                                                                                                                                                                        >  1:638111/1‑239 (MQ=255)

                                                                       cGGGCAGGGCGTGCATAGTCCCTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGtt                                                                                                                                                                        <  1:8206/237‑1 (MQ=255)

                                                                          gCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAGCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCt                                                                                                                                                                    <  1:1039746/238‑1 (MQ=255)

                                                                          gCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGGAAACGCGAAATGACCTTAAACA‑GCAAC‑CTGCCCCATGCCAGTTTTTTAGGTGATTTTCTT‑CCCCCCACAAAGG‑TA‑GGTACCCAGTACGGTTTTCCGGCGCACCGGCAAaccaaacggttttgtt                                                                                                                                                                   >  2:91741/1‑227 (MQ=255)

                                                                            aGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCtaa                                                                                                                                                                  <  1:207921/238‑1 (MQ=255)

                                                                             gggCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑CCAAC‑TGGCCGCATCCCCGTTTTTTTGGGGATTTTCTA‑CCCGCCAAACGGG‑AAGGTTAGCG‑GGCAGGTTTTTAGGGGCACCGGCAAAACAGCAGGTTTgtgcaaa                                                                                                                                                                >  2:203813/1‑233 (MQ=255)

                                                                                      tatCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTg                                                                                                                                                       <  1:300029/239‑1 (MQ=255)

                                                                                        tCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGGCTCCCAGTTTTTTAGGGGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTACGGCGCCCTGGCAAAGCCGCCGGTTTGCGAATAATCTCGGTGAg                                                                                                                                                     >  1:155269/1‑239 (MQ=255)

                                                                                         cccTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGc                                                                                                                                                    >  1:128618/1‑239 (MQ=255)

                                                                                             ttACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTTTGTGGTGCCCGGCGGGCGTAATAGCGAAATGACCTTAACCA‑TCAAC‑TGGCCGCATCCCCGTTTTTTAGGGGCTTTTCTA‑ACCGCCCAAAGGG‑AA‑GATAAGCAGGCAGGGTTTCTGGGGACCCGGCCAAACAGAGGGTTTGCAAAAAATTTCGGGGAGGGttt                                                                                                                                                >  2:1085648/1‑239 (MQ=255)

                                                                                             ttACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑GA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAACAATCTCGGTGAGCGttt                                                                                                                                                >  1:640434/1‑239 (MQ=255)

                                                                                              tACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGtttt                                                                                                                                               >  1:1009372/1‑239 (MQ=255)

                                                                                              tACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGGGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTTGCAAACGAACTGGTTTGCTAATAATCTCGGTGAGCGtttt                                                                                                                                               >  1:1250308/1‑239 (MQ=255)

                                                                                                  ggAAGGGGAGGGGCGGAAAGAGAACTGGAAGCTGAAATGGGGCCGTAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGGGGGCGTATGAGCGAAATGACCTTAACCC‑GCAAG‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCTCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCtt                                                                                                                                           <  2:783337/238‑1 (MQ=255)

                                                                                                  ggAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCAGTGAGCGTTTTGCtt                                                                                                                                           <  1:376461/238‑1 (MQ=255)

                                                                                                   gAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TAGCCGCATCCCCGTTTTTTTGGTGATTTTCTA‑CCAGCCAAAAGGG‑CA‑GTTACCGAGTCAGGGTTTCAGGCGCCCGGGCCAAACAGCTGGTTTGCTTATAATATCGGTGACGGTTTTGTTTc                                                                                                                                          >  2:1056974/1‑239 (MQ=255)

                                                                                                    aaGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATTGCCCTGAATTTAAGCGGTACCGGAGATCTTTCTGTGGTGCCTGGCGGGGGTAATAGCGAAATGACCTTAAACA‑GCAAC‑CGGCCGCATCCCCATTTCTTAGGGGATTTTCTA‑CCCGCCCAACGGG‑CA‑GTTATCGGGTCAGGGTTTCCGGCGCACCGGCAAAACACCCGGGTTGCCACAAATTCCGGGGGGGGTTTttcttcc                                                                                                                                         >  2:994559/1‑238 (MQ=255)

                                                                                                    aaGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCa                                                                                                                                         >  1:770596/1‑239 (MQ=255)

                                                                                                    aaGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGGAATAGCGAAATGACCTTAACCA‑GCAAA‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTTTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCCGGTTTTCAGGGCCACCGGCAAAGGAGCCGGTTTTCTAAAAATTCCGGTGGGGGgtttgtttaa                                                                                                                                         >  2:1031020/1‑237 (MQ=255)

                                                                                                     aGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCa                                                                                                                                         <  1:221098/238‑1 (MQ=255)

                                                                                                          aCTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGtt                                                                                                                                                                        >  1:1015860/1‑202 (MQ=255)

                                                                                                              gCGGAAGCAGCACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGCTCTTTCTGTGGTGCCTCCCGGGCGTAATAGCGCAATGACCTTAACCA‑GCAAC‑CGGCGGCATCCCAGTTTTTTTGGTGTTTTTCTC‑CCAGCCAAACGGG‑CA‑GTTAGCGAGTCAGGTTTTCACGCGCACCGGCAAAGCAGCTGGCTTGCTTATAATCTCGTCGACCGCTTCGCTTCatgacatcgt                                                                                                                                >  2:307516/1‑229 (MQ=255)

                                                                                                                        aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGT‑AA‑GTTAGCGAGTCCGGTTTTCAGGCGCACTGGCAAACCAGCGGGTTTGCTAATAATCTCCGGGAGGCGTTTGCTTCTGGCAATGATACCGGCTGGGa                                                                                                                     >  2:410449/1‑239 (MQ=255)

                                                                                                                        aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCCCTGGCAAAGCAGCTGGTTTGCCAATAAGCCCGGTGGGCGTTTTGGTTCAGGAAAGGAAACCGCTTgggg                                                                                                                     >  2:140252/1‑238 (MQ=255)

                                                                                                                        aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGACCGTTTTGCTTCAGGCAATGATACCGGCTGGGa                                                                                                                     >  1:610348/1‑239 (MQ=255)

                                                                                                                        aaCCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTAGCTAATAATCTCGGTGAGCGCTTTGCTTCAGGCAATGAGAACGGCTGGGa                                                                                                                     >  1:149573/1‑239 (MQ=255)

                                                                                                                                  tGAATATGGCCCTGAATTTAAGCGGTACCGGGGATCTTTCTGTGGTTCCTGGCGGGCGTAATCGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCAACCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCg                                                                                                           <  1:1366572/239‑1 (MQ=255)

                                                                                                                                  tGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTcc                                                                                                             <  1:818321/237‑1 (MQ=255)

                                                                                                                                   taaTAGGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCGGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCg                                                                                                           <  1:140252/237‑1 (MQ=255)

                                                                                                                                     atatGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCCAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCCATGGCCAAGCAGCTGGTTTGCTAATAATCTCGGTGGGCGTTTTGCTTCAGGAAATGATACCGGCTGGGAAAAATTccacgcg                                                                                                        >  2:281909/1‑234 (MQ=255)

                                                                                                                                     atatGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCCCTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGGGGGCGTTTTGTTTCAGGCAAGGATACCGGCTGGGAAAAATTTCCgggg                                                                                                        >  2:1109974/1‑237 (MQ=255)

                                                                                                                                           cccTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGCCGGGGAAAACTTTCCGGCGTTTAGc                                                                                                  >  1:298200/1‑239 (MQ=255)

                                                                                                                                             cTGAATTTAAGCGGGACCGGCGATCGTGCTGTGGGGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAAGCACGGTGAGCGTTTTGCTTCAGGCAAGGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGt                                                                                                <  1:220095/239‑1 (MQ=255)

                                                                                                                                              tGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCgg                                                                                                >  1:868119/1‑237 (MQ=255)

                                                                                                                                              tGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGt                                                                                                >  1:1189664/1‑238 (MQ=255)

                                                                                                                                              tGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGAGt                                                                                                >  1:26546/1‑238 (MQ=255)

                                                                                                                                              tGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGCCAATGATACAGGCTGGGAAAACTTCCCGGCGTTTAGAGTc                                                                                               >  1:189022/1‑239 (MQ=255)

                                                                                                                                              tGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GGTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTTCCGGCGTTTGGCGTc                                                                                               >  1:689074/1‑239 (MQ=255)

                                                                                                                                                aaTTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAACCGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCCGGCAAGAGTACCGGCTGGGGAAAATTCCCCGCGGTTAGCGTCga                                                                                             >  2:386445/1‑238 (MQ=255)

                                                                                                                                                          gTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGGGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTTGGCGTTTTGTTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGGGttttgtgtgcccgtaacaac                                                                                    >  2:282655/1‑221 (MQ=255)

                                                                                                                                                            aCCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTTCCGGCGTTTAGCGTCGCAGTAACGACGcc                                                                                 >  1:977557/1‑239 (MQ=255)

                                                                                                                                                              cGGCGTTTTTTCTGGGGTGCCTGGGGGGCGAAATAGCAAAAGGACCTTACCCA‑GCAAC‑GGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTGGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAg                                                                               <  2:1009372/239‑1 (MQ=255)

                                                                                                                                                               ggCGATCTTTCTGGGGTGCCTGGCGGGCGTAATAGCGAAATGCCCTTAACCA‑GCAAC‑TGGCCGGATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCTCAGTAACGACGCCAGc                                                                              <  1:985276/239‑1 (MQ=255)

                                                                                                                                                                  gtttgttCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GAAAA‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCGCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCg                                                                            <  2:640434/233‑1 (MQ=255)

                                                                                                                                                                      tttCTGTTGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCaa                                                                       <  1:569641/239‑1 (MQ=255)

                                                                                                                                                                      tttCTGGGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCaa                                                                       <  1:1132607/239‑1 (MQ=255)

                                                                                                                                                                         cTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGCTCAATAc                                                                    >  1:1218491/1‑239 (MQ=255)

                                                                                                                                                                          tgtgGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATAcc                                                                   <  1:759394/239‑1 (MQ=255)

                                                                                                                                                                             ggtgCCTGGCGGGCGTAATAGCGAGATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAAt                                                                >  1:515504/1‑239 (MQ=255)

                                                                                                                                                                              gtgCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGGAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAAtt                                                               >  1:1348170/1‑239 (MQ=255)

                                                                                                                                                                               tgCCTGGCGGGCGTAATAGCGAAATGGCCTTAACCA‑GCAAC‑TGGCCGGATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAAtt                                                               <  1:399170/238‑1 (MQ=255)

                                                                                                                                                                               tgCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGTTCAATACCAATTa                                                              >  1:583706/1‑239 (MQ=255)

                                                                                                                                                                                gCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCaaa                                                                                                                                                                                                                            <  1:510151/82‑1 (MQ=255)

                                                                                                                                                                                gCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCaaa                                                                                                                                                                                                                            >  2:510151/1‑82 (MQ=255)

                                                                                                                                                                                gCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTTGGAAAACTTCCCGGCGTTTAGCGTCGCCGTAACGACGCCAGCCGCTCCATACCAATTa                                                              >  1:455643/1‑238 (MQ=255)

                                                                                                                                                                                 ccTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGCCGTTTAGCGTGCCCGTAAAGACGCCAGCCGATCAATACCAATTa                                                              >  1:188969/1‑237 (MQ=255)

                                                                                                                                                                                    ggcgggcgTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATCATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCCGCGTTTAGGGGCGCAGTAACGCCCCCAGcacgtaaaaaaaaaaaaaataa                                                         >  2:213086/1‑218 (MQ=255)

                                                                                                                                                                                    ggcgggcgTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTTGCGTCGCAGGAAAGACGCCAGCCGAACAAAAACaaaaaaaaga                                                         >  2:1211238/1‑231 (MQ=255)

                                                                                                                                                                                       gggggTAATTGGGAAATGCCCTTAACCA‑GCAAT‑TGGCCGGATCCCGGTTTTTTGGGTGATTTTTTA‑CCAGCCAAACGGG‑AA‑TTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCGGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGAcc                                                       <  2:298200/238‑1 (MQ=255)

                                                                                                                                                                                              aTAGGGAAATGACCTTACTCA‑TGATC‑TGGCCGCTTCCCATGTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTGGCGAGTCAGGTTTTCAGGCGCACTGGAAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACt                                               <  2:608365/239‑1 (MQ=255)

                                                                                                                                                                                                   gAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTa                                          <  1:1307159/239‑1 (MQ=255)

                                                                                                                                                                                                    aaagccccTTAGCCT‑GCCAC‑GGGCCCGACCCCAGTTTCTTGTGTGATTCTCAA‑CCAGCCAAACGGG‑AA‑GTTAGGTGGTTAGGTTTTCAGGCACACTGGCAAAGCAGCGGGTTTGCTAATAATCCCGGTGAGCGTTGTGCTTCAGGCAATGATACCGGCTGGGAAAACCTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTa                                          <  2:516174/232‑1 (MQ=255)

                                                                                                                                                                                                     aattgccTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACg                                        <  1:707991/234‑1 (MQ=255)

                                                                                                                                                                                                     aaTGCCCTTAACCA‑GAAAA‑TGGCCGCAACCCAGTTTTTTGGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTGTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACg                                        <  1:1239981/239‑1 (MQ=255)

                                                                                                                                                                                                        gcccTTAACCA‑GAATG‑TGGCCCGTTCCAAGTTTTTTGGGGTATTTGTTA‑CCAGCCAAACGGG‑AA‑TTAAGCGGGTCAGGTTTTCGGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATCCCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCa                                     <  2:474937/237‑1 (MQ=255)

                                                                                                                                                                                                          ccTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACg                                        <  1:1067014/234‑1 (MQ=255)

                                                                                                                                                                                                           cTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTTGCAGTAACGATGCCCGCCGATCACTACCAATTAACTGACCGGGCGGCCAAGTCCGCCac                                    >  2:1209189/1‑236 (MQ=255)

                                                                                                                                                                                                             tAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTctg                                <  1:1109974/239‑1 (MQ=255)

                                                                                                                                                                                                                     cAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTGTCAGTAACGACGCCCGCCGATCAATACCAATTAAATGACCGGGGGACTAAAAACCCCCTTctctttgtg                          >  2:1067014/1‑231 (MQ=255)

                                                                                                                                                                                                                      aaC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAAAGGCGCCAGCCGGTCAATACCAATTAAATGACCGGGGGCCTAAGAACCCCaattcgcggtttgc                        >  2:1111113/1‑226 (MQ=255)

                                                                                                                                                                                                                       aC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCa                       <  1:57586/239‑1 (MQ=255)

                                                                                                                                                                                                                       aC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCa                       <  1:1060750/239‑1 (MQ=255)

                                                                                                                                                                                                                       aC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCCA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCa                       <  1:1056974/239‑1 (MQ=255)

                                                                                                                                                                                                                       aC‑TGGCCGAATCCCAGTTTTTTAGGTGATTTTCGA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGAAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCa                       <  1:1094443/239‑1 (MQ=255)

                                                                                                                                                                                                                        tttGGCCCCCTCCTCGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGTGTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCa                       <  2:1113466/237‑1 (MQ=255)

                                                                                                                                                                                                                          tGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACt                     <  1:246305/239‑1 (MQ=255)

                                                                                                                                                                                                                          gggCCGCTTTCCTTTTTTTTTGGTGTTTTTCTA‑CCAGCCAAAGGGG‑AA‑GTTAGCGATTCAGGTTTTCAGGCGCCCTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGGGCGTTTTGCTTCAGGAAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACt                     <  2:588655/238‑1 (MQ=255)

                                                                                                                                                                                                                                             tAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCCGTAACGACGCCAGCCGAGCAATACCCATTAACTGACCGGGCGGCCAAGGACCCCATCCTGCCGATTGCCCTGACTTTTTTGGGGTTCttt  >  2:953176/1‑239 (MQ=255)

GGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCTGAATATGGCCCTGAATTTAAGCGGTACCGGCGATCTTTCTGTGGTGCCTGGCGGGCGTAATAGCGAAATGACCTTAACCA‑GCAAC‑TGGCCGCATCCCAGTTTTTTAGGTGATTTTCTA‑CCAGCCAAACGGG‑AA‑GTTAGCGAGTCAGGTTTTCAGGCGCACTGGCAAAGCAGCTGGTTTGCTAATAATCTCGGTGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATTGCACTGACTTTTATGGCGTTCTTT  >  NC_000913/4638638‑4639111

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