BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	388,772	T→C	10.2%	E319G (GAG→GGG)	hemB ←	5‑aminolevulinate dehydratase (porphobilinogen synthase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	388,772	0	T	C	10.2%	311.0 / 2.5	112	E319G (GAG→GGG)	hemB	5‑aminolevulinate dehydratase (porphobilinogen synthase)
Reads supporting (aligned to +/- strand): ref base T (71/26); new base C (2/13); total (73/39)
Fisher's exact test for biased strand distribution p-value = 1.43e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

ACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGGCAGTCT  >  NC_000913/388558‑389005

ccggtGCCAAGGTCGATTTCCGGCCACAGCGGCGGATAATGCATCGGGGCACTATCCTTGGGGATAACCCGTTTTATCGGGCGGG‑GTAATCCGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGc                                                                                                                                                                                                                     <  1:381820/235‑1 (MQ=255)

    tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑CAATCCAGCGCAAAGTAGCTgta                                                                                                                                                                                                                 >  2:165470/1‑237 (MQ=255)

    tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGaa                                                                                                                                                                                                                 >  2:657058/1‑239 (MQ=255)

    tGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCAAGCGCAAAAAATCTGaa                                                                                                                                                                                                                 >  1:897280/1‑239 (MQ=255)

        aaTGTCGATTGCCTGCCAAAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTAGCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATc                                                                                                                                                                                                             <  2:451392/239‑1 (MQ=255)

         aTGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGGAAC‑AAAATTAACGCAG‑AATCTTCT‑TTTCAGCC‑AAATCCACCGCAAAGTAGCTGAAAAACa                                                                                                                                                                                                            >  2:1262861/1‑239 (MQ=255)

         aTGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACTCAG‑AATCTCCT‑TCTCAGCC‑AAATCCCACGAAAAGTAGCTGAAAATCa                                                                                                                                                                                                            >  2:464606/1‑239 (MQ=255)

         aTGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAAtaa                                                                                                                                                                                                            >  2:107698/1‑237 (MQ=255)

         aTGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCa                                                                                                                                                                                                            >  2:1157810/1‑239 (MQ=255)

         aTGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGAGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGAAGCTGAAAAtaa                                                                                                                                                                                                            >  2:270372/1‑237 (MQ=255)

            ccGTTTACCTGCCACAGCGGCGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTg                                                                                                                                                                                                                   <  2:548149/229‑1 (MQ=255)

                     ggCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTAGCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATccgcacccg                                                                                                                                                                                                <  2:1136904/238‑1 (MQ=255)

                      gCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGCGAG‑GGTTAGGGTGAGGGGGCGGTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAACC‑CAAATTAACCCAG‑CATCTTCT‑TCTCACCA‑CAATCCAGCGCACAGTAGCTGAAAATCAGGTccccccccg                                                                                                                                                                                                >  1:895647/1‑238 (MQ=255)

                      gCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCACCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTTCGAGAGGAAAC‑AAAGTTAAAGCAA‑AATCTTCT‑TCTCAGCA‑CAACCAACCGCCAACAAGCTGaaaaaaagaacagctctctg                                                                                                                                                                                               >  1:548149/1‑223 (MQ=255)

                                cGGATAATGGAGGCGGCGGCGATTCTTGGGGAAAAACCGGTTTATCGGGCGGG‑GTAATACGGGAGGGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGAGGCGTTCACCGTACTTTCAACAGG‑TTAACTACCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAAt                                                                                                                                                                                     <  1:260748/238‑1 (MQ=255)

                                cGGATAAGGCAGCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GGAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAAt                                                                                                                                                                                     <  2:478667/238‑1 (MQ=255)

                                      gtGCAGCGGGGGACGATTCGTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGGGAGTGGACGGTCCCCTCGCCCCTGTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCt                                                                                                                                                                              <  1:1092526/238‑1 (MQ=255)

                                             gggCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCtt                                                                                                                                                                        >  2:827681/1‑238 (MQ=255)

                                                    ggtCCTTGGGGATAAACCGTTTTATCGGGGGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCGTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCa                                                                                                                                                                 <  2:711514/236‑1 (MQ=255)

                                                    gaTCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGAAACTAAAAAACACATCCCCACCCCC‑ACGCTTAACCGAACCTAAGCTATTGAGCaa                                                                                                                                                                >  1:820246/1‑238 (MQ=255)

                                                             gcttAAACCGTGTGATCGGGTGGGTGTATTGCGAG‑GTGGACGTGCCCCTCCCCCCCTTGGGGTGAG‑GGTTAGGGGTAGGGGGCCTTCACCGTACTTTCAAAAAG‑TTAACTCCCCCTGTCTGAGAGGTAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGCCTTtcgc                                                                                                                                                       <  1:256324/236‑3 (MQ=255)

                                                               aTAACCCTTTTTATTGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGATAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTctct                                                                                                                                                      <  1:646931/238‑1 (MQ=255)

                                                                tAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGAAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGGGCACGCCTTTCTCtt                                                                                                                                                     >  2:1238942/1‑238 (MQ=255)

                                                                     cGTTTTGTCGGGGGGG‑GTAATACGGGAGTGGACGGTCCCCTCGCCCCCTTGGGGGGAG‑GGTTAGGGGGAGGGGGCCTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCt‑a                                                                                                                                              <  1:1218852/239‑1 (MQ=255)

                                                                     cGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCt‑a                                                                                                                                              <  2:37065/239‑1 (MQ=255)

                                                                       ttttATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCt‑ata                                                                                                                                            >  2:164046/1‑239 (MQ=255)

                                                                        tttATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAACCTTTCGAGCACGACTTTCTCTTCATCt‑atac                                                                                                                                           >  2:1317498/1‑238 (MQ=255)

                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGCTCCGCACCCGC‑ACGCTTCATCGAACCTAAACTTTCGAGCACGACTTTCTCTTCATCA‑ATAGCACCCCCCAGCGCGGAGAACta                                                                                                                     >  1:1304327/1‑238 (MQ=255)

                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACtt                                                                                                                     >  2:1266004/1‑239 (MQ=255)

                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCCCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACtt                                                                                                                     >  2:587222/1‑239 (MQ=255)

                                                                                               gagTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTCCTCTTCTCCT‑ATAGAACCCCCCAGCGCGCCGaaata                                                                                                                     >  1:727711/1‑236 (MQ=255)

                                                                                                  ggggcGGCCCCCCCCCCCCTGTGGGGAGGG‑TGGTAGGGTGAGGGGGCGTTCCCCGTTTTTTAAACAGG‑TTAACTCCCCCTTTCAGAGAGGAAAA‑AAAATTACCGCAG‑AACCTTCT‑CCCCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGATTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTaa                                                                                                                   <  1:1262861/234‑1 (MQ=255)

                                                                                                    gACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTTTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCCAAGTAGCTGCAAATCAGATCCGCCCCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCCCGACTTCTCCTTCATCT‑ATAACACACCACACCGCCGCGAACtataaca                                                                                                                >  1:535594/1‑233 (MQ=255)

                                                                                                    gACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCCCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAatc                                                                                                                 >  2:274948/1‑238 (MQ=255)

                                                                                                    gACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑GTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCCGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCCC‑ACGCTTAATCGAACCTAAGCTTTCGAGAACGACTTTCTCTTCACCT‑AGAGCACCCGCCAGCGCGGCCAAATTAata                                                                                                                 >  1:1224649/1‑237 (MQ=255)

                                                                                                        gTCCCCCCCCCCCAGGGGGGAGAG‑GTTGAGGGTAGGGGGACCTCCACAGTAC‑TTCAAAAGGTTTAACTCCCCC‑TTCTGGAAAGAAACAAAAAGTAAC‑CAGAAACCTTCTCTCCCA‑CA‑AAAACCAGGCAAAATTAGC‑GAAAAACAAGTCCCCCCCCCCAACGCTTACTCGAACCAAGGCTTTCGAGCACGCTTTTCTCTTCATCA‑ATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGc                                                                                                            <  1:1281041/239‑1 (MQ=255)

                                                                                                         tCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTc                                                                                                                                                    <  2:197313/199‑1 (MQ=255)

                                                                                                         tCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGCCGAACTTAATCATCGCa                                                                                                           >  2:1093252/1‑239 (MQ=255)

                                                                                                            ccTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATa                                                                                                         >  2:929448/1‑238 (MQ=255)

                                                                                                            ccTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCAACT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATa                                                                                                         >  2:716796/1‑238 (MQ=255)

                                                                                                                   cccTTGGGGGGGGTGGTTAGGGTGAGGGGGCTTACCCCTTCTTTCCAAACGG‑TTAACTCCCCCTTTCTGAAAGGAAAC‑AAAATTAACGCAG‑AACCTGCT‑TCCCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCCC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCg                                                                                                  <  1:164046/239‑1 (MQ=255)

                                                                                                                   cccTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AACATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCCGATCCGCACCCGC‑ACGCGTAATCGACCCTAAGCTTTCGAGCACGACTTTCTCCTCATCT‑ATAGCa                                                                                                                                         >  1:197313/1‑199 (MQ=255)

                                                                                                                    ccTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACCCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCt                                                                                                <  2:320738/239‑1 (MQ=255)

                                                                                                                       tGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAAt                                                                                                                                                                                     >  1:1174396/1‑152 (MQ=255)

                                                                                                                       tGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAAt                                                                                                                                                                                     <  2:1174396/152‑1 (MQ=255)

                                                                                                                              ag‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACCTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTAAATCATCGCATACTcccccccaaccaaaaac                                                                                      >  1:1017532/1‑226 (MQ=255)

                                                                                                                               g‑ggTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCCTACTCACCGCTCACCTGATAc                                                                                      >  2:284751/1‑238 (MQ=255)

                                                                                                                                    tAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGAACGACTTTCTCTTC‑TCTCATAGCCCCCCCCAGCGCGGCGCACTTAATCATCCCAAACCCCCCGCCCACCAGAAACGCGcc                                                                                 >  1:710688/1‑239 (MQ=255)

                                                                                                                                      gggTGAGTGGTCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCCCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTCCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCCCCAACCCGGCGAACTTAACAATCGCCAACTCACCGCCAACCTGATACGCCCCag                                                                               >  1:133616/1‑238 (MQ=255)

                                                                                                                                      gggTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGACCACGACTTTCTCTTCATCC‑ATAGACCCCCCCAACGCGGCGAACTAAATCATCCCAAACCCACCgtccccccaaacacgccca                                                                                >  1:702141/1‑220 (MQ=255)

                                                                                                                                          gcgggggCCTTCCCCTTTCTTTCAATATG‑TTACCTCCCCCTTTTTAAGAGGAAAC‑AAAATTAACACAG‑AATCTCCT‑CCCCAGCC‑AAACCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAACCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCgg                                                                           <  1:587222/237‑1 (MQ=255)

                                                                                                                                           cgggcgCTTTCACCGTACTTTCAACAGG‑TTAGCTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGAAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGAACGCCTTTCTCTTCATCT‑ATAGCACCCCCCAGACCGGCGCATTTAACCATGGCCAACTTCCCGCTCACCCGAAAAGCGCCACTCGGc                                                                          >  1:170825/2‑239 (MQ=255)

                                                                                                                                           aGGGGGCGTTCCCCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCt                                                                                                <  2:307147/217‑1 (MQ=255)

                                                                                                                                           aGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGCAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCACAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACAACTTTCTCTTCATCT‑ATATCACCCGCACGCGCTGCGCACTTAATCATCGCATACCCACCGCTCACCTCATACCCCCCAACCGGc                                                                          >  1:7336/1‑239 (MQ=255)

                                                                                                                                           aGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGCGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTACTCATCCCATACTAACCGCTCACCTGATACGCGCCAACCGGc                                                                          >  1:915026/1‑239 (MQ=255)

                                                                                                                                           aGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGACCTTAATAATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGc                                                                          >  2:888838/1‑239 (MQ=255)

                                                                                                                                           aGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCCCCAATCGGc                                                                          >  2:915241/1‑239 (MQ=255)

                                                                                                                                           aGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCCCGGCGAACTTAATCATCCCCTACCCCCCGCCCACCTGaacacccccaacagcc                                                                          >  1:513907/1‑224 (MQ=255)

                                                                                                                                           aGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCACCT‑ATAGCACCCGCAAGCGGGGCGACCTTAATCCTCGCATACTCCCCGCTCACCTCAAACCCGCCAAACgcc                                                                          >  1:613842/1‑237 (MQ=255)

                                                                                                                                           aGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGAACCCGA‑ACGCTTAATCGAACCTAAGCTTTTGAGCCCGACTTTCTCTTCATCT‑ATACCCCCCCCCACCCCGGCGCAACTAAACATCACATAATCCCCACTCCCCTCACCCGCCCCAAACgac                                                                          >  1:668852/1‑237 (MQ=255)

                                                                                                                                            gggggCTTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCCC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGc                                                                          <  2:76448/238‑1 (MQ=255)

                                                                                                                                            gggggCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCa                                                                         >  2:260344/1‑239 (MQ=255)

                                                                                                                                             ggggCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAGAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCCCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGAACCGCGCCAACCGGCaa                                                                        >  1:1226711/1‑239 (MQ=255)

                                                                                                                                                  tttCCCCGTTTTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCA‑AAATCCAGCGCAAAGTAGCTGAAATCAAGATCCGCCCCGCC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCCCTCCCTCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAg                                                                   <  1:657058/238‑1 (MQ=255)

                                                                                                                                                   ttCACCGTACTTTAAACCGG‑TTACCCCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAA‑AATCTTCT‑TCCCAGCC‑AAATCCAGCGCAAAGGAGCTGAAAATCAGATCCGCCCCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAg                                                                   <  2:1226711/238‑1 (MQ=255)

                                                                                                                                                     cACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGCACTTAATCCTCGCATACTCACCGCCCACCTGAAAcgcgttcactttaaaaaatctga                                                                >  1:307147/1‑220 (MQ=255)

                                                                                                                                                          tACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAGTCCAGCTCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGCACCTAAGCTTTCAAGCAAGACTTTCTCTTCATCT‑ATAGCACCCCCCAGCGCGCCGACCTTAATCATCGCATACTCACCCCTCCCCTGATCACCGCCAAACGCCAAATCAGTACgccaa                                                           >  1:475166/1‑235 (MQ=255)

                                                                                                                                                           aCTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCCCCAGCGCGGCGAACTTAACCATCGCATACTCACCCCCCACCTGATACCCGCCacacgcacatccagaccggtca                                                           >  1:1298962/1‑217 (MQ=255)

                                                                                                                                                            cTTTCAACAGT‑TTAACTCCCCCTTTCTGAAAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACg                                                         <  2:1224649/239‑1 (MQ=255)

                                                                                                                                                            cTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACg                                                         <  2:897280/239‑1 (MQ=255)

                                                                                                                                                            atttCAACCGG‑TTAACCCCCCCTTTCTGAAAGGAAAC‑AAAATTAACGAAA‑AACCTTCT‑TCCCAGCC‑AAATCCAGCCGAAAGTAGCTGAAAATCAGACCCGACCCCGC‑ACGCTTAATCGAACCTACGCTTTCGAGCACGACTTTCTCTTAATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACg                                                         <  2:1173161/238‑1 (MQ=255)

                                                                                                                                                               tCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCa                                                       >  2:934289/1‑238 (MQ=255)

                                                                                                                                                                  aCAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCCACCCCGATACGCGCAAACCGGCAATCGCGTACGTTCACGCAACTc                                                   >  1:287841/1‑239 (MQ=255)

                                                                                                                                                                  aCAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGACCGGCGAACTTAATCATCCCCTACTCACCGACTACCTCAAACGCGCCAATCCGAAATTCAGAACGTTCACCCAGccc                                                   >  1:406936/1‑237 (MQ=255)

                                                                                                                                                                  aCAGG‑TTAACTCCCCCTTTCTCAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCAGCAAATCATATCTGCACCCGC‑ACTCTTAATCGAACCTAAACTTTCACGCAGGACTTTCTCTTCATCT‑ATACCCCCCCCCAACTCCACGACCTAAACCACCGCCTACTCACCCCCCACCCAATACGCGCCAACCGGCAATTCAAACCGTACACGccaca                                                    >  1:572773/1‑234 (MQ=255)

                                                                                                                                                                     gg‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacg                                                >  2:65221/1‑239 (MQ=255)

                                                                                                                                                                      g‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGCCTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCCTCGCATACTCACCCGCCACCTGACACGCGCCAATCGGCAATTCAATACGTTAAcccaccccaccc                                               >  1:1254785/1‑237 (MQ=255)

                                                                                                                                                                      g‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCCCCTGAAACGCGCCAACCGGCAATTAAGTACGTTCACGCAGccacaga                                               >  1:693651/1‑233 (MQ=255)

                                                                                                                                                                      g‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCCCTCACCTGATACGCGCCAACCGGCAATTCAGAACGTTCACGCAGCTcccg                                                >  2:566193/1‑238 (MQ=255)

                                                                                                                                                                          tccTCCCCCTTTAGGAGAGGAAAC‑AAAATAACCCAAG‑AATCTTCC‑TCCCAACC‑AAATCCCGCGCAAAGTAGCTGAAAACCAGACCCCCACCCGC‑ACTCTTCACCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCCTCT‑ACAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcacgcacg                                            <  1:1317498/237‑1 (MQ=255)

                                                                                                                                                                          aaCTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTtc                                                                                                                                                         <  1:1242042/131‑1 (MQ=255)

                                                                                                                                                                          aaCTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTtc                                                                                                                                                         >  2:1242042/1‑131 (MQ=255)

                                                                                                                                                                          aaCTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTcccgcacg                                            >  2:1029913/1‑239 (MQ=255)

                                                                                                                                                                              cccccTTTTTGAAAGGAAAA‑AAAACTAACGCAA‑AATCTTCT‑TCCCCGCC‑AAACCCAGCGCAAAGTAGCGGAAAATCAGATCCCCCCCCGC‑ACCCTTAACCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATGGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGt                                        <  1:1093252/239‑1 (MQ=255)

                                                                                                                                                                               ccccTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTc                                       >  2:78168/1‑239 (MQ=255)

                                                                                                                                                                                 ccTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGa                                     >  2:70195/1‑239 (MQ=255)

                                                                                                                                                                                 ccTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATAATCGCATCCTCACCGCTCACCTGATACGCCCCAATCGGCAATTCAGTACGTTcccccacccaccgcccaaggccca                                     >  1:851769/1‑216 (MQ=255)

                                                                                                                                                                                     tCTGAGAGGAAAC‑AAAATTAACGCAA‑AATCCTCC‑TCCCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGACCCCCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTa                                 <  2:1298962/239‑1 (MQ=255)

                                                                                                                                                                                           agGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTcc                           >  2:342522/1‑239 (MQ=255)

                                                                                                                                                                                           agGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCCCTCACCTGATACGCGCCAATCGGCAATTCAGGACGTTCACGCAGCCCCCGCACGAAGTCGAGGTACGCTcc                           >  1:1045522/1‑239 (MQ=255)

                                                                                                                                                                                             gAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGCAACGCGCCAATCGGCAATTCCGTACGTTACCGCAGCTcacgcaccaggccaaggaaccccccac                         >  1:193447/1‑219 (MQ=255)

                                                                                                                                                                                               aaC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAAt                                                                              >  1:690774/1‑184 (MQ=255)

                                                                                                                                                                                               aaC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCcagca                       >  2:15736/1‑239 (MQ=255)

                                                                                                                                                                                               aaC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAACCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAAt                                                                              <  2:690774/184‑1 (MQ=255)

                                                                                                                                                                                                     aaTTAACGCAA‑AACTTTTT‑TCCCAGCC‑AAACCCAGCGCAAAGTAGCTAAAAACCAGATCCGCCCCCCC‑CCGCTTAACCGAACCAAAGCTTTCGAGCACGACTTTCTCTTCATCT‑AAAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGttt                  <  2:577967/239‑1 (MQ=255)

                                                                                                                                                                                                     aaTTAAACCAA‑AATCTTCT‑CCCGCACC‑AAATCCAGCGAAAAGGAGTCGAAAACCAG‑CCCCCACCCCCAACGCTTAACCGAACCTAAGCTTTCGAACCCGACTTTCTCTTCAACT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGttt                  <  1:716796/239‑1 (MQ=255)

                                                                                                                                                                                                       ttAACGCAC‑AA‑CTTCC‑TCCCCGCCAAAATCAAGCGAAAAGTAGCTAAAAACCAAGCCCGCCCCCGC‑ACGCTTAATCGACCCTAGGCTTCCGAGCCCAACTTTCCCTTACTCT‑ATGGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTaa                <  1:827681/239‑1 (MQ=255)

                                                                                                                                                                                                         aaCGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAAcc              <  2:193447/239‑1 (MQ=255)

                                                                                                                                                                                                                 aaTCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGAACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCCCCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGGc       <  2:702141/239‑1 (MQ=255)

                                                                                                                                                                                                                      tct‑tccCAGCC‑AAATCCGGGACAAAGTAGCTGAAAATCAAATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCCCGACTTTCCCTTCATCT‑ATATCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGGCAGTCt  <  2:619894/239‑1 (MQ=255)

ACTATGCCAATGCCGATTACCTGCCACAGCGACGGATAATGCATCGGGCGACGATCCTTGGGGATAAACCGTTTTATCGGGCGGG‑GTAATACGAGAGTGGACGGTCCCCTCGCCCCCTTGGGGAGAG‑GGTTAGGGTGAGGGGGCGTTCACCGTACTTTCAACAGG‑TTAACTCCCCCTTTCTGAGAGGAAAC‑AAAATTAACGCAG‑AATCTTCT‑TCTCAGCC‑AAATCCAGCGCAAAGTAGCTGAAAATCAGATCCGCACCCGC‑ACGCTTAATCGAACCTAAGCTTTCGAGCACGACTTTCTCTTCATCT‑ATAGCACCCGCCAGCGCGGCGAACTTAATCATCGCATACTCACCGCTCACCTGATACGCGCCAATCGGCAATTCAGTACGTTCACGCAGCTCACGCACGATGTCGAGGTACGCTCCAGCAGGTTTAACCATCAGGCAGTCT  >  NC_000913/388558‑389005

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