BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,360	T→A	11.2%	Y36F (TAT→TTT)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,360	0	T	A	11.2%	259.0 / 6.9	125	Y36F (TAT→TTT)	focA	formate channel
Reads supporting (aligned to +/- strand): ref base T (42/64); new base A (3/12); total (46/79)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGT  >  NC_000913/954126‑954595

caatttctCGCCAACTGACACCAGGTGATGCACCCACTCGCCTTAGCAACAACAATCAAGCCGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGa                                                                                                                                                                                                                                            <  1:1456916/235‑1 (MQ=255)

   tttttCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGaaaa                                                                                                                                                                                                                                         >  1:821074/1‑237 (MQ=255)

     tttCGCCAATTGACCCCAGGTGATGCGCCCACTTGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCCCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGACGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGTTAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TGGCGA‑T‑AGATAAAGTC‑t                                                                                                                                                                                                                                     >  1:1681429/1‑239 (MQ=255)

         gCCAACTGACACCAGGTGATGCGCCCACTCGCCTTAGGAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAg                                                                                                                                                                                                                                 <  1:572687/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      <  1:650997/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      <  1:646554/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      <  1:193211/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      <  1:1252711/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      <  2:1681429/239‑1 (MQ=255)

                    ccAGGTGAGGCGCCCACTCGCCTTAGCAACAACAAGCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      <  1:349674/239‑1 (MQ=255)

                     cAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGGTGGTTCGt                                                                                                                                                                                                                     >  1:180771/1‑239 (MQ=255)

                          gATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCt                                                                                                                                                                                                                 <  2:819399/238‑1 (MQ=255)

                           aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑CG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGGC‑TTTAGCGGATGTTTCGTTGCttt                                                                                                                                                                                                               >  2:541617/1‑239 (MQ=255)

                           aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCtt                                                                                                                                                                                                                <  1:186840/238‑1 (MQ=255)

                           aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAAAAAGTC‑TTTAGCGGATGTTTCGTTGCttt                                                                                                                                                                                                               >  1:746397/1‑239 (MQ=255)

                            tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑CTAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTta                                                                                                                                                                                                              >  2:1272228/1‑239 (MQ=255)

                            tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTta                                                                                                                                                                                                              >  2:1415151/1‑239 (MQ=255)

                                cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      <  1:403230/227‑1 (MQ=255)

                                cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCg                                                                                                                                                                                                                      >  2:403230/1‑227 (MQ=255)

                                cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAAGC‑TTAAGCGGATGGTTCGTTGCTTTATA‑Ga                                                                                                                                                                                                         >  1:1452710/1‑239 (MQ=255)

                                   aCTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑Gaca                                                                                                                                                                                                       >  2:1325284/1‑238 (MQ=255)

                                           tAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCctct                                                                                                                                                                                              >  2:193112/1‑239 (MQ=255)

                                             gCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGaa                                                                                                                                                                                                                                           >  1:297812/1‑194 (MQ=255)

                                             gCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGaa                                                                                                                                                                                                                                           <  2:297812/194‑1 (MQ=255)

                                               aaaaaaaaacAACACGGGGGAAGTAAAGAGAGCGGCTCCGAAGACAACAAAAAGAAACAGCCCCAGAGAAAAGAAAATG‑CCGCCAACCAGTTTTGCCAGGCCGAAGGGTAGTGGGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAAGGCGATTG‑AG‑ATGAAAACCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGGGGTTTTTTCGTTGCTTTATA‑GACCCCCGCCTCTTCgg                                                                                                                                                                                          <  1:739325/230‑1 (MQ=255)

                                                 caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCCGc                                                                                                                                                                                         >  2:1638621/1‑238 (MQ=255)

                                                   acaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCgggc                                                                                                                                                                                        >  2:789160/1‑235 (MQ=255)

                                                       tCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCC‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTg                                                                                                                                                                                  >  2:1355827/1‑239 (MQ=255)

                                                           cccGGTGAAAAAAAAGAGACCGGCCCCGCAAAAAACAAAAAGAAGCACCCCCGGAGAGAAGAAAATG‑CCGCCACCCAGTTTTCCCATGCAGAGGGGATTTGTGCCGGGGCCA‑GTGGTTGTTGTAATAT‑AGAAGAGGAAGGCGATTG‑AG‑ATAAAACCCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTGGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCa                                                                                                                                                                              <  1:1620216/237‑1 (MQ=255)

                                                            acGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTTGCCAt                                                                                                                                                                             >  2:1560654/1‑239 (MQ=255)

                                                                     gTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGGC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GAAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg                                                                                                                                                                    >  2:792278/1‑239 (MQ=255)

                                                                       aaaaagacgggtccgcaaaaaaacaCAAAAAACCAGCCCCAGGAAGAAGCAATGG‑CCCCCATCCAGTTTGGCTAGGCCGAAGGCTTGTGGGCCTGGGCCA‑GTGGTGGGTAGTAAAT‑AAAAGCCAAATGCGATTG‑AG‑AGGAAACCACCGGCGGTAA‑CCGCCA‑G‑ATAGAAATTC‑TGAAGGGGATTTTTCGTTGTTTTAAA‑AACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg                                                                                                                                                                    <  1:541617/216‑1 (MQ=255)

                                                                            gaTCGGCACCGCAGAAAACAAAAAGAATCACCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAGGGGTATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AAAAGACGAATGCGATTG‑AG‑ATAAAAACCCCGGCGGTAA‑GCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaa                                                                                                                                                              <  1:1328079/238‑1 (MQ=255)

                                                                             aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCCCTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                            >  2:1617028/1‑239 (MQ=255)

                                                                             aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TAGCCA‑G‑ATAGAAAGTC‑TTAAGCGGGTGTTTCGTTTCTTTTTA‑GACACCCCCCCCTTCCGCCCCTTTGGGCCTTTCAGCCGGAAGTaaaa                                                                                                                                                            >  1:826810/1‑239 (MQ=255)

                                                                               cGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTTGCCATTGCAGCCAGAAGTaaaaaa                                                                                                                                                          >  1:1114152/1‑237 (MQ=255)

                                                                                       cAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAgg                                                                                                                                                  >  1:1334059/1‑239 (MQ=255)

                                                                                         agcaacaCAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGa                                                                                                                                                                                                                                            <  2:407877/147‑1 (MQ=255)

                                                                                         agcaacaCAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TAGCCA‑G‑ATAGa                                                                                                                                                                                                                                            >  1:407877/3‑149 (MQ=255)

                                                                                            aacaCAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAATC‑TTAAGCGGGTGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAAGGTTg                                                                                                                                              >  2:287540/1‑238 (MQ=255)

                                                                                                 aaaGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCACCCGTTTTGGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTAATAT‑AGAAGACGAAGGCGATGG‑AG‑ATAAAACCACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCt                                                                                                                                        <  2:640238/239‑1 (MQ=255)

                                                                                                 aaaGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATAAAAACACCGGGGGTAA‑TCGCAA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCt                                                                                                                                        <  2:180771/239‑1 (MQ=255)

                                                                                                   aacgtcgccacgagaggaagacAAATC‑CCGCGAACCTGGTTGGCAATGCAGAAGGACATGGGGCAGGCGCTG‑GTGGGTGCTGTGATAT‑AGAAGACGACTGCTATCG‑AA‑ATAAAAACCCCGGCGTAAA‑TCGCCA‑G‑AAAGAAAGTC‑TTAAGCGGATTTTTCTTTGTTTTAAA‑AACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGCGTTGTCAGCttt                                                                                                                                      <  1:1236827/218‑1 (MQ=255)

                                                                                                   aGAACCAGCCCAAGAGAGAAGCAAATG‑CCGCCAACCATTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACCCCGGCGGTAA‑CCGCCA‑G‑ATAAAAAGTC‑TTAAGCGGATTTTTCTTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCttt                                                                                                                                      <  2:1199886/239‑1 (MQ=255)

                                                                                                            cccGGAAAGAAAAAAAGG‑CGGCTACCCTGTTTTGCCAGGCCGAGGGGCGTGGGCGATGAACGA‑GTGGTTGTGGTGAAAA‑GAAAGCAGAATCGGATTG‑AG‑AAGAAAGCCCGGGGGGTAA‑TCCCCA‑G‑ATAAAATGTA‑TTAAGCGGATGTTTTTTTGTTTTAAA‑AACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACt                                                                                                                                <  1:193112/236‑1 (MQ=255)

                                                                                                               agagagAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TACGCA‑G‑ATAGAAAGTC‑CTTAGCGGATGTTTCCTTGCTTTATC‑AACACCCGCCTCGTCGGCCACTTTGGCCCTTGCAGCAGGCAGACAAAAGTTAAAAGGGTTGTCAGCTTTCTAactcactc                                                                                                                           >  1:1633141/1‑238 (MQ=255)

                                                                                                                 gaggaAGCAAAGG‑CCGCCAACCAGTTTTGCCATGCGGAAGGGCATTGTGCCTGTGCCA‑GGGGTTGCTGTGATAT‑AAAAGACGAATGCGATTG‑AG‑ATAAAACCCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GCCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc                                                                                                                         <  2:737959/235‑1 (MQ=255)

                                                                                                                  ggaaaagaaaTGCCCGCAAACCGTTTTTGCCAGGCGGAAGGGCATTGTTCGTGGGCAA‑GTGGTTGTTGTGATAA‑AGAAGAAGAATGCGA‑TG‑AGAATAAAACCCCCGGCGGTAA‑TCGCAA‑G‑ATAGAAATTC‑TAAAGCGGAGTTTTCGTTGCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActc                                                                                                                           <  1:1089205/229‑1 (MQ=255)

                                                                                                                     aaGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGCTGTTTCGTTGCTTTATA‑AACACCCGGCTCTTCGGCCACTTTGGCCATTGCACGAGGAAGGAAAAGATCAAAAGGGTTGTCTGCTTTCACACCAACTCTCTCttt                                                                                                                    >  1:920741/1‑239 (MQ=255)

                                                                                                                     aaGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGa                                                                                                                                                          >  2:12598/1‑201 (MQ=255)

                                                                                                                         aaaTG‑CCGCCAACCAGTTTTGCCAGGCGGAGGGGCATTGTGCCGGGTCCA‑GGGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATAAAACCCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TAAGGCGTATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTttat                                                                                                                  <  1:1299285/237‑1 (MQ=255)

                                                                                                                                 gcaacccATTTTTGCCATGCCGAAGGGAATTGTGCCTGTGCCG‑GTGGTGGCTGTGATAT‑AGAAGACGAAGGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGc                                                                                                         <  2:735529/234‑1 (MQ=255)

                                                                                                                                    aaacttttttGACATGCAGAAGGGTATTGTGCCGGTGCCA‑GGGGTTGGTGTAAAAA‑AAAAAAAGAAGGAGAGTG‑AG‑ATAAAACCCCCGGGTGTAA‑CCGCCA‑G‑ATAAAAATTT‑TTAAGCGGTTGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                      <  1:1617028/233‑1 (MQ=255)

                                                                                                                                    aaCCATTTTTGCCAGGCGGAAGGGCTTTGTGCCTGGGCCA‑GTGGTTGCTGTGATAT‑AAAAGACAAATGCGATTG‑AG‑ATAAAAACACCGGCGGTAA‑TCGCCA‑A‑TAAGAAAGTC‑TTACGCGGATTTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                      <  2:240088/239‑1 (MQ=255)

                                                                                                                                          ttttGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTT‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCCTCGGCCACCTTTGCCATTGCAGAAGGAAGTAAAAGATCACAAGGGTTGTTAGCTTTcac                                                                                                                                   >  2:1441567/1‑204 (MQ=255)

                                                                                                                                           tttGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAACCCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  2:497753/237‑1 (MQ=255)

                                                                                                                                           ttggCCAGGCAGAGGGGAATTGTGCCGGTGCAA‑GTGGTGGCTGTGATAA‑AGAAGACGAAGGCAATTG‑AG‑ATAAAAACCCCGGCGGTAA‑CCGCCA‑G‑AAAAAAAGTC‑TTAAGCGGATTTTTCTTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  2:821074/234‑1 (MQ=255)

                                                                                                                                           gtgccaatgcggaggggattggggcctgtgccA‑GTGGTGGCTGGGATAA‑AGAAGACAAATGCGATGG‑AG‑AGAAAACCCCCGGGGGTAA‑CCGCCA‑G‑AAAGAAAGTT‑TAAAGCGGATTTTTCGTTTTTTTATA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  1:1415151/214‑1 (MQ=255)

                                                                                                                                            ttGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑A‑ATAGAAAAGC‑TTAAGCGGATGTTTCGTTGCTTTATT‑GACACCCGCCTCTTCGGCCACTTTGGCCCTTGCAGCAGGAAGGAAAAGATCAAAAGGGTTGTCAGCTTTCCCCCTAACTCTCCCTTTATTAAGTCGGCGACGAGATACt                                                                                              >  1:882305/1‑239 (MQ=255)

                                                                                                                                                ttgcccaaaGGCCATTGTCCCTG‑GCCATTTGTTTGTTGTGATAT‑AAAAGACAAGTGCGTTAG‑AG‑ATGAAAACACCGGGGATAA‑CCGCCA‑G‑ATAGAAAGTT‑TTAATGTGATGTTTCTTTGTTTTAAA‑GCAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                          <  2:1287763/232‑1 (MQ=255)

                                                                                                                                                cATGCGGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATAAAAACCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TAAAGCGGATGTTTCGTTGTTTTATA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                          <  2:1114152/239‑1 (MQ=255)

                                                                                                                                                    ccGAAGGGCATTGTGCCTGTGCAA‑GTGGTTGCTGGGATAT‑AGAAGACGAAGGCGATTG‑AG‑ATGAAAACCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAGGCGGATGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                                                                       <  2:746397/238‑1 (MQ=255)

                                                                                                                                                            ctttGTGCCTGTCCCG‑GTGGTTGCTGTGATTG‑AAAAGCCGAATGCGATTG‑AG‑ATGAAAACCCTGGGGGTAA‑GCGCAA‑G‑ATAAAAATTC‑TTAAGGGGTTTTTTCGTTGTTTAAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg                                                                                                                                                                    <  1:1686316/151‑1 (MQ=255)

                                                                                                                                                            cATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGa                                                                              >  2:1663671/1‑239 (MQ=255)

                                                                                                                                                               gtggcttgcgccA‑GTGGTGGCTG‑GAAATGAGAAGACGGAGGCGTTGA‑AG‑AGAAAAACCCCGGCGGGAA‑TCGCCA‑G‑ATAGAAATTC‑TTAGGCGGGTGTTTCGTTTTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                           <  1:287540/230‑1 (MQ=255)

                                                                                                                                                                gttgccgtggcgatggtttGCTGTGATAT‑AGAAGA‑GAAAGGGGTTGAAG‑AGAAAACCTCGGGCGTTAA‑TCGCAA‑T‑AAAAAAATCC‑TTAAGCGGATGTTTGGTTTTTTAATA‑GAGCTCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                           <  1:859454/222‑1 (MQ=255)

                                                                                                                                                                gggaaatgtgtaaaagagaaaggtgATAA‑AGAAGACGAATGGGATTG‑AG‑ATAAAAACCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TGAAGCGGATGTTTCTTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTCATTAAGTCGGCTACGAGATACTAACAAAGCATTATAGATga                                                                           <  2:1305198/217‑1 (MQ=255)

                                                                                                                                                                    ctgtgccA‑GTGTTTGCTGTAATAT‑AAAAGACAAATGCGATTG‑AG‑ATGAAACCACCGGCGTTAA‑TCGCAA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                       <  2:191841/238‑1 (MQ=255)

                                                                                                                                                                     cggtgcgggggttttggtgaataa‑AGAAAACGAAGCCTATTG‑GA‑AAGAAAACACCGG‑GGTAA‑CTGCCATA‑ATAG‑AAGTCTTTAGGGGGATTTTTCTTTGATTAATG‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                       <  1:1397905/218‑1 (MQ=255)

                                                                                                                                                                        gccA‑GTGGTTGTTGTGAAAT‑AGAAGACGAAGGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGat                                                                  <  2:1160712/239‑1 (MQ=255)

                                                                                                                                                                         ccA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAACGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTTACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAACAATATAAAAGAGAAATTtata                                                                 >  1:107315/1‑239 (MQ=255)

                                                                                                                                                                          ccatgggcgcgggtgAGAG‑AGAAGAAGAAGGGGGTTG‑AG‑AAAAAAACACCGGCGGTAG‑TCGCCA‑G‑AGAGAATGTA‑TGAAGCGGTTGTGTTGGTGTTTAATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  1:411127/227‑1 (MQ=255)

                                                                                                                                                                          cA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCAATTG‑AG‑ATGAAACCACCGGCGGTAA‑CCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  2:1452710/238‑1 (MQ=255)

                                                                                                                                                                          cA‑GTGGTTGCTGTGATAT‑AGAAGACGAAGGCGATTG‑AG‑ATAAAAACCCCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  2:1680382/238‑1 (MQ=255)

                                                                                                                                                                          cA‑GTGGTTGCGGGGATAT‑AAAAGACAAATGCGATGG‑AG‑AGAAAAACCCGGGCGGTAA‑CCGCCA‑G‑ATAGAAAGTC‑TAAAGCGGATGTTTCGTTGTTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  2:868158/238‑1 (MQ=255)

                                                                                                                                                                          cA‑GGGGTTGTTGTGATAA‑AAAAGACGAAGGCGATTG‑AA‑ATAAAAACACCGGCGTTAA‑TCGCAA‑G‑ATAGAAATTC‑TTAAGGGGTTGTTTCGTGGCTTAATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGatat                                                                <  1:15985/239‑1 (MQ=255)

                                                                                                                                                                           gagggttttatgtgATGT‑AAAAGACGAATGCGATTG‑AG‑ATAAAAACCCCGGCGTAAA‑CCGCCA‑G‑AAGGAAAGTC‑TTAAGCGAGTGTTTCGTTGCTTTATA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  1:789160/228‑1 (MQ=255)

                                                                                                                                                                              tGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGAAGTTTCGTTGCTTTATA‑GACACCCGCCTCTTTGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTTTTAAGACGGCGACGAGGAACTAAAAAAGCATTATAGATGAGAAATTGAtttaaa                                                             >  1:1322375/1‑237 (MQ=255)

                                                                                                                                                                                      tttgaaaaaaaagacaagggCGATTG‑AG‑ATGAAAACCCCGGGGGTAA‑TGCCTA‑T‑ATAGAAAGTC‑GAAAGCTGATGTTCCGTTGCTTTATA‑TGCTCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCatat                                                       <  2:1633141/219‑1 (MQ=255)

                                                                                                                                                                                           aaaaaaagtggagtgtgATGG‑AG‑AAGAAAACACCGGGGGTGA‑‑CACCA‑GTATGAAAAGTT‑TGAAGCGTTGGTTTCTTTCCTTTAAA‑CACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  1:1219004/223‑1 (MQ=255)

                                                                                                                                                                                                  gaCGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTta                                          >  1:1069068/1‑239 (MQ=255)

                                                                                                                                                                                                  gaCGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCCCTTCGGCCACTTTTGCCATTGCAGCAGGGAGTAAAAGATCAAACGGGTTGTCAGCTTTCACACTAACTCTCCCTTTATTAAGTCGCCGCAGGGATACTAACAAAGCATTATAGAAGAGAAATTGaaatttgtctcttttacacatcttc                                           >  1:1305198/1‑217 (MQ=255)

                                                                                                                                                                                                  gaCGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑AAAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGGTCAAAAGGGTTGTCAGATTTCACACTAACTCTCTCTTTATTAAATCGGCGACGAGATAATAACAAAGAATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTta                                          >  2:147667/1‑239 (MQ=255)

                                                                                                                                                                                                    gaagggggagtG‑AG‑ATG‑AAACACCGGCGGGAAGTCGCCA‑G‑AAAGAAAGTC‑TTAAGGGGA‑GTTTGCGTGTTTTAAAGGCCACCCGCCTCTGCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATATCGCCTGGCTta                                          <  1:1682161/227‑1 (MQ=255)

                                                                                                                                                                                                       aTGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGGTGTTTCGTTGATTTATA‑AACACCCGCCTCTTCGGCCAATTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACCCTCCCTTTATTAAGTCGGGGACGAGATACTAACAAAGCATTATAGCTGAGAAATTGATTTAGATCATATCTTGCCTGGCATATAgg                                      >  1:357890/1‑238 (MQ=255)

                                                                                                                                                                                                       aTGCGATTG‑AG‑ATGAAAACACCGGCGGAAA‑TTGCCA‑G‑ATCGCAAGAC‑TTAAGTGGTTGATTCCTTGCTTTATA‑ACCACCCGACTCTTCGGCCACTTTTGTCCTTGGATAAGGAAATACACAATCAAAGGGGTTTTTACAGTTACCACTTAAACTCTCTTTATTAAAACGGCGACGCGAAACTAACAAAGCATAATTACTGTCAAATTGGTATAGATCCTAtctctattctacatctta                                      >  1:926996/1‑224 (MQ=255)

                                                                                                                                                                                                              ttggagagaaaaaCCCCGGGGGTAA‑GCGCAA‑G‑AAAAAAAGTC‑TTAAGAGGATGTTTCGTTGTTTAAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGt                                 <  2:1388968/229‑1 (MQ=255)

                                                                                                                                                                                                               g‑ag‑aTGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACt                             <  2:1069068/239‑1 (MQ=255)

                                                                                                                                                                                                                     taaaacccccGGCTGTAA‑TCGCCA‑G‑ATTAAATGTC‑TTAAGCGGTTGTTTTTTTGTTTTATA‑GCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                         <  2:301342/231‑1 (MQ=255)

                                                                                                                                                                                                                     taaaaaCCCCGGCGTTAA‑CCGCCG‑G‑ATAGAAAGTC‑GAAAGCGGATGTTTCGTTTCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                         <  2:63770/237‑1 (MQ=255)

                                                                                                                                                                                                                     taaaaaCACCGGCGGTAA‑TGGCAA‑G‑AAAAAAAGTC‑TGAAGAGGATTTTTCGTTTTTTTATA‑AACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACGCGCa                         <  2:1129532/237‑1 (MQ=255)

                                                                                                                                                                                                                     tGAAACCCCCGGGGGTAG‑CCGCGA‑G‑AGAGAAAGTC‑TTAAGTGGGTTTTGGTTTGCTTAAAA‑GACGTCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                         <  2:357890/239‑1 (MQ=255)

                                                                                                                                                                                                                      aaaacCACCGGCGGAAA‑TCGCAA‑G‑ATAGAAAGTC‑TAAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                        <  2:882305/238‑1 (MQ=255)

                                                                                                                                                                                                                      aaaaaCCCCGGCGGTAA‑TCCCCG‑G‑ATAGAAAGTT‑TTAAGCGGATGTTTCGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                        <  2:1322375/238‑1 (MQ=255)

                                                                                                                                                                                                                              cgtgggTAA‑TCGCAG‑A‑ATAGAAAGTC‑TGAGGCGGATTTTTCGTGGCTTTATA‑CACCACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTa                <  1:1189752/235‑1 (MQ=255)

                                                                                                                                                                                                                                      a‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAAt         >  2:570452/1‑238 (MQ=255)

                                                                                                                                                                                                                                       acccccA‑G‑ATAAAAAGTG‑TAAAGGGGATGTTTGTTTTCTATAAA‑AACAACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa        <  1:1355827/235‑1 (MQ=255)

                                                                                                                                                                                                                                          tgcacg‑ATGGAAAGTT‑TAAGGCGGATGTTTCGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAGGTAAAAGATCAAAAGGGTTGTCAGCTTTCTCACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACg      <  2:875911/234‑1 (MQ=255)

                                                                                                                                                                                                                                               aaaTAGAAGTTC‑TTAAGGGGATGTTTCGTTGCTTAAGA‑GCAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAg   <  1:792278/236‑1 (MQ=255)

                                                                                                                                                                                                                                                  tAGAAAGTC‑TTAGCCGGTTGTTTCGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt  <  2:1411825/236‑1 (MQ=255)

CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTG‑AG‑ATGAAAACACCGGCGGTAA‑TCGCCA‑G‑ATAGAAAGTC‑TTAAGCGGATGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGT  >  NC_000913/954126‑954595

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