BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,107,591	C→G	7.9%	A14G (GCG→GGG)	pfkA →	6‑phosphofructokinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,107,591	0	C	G	7.9%	398.0 / 5.2	163	A14G (GCG→GGG)	pfkA	6‑phosphofructokinase I
Reads supporting (aligned to +/- strand): ref base C (60/80); new base G (11/3); total (79/84)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 172 total aligned reads displayed.

GGTCTATGCTTTCATAATCAGTATAAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTG  >  NC_000913/4107354‑4107820

gggttgttgtttCATAACCAGTAAAAAAGGGGGCC‑CACCCCCCGTTTTTGTGAAAAAATTCCTCCATGTTGCCTTACCTGAATCAATTCAGCAGGAAGTGATTGTTCTACTAGTTGCACATTCGTTGGATCGCTTCGATGTGCAAGAAGGCTTCCGGCAACAGATTTCAGTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑gcg                                                                                                                                                                                                                                       <  1:723833/230‑1 (MQ=255)

 ttCTTTGCTTTAAAAATCGGTAAAAAGGGGGGCCAAGTCCGCATTTGGTTT‑ATAAATTACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCAGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑gcgc                                                                                                                                                                                                                                      <  2:838737/238‑1 (MQ=255)

  tCTATGCTTTCATAATCAGTATAAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑gcgc                                                                                                                                                                                                                                      <  2:923184/238‑1 (MQ=255)

  tCTATGCTTTCATAATCAGTATAAAAGAGAGCAAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑gcg                                                                                                                                                                                                                                       <  2:832106/237‑1 (MQ=255)

   ttttgttttAATAACCAGTTAAAAAGAGAGCCCGACCCGCTTTTTGGGT‑ATAAAATTCCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGAGTGTTATACTATTTGCAAATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCa                                                                                                                                                                                                                                    <  1:460574/233‑1 (MQ=255)

   cTATGCTTTCATAATCAGTATAAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGCCTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGCCATGATTAAGAAAATCGGGGGGTTGACAAACGGC‑GGGGAT‑GTGccc                                                                                                                                                                                                                                    >  2:1276346/1‑238 (MQ=255)

   cTATGCTTTCATAATCAGTATAAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTAATTTTGCATTCCAAAGTTCAGAGGAATCCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTTAT‑GCGccg                                                                                                                                                                                                                                    >  2:1047340/1‑238 (MQ=255)

    tATGCTTTCATAATCAGTATAAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTAAGAGGTAGTCATTGTTAAGAAAATCGGTGTGTTGACAACCGGC‑GGGGAT‑GGGACAg                                                                                                                                                                                                                                   >  1:239836/1‑239 (MQ=255)

     aTGCTTTCATAATCAGTATAAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGCCAAGCGGC‑GGTGAT‑GGGCCggg                                                                                                                                                                                                                                  >  2:329080/1‑239 (MQ=255)

            cATACTCAGCATAAAAGAGAGCCCGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTAGCCTGACCTGAATCAATTCAGCAGGAAGAGATTGTTCTACTATTTGCACATTCGTTGGATAACTTCGATGGGCAAGAAGCCTTCCGGCACCAGATTTCTTTTTGCTTTCCAAGGCTAGGAGGGAGCCCTGATTAAGGACATCGTGGGGTTTACACCCTGA‑GGGGAG‑GCGCCAGGGATGaaa                                                                                                                                                                                                                           >  1:138763/1‑238 (MQ=255)

            aaTAATCAGTATAAAAGAGAGCCACACCCGCATTTTGTGT‑CTAAATTACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTTTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAAc                                                                                                                                                                                                                           <  1:745745/238‑1 (MQ=255)

                     tataAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATAGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCCGGGATGAACGCCGCAAtt                                                                                                                                                                                                                  >  2:1246302/1‑239 (MQ=255)

                     tataAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGCCATGATTAAGAAAATCGTTGTTTTTACAAGCGGC‑GGGGAT‑GCGCCAGGCATGAAAGCCGGAAt                                                                                                                                                                                                                   >  1:1198800/1‑238 (MQ=255)

                           agaggcccgtccccgatttttGTTT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTTATTGTTATTCTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGCCATGATTAAGAAAAGCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCggg                                                                                                                                                                                                            <  1:487304/222‑1 (MQ=255)

                               agccagccCCGCATTTTGAGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTCTTATCCTATTTGCACCTTCGTTGGATCACTTCGATGTGCAAGATGACTTCCGGTAACAGCTTTCATTTTGCATTCCAAAGTTCACAGGTAGTAATGATTCAGAAAACCGGGGTGTTGACAAGCGGC‑GGGGAT‑GCGCCGGGGATACAAGCCGCAACTACCGGGgtgt                                                                                                                                                                                                        >  1:995319/1‑237 (MQ=255)

                               agCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCACCAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GGGCCGGGGCTGAAAGCCCCAATTCCCGGGgttg                                                                                                                                                                                                        >  2:88595/1‑239 (MQ=255)

                               agCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCCGGCATGAACGCCGCAATTCGCGGGggtg                                                                                                                                                                                                        >  2:108167/1‑239 (MQ=255)

                               agCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCTCGGGgtgg                                                                                                                                                                                                        >  2:69633/1‑237 (MQ=255)

                               agCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGGGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGgtgg                                                                                                                                                                                                        >  2:418984/1‑237 (MQ=255)

                               agCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGGAGTCATGATTAAGAAAATCGGTGGGTTGCCAAGCGGC‑GGTGAT‑GTGCCAAGCATGAACGCCCCAATTTCGGGGgttt                                                                                                                                                                                                        >  2:1285502/1‑238 (MQ=255)

                               agCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGAAGTCATGATTAAGAAAAACGGTGTGTTGACAAGCGGA‑GGGGAT‑GCGACAAGCAATAAAGTCGCAAATCggggggggg                                                                                                                                                                                                        >  1:223895/1‑236 (MQ=255)

                               agCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGGGCAAGAAGACTTCCGGCCACAGATTTCGTTTTGCATTCCAAAGTTCAGTGGAAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGCCGGGGGAT‑G‑GCCAGGAAAGAAAGCCACAATTGGAggggggg                                                                                                                                                                                                        >  1:1106863/1‑236 (MQ=255)

                                 cctaccccgttttttGGGT‑TTAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATTCTATTTGCACATTCGTTGGATTACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGgttgt                                                                                                                                                                                                       <  1:1285502/227‑1 (MQ=255)

                                 ccgggcccGCATTTTGTGT‑AAAAAATTCCGCCAATTGGGCGGACCTGAATCAATTCAGCAGGAAGTGATTGTTATTATAAGTGCAAATTCGTTGGAACACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTGCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGgttgt                                                                                                                                                                                                       <  1:373087/233‑1 (MQ=255)

                                 cccgACCCGTATTTTGGTT‑AAAAAATACCTCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGGTATACTATTGGCACATTCGTTGGATCACTTCGATGGGCAAGAAGACTTCCGGCAACAGATTGCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGgttgt                                                                                                                                                                                                       <  1:108167/235‑1 (MQ=255)

                                 ccAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑TCGCCAGGCCTGAAAGCCCCCATTCCCCGGgtgtgt                                                                                                                                                                                                      >  2:968407/1‑235 (MQ=255)

                                     aCCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGGGAT‑GCGCCAGGCATGAACGCCGCAATTCTCGGGGTTgttcgg                                                                                                                                                                                                   >  2:484527/1‑237 (MQ=255)

                                              ttGTGT‑ATAAGATACCGCCATTTGGCCTGACCTTAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGa                                                                                                                                                                                         <  2:1158333/239‑1 (MQ=255)

                                                 tgt‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGGGATTGTTATACTATTTGCCCATTCGTTGGATCACTTCGATGGGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCACAGTTCAGAGGTATGCATGATTAAGAAAATCGTTGTGTTGACAAACGGC‑GGTGGTAGCG‑CAGGCACGAAAGACGGAATTCGCggggggtatagtgctcgcccggaca                                                                                                                                                                                      >  1:1173945/1‑218 (MQ=255)

                                                  tt‑ataAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGa                                                                                                                                                                                     <  2:451694/238‑1 (MQ=255)

                                                  gt‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGa                                                                                                                                                                                     <  1:972978/239‑1 (MQ=255)

                                                  gt‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGa                                                                                                                                                                                     <  2:499263/239‑1 (MQ=255)

                                                   ttaaaaaaatccGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGa                                                                                                                                                                                     <  2:662454/229‑1 (MQ=255)

                                                   ttaaaaaaTTCCGCCTTTTGGCCTTACCTGGATGCATTCAGCAGGAAGTTATTGTTAAGCTATTTGCACATTCGTTGGATCGCTTCGAGGTGCAAGAAGACTGCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGATGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCTCTGACAGa                                                                                                                                                                                     <  1:1276346/235‑1 (MQ=255)

                                                   t‑ataAAATACCGACATTTGGCCTGACCTGAATCAATTCAGCAGGATGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGTTGTGCAAGAAGACTTCAGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGAAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GTGGCT‑GGGCCAGGCATGAACGGAGCAATTCCCGCGGTAGTTCGTTGTGCGCTGACAGaa                                                                                                                                                                                    >  2:1094031/1‑239 (MQ=255)

                                                           tACCGCCATTTGGCCTGAACTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTgg                                                                                                                                                                             <  2:719491/239‑1 (MQ=255)

                                                              cGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCCCGGGGTTGTTCGTTCTGCGCCGACGAAAGGTCTGGAAg                                                                                                                                                                          >  2:894396/1‑239 (MQ=255)

                                                              cGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCCCACTTCCCGGGGGTGTTCTTTCTGCGCTGCCAGAAGGGCTGGagg                                                                                                                                                                          >  2:768759/1‑237 (MQ=255)

                                                                ccATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTa                                                                                                                                                                        <  2:223895/239‑1 (MQ=255)

                                                                   tttGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCCATTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATg                                                                                                                                                                     <  2:286627/239‑1 (MQ=255)

                                                                           gccctgaATCAATTTAGCAGGAAGTGATTTTTATACTAGTTGCACATTCGTTGGATCACTGCGATGTGAAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGGCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTa                                                                                                                                                             <  1:610399/237‑1 (MQ=255)

                                                                           gacctgaATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GGGCCAGGCATGAACGCCGCAATTCGCGGGGTGGTTCGTTCTGCGCTGACAGAAGGACTGGAAGTAATGGGTATTTa                                                                                                                                                             >  2:940623/1‑239 (MQ=255)

                                                                            acctgaATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑TCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTTCGCTGGCAGAAGGGCTGGGAGTAATGGGGAttttt                                                                                                                                                            >  2:70723/1‑237 (MQ=255)

                                                                            acctgaATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTTCCGGGGTTGTTCGTTCTGCGCTGACAGAAGGGCTGGGAGTAATGGGGATTTAt                                                                                                                                                            >  2:458648/1‑239 (MQ=255)

                                                                             actgaATCAATTCAGCAGGAAGTTATTGTTATACTATTTGCACATTCGTTGGATCCCTTCGATGTGCAAGAAGACTTCCGGCAACAGATTCCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTAt                                                                                                                                                            <  2:1040231/237‑1 (MQ=255)

                                                                               ttgattaaattagaagGTAGTGATTGTTATACTATTTGCAAATTCGTTGGATCACTTCGATGTGCAAGAAGGTTTCCGGCAACAGATTTCATTTTGCATGCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGTGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGAc                                                                                                                                                         <  1:1094059/225‑1 (MQ=255)

                                                                                 aaTCAATTCAGCAGGATGTTGTTGTTATACTTTTTGCACATTCGTTGGATCACTTCGATGTGCAAGGAGGCTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACgg                                                                                                                                                       <  1:880502/239‑1 (MQ=255)

                                                                                 aaTCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGGGAC‑GAGCCAGGCATTAACGCCGCAAATTGCGGGGGTGTTTGTTCTCCGCTgtgaaaagagctgcga                                                                                                                                                                           >  1:839034/1‑204 (MQ=255)

                                                                                      aTTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCCATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGAAATGGGGATTTTTTACGGcca                                                                                                                                                    >  2:516268/1‑235 (MQ=255)

                                                                                      aTTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGGTCGTTTTGCGCTGACAGAAGGTCTGGAAGGAATGGGGATTTTTGACGGCTa                                                                                                                                                    >  2:1002860/1‑237 (MQ=255)

                                                                                       ttCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑gcgc                                                                                                                                                                                                                                      <  2:467425/154‑1 (MQ=255)

                                                                                            gtggAAGTGATTTTTTTTCTATTTGCACATTCGTTGGATCACTTCGATGTTCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGGGCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTc                                                                                                                                            <  1:1204642/237‑1 (MQ=255)

                                                                                            gagGAAGTGATTTTTATACTATTTGCACATGTGTTGGATCACTTCGATGGTCAATAAGGCGTCCGGCAACAGATTTCAGTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAAGCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTc                                                                                                                                            <  1:287146/238‑1 (MQ=255)

                                                                                            cgtgtaGTGATTTGTTTTCTTTTTGCACTTGTGTTTGATCCCTTCGCTGTGAAAGAAGTCTGCCGGCAACAGAATTCATTTTGCATTCCAAAGTTCAGAGGGAGGTATGAGTAAGAAAATCGGTGTGTTGACAAGTGGC‑GGTGAT‑GCTCCAGGCATGAACGCAGGAATTCGCGGGGTTGTTCGTTCTGCGCTGGCAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTAGCTGGGTc                                                                                                                                            <  1:633538/234‑1 (MQ=255)

                                                                                             agggaggtattattATACTATTTGCCCATTCGTTTGATCACTTCGATGTGCAAGAAGGCTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGACATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTc                                                                                                                                            <  2:688121/230‑1 (MQ=255)

                                                                                             agGATGGGATTTTTATACTTTTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACATCCGGCAACAGATTTCATTTTGCATTCCAAAGGTCAGGGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTc                                                                                                                                            <  1:940623/238‑1 (MQ=255)

                                                                                                aaGTGATTTTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTc                                                                                                                                            <  2:1128830/235‑1 (MQ=255)

                                                                                                 aGTGATTTTTATACTATTTGCACTTTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTgcg                                                                                                                                                                                             <  2:1258611/185‑1 (MQ=255)

                                                                                                    gATTGTTATATTTTTTGCACATTTGTTGGATCACTGCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTGGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGaa                                                                                                                                    <  1:1002860/239‑1 (MQ=255)

                                                                                                         tttttctttttGCCCTTTCGTTGGATTCCTTCGATGTGCAAGAAGGCTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTTGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCg                                                                                                                               <  1:491809/231‑1 (MQ=255)

                                                                                                           atgcTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTa                                                                                                                             <  2:186070/236‑1 (MQ=255)

                                                                                                              ctttttGCCAATTCGGTGGCTTCCTTCGAGGTGGAAGAATGCCTCCGGCAAAAGAGTTCATTTTTGATTCCAAAGGTCAGAGGTAGTCATGATTAAGAAAATAGGTGTGTGTACGAGGGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTGTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATgg                                                                                                                          <  1:1246302/236‑1 (MQ=255)

                                                                                                              ctttttGCACATTCGTTGGATCACTTCGAAGTGCAAGAAGACTTCCGGCAACAGATTTAATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATgg                                                                                                                          <  2:883594/236‑1 (MQ=255)

                                                                                                              cTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACgg                                                                                                                                                       <  2:147558/210‑1 (MQ=255)

                                                                                                              cTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATgg                                                                                                                          <  2:1299865/239‑1 (MQ=255)

                                                                                                               tATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGGGAA‑GCGCCAGGCATGAAAGCCGCAATTCGCGGGGTTGTTGGGTCTGCGCTGACGGAAGGGCTGGGAGCAATGGGGCATTATGAGGGCTATCTGGGGTTGTGGTAAGACCGAATGGt                                                                                                                         >  1:938922/1‑239 (MQ=255)

                                                                                                                aTTTGCCAATTCGTTGGATCCCTTCGATGTGCAAGAAGGCTTCCGGCAACCGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAAGCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTa                                                                                                                        <  2:239836/239‑1 (MQ=255)

                                                                                                                 tttGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAGCAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGCT‑GGGCCAGGCATGCACGCCGCAATTTGCGGGGGTGTTCGTTCTGCGCTGACAGAGGGGCTGGGACGAATGGGGATTTATGACGGCTATCTGGGGGTGTagatagaaggagagggaa                                                                                                                       >  1:1326824/1‑222 (MQ=255)

                                                                                                                 tttGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GTGCCAGGCATGAACGCCGCAAATTGCGGGGGTGTTCGGTTTGGGCTGGAAGAAGGTCTGGAAGGAATGGGGATTTTTGGGGGCTATTCGGGGCGGTATGAAGGCCgcagtgga                                                                                                                        >  1:277503/1‑231 (MQ=255)

                                                                                                                 tttGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GTGCCAGGCATGAAAGCCGCAATTCGCGGGGGTGTTCGTTCTGCGCTGACAGAAGGTATGGAAGAAATGGGGATTTATGAGGGCTATCTTGGTCTGTATGAAGACCGTATGGGAc                                                                                                                       >  2:82626/1‑239 (MQ=255)

                                                                                                                   tGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAAAGCCGCAAtttg                                                                                                                                                                                                                >  1:781843/1‑146 (MQ=255)

                                                                                                                    gacattttgTTGGATCACTTCGATGTGCAAGAAGACTTCCGGAAAAAGATTTCAGTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGc                                                                                                                    <  1:69633/231‑1 (MQ=255)

                                                                                                                      acaTTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTa                                                                                                                  <  2:568805/239‑1 (MQ=255)

                                                                                                                        aTTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGa                                                                                                                <  2:950250/239‑1 (MQ=255)

                                                                                                                         ttCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGa                                                                                                                >  2:426024/1‑238 (MQ=255)

                                                                                                                         ttCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAAGCCGGATGGTACAGCTTGa                                                                                                                >  2:721307/1‑238 (MQ=255)

                                                                                                                         ttCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGACGCAATTCGCGGGGTTGTTCGTTCTGCGCTGCCCGAAGGTCCGGAAGTAATGGGTATTTATTCCGGCGATCTGGGCGTGTATGAAGACCGTATGGGACAGCCAGa                                                                                                                >  2:376703/1‑238 (MQ=255)

                                                                                                                                  tctcTTCTATGTGCCAGAAGGCTGCCGGCATCAGATTTCAGTTTGCCTTCAAAAGGTCAGAGGTATGCGTGATTAAGAAAAGCGGTGGTTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGACGCAAGTCGTGGTGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑g                                                                                                      <  1:272874/235‑1 (MQ=255)

                                                                                                                                  tCACTTCGATTGTCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑g                                                                                                      <  2:293655/238‑1 (MQ=255)

                                                                                                                                    aCTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGGATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑CCgg                                                                                                   >  2:976692/1‑238 (MQ=255)

                                                                                                                                         cttgtgAAAAAAGATTGCCGGCAAAAGATTTGATTTTGCATTCCAAAGTTCAGAGGTAGTCATGTTTAAGAAAATCGGTTTGTTGGCAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTg                                                                                              <  1:836014/237‑1 (MQ=255)

                                                                                                                                           tgtgCAAGAAGGCTTCCGGCAGCAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCAGGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGAc                                                                                            <  2:1326824/239‑1 (MQ=255)

                                                                                                                                             tgCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAAAGCCGCAATTCGCGGGGGTTGTCGTTCGGCGCTGTCGGAAGGGCTGGGAGTAAGGGGGATTTTTGACGGCCTCCTGGGCGTGTATGAAAACCGTATGGGACGCGTAAACCGTTAAA‑gggggctcggaaa                                                                                          >  1:817080/1‑227 (MQ=255)

                                                                                                                                                       cTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGGCAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGGCTGTTTGAAGACCCTATGGTACAGCTAGACCGTTACA‑GCGGGTCTGACATGATCAACCGt                                                                                >  2:597311/1‑239 (MQ=255)

                                                                                                                                                       cTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGGAATGGGGATTTATGACGGCTATCTGGGGCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑CCGTGTCTGACATGATCAACCGt                                                                                >  2:1063882/1‑239 (MQ=255)

                                                                                                                                                       cTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAAAGCCGCAAATCGCGGGGGTGTTCGTTCTGCGCCGGCCGAAGGTCTGGGAGGCATGGGTATTTTTGATGGTGTTCTGGGGCGAGATTGAGACCGAATGGTACAGCTAGACCGTTAAA‑GAGGGTGTGAAATGAGAAACCGt                                                                                >  1:586432/1‑239 (MQ=255)

                                                                                                                                                        ttCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCCATTCGCGGGGGTGTTCGGTCTGCGCTGCCCAACAGGATGGGAGTAATGGGGATTTTTGACGGCTTTCTGGGTCTGTATGAAGACCGCATGGGACAGATAGACCGTTACACGCG‑GTCTGAAATGaataaacagg                                                                               >  1:903009/1‑230 (MQ=255)

                                                                                                                                                          ccGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGAATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTTCA‑GCGGGTCTGACATGATCAACCGGGGc                                                                             >  2:655207/1‑239 (MQ=255)

                                                                                                                                                               aaCATATTTCATTTTGCATTCCAAAGTTCGGGGGTAGTCATGATTAAGAAAAGCGGGTTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTa                                                                         <  1:581814/238‑1 (MQ=255)

                                                                                                                                                                     ttttcttttGCTTTCCAAAGTTCGGGGGTGTTTATGATTAAGGAAATCGGTGTGTCGACAAGCGGC‑GGTGAT‑GAGCCAGGCATGGACGCCGTAATTCGCGGGGTTTGTCTTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTACGTTCCt                                                                  <  1:652729/234‑1 (MQ=255)

                                                                                                                                                                      ttCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGGTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGt                                                                                                                              >  2:1203879/1‑179 (MQ=255)

                                                                                                                                                                          ttttGCCTTCCAAAGTTCAGAGGTAGGCTTGATTAAGAAAAGCGTGTTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTACGTTCCTCGGtt                                                             <  1:82626/239‑1 (MQ=255)

                                                                                                                                                                                 ttCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGcgcg                                                      <  2:817080/239‑1 (MQ=255)

                                                                                                                                                                                      aaGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTcc                                                 <  2:583969/239‑1 (MQ=255)

                                                                                                                                                                                              agGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAAAGCCGCAATTCGCGGGGGTGTTCGTTCTGCGCTGGCCAAAGGACTGGGAGGAATGGGGATTTATGACGGCCATCTGGGGCTGTATGGAGACCGTATGGGACAGCTAGACCGTTACA‑GCGGGTCTGACATGATCAAACGAGGCGCGACGGTCCCCGGTTCCGCGCGGTTCCCgggaat                                          >  1:345701/1‑234 (MQ=255)

                                                                                                                                                                                                   gTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAAt                                                                                                                                                                      >  2:1101012/1‑110 (MQ=255)

                                                                                                                                                                                                   gTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTa                                                                                                          <  2:199849/170‑1 (MQ=255)

                                                                                                                                                                                                   gTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGGGAT‑GTGCCAGGGCTGAACGCCGCAATTGGAGGGGTTGGTTGTTCTGCGGTGGAAGAAGGTCTGGAAGTAAAGGGTGTTTATGACTGTTCTCTGGGTGTGTATGAAGACCGTATGGTAAATCTAGa                                                                                                                >  1:735097/1‑164 (MQ=255)

                                                                                                                                                                                                                  aTCGGTGTGTTGACAAGCGGC‑GGCGAT‑GCGCCAGGCATGAAAGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGGTTTTTTGACGGCTATCTGGGGCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACAAGATCAACCGGGGCGGTACGATTCCCCGTTCGTGGGGTTTCCCGGAAATCCGCTGCGAGaaaacaacccc                     >  1:417156/1‑230 (MQ=255)

                                                                                                                                                                                                                                     gc‑ggTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCCCGACGAGAAAATCCCCGCCTGGGCTATCGCAAACCTg  >  2:413205/1‑239 (MQ=255)

GGTCTATGCTTTCATAATCAGTATAAAAGAGAGCCAGACCCGCATTTTGTGT‑ATAAAATACCGCCATTTGGCCTGACCTGAATCAATTCAGCAGGAAGTGATTGTTATACTATTTGCACATTCGTTGGATCACTTCGATGTGCAAGAAGACTTCCGGCAACAGATTTCATTTTGCATTCCAAAGTTCAGAGGTAGTCATGATTAAGAAAATCGGTGTGTTGACAAGCGGC‑GGTGAT‑GCGCCAGGCATGAACGCCGCAATTCGCGGGGTTGTTCGTTCTGCGCTGACAGAAGGTCTGGAAGTAATGGGTATTTATGACGGCTATCTGGGTCTGTATGAAGACCGTATGGTACAGCTAGACCGTTACA‑GCGTGTCTGACATGATCAACCGTGGCGGTACGTTCCTCGGTTCTGCGCGTTTCCCGGAATTCCGCGACGAGAACATCCGCGCCGTGGCTATCGAAAACCTG  >  NC_000913/4107354‑4107820

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