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
A1 F35 I0 R1 338 60.8 3106128 78.2% 2428992 142.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,780 A→G 100% E244E (GAA→GAG stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,7800AG98.2% 6753.9 / 32.6 2377E244E (GAA→GAGstfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base A (20/25);  new base G (1131/1194);  total (1152/1225)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2411 total aligned reads displayed.

GCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCG  >  NC_000913/1429634‑1429927
                                                                                                                                                  |                                                                                                                                                   
 cggcggCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGca                                                                                                                                                  <  5:355011/149‑2 (MQ=255)
      gcCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCaaa                                                                                                                                             <  2:386287/149‑1 (MQ=255)
        caccaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAat                                                                                                                                           >  3:41850/1‑149 (MQ=255)
        caccaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCGCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGcag                                                                                                                                                 >  6:261449/1‑141 (MQ=255)
         accaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatc                                                                                                                                          <  6:166791/149‑1 (MQ=255)
         accaGTGCCGGTGCGGCGAAAACGGCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatc                                                                                                                                          <  5:325517/149‑1 (MQ=255)
          ccaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatca                                                                                                                                         >  1:100362/1‑149 (MQ=255)
           caGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatcat                                                                                                                                        <  4:294813/149‑1 (MQ=255)
             gTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatcat                                                                                                                                        >  8:249164/1‑147 (MQ=255)
               gCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAg                                                                                                                                     >  8:356826/1‑148 (MQ=255)
                ccGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGaa                                                                                                                                   <  7:241756/149‑1 (MQ=255)
                    tGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatcatc                                                                                                                                       >  8:103990/1‑141 (MQ=255)
                        gcgAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAaacgaacg                                                                                                                            >  1:337816/1‑148 (MQ=255)
                        gcgAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGc                                                                                                                           <  1:18222/149‑1 (MQ=255)
                         cgAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGca                                                                                                                          <  7:74354/149‑1 (MQ=255)
                          gAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGCGATGCGGTGGCCTCAAAAGAGGCAGCAAAATAATCAGAAACGACCGcat                                                                                                                         >  2:108903/1‑149 (MQ=255)
                           aaaaCGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCCCCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGcat                                                                                                                         >  6:315575/1‑148 (MQ=255)
                           aaaaCGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGcatc                                                                                                                        <  4:51565/149‑1 (MQ=255)
                            aaaCGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGcatca                                                                                                                       <  3:352177/149‑1 (MQ=255)
                             aaCGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGcatcat                                                                                                                      <  7:326127/149‑1 (MQ=255)
                                gTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCaa                                                                                                                   <  1:61454/149‑1 (MQ=255)
                                 tCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAg                                                                                                                  >  1:310095/1‑149 (MQ=255)
                                    gAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAgtgc                                                                                                               <  8:278617/149‑1 (MQ=255)
                                    gAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAgtgc                                                                                                               <  3:164609/149‑1 (MQ=255)
                                    gAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCAGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAaacgaac                                                                                                                             >  8:165525/1‑135 (MQ=255)
                                     aaaCGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGcatca                                                                                                                       >  6:287378/1‑140 (MQ=255)
                                     aaaCGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGcc                                                                                                              >  3:283702/1‑149 (MQ=255)
                                     aaaCGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGcc                                                                                                              >  2:59479/1‑149 (MQ=255)
                                      aaCGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAaacgaac                                                                                                                             >  4:209923/1‑133 (MQ=255)
                                       aCGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGccgg                                                                                                            <  1:149399/149‑3 (MQ=255)
                                              ctgcAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGc                                                                                                      <  5:250566/145‑1 (MQ=255)
                                              ctgcAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGc                                                                                                      <  2:310095/145‑1 (MQ=255)
                                                  aGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCtt                                                                                                 >  5:81170/1‑149 (MQ=255)
                                                      tCACAACCATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGCAAGCGTCAGAGGCCGCCACTTCAGCACGAGCTGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGGGCAGCTTCCTc                                                                                             <  5:226698/149‑1 (MQ=255)
                                                      tCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTc                                                                                             <  4:145410/149‑1 (MQ=255)
                                                          aaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGc                                                                                          >  8:377012/2‑148 (MQ=255)
                                                           aCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGcca                                                                                        >  5:111801/1‑147 (MQ=255)
                                                             aaTCAGCCGCCACGTCTGCCTCCACCGCGGCCTCGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcg                                                                                      >  6:141381/1‑149 (MQ=255)
                                                             aaTCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcg                                                                                      >  4:265055/1‑149 (MQ=255)
                                                               tCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggc                                                                                    >  7:343563/1‑149 (MQ=255)
                                                               tCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggc                                                                                    >  2:287505/1‑149 (MQ=255)
                                                                  gccgcCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCaa                                                                                        <  8:310828/139‑1 (MQ=255)
                                                                    cgcCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGCACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggcgg                                                                                  >  6:97570/2‑146 (MQ=255)
                                                                          gtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAAtt                                                                         >  3:132516/2‑149 (MQ=255)
                                                                          gtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAAtt                                                                         <  2:36241/148‑1 (MQ=255)
                                                                          gtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTCGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCACCGGCGGCAGAAAAtt                                                                         >  5:172840/2‑149 (MQ=255)
                                                                            cTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcgg                                                                                     <  2:162808/135‑1 (MQ=255)
                                                                            cTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTct                                                                       >  6:224104/1‑148 (MQ=255)
                                                                              gcctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                                                                     >  7:272998/4‑149 (MQ=255)
                                                                               cctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                                                                     <  5:149827/146‑1 (MQ=255)
                                                                               cctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                                                                     <  5:38875/146‑1 (MQ=255)
                                                                                ctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAaacgaac                                                                                                                             >  1:390906/2‑91 (MQ=255)
                                                                                ctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTct                                                                       >  4:317424/2‑144 (MQ=255)
                                                                                ctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgcca                                                                   <  1:8626/148‑3 (MQ=255)
                                                                                 tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                  >  5:304409/1‑146 (MQ=255)
                                                                                  ccACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagg                                                                 >  6:191181/1‑145 (MQ=255)
                                                                                   cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACAGCGGCAGAAAATTCtgccaggg                                                                <  4:173303/149‑6 (MQ=255)
                                                                                    aCCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggc                                                               >  1:279253/1‑149 (MQ=255)
                                                                                             aCGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcggc                                                            >  3:84242/1‑143 (MQ=255)
                                                                                               gAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGACAATTCTGCCAGGGCGGCAAAAACGt                                                    >  6:110062/1‑149 (MQ=255)
                                                                                               gAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGt                                                    >  8:47738/1‑149 (MQ=255)
                                                                                                aaagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTc                                                   <  7:191643/143‑1 (MQ=255)
                                                                                                 aagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTcc                                                  <  7:128793/144‑1 (MQ=255)
                                                                                                 aagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGGcc                                                  <  1:126271/144‑1 (MQ=255)
                                                                                                  agcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCg                                                 <  7:300512/145‑1 (MQ=255)
                                                                                                      tCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGAc                                             >  5:203579/1‑149 (MQ=255)
                                                                                                                       agcacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGACCGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGgcatct                             <  1:363185/141‑1 (MQ=255)
                                                                                                                       agcacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAACACGTCCGAGACGCAGGCCAGGtcatct                             >  4:348160/8‑148 (MQ=255)
                                                                                                                       agcacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGtcatct                             <  2:902/141‑1 (MQ=255)
                                                                                                                       agcacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGtcatct                             >  8:133247/8‑148 (MQ=255)
                                                                                                                        gcacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGa                           <  6:334025/143‑1 (MQ=255)
                                                                                                                         cacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCACCGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGa                           >  3:309334/6‑148 (MQ=255)
                                                                                                                         cacgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGa                           >  8:336896/6‑148 (MQ=255)
                                                                                                                          acgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                         <  7:39918/145‑1 (MQ=255)
                                                                                                                          acgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                         >  7:15579/5‑149 (MQ=255)
                                                                                                                          acgagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCAACAGAAACGAACGCATCTTCTAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAACTTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                         <  4:235066/145‑1 (MQ=255)
                                                                                                                            gagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGaat                                         <  2:216208/129‑2 (MQ=255)
                                                                                                                            gagATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGaat                                         >  1:216208/3‑130 (MQ=255)
                                                                                                                             ggATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACgg                      >  6:48347/1‑149 (MQ=255)
                                                                                                                                 gcggtgGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACagcagc                  >  2:340615/1‑149 (MQ=255)
                                                                                                                                 gcggtgGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACagcagc                  >  5:279790/1‑149 (MQ=255)
                                                                                                                                     tggcCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaa              >  2:196512/2‑147 (MQ=255)
                                                                                                                                       gcCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg            >  8:65848/1‑148 (MQ=255)
                                                                                                                                       gcCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg            <  3:225626/149‑2 (MQ=255)
                                                                                                                                        ccTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg            >  8:19372/1‑147 (MQ=255)
                                                                                                                                         cTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctg                                                                      >  2:255841/1‑87 (MQ=255)
                                                                                                                                         cTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTcc                                                  >  7:37982/1‑109 (MQ=255)
                                                                                                                                         cTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg            >  5:55403/1‑146 (MQ=255)
                                                                                                                                         cTCAAAAGAGGCAGCAAAAGCATCAGCAACGAACGCATCAACAAGTGCCGGTCGTGCAGCTTCCTCGCCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACgga          <  4:89130/149‑1 (MQ=255)
                                                                                                                                          tCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACggag         >  3:255823/1‑149 (MQ=255)
                                                                                                                                           cAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAACTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAAAACTGCCAGGGCGGAAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgc        >  3:329835/1‑149 (MQ=255)
                                                                                                                                            aaaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGATACAGCAGCGGAACGGAgcg       >  1:252806/1‑149 (MQ=255)
                                                                                                                                            aaaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcg       >  3:19133/1‑149 (MQ=255)
                                                                                                                                            aaaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGAAGCGGAACggag         <  6:77251/147‑1 (MQ=255)
                                                                                                                                             aaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcgc      <  7:221927/149‑1 (MQ=255)
                                                                                                                                             aaaGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGGCCGAGACGCATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcgc      >  2:131315/1‑149 (MQ=255)
                                                                                                                                              aagaggCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACggag         >  4:110368/1‑145 (MQ=255)
                                                                                                                                              aagaggCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGcc     <  8:228735/149‑1 (MQ=255)
                                                                                                                                              aagaggCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGcc     >  6:362981/1‑149 (MQ=255)
                                                                                                                                              aagaggCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGcc     >  3:107264/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCG  >  NC_000913/1429634‑1429927

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ 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 Only 200 of 1365 total aligned reads displayed.
AAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGG  >  NC_000913/1429630‑1429928
                                                                                                                                                      |                                                                                                                                                    
  ACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGC                                                                                                                                                    <  NB501016:230:HFVJLAFXY:4:11507:11732:7677/151‑1 (MQ=60)
  ACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGC                                                                                                                                                    <  NB501016:230:HFVJLAFXY:2:11104:19351:17218/151‑1 (MQ=60)
  ACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGC                                                                                                                                                    <  NB501016:230:HFVJLAFXY:1:21110:3792:15058/151‑1 (MQ=60)
        CGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAA                                                                                                                                               >  NB501016:230:HFVJLAFXY:3:11409:13754:12078/1‑150 (MQ=60)
        CGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAA                                                                                                                                               >  NB501016:230:HFVJLAFXY:3:21512:18288:5504/1‑150 (MQ=60)
         GGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACACCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAA                                                                                                                                              >  NB501016:230:HFVJLAFXY:4:21408:26736:1270/1‑150 (MQ=50)
         GGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAA                                                                                                                                             <  NB501016:230:HFVJLAFXY:4:21604:7087:18508/151‑1 (MQ=58)
          GCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAAT                                                                                                                                            <  NB501016:230:HFVJLAFXY:2:11204:22799:14420/151‑1 (MQ=56)
           CCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATC                                                                                                                                           >  NB501016:230:HFVJLAFXY:2:11312:3814:15928/1‑151 (MQ=54)
            CACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCA                                                                                                                                          <  NB501016:230:HFVJLAFXY:1:11207:18240:6058/151‑1 (MQ=54)
            CACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATC                                                                                                                                           >  NB501016:230:HFVJLAFXY:1:21307:14465:4830/1‑150 (MQ=53)
             ACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCAT                                                                                                                                         >  NB501016:230:HFVJLAFXY:1:11206:13388:10386/1‑151 (MQ=54)
              CCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATC                                                                                                                                        >  NB501016:230:HFVJLAFXY:1:21206:26557:14305/1‑151 (MQ=54)
               CAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCA                                                                                                                                       <  NB501016:230:HFVJLAFXY:4:21402:6259:1890/151‑1 (MQ=54)
               CAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCA                                                                                                                                       <  NB501016:230:HFVJLAFXY:3:11606:23645:6234/151‑1 (MQ=54)
               CAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCA                                                                                                                                       <  NB501016:230:HFVJLAFXY:4:11605:1451:9752/151‑1 (MQ=54)
                AGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAG                                                                                                                                      <  NB501016:230:HFVJLAFXY:2:11303:3033:3385/151‑1 (MQ=60)
                AGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAG                                                                                                                                      <  NB501016:230:HFVJLAFXY:2:11301:19208:13602/151‑1 (MQ=54)
                AGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAG                                                                                                                                      <  NB501016:230:HFVJLAFXY:1:21107:14102:3700/151‑1 (MQ=54)
                 GTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATC                                                                                                                                        >  NB501016:230:HFVJLAFXY:4:21507:9830:19013/1‑148 (MQ=50)
                   GCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCCGCAAAATCATCAGAA                                                                                                                                    >  NB501016:230:HFVJLAFXY:4:11505:8630:20305/1‑150 (MQ=53)
                   GCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAA                                                                                                                                   <  NB501016:230:HFVJLAFXY:4:21407:20103:12106/151‑1 (MQ=54)
                   GCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAA                                                                                                                                    <  NB501016:230:HFVJLAFXY:4:11502:3736:2882/150‑1 (MQ=53)
                   GCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGC                                                                                                                                                    <  NB501016:230:HFVJLAFXY:3:11410:18722:2558/134‑1 (MQ=46)
                    CCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAAC                                                                                                                                  <  NB501016:230:HFVJLAFXY:2:11207:15491:1115/151‑1 (MQ=54)
                         GCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACG                                                                                                                             <  NB501016:230:HFVJLAFXY:4:11409:11411:19104/151‑1 (MQ=54)
                         GCGGCGAAAACGTCAGAAACGAATGCCTCAGCGTCACAACAATCAGCCGCCAAGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACG                                                                                                                             <  NB501016:230:HFVJLAFXY:1:21201:7733:3240/151‑1 (MQ=47)
                          CGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGC                                                                                                                            <  NB501016:230:HFVJLAFXY:2:21106:14968:9535/151‑1 (MQ=54)
                           GGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGC                                                                                                                            >  NB501016:230:HFVJLAFXY:1:11208:11506:10126/1‑150 (MQ=53)
                           GGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGC                                                                                                                            >  NB501016:230:HFVJLAFXY:2:21103:13604:14399/1‑150 (MQ=53)
                                   CGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAG                                                                                                                   <  NB501016:230:HFVJLAFXY:4:11409:7638:15669/151‑1 (MQ=54)
                                   CGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCA                                                                                                                     >  NB501016:230:HFVJLAFXY:4:11406:18829:13831/1‑149 (MQ=52)
                                   CGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCA                                                                                                                     >  NB501016:230:HFVJLAFXY:4:21601:10135:18256/1‑149 (MQ=52)
                                      CAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGTCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTG                                                                                                                 <  NB501016:230:HFVJLAFXY:4:21407:4731:19211/150‑1 (MQ=42)
                                      CAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGC                                                                                                                >  NB501016:230:HFVJLAFXY:3:11602:14498:13924/1‑151 (MQ=50)
                                      CAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGC                                                                                                                >  NB501016:230:HFVJLAFXY:3:21602:23723:5248/1‑151 (MQ=50)
                                      CAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGC                                                                                                                <  NB501016:230:HFVJLAFXY:3:21510:22005:12954/151‑1 (MQ=50)
                                      CAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGC                                                                                                                <  NB501016:230:HFVJLAFXY:4:11510:16058:9604/151‑1 (MQ=50)
                                       AGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCC                                                                                                               >  NB501016:230:HFVJLAFXY:2:11108:1762:6075/1‑151 (MQ=48)
                                       AGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCC                                                                                                               >  NB501016:230:HFVJLAFXY:4:21510:13348:13878/1‑151 (MQ=48)
                                        GAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGtgccg                                                                                                              >  NB501016:230:HFVJLAFXY:2:21208:24061:13707/1‑146 (MQ=48)
                                        GAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGtgccg                                                                                                              <  NB501016:230:HFVJLAFXY:2:11209:12098:19266/151‑6 (MQ=48)
                                        GAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCAT                                                                                                                       <  NB501016:230:HFVJLAFXY:4:21610:11789:8893/142‑1 (MQ=42)
                                          AACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGT                                                                                                            <  NB501016:230:HFVJLAFXY:2:21207:14630:13227/151‑1 (MQ=43)
                                            CGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGtgccggtcg                                                                                                          <  NB501016:230:HFVJLAFXY:1:11101:22220:5625/151‑10 (MQ=42)
                                             GAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGtgccggtcgt                                                                                                         >  NB501016:230:HFVJLAFXY:3:11408:15202:14460/1‑141 (MQ=41)
                                               ATGCTGCAGCGGCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGC                                                                                                       <  NB501016:230:HFVJLAFXY:3:21401:26582:9433/151‑1 (MQ=29)
                                                 GCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAG                                                                                                     <  NB501016:230:HFVJLAFXY:1:21106:10643:10867/151‑1 (MQ=35)
                                                 GCTGCAGCGGCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCA                                                                                                      <  NB501016:230:HFVJLAFXY:3:11603:24110:13318/150‑1 (MQ=28)
                                                   TGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCT                                                                                                   >  NB501016:230:HFVJLAFXY:1:21304:1777:12563/1‑151 (MQ=36)
                                                    GCAGCGTCACAACAATCAGCCGCCACGTCTGCCACCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTT                                                                                                  >  NB501016:230:HFVJLAFXY:4:21410:3866:19695/1‑151 (MQ=31)
                                                        CGTCACAACAATCAGCCGCGACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATTCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTC                                                                                              <  NB501016:230:HFVJLAFXY:3:21408:16695:16703/151‑1 (MQ=29)
                                                        CGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAAC                                                                                                                              >  NB501016:230:HFVJLAFXY:1:11108:13265:3594/1‑119 (MQ=19)
                                                         GTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCG                                                                                             >  NB501016:230:HFVJLAFXY:3:11612:19247:15049/1‑151 (MQ=39)
                                                         GTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCA                                                                                                                        >  NB501016:230:HFVJLAFXY:4:21404:15451:10124/1‑124 (MQ=26)
                                                          TCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAACGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGG                                                                                            >  NB501016:230:HFVJLAFXY:2:11304:6889:10319/1‑151 (MQ=39)
                                                           CACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGC                                                                                           <  NB501016:230:HFVJLAFXY:2:11309:21507:14912/151‑1 (MQ=40)
                                                           CACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGC                                                                                           >  NB501016:230:HFVJLAFXY:3:21503:4837:12024/1‑151 (MQ=40)
                                                           CACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGG                                                                                            <  NB501016:230:HFVJLAFXY:1:11211:23971:18430/150‑1 (MQ=39)
                                                           CACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGG                                                                                            <  NB501016:230:HFVJLAFXY:1:11203:6585:10575/150‑1 (MQ=39)
                                                            ACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGC                                                                                           >  NB501016:230:HFVJLAFXY:2:21211:15913:9843/1‑150 (MQ=41)
                                                             CAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCA                                                                                          >  NB501016:230:HFVJLAFXY:3:21411:21329:10891/1‑150 (MQ=42)
                                                             CAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCA                                                                                          >  NB501016:230:HFVJLAFXY:3:11509:22700:7431/1‑150 (MQ=42)
                                                             CAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGC                                                                                           >  NB501016:230:HFVJLAFXY:4:21611:10109:7478/1‑149 (MQ=41)
                                                               ACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACG                                                                                       >  NB501016:230:HFVJLAFXY:1:21109:16959:17171/1‑151 (MQ=47)
                                                               ACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACG                                                                                       >  NB501016:230:HFVJLAFXY:1:21203:12261:14999/1‑151 (MQ=47)
                                                               ACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGACGGACGGGCAGCTTCCTCGGCAACG                                                                                       >  NB501016:230:HFVJLAFXY:1:11204:21711:9329/1‑151 (MQ=28)
                                                               ACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATC                                                                                                                      <  NB501016:230:HFVJLAFXY:4:21608:19212:11756/120‑1 (MQ=28)
                                                                 AATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGC                                                                                     <  NB501016:230:HFVJLAFXY:2:21312:2296:14239/151‑1 (MQ=47)
                                                                  ATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATagatcggaagcgcacactctgaactccagtcaccg                                                                                    >  NB501016:230:HFVJLAFXY:1:21204:9487:12210/1‑116 (MQ=24)
                                                                      gccgccacgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAG                                                                                <  NB501016:230:HFVJLAFXY:2:21304:25195:9398/91‑1 (MQ=60)
                                                                      gccgccacgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAG                                                                                <  NB501016:230:HFVJLAFXY:4:21602:10000:6074/91‑1 (MQ=60)
                                                                      gccgccacgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCC                                                                                                               <  NB501016:230:HFVJLAFXY:3:11404:2876:4567/60‑1 (MQ=55)
                                                                        cgccacgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAG                                                                                                     >  NB501016:230:HFVJLAFXY:2:11201:10489:13940/59‑128 (MQ=47)
                                                                         GCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAA                                                                              >  NB501016:230:HFVJLAFXY:1:11312:26362:10384/1‑150 (MQ=55)
                                                                          CCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCCTCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAA                                                                             >  NB501016:230:HFVJLAFXY:3:21406:14948:7694/1‑150 (MQ=47)
                                                                           cacgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAAT                                                                           >  NB501016:230:HFVJLAFXY:2:11101:18287:12118/56‑151 (MQ=60)
                                                                           cacgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAAT                                                                           <  NB501016:230:HFVJLAFXY:1:11106:4780:3164/96‑1 (MQ=60)
                                                                             cgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgagatgcggtcGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCACCGGCGGCAGAAAATTC                                                                         >  NB501016:230:HFVJLAFXY:3:11608:20154:20338/63‑151 (MQ=60)
                                                                             cgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTTCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTC                                                                         >  NB501016:230:HFVJLAFXY:3:21507:26401:17643/54‑151 (MQ=60)
                                                                             cgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTC                                                                         >  NB501016:230:HFVJLAFXY:2:21308:20755:14779/54‑151 (MQ=60)
                                                                             cgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTC                                                                         >  NB501016:230:HFVJLAFXY:4:11412:8176:17519/54‑151 (MQ=60)
                                                                             cgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTC                                                                         >  NB501016:230:HFVJLAFXY:4:11405:16390:18274/54‑151 (MQ=60)
                                                                             cgtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTC                                                                         >  NB501016:230:HFVJLAFXY:3:21504:7684:13166/54‑151 (MQ=60)
                                                                             CGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACCTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTC                                                                         <  NB501016:230:HFVJLAFXY:1:21308:20481:16726/151‑1 (MQ=59)
                                                                              gtctgcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCT                                                                        <  NB501016:230:HFVJLAFXY:2:21302:15923:15960/99‑1 (MQ=60)
                                                                               gctgcctccaccgcggccacgaaagcggcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTG                                                                       <  NB501016:230:HFVJLAFXY:1:11307:11848:8692/100‑1 (MQ=60)
                                                                                  gcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGagatcggaagagcacacggtctg                                                                     >  NB501016:230:HFVJLAFXY:4:21401:2410:1302/49‑127 (MQ=60)
                                                                                  gcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgcca                                                                    <  NB501016:230:HFVJLAFXY:3:11506:15209:12819/103‑6 (MQ=60)
                                                                                  gcctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTG                                                                       <  NB501016:230:HFVJLAFXY:3:11607:24763:6362/100‑1 (MQ=60)
                                                                                   cctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                   <  NB501016:230:HFVJLAFXY:1:11103:17341:11070/104‑7 (MQ=60)
                                                                                   cctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                   <  NB501016:230:HFVJLAFXY:2:21104:11619:11644/104‑7 (MQ=60)
                                                                                   cctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTG                                                                       <  NB501016:230:HFVJLAFXY:2:21303:16713:7891/100‑1 (MQ=60)
                                                                                   ccaccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgcc                                                                     >  NB501016:230:HFVJLAFXY:4:11401:6808:8390/48‑145 (MQ=60)
                                                                                    ctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagg                                                                  >  NB501016:230:HFVJLAFXY:1:11210:22350:8610/47‑144 (MQ=60)
                                                                                    ctccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagg                                                                  >  NB501016:230:HFVJLAFXY:1:21305:4564:7862/47‑144 (MQ=60)
                                                                                     tccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                                                                 >  NB501016:230:HFVJLAFXY:1:21301:18381:14012/46‑143 (MQ=60)
                                                                                      ccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAACAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCCACGGCGGCAGAAAATTCtgccagggc                                                                >  NB501016:230:HFVJLAFXY:2:11209:24170:11247/45‑142 (MQ=60)
                                                                                      ccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGaattttctgccgcc                                                                  <  NB501016:230:HFVJLAFXY:3:11504:10441:6163/105‑15 (MQ=60)
                                                                                      ccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagggc                                                                >  NB501016:230:HFVJLAFXY:4:21504:8343:13223/45‑142 (MQ=60)
                                                                                      ccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagggc                                                                >  NB501016:230:HFVJLAFXY:4:11610:12514:18346/45‑142 (MQ=60)
                                                                                      ccaccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagggc                                                                <  NB501016:230:HFVJLAFXY:3:11609:5915:16219/107‑10 (MQ=60)
                                                                                         ccgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGC                                                             <  NB501016:230:HFVJLAFXY:1:21309:18035:18616/110‑1 (MQ=60)
                                                                                          cgcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCA                                                            <  NB501016:230:HFVJLAFXY:1:11305:1879:10292/111‑1 (MQ=60)
                                                                                           gcggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTCCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCA                                                            >  NB501016:230:HFVJLAFXY:3:21604:5146:18473/40‑150 (MQ=60)
                                                                                            cggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAA                                                           >  NB501016:230:HFVJLAFXY:1:11102:19761:12788/39‑150 (MQ=60)
                                                                                             ggccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAA                                                          >  NB501016:230:HFVJLAFXY:2:21201:7608:16670/38‑150 (MQ=60)
                                                                                              gccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAA                                                        <  NB501016:230:HFVJLAFXY:2:21204:14703:19627/115‑1 (MQ=60)
                                                                                              gccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTATGCCAGGGCGGCAAAA                                                         >  NB501016:230:HFVJLAFXY:3:11507:4344:1849/37‑150 (MQ=60)
                                                                                               ccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAAC                                                       <  NB501016:230:HFVJLAFXY:1:21210:7463:12556/116‑1 (MQ=60)
                                                                                               ccacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAAC                                                       >  NB501016:230:HFVJLAFXY:2:21212:20225:16927/36‑151 (MQ=60)
                                                                                                cacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCA                                                            >  NB501016:230:HFVJLAFXY:2:11205:19151:16341/35‑145 (MQ=60)
                                                                                                cacgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCC                                                                                                >  NB501016:230:HFVJLAFXY:4:11611:9412:3329/35‑109 (MQ=54)
                                                                                                 acgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGACACGCACGCATCATCAAGTGCCGGTCGTGCCGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGT                                                     >  NB501016:230:HFVJLAFXY:4:21501:15916:13312/34‑151 (MQ=51)
                                                                                                  cgaaagcgtcagaggccggcacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCCTCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTC                                                    <  NB501016:230:HFVJLAFXY:1:21203:3259:12303/119‑1 (MQ=60)
                                                                                                  cgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCCACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTC                                                    >  NB501016:230:HFVJLAFXY:3:11502:9650:5410/33‑151 (MQ=60)
                                                                                                  cgaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTC                                                    >  NB501016:230:HFVJLAFXY:3:21412:2280:7311/33‑151 (MQ=60)
                                                                                                   gaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCC                                                   <  NB501016:230:HFVJLAFXY:4:21509:11226:6695/120‑1 (MQ=60)
                                                                                                   gaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCC                                                   >  NB501016:230:HFVJLAFXY:3:11510:19572:14371/32‑151 (MQ=60)
                                                                                                   gaaagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACCCATCATCAAGTGCCGGTCGTGGAGCTTCCTCGGCAACGGCTGCAGAAAATTCTGCCAGGGCGGCAAAATCGTCC                                                   <  NB501016:230:HFVJLAFXY:4:11504:21714:6899/120‑1 (MQ=45)
                                                                                                     aagcgtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGA                                                 <  NB501016:230:HFVJLAFXY:4:21406:20700:5004/122‑1 (MQ=60)
                                                                                                      ctcttttgaggccaccgcatctcgtgctgaagtGGCGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCG                                                  >  NB501016:230:HFVJLAFXY:1:21208:21101:7298/34‑149 (MQ=60)
                                                                                                         gtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACG                                             <  NB501016:230:HFVJLAFXY:3:21601:19480:16221/126‑1 (MQ=60)
                                                                                                         gtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACG                                             <  NB501016:230:HFVJLAFXY:1:11202:24869:9761/126‑1 (MQ=60)
                                                                                                         gtcagaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAG                                                >  NB501016:230:HFVJLAFXY:3:21407:15585:19421/26‑148 (MQ=60)
                                                                                                             gaggccgccacttcagcacgaggTGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATG                                         <  NB501016:230:HFVJLAFXY:4:11411:3420:7087/128‑1 (MQ=60)
                                                                                                             gaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATG                                         >  NB501016:230:HFVJLAFXY:4:21606:13570:4954/22‑151 (MQ=60)
                                                                                                             gaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATG                                         <  NB501016:230:HFVJLAFXY:3:21606:25805:6810/130‑1 (MQ=60)
                                                                                                             gaggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGACGCTGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATG                                         <  NB501016:230:HFVJLAFXY:4:11609:12971:7856/130‑1 (MQ=60)
                                                                                                               ggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCC                                       <  NB501016:230:HFVJLAFXY:1:11308:21162:5961/132‑1 (MQ=60)
                                                                                                               ggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCC                                       >  NB501016:230:HFVJLAFXY:2:11109:22444:12407/20‑151 (MQ=60)
                                                                                                               ggccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCC                                       <  NB501016:230:HFVJLAFXY:1:21310:5682:16142/132‑1 (MQ=60)
                                                                                                                gccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGTGCGGCAAAAACGTCCGAGACGAATGCCA                                      <  NB501016:230:HFVJLAFXY:3:21612:2807:3709/133‑1 (MQ=60)
                                                                                                                 ccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAG                                     <  NB501016:230:HFVJLAFXY:4:21405:11627:9939/134‑1 (MQ=60)
                                                                                                                 ccgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCC                                       >  NB501016:230:HFVJLAFXY:1:11105:18843:17533/18‑149 (MQ=60)
                                                                                                                  cgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAACCGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGG                                    <  NB501016:230:HFVJLAFXY:1:21104:3174:6312/135‑1 (MQ=60)
                                                                                                                  cgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGTCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGG                                    <  NB501016:230:HFVJLAFXY:3:21403:10215:4875/135‑1 (MQ=60)
                                                                                                                  cgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGG                                    >  NB501016:230:HFVJLAFXY:4:21512:14726:11357/17‑151 (MQ=60)
                                                                                                                  cgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGG                                    >  NB501016:230:HFVJLAFXY:2:21305:26028:9474/17‑151 (MQ=60)
                                                                                                                  cgccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGG                                    <  NB501016:230:HFVJLAFXY:4:11508:13186:2967/135‑1 (MQ=60)
                                                                                                                   gccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACTCATCATCAAGTGCCTGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCTGGGCGGCAAAAACGTCTGAGACGAATGCCAGGT                                   <  NB501016:230:HFVJLAFXY:2:11210:24649:8023/136‑1 (MQ=60)
                                                                                                                   gccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGT                                   <  NB501016:230:HFVJLAFXY:3:21603:20969:8722/136‑1 (MQ=60)
                                                                                                                    ccacttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTC                                  <  NB501016:230:HFVJLAFXY:3:11411:7995:10198/137‑1 (MQ=60)
                                                                                                                       cttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGCACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATC                               >  NB501016:230:HFVJLAFXY:2:21309:18272:12184/12‑151 (MQ=60)
                                                                                                                       cttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATC                               >  NB501016:230:HFVJLAFXY:3:21506:2438:14219/12‑151 (MQ=60)
                                                                                                                       cttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATC                               >  NB501016:230:HFVJLAFXY:1:21211:8899:11610/12‑151 (MQ=60)
                                                                                                                       cttcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATC                               >  NB501016:230:HFVJLAFXY:3:21605:24285:17342/12‑151 (MQ=60)
                                                                                                                         tcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCT                              <  NB501016:230:HFVJLAFXY:2:21112:16168:7053/141‑1 (MQ=60)
                                                                                                                         tcagcacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATC                               <  NB501016:230:HFVJLAFXY:2:21301:3178:2238/140‑1 (MQ=60)
                                                                                                                          cagcacgagatgaggtGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCATGTGCCGGTCGTGCAGCTTCCTCGTCAACTGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTG                             <  NB501016:230:HFVJLAFXY:4:11602:9303:14138/134‑1 (MQ=60)
                                                                                                                            GCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAG                                     >  NB501016:230:HFVJLAFXY:1:11111:1490:12603/1‑140 (MQ=60)
                                                                                                                             cacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAC                         >  NB501016:230:HFVJLAFXY:3:21608:22627:8616/6‑151 (MQ=60)
                                                                                                                             cacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAC                         <  NB501016:230:HFVJLAFXY:3:11406:5449:9505/146‑1 (MQ=60)
                                                                                                                             cacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAC                         <  NB501016:230:HFVJLAFXY:1:11302:10152:3654/146‑1 (MQ=60)
                                                                                                                             cacgaGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAAGTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAC                         <  NB501016:230:HFVJLAFXY:3:21611:1461:3845/146‑1 (MQ=60)
                                                                                                                               CGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAC                         >  NB501016:230:HFVJLAFXY:2:11307:8778:3354/1‑149 (MQ=60)
                                                                                                                                GAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGC                      <  NB501016:230:HFVJLAFXY:1:11212:25959:11182/151‑1 (MQ=60)
                                                                                                                                 AGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAACACGTCCGAGACGAATGCCAGGTCATCTGAAACAGC                      >  NB501016:230:HFVJLAFXY:1:11309:2751:19206/1‑150 (MQ=60)
                                                                                                                                 AGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGC                      >  NB501016:230:HFVJLAFXY:4:11603:4984:8333/1‑150 (MQ=60)
                                                                                                                                 AGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTG                                                                                                        >  NB501016:230:HFVJLAFXY:3:11501:2905:3604/1‑68 (MQ=48)
                                                                                                                                  GATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAG                    <  NB501016:230:HFVJLAFXY:3:11405:24928:16958/151‑1 (MQ=60)
                                                                                                                                  GATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGC                      >  NB501016:230:HFVJLAFXY:4:11407:3001:1634/1‑149 (MQ=60)
                                                                                                                                   ATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGC                   <  NB501016:230:HFVJLAFXY:4:21612:18698:10790/151‑1 (MQ=60)
                                                                                                                                    TGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGG                                                              >  NB501016:230:HFVJLAFXY:1:11303:17142:2197/1‑107 (MQ=60)
                                                                                                                                    CGCGGTGGCCTCAAAAGAGGCATCAAAATAATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGCCATAAATTTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCG                  <  NB501016:230:HFVJLAFXY:1:21312:11246:4899/151‑1 (MQ=53)
                                                                                                                                     GCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCT                  <  NB501016:230:HFVJLAFXY:4:21502:21008:5707/150‑1 (MQ=60)
                                                                                                                                      CGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAACAACGTCCGAGACGCATGCCAGGTCATCTGAAACAGCAGCGG                 >  NB501016:230:HFVJLAFXY:3:11403:10885:15729/1‑150 (MQ=60)
                                                                                                                                      CGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGG                 >  NB501016:230:HFVJLAFXY:4:21409:20805:5422/1‑150 (MQ=60)
                                                                                                                                      CGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAC                         >  NB501016:230:HFVJLAFXY:4:11602:11948:11116/1‑142 (MQ=60)
                                                                                                                                        GTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAC              >  NB501016:230:HFVJLAFXY:3:11505:10533:8523/1‑151 (MQ=60)
                                                                                                                                        GTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACA                        <  NB501016:230:HFVJLAFXY:4:21603:10969:4817/141‑1 (MQ=60)
                                                                                                                                          GGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            <  NB501016:230:HFVJLAFXY:2:21108:15910:15994/151‑6 (MQ=60)
                                                                                                                                          GGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            >  NB501016:230:HFVJLAFXY:1:11109:22219:5070/1‑146 (MQ=60)
                                                                                                                                          GGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            <  NB501016:230:HFVJLAFXY:2:11110:2580:14552/151‑6 (MQ=60)
                                                                                                                                          GGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            >  NB501016:230:HFVJLAFXY:3:11402:23379:8669/1‑146 (MQ=60)
                                                                                                                                            CCTCAAAAGCGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACCGCAGCGGaccgc            >  NB501016:230:HFVJLAFXY:3:21508:22158:2768/1‑144 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGACCGCATCAACAAGTGCCGGACGTGCAGCTTCCTAGGCACCGGCGGCAGAAAAATCAGCCAGGGCGGCAAAAACGTCCGAcacgagtcgcaggtcatctgatacagcagcgcaacgg            >  NB501016:230:HFVJLAFXY:3:11511:9730:7057/1‑112 (MQ=25)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgcc                                                                     >  NB501016:230:HFVJLAFXY:3:21404:3520:17533/1‑88 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            >  NB501016:230:HFVJLAFXY:1:21101:2956:10092/1‑144 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            >  NB501016:230:HFVJLAFXY:2:21111:11909:6370/1‑144 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            >  NB501016:230:HFVJLAFXY:1:21303:24777:17468/1‑144 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            >  NB501016:230:HFVJLAFXY:3:11411:11882:13985/1‑144 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGaacgg            >  NB501016:230:HFVJLAFXY:3:11604:8361:1775/1‑144 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAG          <  NB501016:230:HFVJLAFXY:4:11403:7858:1448/151‑1 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAG          <  NB501016:230:HFVJLAFXY:2:11103:10542:8156/151‑1 (MQ=60)
                                                                                                                                            CCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCT                              >  NB501016:230:HFVJLAFXY:2:21310:8493:1885/1‑131 (MQ=60)
                                                                                                                                             CTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAACCAGCAGCGGAACGGAGC         >  NB501016:230:HFVJLAFXY:3:21410:10941:7039/1‑151 (MQ=60)
                                                                                                                                             CTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGC         >  NB501016:230:HFVJLAFXY:4:21506:10841:5602/1‑151 (MQ=60)
                                                                                                                                             CTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGC         >  NB501016:230:HFVJLAFXY:2:11107:10696:19948/1‑151 (MQ=60)
                                                                                                                                             CTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGC         >  NB501016:230:HFVJLAFXY:2:21307:1057:8165/1‑151 (MQ=60)
                                                                                                                                             CTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGC         >  NB501016:230:HFVJLAFXY:2:11306:3321:10617/1‑151 (MQ=60)
                                                                                                                                               CAAAAGAGGCAGCAAATTCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGC       >  NB501016:230:HFVJLAFXY:1:21112:15683:17804/1‑151 (MQ=60)
                                                                                                                                               CAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGC       >  NB501016:230:HFVJLAFXY:1:21103:10838:19849/1‑151 (MQ=60)
                                                                                                                                               CAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGC       >  NB501016:230:HFVJLAFXY:4:11511:4488:5873/1‑151 (MQ=60)
                                                                                                                                               CAAAAGAGGCAGCAAAATCATAAGAAACGAACGCATCATCAACTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTCCCAGTGCGCCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGC       >  NB501016:230:HFVJLAFXY:2:11211:8753:2961/1‑151 (MQ=47)
                                                                                                                                                 AAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCC                                                   >  NB501016:230:HFVJLAFXY:1:11112:2418:10192/1‑105 (MQ=60)
                                                                                                                                                  AAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTC    >  NB501016:230:HFVJLAFXY:2:21110:15279:15877/1‑151 (MQ=60)
                                                                                                                                                   AGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGATACTAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCT   <  NB501016:230:HFVJLAFXY:2:11311:17824:3288/151‑1 (MQ=60)
                                                                                                                                                   AGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATG                                         >  NB501016:230:HFVJLAFXY:4:11607:6902:17505/1‑113 (MQ=60)
                                                                                                                                                    GAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTG  <  NB501016:230:HFVJLAFXY:2:21205:12914:10373/151‑1 (MQ=60)
                                                                                                                                                      |                                                                                                                                                    
AAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGG  >  NC_000913/1429630‑1429928

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