BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,360	T→A	9.7%	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	9.7%	247.3 / 4.8	114	Y36F (TAT→TTT)	focA	formate channel
Reads supporting (aligned to +/- strand): ref base T (40/59); new base A (3/9); total (44/70)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGT  >  NC_000913/954126‑954595

caatttctCGCCAACTGACACCAGGTGATGCACCCACTCGCCTTAGCAACAACAATCAAGCCGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑a                                                                                                                                                                                                                                            <  2:1456916/235‑1 (MQ=255)

   tttttCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑aaaa                                                                                                                                                                                                                                         >  2:821074/1‑237 (MQ=255)

     tttCGCCAATTGACCCCAGGTGATGCGCCCACTTGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCCCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGACGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGTTAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TGGCGATAGAT‑AAAGTCt                                                                                                                                                                                                                                      >  2:1681429/1‑239 (MQ=255)

        cGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAATCTTaa                                                                                                                                                                                                                                   >  1:1299285/1‑239 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  2:1252711/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  1:1681429/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  2:193211/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  2:646554/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  2:650997/239‑1 (MQ=255)

                    ccAGGTGATGCGCCCACGCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  2:1402692/239‑1 (MQ=255)

                    ccAGGTGAGGCGCCCACTCGCCTTAGCAACAACAAGCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  2:349674/239‑1 (MQ=255)

                     cAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGG‑TGGTTCGt                                                                                                                                                                                                                     >  2:180771/1‑239 (MQ=255)

                        gTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TTGCCCGAAAG‑CAAACCTTAAGCGGA‑TGTTTCGtt                                                                                                                                                                                                                    >  2:737959/1‑237 (MQ=255)

                          gATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCt                                                                                                                                                                                                                 <  1:819399/238‑1 (MQ=255)

                           aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGCG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGGCTTTAGCGGA‑TGTTTCGTTGCttt                                                                                                                                                                                                               >  1:541617/1‑239 (MQ=255)

                           aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCtt                                                                                                                                                                                                                <  2:186840/238‑1 (MQ=255)

                           aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAA‑AAAGTCTTTAGCGGA‑TGTTTCGTTGCttt                                                                                                                                                                                                               >  2:746397/1‑239 (MQ=255)

                            tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGCTAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTta                                                                                                                                                                                                              >  1:1272228/1‑239 (MQ=255)

                            tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTta                                                                                                                                                                                                              >  1:1415151/1‑239 (MQ=255)

                                cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCg                                                                                                                                                                                                                      <  2:403230/227‑1 (MQ=255)

                                cccACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAAGCTTAAGCGGA‑TGGTTCGTTGCTTTATA‑Ga                                                                                                                                                                                                         >  2:1452710/1‑239 (MQ=255)

                                   aCTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑Gaca                                                                                                                                                                                                       >  1:1325284/1‑238 (MQ=255)

                                           tAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCctct                                                                                                                                                                                              >  1:193112/1‑239 (MQ=255)

                                             gCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑aa                                                                                                                                                                                                                                           <  1:297812/194‑1 (MQ=255)

                                               aaaaaaaaacAACACGGGGGAAGTAAAGAGAGCGGCTCCGAAGACAACAAAAAGAAACAGCCCCAGAGAAAAGAAAATG‑CCGCCAACCAGTTTTGCCAGGCCGAAGGGTAGTGGGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAAGGCGATTGAG‑ATGAAAACCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGGGGT‑TTTTTCGTTGCTTTATA‑GACCCCCGCCTCTTCgg                                                                                                                                                                                          <  2:739325/230‑1 (MQ=255)

                                                 caacaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCCGc                                                                                                                                                                                         >  1:1638621/1‑238 (MQ=255)

                                                   acaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCgggc                                                                                                                                                                                        >  1:789160/1‑235 (MQ=255)

                                                   acaaTCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGAAAG‑AAAGTCTTAAGCGGA‑TGTTTCCTTGCTTTATA‑GGCACCCGCCTCCTCGGCCa                                                                                                                                                                                       >  2:1160712/1‑238 (MQ=255)

                                                       tCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCCGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTg                                                                                                                                                                                  >  1:1355827/1‑239 (MQ=255)

                                                       tCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAATCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GAAACCCGCCTCTTCGGCCACTTTg                                                                                                                                                                                  >  1:1194209/1‑239 (MQ=255)

                                                           cccGGTGAAAAAAAAGAGACCGGCCCCGCAAAAAACAAAAAGAAGCACCCCCGGAGAGAAGAAAATG‑CCGCCACCCAGTTTTCCCATGCAGAGGGGATTTGTGCCGGG‑GCCA‑GTGGTTGTTGTAATAT‑AGAAGAGGAAGGCGATTGAG‑ATAAAACCCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTGGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCa                                                                                                                                                                              <  2:1620216/237‑1 (MQ=255)

                                                            acGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTTGCCAt                                                                                                                                                                             >  1:1560654/1‑239 (MQ=255)

                                                             cGGTGGAAATAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGCCGAATGCGATTGAG‑ATGAAAACCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCAtt                                                                                                                                                                            <  1:1506077/239‑1 (MQ=255)

                                                                            gaTCGGCACCGCAGAAAACAAAAAGAATCACCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAGGGGTATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AAAAGACGAATGCGATTGAG‑ATAAAAACCCCGGCGGTAA‑GCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaa                                                                                                                                                              <  2:1328079/238‑1 (MQ=255)

                                                                             aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCCCTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                            >  1:1617028/1‑239 (MQ=255)

                                                                             aTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TAGCCAGATAG‑AAAGTCTTAAGCGGG‑TGTTTCGTTTCTTTTTA‑GACACCCCCCCCTTCCGCCCCTTTGGGCCTTTCAGCCGGAAGTaaaa                                                                                                                                                            >  2:826810/1‑239 (MQ=255)

                                                                              gcGGCCCCGCAGACAAAAAAAAAAATACGCCCCAGAGAAAAGCAAAGG‑CCGCCAACCAGTTTTGCCAGGCCGAAGGG‑ATTTGGCCCGTGGCCA‑GGGGGTGCTGGAAAAT‑AAAAGACAAATGCGATTGAG‑AAAAAAACCCCGGGGGTAA‑TCGCCAGAAAG‑AAAGTCTAAAGCGGA‑TGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTaaaa                                                                                                                                                            <  2:399693/237‑1 (MQ=255)

                                                                               cGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTTGCCATTGCAGCCAGAAGTaaaaaa                                                                                                                                                          >  2:1114152/1‑237 (MQ=255)

                                                                                     cGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCCCTTTTGCCATTGCAGCAGGAAGTAAAAGATCaaaa                                                                                                                                                    >  1:379034/1‑239 (MQ=255)

                                                                                         agcaacaCAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑a                                                                                                                                                                                                                                            <  1:407877/147‑1 (MQ=255)

                                                                                         agcaacaCAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TAGCCAGATAG‑a                                                                                                                                                                                                                                            >  2:407877/3‑149 (MQ=255)

                                                                                           aacgacaAAGAGCCCGCCCGAGGGAAAAGAAAGTG‑CCCACAACCAGTTTTGACATGCGGAGGGTTATGGTGCATGG‑ACAA‑GTGGTGGCGGTAATAG‑AAAAGAAGAATTCGATTGAA‑AAAAAAACACCGGCGGAAC‑TCCCCAAATAA‑AATGTCTAAAGCGTT‑TGTTTTGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTg                                                                                                                                              <  2:516094/235‑1 (MQ=255)

                                                                                            aacaCAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAATCTTAAGCGGG‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAGATCAAAAAGGTTg                                                                                                                                              >  1:287540/1‑238 (MQ=255)

                                                                                                 aaaGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATAAAAACACCGGGGGTAA‑TCGCAAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCt                                                                                                                                        <  1:180771/239‑1 (MQ=255)

                                                                                                   agacacaccccaAGAGAGAGGCAAGGC‑CCGCAAACCGTTTTTCCCATCCGGAAGGCCTTTTTGCTTGT‑GCAA‑GTGGTGGCGGTAATAT‑AAAAACCGAATGCGATTGAG‑AGAAAAACCCCGGCGGTAA‑TCCCAAGATAG‑AAATTTTTAAGGGGT‑TTTTTCTTTGTTTTATA‑AACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCttt                                                                                                                                      <  2:878063/231‑1 (MQ=255)

                                                                                                           ccccAGAGAAAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCAGAGGGGCATTGGGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACCCCGGCGGTAA‑CCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg                                                                                                                                                                    <  2:12598/201‑1 (MQ=255)

                                                                                                            cccGGAAAGAAAAAAAGG‑CGGCTACCCTGTTTTGCCAGGCCGAGGGGCGTGGGCGATGA‑ACGA‑GTGGTTGTGGTGAAAA‑GAAAGCAGAATCGGATTGAG‑AAGAAAGCCCGGGGGGTAA‑TCCCCAGATAA‑AATGTATTAAGCGGA‑TGTTTTTTTGTTTTAAA‑AACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACt                                                                                                                                <  2:193112/236‑1 (MQ=255)

                                                                                                               agagagAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TACGCAGATAG‑AAAGTCCTTAGCGGA‑TGTTTCCTTGCTTTATC‑AACACCCGCCTCGTCGGCCACTTTGGCCCTTGCAGCAGGCAGACAAAAGTTAAAAGGGTTGTCAGCTTTCTAactcactc                                                                                                                           >  2:1633141/1‑238 (MQ=255)

                                                                                                                 gaggaAGCAAAGG‑CCGCCAACCAGTTTTGCCATGCGGAAGGGCATTGTGCCTGT‑GCCA‑GGGGTTGCTGTGATAT‑AAAAGACGAATGCGATTGAG‑ATAAAACCCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GCCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActctc                                                                                                                         <  1:737959/235‑1 (MQ=255)

                                                                                                                  ggaaaagaaaTGCCCGCAAACCGTTTTTGCCAGGCGGAAGGGCATTGTTCGTGG‑GCAA‑GTGGTTGTTGTGATAA‑AGAAGAAGAATGCGA‑TGAGAATAAAACCCCCGGCGGTAA‑TCGCAAGATAG‑AAATTCTAAAGCGGA‑GTTTTCGTTGCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAActc                                                                                                                           <  2:1089205/229‑1 (MQ=255)

                                                                                                                     aaGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGa                                                                                                                                                          >  1:12598/1‑201 (MQ=255)

                                                                                                                        cAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGGCTTAAGCGGA‑AGTTTAGTTGCTTTATA‑GACACCCGCCTCTTCCGCCACTTTGGGCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACTCTCTCTttatt                                                                                                                 >  1:284585/1‑239 (MQ=255)

                                                                                                                         aaaTG‑CCGCCAACCAGTTTTGCCAGGCGGAGGGGCATTGTGCCGGG‑TCCA‑GGGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATAAAACCCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTAAGGCGTA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTttat                                                                                                                  <  2:1299285/237‑1 (MQ=255)

                                                                                                                                 gcaacccATTTTTGCCATGCCGAAGGGAATTGTGCCTGT‑GCCG‑GTGGTGGCTGTGATAT‑AGAAGACGAAGGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGc                                                                                                         <  1:735529/234‑1 (MQ=255)

                                                                                                                                          ttttGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTTTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCCTCGGCCACCTTTGCCATTGCAGAAGGAAGTAAAAGATCACAAGGGTTGTTAGCTTTcac                                                                                                                                   >  1:1441567/1‑204 (MQ=255)

                                                                                                                                           tttGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAACCCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  1:497753/237‑1 (MQ=255)

                                                                                                                                           ttggCCAGGCAGAGGGGAATTGTGCCGGT‑GCAA‑GTGGTGGCTGTGATAA‑AGAAGACGAAGGCAATTGAG‑ATAAAAACCCCGGCGGTAA‑CCGCCAGAAAA‑AAAGTCTTAAGCGGA‑TTTTTCTTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGAt                                                                                                 <  1:821074/234‑1 (MQ=255)

                                                                                                                                            ttGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAAATAG‑AAAAGCTTAAGCGGA‑TGTTTCGTTGCTTTATT‑GACACCCGCCTCTTCGGCCACTTTGGCCCTTGCAGCAGGAAGGAAAAGATCAAAAGGGTTGTCAGCTTTCCCCCTAACTCTCCCTTTATTAAGTCGGCGACGAGATACt                                                                                              >  2:882305/1‑239 (MQ=255)

                                                                                                                                               ccATGCCGAAGGGCATTGTGCCGGT‑GCCA‑GGGGTTGCTGTGATAA‑AGAAGACGAATGCGATTGAG‑ATGAAAACCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTAAAGCGGA‑TGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTaa                                                                                            <  1:826810/238‑1 (MQ=255)

                                                                                                                                                ttgcccaaaGGCCATTGTCCCTG‑‑GCCATTTGTTTGTTGTGATAT‑AAAAGACAAGTGCGTTAGAG‑ATGAAAACACCGGGGATAA‑CCGCCAGATAG‑AAAGTTTTAATGTGA‑TGTTTCTTTGTTTTAAA‑GCAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                          <  1:1287763/232‑1 (MQ=255)

                                                                                                                                                cATGCGGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATAAAAACCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTAAAGCGGA‑TGTTTCGTTGTTTTATA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                          <  1:1114152/239‑1 (MQ=255)

                                                                                                                                                    ccGAAGGGCATTGTGCCTGT‑GCAA‑GTGGTTGCTGGGATAT‑AGAAGACGAAGGCGATTGAG‑ATGAAAACCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAGGCGGA‑TGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                                                                       <  1:746397/238‑1 (MQ=255)

                                                                                                                                                            ctttGTGCCTGT‑CCCG‑GTGGTTGCTGTGATTG‑AAAAGCCGAATGCGATTGAG‑ATGAAAACCCTGGGGGTAA‑GCGCAAGATAA‑AAATTCTTAAGGGGT‑TTTTTCGTTGTTTAAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAgg                                                                                                                                                                    <  2:1686316/151‑1 (MQ=255)

                                                                                                                                                                gtggcttgcgccA‑GTGGTGGCTG‑GAAATGAGAAGACGGAGGCGTTGAAG‑AGAAAAACCCCGGCGGGAA‑TCGCCAGATAG‑AAATTCTTAGGCGGG‑TGTTTCGTTTTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                           <  2:287540/230‑1 (MQ=255)

                                                                                                                                                                    ctgt‑gccA‑GTGTTTGCTGTAATAT‑AAAAGACAAATGCGATTGAG‑ATGAAACCACCGGCGTTAA‑TCGCAAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                       <  1:191841/238‑1 (MQ=255)

                                                                                                                                                                     ggtggccT‑GTGGTGGTTTTGATAA‑AAAAGACGAATGCGATTGAG‑ATGAAACCCCCGGCGGTAC‑TCGCCAGATAG‑AAAGTCTAAAGGGGA‑TGTCTCGTTGCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAAt                                                                      <  1:920741/235‑1 (MQ=255)

                                                                                                                                                                      gt‑gccA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGACGCAGGAAGTAAAAGa                                                                                                                                                          >  1:1686316/1‑153 (MQ=255)

                                                                                                                                                                         gccA‑GTGGTTGTTGTGAAAT‑AGAAGACGAAGGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGat                                                                  <  1:1160712/239‑1 (MQ=255)

                                                                                                                                                                          ccA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAACGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTTACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAACAATATAAAAGAGAAATTtata                                                                 >  2:107315/1‑239 (MQ=255)

                                                                                                                                                                           ctgttgttggCTGTGATAT‑AGAAGGGGAATGCGATTGAG‑AGGAACACACCGCTGGTAA‑TGGCCAGTTAG‑AATTTCTTAGGTGGA‑TGTTTCGTTGCTTAATA‑GACTCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  2:1588772/230‑1 (MQ=255)

                                                                                                                                                                           cA‑GTGTTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTTTTAAGCGGA‑TGTTTCGTTGTTTTATA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  1:647106/238‑1 (MQ=255)

                                                                                                                                                                           cA‑GTGGTTGCTGTGATAT‑AGAAGACGAAGGCGATTGAG‑ATAAAAACCCCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  1:1680382/238‑1 (MQ=255)

                                                                                                                                                                           cA‑GTGGTTGCGGGGATAT‑AAAAGACAAATGCGATGGAG‑AGAAAAACCCGGGCGGTAA‑CCGCCAGATAG‑AAAGTCTAAAGCGGA‑TGTTTCGTTGTTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  1:868158/238‑1 (MQ=255)

                                                                                                                                                                           cA‑GGGGTTGTTGTGATAA‑AAAAGACGAAGGCGATTGAA‑ATAAAAACACCGGCGTTAA‑TCGCAAGATAG‑AAATTCTTAAGGGGT‑TGTTTCGTGGCTTAATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGatat                                                                <  2:15985/239‑1 (MQ=255)

                                                                                                                                                                            gagggttttatgtgATGT‑AAAAGACGAATGCGATTGAG‑ATAAAAACCCCGGCGTAAA‑CCGCCAGAAGG‑AAAGTCTTAAGCGAG‑TGTTTCGTTGCTTTATA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                                 <  2:789160/228‑1 (MQ=255)

                                                                                                                                                                               tGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑AGTTTCGTTGCTTTATA‑GACACCCGCCTCTTTGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTTTTAAGACGGCGACGAGGAACTAAAAAAGCATTATAGATGAGAAATTGAtttaaa                                                             >  2:1322375/1‑237 (MQ=255)

                                                                                                                                                                                           tat‑aAAAGACGAAGGCGATTGAA‑ATAAAAACCCCGGCGGTAA‑CCGACAGATAG‑AAAGTTTAAGGCGGA‑TTTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGcc                                                 <  2:1663671/239‑1 (MQ=255)

                                                                                                                                                                                           aaaaagaagtCGAAGGCGATTGAG‑ATGAAAACCCGGGCGTTAA‑TC‑CCAAAAAGAAAAGTTTAAACCGGT‑TTTTGCGTTGTTTTAGA‑GCCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGc                                                  <  2:284585/235‑1 (MQ=255)

                                                                                                                                                                                             atggaagaCGAATGCGATGGAG‑AGAAAAACCCGGCCGGTAA‑TCGCCAGATAG‑AAAGTCTAAGGCGGT‑TGTTTCTTTGTTTTATG‑ACCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCt                                                <  1:506524/236‑1 (MQ=255)

                                                                                                                                                                                                   gaCGAATGCGCTTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑ACAACCTTAAGCGGA‑TGTTTCTTTGCTTTATA‑GACCCCCGCCTCTTCGGCCCCTTTGGCCATTTCAGCAGGAAGTAAAAATTTAAAAGGGTTGTCAGCTTTCCCACTAACTCTCTCTTTTTTAAGTCGGCGACGCGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCACGCCTGGCt                                            >  1:751079/1‑237 (MQ=255)

                                                                                                                                                                                                   gaCGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTta                                          >  2:1069068/1‑239 (MQ=255)

                                                                                                                                                                                                   gaCGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGAAAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGGTCAAAAGGGTTGTCAGATTTCACACTAACTCTCTCTTTATTAAATCGGCGACGAGATAATAACAAAGAATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTta                                          >  1:147667/1‑239 (MQ=255)

                                                                                                                                                                                                     gaagggggagtGAG‑ATG‑AAACACCGGCGGGAAGTCGCCAGAAAG‑AAAGTCTTAAGGGGA‑‑GTTTGCGTGTTTTAAAGGCCACCCGCCTCTGCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATATCGCCTGGCTta                                          <  2:1682161/227‑1 (MQ=255)

                                                                                                                                                                                                        aTGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGG‑TGTTTCGTTGATTTATA‑AACACCCGCCTCTTCGGCCAATTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAACCCTCCCTTTATTAAGTCGGGGACGAGATACTAACAAAGCATTATAGCTGAGAAATTGATTTAGATCATATCTTGCCTGGCATATAgg                                      >  2:357890/1‑238 (MQ=255)

                                                                                                                                                                                                              ttggagagaaaaaCCCCGGGGGTAA‑GCGCAAGAAAA‑AAAGTCTTAAGAGGA‑TGTTTCGTTGTTTAAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGt                                 <  1:1388968/229‑1 (MQ=255)

                                                                                                                                                                                                                gag‑aTGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGTTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACt                             <  1:1069068/239‑1 (MQ=255)

                                                                                                                                                                                                                     taaaacccccGGCTGTAA‑TCGCCAGATTA‑AATGTCTTAAGCGGT‑TGTTTTTTTGTTTTATA‑GCCCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                         <  1:301342/231‑1 (MQ=255)

                                                                                                                                                                                                                     taaaaaCCCCGGCGTTAA‑CCGCCGGATAG‑AAAGTCGAAAGCGGA‑TGTTTCGTTTCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                         <  1:63770/237‑1 (MQ=255)

                                                                                                                                                                                                                     taaaaaCACCGGCGGTAA‑TGGCAAGAAAA‑AAAGTCTGAAGAGGA‑TTTTTCGTTTTTTTATA‑AACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACGCGCa                         <  1:1129532/237‑1 (MQ=255)

                                                                                                                                                                                                                     tGAAACCCCCGGGGGTAG‑CCGCGAGAGAG‑AAAGTCTTAAGTGGG‑TTTTGGTTTGCTTAAAA‑GACGTCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                         <  1:357890/239‑1 (MQ=255)

                                                                                                                                                                                                                     agAAAACCCGGGGGTTAA‑CCCCAAGATAG‑AAA‑TCTTTAGGGGATTTTTTCGTTGTTTTAAA‑CACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCa                         <  2:1325284/238‑1 (MQ=255)

                                                                                                                                                                                                                      aaaaaCCCCGGCGGTAA‑TCCCCGGATAG‑AAAGTTTTAAGCGGA‑TGTTTCGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                        <  1:1322375/238‑1 (MQ=255)

                                                                                                                                                                                                                                      a‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAAt         >  1:570452/1‑238 (MQ=255)

                                                                                                                                                                                                                                       agcGCCAGAAAG‑AAATTCTTAAGCGGT‑TGTTTGTTTACTTTGTA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa        <  1:1040020/237‑1 (MQ=255)

                                                                                                                                                                                                                                       acccccAGATAA‑AAAGTGTAAAGGGGA‑TGTTTGTTTTCTATAAA‑AACAACCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATa        <  2:1355827/235‑1 (MQ=255)

                                                                                                                                                                                                                                         tgcacgATGG‑AAAGTTTAAGGCGGA‑TGTTTCGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAGGTAAAAGATCAAAAGGGTTGTCAGCTTTCTCACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACg      <  1:875911/234‑1 (MQ=255)

                                                                                                                                                                                                                                            cAGATAA‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTAAA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAg   <  1:1649144/239‑1 (MQ=255)

                                                                                                                                                                                                                                             agaaag‑aaATTCGTAAGCGGT‑TGTTGCGTTGCTTTATA‑GACACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt  <  1:1418791/239‑1 (MQ=255)

                                                                                                                                                                                                                                             aaaTAG‑AAGTTCTTAAGGGGA‑TGTTTCGTTGCTTAAGA‑GCAACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAg   <  2:792278/236‑1 (MQ=255)

                                                                                                                                                                                                                                                tAG‑AAAGTCTTAGCCGGT‑TGTTTCGTTGCTTTAAA‑GACCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACGGAGt  <  1:1411825/236‑1 (MQ=255)

CAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATG‑CCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGT‑GCCA‑GTGGTTGCTGTGATAT‑AGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGGTAA‑TCGCCAGATAG‑AAAGTCTTAAGCGGA‑TGTTTCGTTGCTTTATA‑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: `‑`