Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A10 F13 I0 R1 150 9496.0 7621532 98.4% 7499587 136.2

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NZ_CP022686 2,502,553 A→C 41.2% I539S (ATC→AGC)  selB ← selenocysteine‑specific translation elongation factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0226862,502,5530AC41.2% 58.0 / 189.7 178I539S (ATC→AGC) selBselenocysteine‑specific translation elongation factor
Reads supporting (aligned to +/- strand):  ref base A (54/51);  new base C (37/36);  total (91/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 178 total aligned reads displayed.

TGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGA  >  NZ_CP022686/2502415‑2502688
                                                                                                                                          |                                                                                                                                       
tGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGAt                                                                                                                                        <  6:795955/140‑1 (MQ=255)
tGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGAt                                                                                                                                        >  6:888907/1‑140 (MQ=255)
tGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGAt                                                                                                                                        <  4:899663/140‑1 (MQ=255)
  aaTTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCg                                                                                                                                      <  7:163158/140‑1 (MQ=255)
   aTTGGCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACACTCCGATCGTTACGGTAATAACGATCTTTGACGCTCGc                                                                                                                                     <  5:175460/140‑1 (MQ=255)
    ttGCCAGCTTGCGGCCTACGACTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCa                                                                                                                                    <  1:845635/140‑1 (MQ=255)
     tGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAg                                                                                                                                   <  3:788667/140‑1 (MQ=255)
      gCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGt                                                                                                                                  >  1:112321/1‑140 (MQ=255)
        cAGCTTGCGGCCTACGCATAAGCGATCGCGGAAATCCGCCTCGCAGGGTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTaa                                                                                                                                <  4:233682/140‑1 (MQ=255)
         aGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAAt                                                                                                                               >  4:20359/1‑140 (MQ=255)
          gCTTGCGGCCTACGCCTAAGCGATCTCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACTCTCGCAGTAAtt                                                                                                                              >  1:889141/1‑140 (MQ=255)
            ttGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGAt                                                                                                                                        <  5:376342/128‑1 (MQ=255)
             tGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATc                                                                                                                           <  7:201213/140‑1 (MQ=255)
             tGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATc                                                                                                                           <  7:537114/140‑1 (MQ=255)
             tGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATc                                                                                                                           <  3:403364/140‑1 (MQ=255)
                     aCGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTg                                                                                                                   >  5:241275/1‑140 (MQ=255)
                     aCGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTg                                                                                                                   >  2:442685/1‑140 (MQ=255)
                     aCGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGATGAACCGCATTCCTGATCGAGATAGCGGATCAAATTGGCAAATACTACAATCCGATCGTTACGGTAATAACGATCTTAGACGATCGCAGAAATTATCCACTGtag                                                                                                                   >  5:4545/1‑138 (MQ=255)
                      cGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTgg                                                                                                                  >  2:393169/1‑140 (MQ=255)
                       gCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTggg                                                                                                                 <  3:763373/140‑1 (MQ=255)
                           aaGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGgccgc                                                                                                             <  5:306930/140‑1 (MQ=255)
                            aGCGCTCGCGGAAATCCGCCGCGCAGGTTGAACAGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTTACGCTCGCAGTAAGTATCCCCTGTTTGtccgcc                                                                                                            >  6:843361/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGCACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGCCAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            >  1:879019/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGCCGCTCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            >  2:843665/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            >  2:27017/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            >  7:44988/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            >  3:940036/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            >  3:669835/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            <  2:161003/140‑1 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGgccgcc                                                                                                            <  5:382203/140‑1 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGAAATTATCCCCTGTTTGgccgcc                                                                                                            >  8:889555/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGGAATAACGATCTTTGACGATCGCTGTAATTATCCCCTGTTGGgccgcc                                                                                                            >  5:2471/1‑140 (MQ=255)
                            aGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGAAAATCCGAACGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTAACCCCTGTTGGgccgcc                                                                                                            >  3:610363/1‑140 (MQ=255)
                                aTCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGc                                                                                                                                     <  8:332223/111‑1 (MQ=255)
                                aTCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGc                                                                                                        >  4:413299/1‑140 (MQ=255)
                                aTCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGc                                                                                                        >  7:357227/1‑140 (MQ=255)
                                aTCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGc                                                                                                        >  7:406011/1‑140 (MQ=255)
                                aTCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATc                                                                                                                           <  2:660611/121‑1 (MQ=255)
                                aTCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGc                                                                                                        >  8:457829/1‑140 (MQ=255)
                                 tCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCg                                                                                                       >  1:667610/1‑140 (MQ=255)
                                 tCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCg                                                                                                       <  2:197113/140‑1 (MQ=255)
                                 tCGCGGAAATCAGCCGCGCTGGTTGATCCGCATTCCTGATCTAGATCGCGGATCATATTGGCAAATTCTACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCg                                                                                                       >  8:779163/1‑140 (MQ=255)
                                  cgcgGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCg                                                                                                                                      >  7:145842/1‑108 (MQ=255)
                                  cgcgGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCg                                                                                                                                      <  8:145842/108‑1 (MQ=255)
                                  cgcgGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGt                                                                                                      >  6:475476/1‑140 (MQ=255)
                                  cgcgGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGt                                                                                                      >  3:279401/1‑140 (MQ=255)
                                       aaaTCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGt                                                                                                 <  1:699522/140‑1 (MQ=255)
                                        actCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTGGGGCCGCCTGGCGTAGCGTc                                                                                                <  6:238188/138‑1 (MQ=255)
                                         aTCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCTGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCa                                                                                               <  3:395172/140‑1 (MQ=255)
                                             gccgcGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGgcg                                                                                           <  8:304751/140‑1 (MQ=255)
                                              ccgcGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGgcgc                                                                                          >  2:623708/1‑140 (MQ=255)
                                               cgcgCAGGTTGAACCGCATTCCTGAGCGAGATCTCGGATCATATTGGCAAATTCGACACTCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCa                                                                                         <  5:847518/140‑1 (MQ=255)
                                                 cgcAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCAtt                                                                                       <  1:205857/140‑1 (MQ=255)
                                                      gTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTg                                                                                  <  3:736952/140‑1 (MQ=255)
                                                           aCCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCtcgt                                                                             <  3:135484/140‑1 (MQ=255)
                                                           aCCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCtcgt                                                                             >  1:926312/1‑140 (MQ=255)
                                                                aTTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGtt                                                                        >  6:656567/1‑140 (MQ=255)
                                                                aTTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTAGTCGGtt                                                                        >  5:715750/1‑140 (MQ=255)
                                                                    cTGATCGAGATCGCGGCTCATATTGGCAACTTCGACAATCCGAGCGTTACGGTCATAACGATCTTTGACGATCGCAGTAATTCTCCCCTGTAGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCg                                                                    <  1:600377/140‑1 (MQ=255)
                                                                    cTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGAAGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCg                                                                    <  5:12633/140‑1 (MQ=255)
                                                                       aTCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGAGCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGtcg                                                                 >  4:532627/1‑139 (MQ=255)
                                                                         cGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTctct                                                               >  3:163883/1‑140 (MQ=255)
                                                                            gaTCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGgccgc                                                                                                             <  4:714007/91‑1 (MQ=255)
                                                                             aTCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTg                                                           <  5:5700/140‑1 (MQ=255)
                                                                               cgcgGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGcc                                                         <  8:554949/140‑1 (MQ=255)
                                                                                 cgGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCGTTTGCCAg                                                       >  4:631451/1‑140 (MQ=255)
                                                                                      cATATTGGCAAATTCGACAATCCGATCGTTAAGGTAATAACGATCTTTGACGATCGCAGTAAATATCCCCTGTTGGGCCGCCTGGCGTAGCTTCAGGCGCATTGCCTGCTCGTCGGTGCCCGTGTATTTTGCCAGGTccc                                                  >  5:639119/1‑138 (MQ=255)
                                                                                         tttGGCAAATTCGACAATCCGATCGTTACGGTAATTACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGca                                               <  1:124948/139‑1 (MQ=255)
                                                                                         tgtGGCAAATTCGACAATACGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGTCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGca                                               <  3:348859/138‑1 (MQ=255)
                                                                                         atTGGCAAATTCGACAATCCGATCGTTACGTTAATAACTATCTTTGACGCTCGCAGTAATTATCCCCTGATGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGca                                               <  1:616129/140‑1 (MQ=255)
                                                                                         atTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGca                                               <  3:20359/140‑1 (MQ=255)
                                                                                         atTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGGCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGca                                               <  2:576284/140‑1 (MQ=255)
                                                                                          ttGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGcac                                              <  2:885648/140‑1 (MQ=255)
                                                                                           tGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGGCGGTTCCCGTCTCTTTTGCCAGGTCACGcacc                                             <  1:421421/140‑1 (MQ=255)
                                                                                              cAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGaacccac                                          >  1:680865/1‑140 (MQ=255)
                                                                                                 aTTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGACGGTTCCCGTCTCTTTTGCCAGGTCACGCACccaccaa                                       >  5:88386/1‑139 (MQ=255)
                                                                                                  ttCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTATGCCAGGTCACGCACCCACCATg                                      >  1:819417/1‑140 (MQ=255)
                                                                                                  ttCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTAGCGTTTGCGAGGTCACGCACCCACCATg                                      >  1:938837/1‑140 (MQ=255)
                                                                                                     gaaaaTCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGTTCACGCACCCACCATGGCt                                   <  2:744399/137‑1 (MQ=255)
                                                                                                     gACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCTTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCt                                   <  4:345388/140‑1 (MQ=255)
                                                                                                     gACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCt                                   <  4:794837/140‑1 (MQ=255)
                                                                                                     gACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCt                                   <  8:273869/140‑1 (MQ=255)
                                                                                                     gACAATCCGATCGTTACGGTAATAACGAGCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCt                                   <  4:898422/140‑1 (MQ=255)
                                                                                                        aaTCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGt                                <  5:403326/140‑1 (MQ=255)
                                                                                                             gATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCa                           >  1:778820/1‑140 (MQ=255)
                                                                                                             gATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCGCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCa                           >  4:670810/1‑140 (MQ=255)
                                                                                                                  ttACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAg                      >  6:5682/1‑140 (MQ=255)
                                                                                                                  ttACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAg                      <  7:469640/140‑1 (MQ=255)
                                                                                                                      ggTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCATGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCCCCAAACAGCGGc                  >  3:575687/1‑140 (MQ=255)
                                                                                                                        taataaCGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGctc                <  3:389697/140‑1 (MQ=255)
                                                                                                                          ataaCGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACCGCGGCGCTg              >  7:65266/1‑140 (MQ=255)
                                                                                                                          ataaCGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCCCCAAACAGCGGCTCTg              >  4:427352/1‑140 (MQ=255)
                                                                                                                          ataaCGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTg              >  5:229335/1‑140 (MQ=255)
                                                                                                                            aaCGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCt            <  6:229335/140‑1 (MQ=255)
                                                                                                                                aTCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCACGCGCATTGCCTGCTCGACGGTTCCCGTCACAGTTGCCAGGGCACGCACCCACCATCGCTCGGCACCACACAGCAGCTCTGCttttt        >  8:84322/1‑140 (MQ=255)
                                                                                                                                aTCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTATGCtgttt        >  4:832454/1‑140 (MQ=255)
                                                                                                                                  cTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTgg      >  5:825238/1‑139 (MQ=255)
                                                                                                                                      gACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGa  <  2:907708/140‑1 (MQ=255)
                                                                                                                                      gACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCGCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGa  <  6:324940/140‑1 (MQ=255)
                                                                                                                                      gACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGGCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGa  <  7:779163/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGA  >  NZ_CP022686/2502415‑2502688

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATTCCAGAATCTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCCTGCTGCTC  >  NZ_CP022686/2502404‑2502701
                                                                                                                                                     |                                                                                                                                                    
ATTCCAGAATCTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCT                                                                                                                                                     >  NB501016:372:H5LTGBGXK:1:21302:16500:11014/1‑151 (MQ=60)
ATTCCAGAATCTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGAT                                                                                                                                                     >  NB501016:372:H5LTGBGXK:1:22103:9825:9934/1‑151 (MQ=60)
   CCAGAATCTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGC                                                                                                                                                  >  NB501016:372:H5LTGBGXK:3:12507:21839:7868/1‑151 (MQ=60)
      GAATCTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGT                                                                                                                                               >  NB501016:372:H5LTGBGXK:1:11301:3126:5567/1‑151 (MQ=60)
          CTGAATTGCCAGCTTGCGGCCTCCGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATT                                                                                                                                           >  NB501016:372:H5LTGBGXK:1:23212:6495:13969/1‑151 (MQ=60)
          CTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCTCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACTCTCGCAGTAATT                                                                                                                                           >  NB501016:372:H5LTGBGXK:1:23202:21976:3918/1‑151 (MQ=60)
          CTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGTTCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATT                                                                                                                                           >  NB501016:372:H5LTGBGXK:2:13205:16076:11919/1‑151 (MQ=60)
          CTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATT                                                                                                                                           >  NB501016:372:H5LTGBGXK:1:22106:7841:14629/1‑151 (MQ=60)
           TTAATTGCCAGCTTGCGTCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCTGTTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGTCAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTA                                                                                                                                          <  NB501016:372:H5LTGBGXK:2:11204:22131:18702/151‑1 (MQ=60)
           TGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTA                                                                                                                                          <  NB501016:372:H5LTGBGXK:1:21111:13842:7034/151‑1 (MQ=60)
           TGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTA                                                                                                                                          <  NB501016:372:H5LTGBGXK:2:11204:22122:18697/151‑1 (MQ=60)
           GGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTA                                                                                                                                          <  NB501016:372:H5LTGBGXK:3:13504:2758:2890/151‑1 (MQ=60)
            GAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTAT                                                                                                                                         >  NB501016:372:H5LTGBGXK:3:13602:18795:19961/1‑151 (MQ=60)
             AATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATC                                                                                                                                        <  NB501016:372:H5LTGBGXK:4:12403:2045:16286/151‑1 (MQ=60)
             AATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATC                                                                                                                                        <  NB501016:372:H5LTGBGXK:2:12102:23833:11322/151‑1 (MQ=60)
              ATTGGCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACACTCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCC                                                                                                                                       <  NB501016:372:H5LTGBGXK:3:12405:24865:8557/151‑1 (MQ=60)
              ATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCC                                                                                                                                       <  NB501016:372:H5LTGBGXK:3:23505:17393:13415/151‑1 (MQ=60)
              ATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCC                                                                                                                                       <  NB501016:372:H5LTGBGXK:3:12408:1762:8540/151‑1 (MQ=60)
               TTGCCAGCTTGCGGCCTACGACTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCC                                                                                                                                      <  NB501016:372:H5LTGBGXK:1:23104:3302:13559/151‑1 (MQ=60)
                TGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCC                                                                                                                                     <  NB501016:372:H5LTGBGXK:2:22308:22289:10855/151‑1 (MQ=60)
                     GCTTGCGGCCTACTCCTAAGCGATCGCGGAAATCCGCCGCGCATGTTGAACCTCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATACGATCGTTACTGTAATAACGATCTTTGACTCTCGCAGTAATTATCCCCTGTTG                                                                                                                                <  NB501016:372:H5LTGBGXK:4:21409:21758:2961/151‑1 (MQ=60)
                     GCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTTCGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTG                                                                                                                                >  NB501016:372:H5LTGBGXK:3:21406:23800:9042/1‑151 (MQ=60)
                     GCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTG                                                                                                                                >  NB501016:372:H5LTGBGXK:3:12509:6668:4817/1‑151 (MQ=60)
                     GCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGATGAACCGCATTCCTGATCGAGATAGCGGATCAAATTGGCAAATACTACAATCCGATCGTTACGGTAATAACGATCTTAGACGATCGCAGAAATTATCCACTGTAG                                                                                                                                >  NB501016:372:H5LTGBGXK:3:11402:9366:1939/1‑151 (MQ=60)
                        TGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGC                                                                                                                             <  NB501016:372:H5LTGBGXK:4:12412:7227:15532/151‑1 (MQ=60)
                        TGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGC                                                                                                                             <  NB501016:372:H5LTGBGXK:4:21508:11564:15056/151‑1 (MQ=60)
                        TGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGC                                                                                                                             <  NB501016:372:H5LTGBGXK:2:13206:14739:18488/151‑1 (MQ=60)
                        TGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGC                                                                                                                             <  NB501016:372:H5LTGBGXK:2:23108:10785:15298/151‑1 (MQ=60)
                        GGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGC                                                                                                                             <  NB501016:372:H5LTGBGXK:2:22209:4846:17844/151‑1 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGCACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGCCAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:1:23111:5034:20003/1‑151 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:2:23306:22638:7627/1‑151 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:3:11508:22579:19227/1‑151 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:4:11412:19175:2506/1‑151 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:2:22102:19083:6569/1‑151 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:2:22105:14766:13164/1‑151 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGGAATAACGATCTTTGACGATCGCTGTAATTATCCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:3:11401:15624:12002/1‑151 (MQ=60)
                            GCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGAAAATCCGAACGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTAACCCCTGTTGGGCCGCC                                                                                                                         >  NB501016:372:H5LTGBGXK:2:21304:15195:17120/1‑151 (MQ=60)
                                ACGCCTAAGCGCTCGCGGAAATCCGCCGCGCAGGTTGAACAGCTTTCGTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGCTCGTTACGGGAATAACGATCTTTGACGATCGCAGTATTTATCCCCTGTTGGGACGGCTGGC                                                                                                                     >  NB501016:372:H5LTGBGXK:3:12606:4720:4788/1‑151 (MQ=60)
                                ACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGCCAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGC                                                                                                                     >  NB501016:372:H5LTGBGXK:4:11502:10483:13293/1‑151 (MQ=60)
                                ACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCtgtctcttatacacatc                                                                                                                     >  NB501016:372:H5LTGBGXK:4:13407:10888:18872/1‑134 (MQ=60)
                                ACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGC                                                                                                                     >  NB501016:372:H5LTGBGXK:1:12208:25429:8264/1‑151 (MQ=60)
                                ACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGC                                                                                                                     >  NB501016:372:H5LTGBGXK:4:13501:16728:12772/1‑151 (MQ=60)
                                ACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGC                                                                                                                     >  NB501016:372:H5LTGBGXK:4:13512:24741:17074/1‑151 (MQ=60)
                                ACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCtgtctct                                                                                                                     >  NB501016:372:H5LTGBGXK:1:21312:22624:15285/1‑144 (MQ=60)
                                 CGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCG                                                                                                                    >  NB501016:372:H5LTGBGXK:1:22102:17062:2378/1‑151 (MQ=60)
                                  GCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGT                                                                                                                   >  NB501016:372:H5LTGBGXK:2:12303:13997:7080/1‑151 (MQ=60)
                                  GCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGT                                                                                                                   <  NB501016:372:H5LTGBGXK:2:22303:10194:4082/151‑1 (MQ=60)
                                  GCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATACGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCtgtctcttatacacatctcc                                                                                                                   >  NB501016:372:H5LTGBGXK:4:11611:26095:10837/1‑131 (MQ=60)
                                      AAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCG                                                                                                               <  NB501016:372:H5LTGBGXK:3:12610:5002:16899/151‑1 (MQ=60)
                                      AAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCG                                                                                                               <  NB501016:372:H5LTGBGXK:1:22312:7877:18643/151‑1 (MQ=60)
                                       AGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGT                                                                                                              <  NB501016:372:H5LTGBGXK:3:22401:10316:9883/151‑1 (MQ=60)
                                       AGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGT                                                                                                              <  NB501016:372:H5LTGBGXK:3:13505:9595:8853/151‑1 (MQ=60)
                                                  AAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTG                                                                                                   <  NB501016:372:H5LTGBGXK:1:22108:5989:10140/151‑1 (MQ=60)
                                                  AAATCCGCCGCGCAGGTGGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTG                                                                                                   <  NB501016:372:H5LTGBGXK:3:13501:17143:7502/151‑1 (MQ=60)
                                                    ATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCTGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCC                                                                                                 <  NB501016:372:H5LTGBGXK:2:13204:26523:15587/151‑1 (MQ=60)
                                                    ATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCC                                                                                                 <  NB501016:372:H5LTGBGXK:4:22501:19263:4348/151‑1 (MQ=60)
                                                    ATCCGCAGCGCAGGTTGCACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCC                                                                                                 <  NB501016:372:H5LTGBGXK:2:13203:1067:15600/151‑1 (MQ=60)
                                                          CGCGGAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATGGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCG                                                                                           <  NB501016:372:H5LTGBGXK:3:11408:15141:18093/151‑1 (MQ=60)
                                                          CGCGCAGGTTGAACCGCATTCCTGAGCGAGATCTCGGATCATATTGGCAAATTCGACACTCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCG                                                                                           <  NB501016:372:H5LTGBGXK:3:23505:6653:6541/151‑1 (MQ=60)
                                                           GCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGT                                                                                          >  NB501016:372:H5LTGBGXK:1:23210:6624:7937/1‑151 (MQ=60)
                                                            CGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTC                                                                                         <  NB501016:372:H5LTGBGXK:1:12109:22088:4841/151‑1 (MQ=60)
                                                                GGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTAGTCGGTT                                                                                     >  NB501016:372:H5LTGBGXK:3:22510:1289:12980/1‑151 (MQ=60)
                                                                 GTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTC                                                                                    <  NB501016:372:H5LTGBGXK:2:22209:19360:1291/151‑1 (MQ=60)
                                                                      ACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTC                                                                               <  NB501016:372:H5LTGBGXK:2:11307:24176:8439/151‑1 (MQ=60)
                                                                         GCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCT                                                                            >  NB501016:372:H5LTGBGXK:2:12101:13851:2109/1‑151 (MQ=60)
                                                                            TTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGctgtctcttatacacatctccgagcccacgagactcaa                                                                         >  NB501016:372:H5LTGBGXK:2:22203:9205:14834/1‑113 (MQ=60)
                                                                             TCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTT                                                                        >  NB501016:372:H5LTGBGXK:2:12111:15687:3295/1‑151 (MQ=60)
                                                                              CCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTATGACGCTCGCAGTAATTATCCCATGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGC                                                                       >  NB501016:372:H5LTGBGXK:4:13604:25427:9891/1‑151 (MQ=60)
                                                                               CTGATCGAGATCGCGGCTCATATTGGCAACTTCGACAATCCGAGCGTTACGGTCATAACGATCTTTGACGATCGCAGTAATTCTCCCCTGTAGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCC                                                                      <  NB501016:372:H5LTGBGXK:1:21211:2875:11595/151‑1 (MQ=60)
                                                                               CTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGAAGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCC                                                                      <  NB501016:372:H5LTGBGXK:3:11404:10939:4877/151‑1 (MQ=60)
                                                                                  ATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGG                                                                   <  NB501016:372:H5LTGBGXK:1:23109:1715:11357/151‑1 (MQ=60)
                                                                                      AGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTAAGGTAATAACGATCTTTGACGATCGCAGTAAATATCCCCTGTTGGGCCGCCTGGCGTAGCTTCAGGCGCATTGCCTGCTCGTCGGTGCCCGTGTATTTTGCCAGGTCCC                                                               >  NB501016:372:H5LTGBGXK:3:22404:7129:4671/1‑151 (MQ=60)
                                                                                        ATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGC                                                             <  NB501016:372:H5LTGBGXK:3:11402:8840:7778/151‑1 (MQ=60)
                                                                                          CGCGGATCATATTGGCAAATTCTACAATCCGATCGTTACGGTTATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCAC                                                           <  NB501016:372:H5LTGBGXK:3:13504:18686:18052/151‑1 (MQ=60)
                                                                                              GATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGAACCCAC                                                       >  NB501016:372:H5LTGBGXK:1:22104:24541:13419/1‑151 (MQ=60)
                                                                                              GATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCAC                                                       <  NB501016:372:H5LTGBGXK:4:13508:16416:5056/151‑1 (MQ=60)
                                                                                                 CATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGACGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCAA                                                    >  NB501016:372:H5LTGBGXK:3:11510:23677:5289/1‑151 (MQ=60)
                                                                                                  ATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTATGCCAGGTCACGCACCCACCATG                                                   >  NB501016:372:H5LTGBGXK:1:22310:23812:20246/1‑151 (MQ=60)
                                                                                                  ATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTAGCGTTTGCGAGGTCACGCACCCACCATG                                                   >  NB501016:372:H5LTGBGXK:1:23301:15706:3126/1‑151 (MQ=60)
                                                                                                   TATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGG                                                  <  NB501016:372:H5LTGBGXK:3:22408:6554:13074/151‑1 (MQ=60)
                                                                                                    tgTGGCAAATTCGACAATACGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGTCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGC                                                 <  NB501016:372:H5LTGBGXK:2:13106:25844:6157/149‑1 (MQ=60)
                                                                                                    TTTGGCAAATTCGACAATCCGATCGTTACGGTAATTACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGC                                                 <  NB501016:372:H5LTGBGXK:1:11303:15779:20039/151‑1 (MQ=60)
                                                                                                    ATTGGCAAATTCGACAATCCGATCGTTACGTTAATAACTATCTTTGACGCTCGCAGTAATTATCCCCTGATGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGC                                                 <  NB501016:372:H5LTGBGXK:1:21303:13555:1490/151‑1 (MQ=60)
                                                                                                    ATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGC                                                 <  NB501016:372:H5LTGBGXK:2:11105:6547:12606/151‑1 (MQ=60)
                                                                                                    ATTGGCAAATTAGACAATCCGATCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGC                                                 <  NB501016:372:H5LTGBGXK:2:13208:14615:5890/151‑1 (MQ=60)
                                                                                                      TGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGGCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTC                                               <  NB501016:372:H5LTGBGXK:1:13207:4087:18974/151‑1 (MQ=60)
                                                                                                             TTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCA                                        >  NB501016:372:H5LTGBGXK:1:22301:25694:19277/1‑151 (MQ=60)
                                                                                                                 ACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACC                                    <  NB501016:372:H5LTGBGXK:1:12104:6677:13727/151‑1 (MQ=60)
                                                                                                                   AATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGC                                  <  NB501016:372:H5LTGBGXK:3:13510:17714:5249/151‑1 (MQ=60)
                                                                                                                   AATCCGATCGTTACGGTAATAACGATCTCTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGC                                  <  NB501016:372:H5LTGBGXK:2:11310:14078:14963/151‑1 (MQ=60)
                                                                                                                    ATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCG                                 >  NB501016:372:H5LTGBGXK:1:22210:10129:6304/1‑151 (MQ=60)
                                                                                                                      CCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCATGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCCCCAAACAGCGGC                               >  NB501016:372:H5LTGBGXK:2:21209:23131:5475/1‑151 (MQ=60)
                                                                                                                          TCGTTACGGTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACCGCGGCGCTG                           >  NB501016:372:H5LTGBGXK:4:11504:20232:15463/1‑151 (MQ=60)
                                                                                                                          TCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTG                           >  NB501016:372:H5LTGBGXK:3:12506:10114:10579/1‑151 (MQ=60)
                                                                                                                             TTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTT                        <  NB501016:372:H5LTGBGXK:4:21403:9591:19319/151‑1 (MQ=60)
                                                                                                                             TTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCGGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCGCGTCGGGTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTT                        <  NB501016:372:H5LTGBGXK:3:11405:11818:3499/151‑1 (MQ=60)
                                                                                                                             TTACGGTAATAACGATCTTTGACGATCGAAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCGCGTCGGTTCCCTTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACCGCGGCTCTGCTT                        <  NB501016:372:H5LTGBGXK:1:13311:5575:5250/151‑1 (MQ=60)
                                                                                                                                  GTAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGC                   >  NB501016:372:H5LTGBGXK:4:12411:7900:1967/1‑151 (MQ=60)
                                                                                                                                  GTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGG                   >  NB501016:372:H5LTGBGXK:3:23412:17040:4692/1‑151 (MQ=60)
                                                                                                                                   TAATAACGATCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCC                  <  NB501016:372:H5LTGBGXK:2:13203:11812:10473/151‑1 (MQ=60)
                                                                                                                                          GATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGC           <  NB501016:372:H5LTGBGXK:1:23201:19166:12636/151‑1 (MQ=60)
                                                                                                                                           ATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCC          >  NB501016:372:H5LTGBGXK:3:21610:5393:9735/1‑151 (MQ=60)
                                                                                                                                            TCTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCGCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCCT         >  NB501016:372:H5LTGBGXK:4:23411:2049:1140/1‑151 (MQ=60)
                                                                                                                                            TCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCCCCCACCATGGCTCGTCACCAAACAGCGGCTCTGCGTTTTGCCAGATGGCCT         >  NB501016:372:H5LTGBGXK:4:12409:12695:9221/1‑151 (MQ=60)
                                                                                                                                             CTTTGACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCCTG        >  NB501016:372:H5LTGBGXK:4:11403:10654:16139/1‑151 (MQ=60)
                                                                                                                                                 GACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGGCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCCTGCTGC    <  NB501016:372:H5LTGBGXK:2:13109:10493:10296/151‑1 (MQ=60)
                                                                                                                                                 GACGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGGCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCCTGCTGC    <  NB501016:372:H5LTGBGXK:4:23406:4952:19218/151‑1 (MQ=60)
                                                                                                                                                   CGCTCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGGCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCCTGCTGCTC  <  NB501016:372:H5LTGBGXK:2:12109:17617:10248/151‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
ATTCCAGAATCTGAATTGCCAGCTTGCGGCCTACGCCTAAGCGATCGCGGAAATCCGCCGCGCAGGTTGAACCGCATTCCTGATCGAGATCGCGGATCATATTGGCAAATTCGACAATCCGATCGTTACGGTAATAACGATCTTTGACGATCGCAGTAATTATCCCCTGTTGGGCCGCCTGGCGTAGCGTCAGGCGCATTGCCTGCTCGTCGGTTCCCGTCTCTTTTGCCAGGTCACGCACCCACCATGGCTCGTCACCAAACAGCGGCTCTGCTTTTTGCCAGATGGCCTGCTGCTC  >  NZ_CP022686/2502404‑2502701

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑