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,423,244 G→C 100% pseudogene (44/444 nt) rzpR → pseudogene, Rac prophage; Bacteriophage Rz lysis protein family;Phage or Prophage Related; putative Rac prophage endopeptidase; completely contained in another CDS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,423,2440GC96.4% 1614.2 / 64.1 1052pseudogene (44/444 nt)rzpRpseudogene, Rac prophage; Bacteriophage Rz lysis protein family;Phage or Prophage Related; putative Rac prophage endopeptidase; completely contained in another CDS
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base C (577/432);  total (599/453)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1078 total aligned reads displayed.

AATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAG  >  NC_000913/1423096‑1423389
                                                                                                                                                    |                                                                                                                                                 
                                                                                                 gcaTCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTgag                                                     <  8:8851/146‑1 (MQ=255)
                                                                                                 ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTTTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACATACAAAGCCCAGCGCGACAAAAAAGCCAGAGAACTGAAGCTGGCGAACGCCGCAATTACTGaaatg                                                  >  3:320232/8‑149 (MQ=11)
                                                                                                 ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                                                  >  1:203707/8‑149 (MQ=11)
                                                                                                 ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                                                  >  1:71891/8‑149 (MQ=11)
                                                                                                 ccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTATGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                                                  >  5:340059/8‑149 (MQ=11)
                                                                                                  cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgc                                                 <  3:116970/143‑1 (MQ=11)
                                                                                                  cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgc                                                 <  4:237535/143‑1 (MQ=11)
                                                                                                  cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                                                  >  8:104908/7‑148 (MQ=11)
                                                                                                  cgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                                                  >  2:289956/7‑148 (MQ=11)
                                                                                                   gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgca                                                <  2:143002/144‑1 (MQ=11)
                                                                                                   gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgc                                                 >  5:380791/6‑148 (MQ=11)
                                                                                                   gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                                                  >  6:319893/6‑147 (MQ=11)
                                                                                                   gcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAAGTACTGacatgca                                                >  2:194008/6‑149 (MQ=11)
                                                                                                    cgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatg                                                  >  7:87046/5‑146 (MQ=11)
                                                                                                      attATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagat                                             <  6:178677/147‑1 (MQ=11)
                                                                                                       ttATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTga                                                      >  2:224948/2‑139 (MQ=14)
                                                                                                       ttATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatg                                            >  6:59644/2‑149 (MQ=11)
                                                                                                         aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGATGCTGGCGAACGCGGCAATTACTGacatgcagatgc                                           <  3:166953/148‑1 (MQ=11)
                                                                                                         aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                                           <  3:188675/148‑1 (MQ=11)
                                                                                                         aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                                           <  7:210720/148‑1 (MQ=11)
                                                                                                         aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                                           <  1:305692/148‑1 (MQ=11)
                                                                                                         aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                                           <  2:58001/148‑1 (MQ=11)
                                                                                                         aTCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCg                                          >  3:9561/1‑149 (MQ=11)
                                                                                                          tctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgc                                                 <  6:2865/141‑1 (MQ=11)
                                                                                                            tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGaa                                                                              <  6:265144/110‑1 (MQ=11)
                                                                                                            tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagatgc                                           >  2:92526/1‑145 (MQ=11)
                                                                                                            tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                                         >  4:7261/1‑147 (MQ=11)
                                                                                                            tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                                         >  7:320605/1‑147 (MQ=11)
                                                                                                            tcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAACCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                                         >  1:41721/1‑147 (MQ=11)
                                                                                                             ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGa                                                                                                   <  1:156664/88‑1 (MQ=11)
                                                                                                             ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAg                                      >  6:7780/1‑149 (MQ=11)
                                                                                                             ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAg                                      >  7:344156/1‑149 (MQ=11)
                                                                                                              cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGt                                         <  6:236505/145‑1 (MQ=11)
                                                                                                              cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                                     <  5:365493/149‑1 (MQ=11)
                                                                                                              cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                                     >  4:147169/1‑149 (MQ=11)
                                                                                                              cGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGc                                     >  1:320074/1‑149 (MQ=11)
                                                                                                               gCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCg                                    >  4:333169/1‑149 (MQ=11)
                                                                                                                 tctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGt                                   >  5:30091/1‑148 (MQ=11)
                                                                                                                 gctGGTGAGCTGCCTCATCGTCTGCCTGTCATGCGCTGTTCATCAATACCGTGCTAACGCCATAACCTACCAAGCCCACCGCGACAACACAGCCAGAGAACTGACGCTGGCGAACGCGGCAATTACAGACATCCAGATGCGTCAGCGTg                                  <  8:73109/148‑1 (MQ=255)
                                                                                                                  ctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGa                                 >  6:340465/1‑149 (MQ=11)
                                                                                                                    ggTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATg                               >  1:286725/1‑149 (MQ=11)
                                                                                                                    ggTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATg                               >  3:337994/1‑149 (MQ=11)
                                                                                                                     gTTATCTGCATCCTCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAACGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                              <  8:350554/149‑1 (MQ=11)
                                                                                                                     gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGATAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGAt                                <  4:220567/147‑1 (MQ=11)
                                                                                                                     gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                              >  1:345148/1‑149 (MQ=11)
                                                                                                                     gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                              >  2:100103/1‑149 (MQ=11)
                                                                                                                     gTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGt                              >  1:255937/1‑149 (MQ=11)
                                                                                                                      ttATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGAAAACGCCATTCCCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGtt                             >  6:132956/1‑149 (MQ=11)
                                                                                                                        aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTGATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCAAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGAtgc                           >  1:106056/1‑149 (MQ=11)
                                                                                                                        aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgc                           >  1:374243/1‑149 (MQ=11)
                                                                                                                        aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgc                           >  8:13792/1‑149 (MQ=11)
                                                                                                                        aTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATACCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgc                           >  6:160143/1‑149 (MQ=11)
                                                                                                                         tCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCgg                                                                >  3:217018/1‑111 (MQ=11)
                                                                                                                         tCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgcagat                                             <  4:300550/130‑1 (MQ=11)
                                                                                                                         tCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct                          >  8:2272/1‑149 (MQ=11)
                                                                                                                         tCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGACAGCGTGATGTtgct                          >  6:100037/1‑149 (MQ=11)
                                                                                                                          cTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctg                         >  2:258387/1‑149 (MQ=11)
                                                                                                                           tGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctgc                        >  7:26540/1‑149 (MQ=11)
                                                                                                                           tGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAATGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgctgc                        >  2:311334/1‑149 (MQ=11)
                                                                                                                            gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCa                                                                                        >  7:13958/1‑84 (MQ=11)
                                                                                                                            gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTgcg                       >  8:199085/1‑149 (MQ=11)
                                                                                                                             catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGacatgc                                                 <  3:295792/122‑1 (MQ=11)
                                                                                                                             catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct                          <  6:92621/145‑1 (MQ=11)
                                                                                                                             catcatCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTtgct                          <  8:389946/145‑1 (MQ=11)
                                                                                                                              atcatcGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCt                     <  3:366722/149‑1 (MQ=11)
                                                                                                                              atcatcGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCt                     >  7:247441/1‑149 (MQ=11)
                                                                                                                                catcgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTc                    >  7:59748/3‑148 (MQ=11)
                                                                                                                                 atcgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGa                  >  3:268057/2‑149 (MQ=11)
                                                                                                                                  tcgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTTGAt                 >  2:342343/1‑149 (MQ=11)
                                                                                                                                  tcgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAt                 >  4:11548/1‑149 (MQ=11)
                                                                                                                                   cgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATg                >  6:205881/1‑149 (MQ=11)
                                                                                                                                    gtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATg                >  5:63194/1‑148 (MQ=11)
                                                                                                                                    gtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGc               >  7:381594/1‑148 (MQ=11)
                                                                                                                                     tcTGCCTGTCATGGGCTGTTAATCATTACCGTGATACCGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCa              <  5:181661/149‑1 (MQ=11)
                                                                                                                                       tGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGCCAGCGTGATGTTGCTGCGCTCGATGCaa             >  8:366294/1‑148 (MQ=11)
                                                                                                                                        gCCTGTCATGGGCTGTTAATCAATACCCTGCTAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGACAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCaaaa           <  7:287120/149‑1 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTTAAAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACTCGGCAATTACTGACATGCAGAAGCGTCAGCGTGATGTTGCAGCGCTCGATGCaaaa           <  5:131800/148‑1 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGATAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCACTCGATTCaaaa           <  2:364512/148‑1 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTAc                                                         <  7:16021/102‑1 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCaaaa           <  5:156798/148‑1 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCaaaa           <  4:27493/148‑1 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCaaaa           >  5:103651/1‑148 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCaa             <  3:68295/146‑1 (MQ=11)
                                                                                                                                         ccTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAAt          <  8:167894/149‑1 (MQ=11)
                                                                                                                                          cTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTGCCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATa         <  3:151679/149‑1 (MQ=11)
                                                                                                                                          cTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGTGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATa         <  5:237517/149‑1 (MQ=11)
                                                                                                                                          cTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATa         >  7:188397/1‑149 (MQ=11)
                                                                                                                                           tGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAt                 <  5:280447/140‑1 (MQ=11)
                                                                                                                                           tGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATa         >  2:393500/1‑148 (MQ=11)
                                                                                                                                           tGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATa         >  4:16663/1‑148 (MQ=11)
                                                                                                                                            gTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATaca       >  7:1217/1‑149 (MQ=11)
                                                                                                                                             tCATGGGCTGTTAATGATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATac        <  8:32766/147‑1 (MQ=11)
                                                                                                                                             tCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACACAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATacac      <  4:196246/149‑1 (MQ=11)
                                                                                                                                             tCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGc                                                               <  8:290882/92‑1 (MQ=11)
                                                                                                                                             tCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATaca       <  8:309867/148‑1 (MQ=11)
                                                                                                                                             gactGGGCTGTTAATCATTCCCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACATAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGGCAGCGTGATGTTGCTGCGCTCGATGCAAAATaca       <  5:211257/145‑1 (MQ=11)
                                                                                                                                               aTGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACg     >  2:170507/1‑148 (MQ=11)
                                                                                                                                               aTGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACg     >  5:321560/1‑148 (MQ=11)
                                                                                                                                               aTGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACg     >  6:365426/1‑148 (MQ=11)
                                                                                                                                                tGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACg     >  4:37089/1‑147 (MQ=11)
                                                                                                                                                    |                                                                                                                                                 
AATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAG  >  NC_000913/1423096‑1423389

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output Only 200 of 1607 total aligned reads displayed.
CAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTT  >  NC_000913/1423095‑1423394
                                                                                                                                                     |                                                                                                                                                      
ctggttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCT                                                                                                                                                       >  NB501016:230:HFVJLAFXY:4:21611:11492:9188/106‑151 (MQ=60)
 tggttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTG                                                                                                                                                      >  NB501016:230:HFVJLAFXY:4:21603:20592:11760/105‑151 (MQ=60)
  ggttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagacgcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTTCATCATCGTCTGCCTGTCATGGGCTG                                                                                                                                                      <  NB501016:230:HFVJLAFXY:3:21608:21899:16995/47‑1 (MQ=60)
  ggttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTGCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTG                                                                                                                                                      <  NB501016:230:HFVJLAFXY:2:21108:8001:14611/47‑1 (MQ=60)
  ggttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGT                                                                                                                                                     <  NB501016:230:HFVJLAFXY:1:21202:15902:10714/48‑1 (MQ=60)
  ggttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGT                                                                                                                                                     <  NB501016:230:HFVJLAFXY:3:11412:2506:12781/48‑1 (MQ=60)
  ggttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacgggcagagagattgatgtatgagcagaggcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCT                                                                                                                                                       <  NB501016:230:HFVJLAFXY:3:11411:21917:11622/46‑1 (MQ=60)
   gttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTT                                                                                                                                                    <  NB501016:230:HFVJLAFXY:1:11108:24450:5232/49‑1 (MQ=60)
   gttatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacgggcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTT                                                                                                                                                    <  NB501016:230:HFVJLAFXY:1:21312:7846:3428/49‑1 (MQ=60)
   gttatggtcagttcgagcataaggctgacagactgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTTCATCATCGTCTGCCTGTCATGGGCTGTT                                                                                                                                                    <  NB501016:230:HFVJLAFXY:2:11109:24665:1889/49‑1 (MQ=60)
     tatggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTA                                                                                                                                                   >  NB501016:230:HFVJLAFXY:1:11212:3144:14573/101‑150 (MQ=60)
      atggtcagttcgagcataaggctaacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAAT                                                                                                                                                 <  NB501016:230:HFVJLAFXY:3:11611:15049:6349/52‑1 (MQ=60)
        ggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCA                                                                                                                                               <  NB501016:230:HFVJLAFXY:3:11512:25894:11627/54‑1 (MQ=60)
        ggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATC                                                                                                                                                <  NB501016:230:HFVJLAFXY:4:11609:6894:10711/53‑1 (MQ=60)
        ggtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagagggcggaacggtcagagagattgatgtatgagcagagtcaacgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATC                                                                                                                                                >  NB501016:230:HFVJLAFXY:3:11507:24070:12367/98‑150 (MQ=60)
        ggtcacttccatgataacgccgtaaggccgaatgcaaaaatcaaaccagcgggccgcaagggcatacaccatgcaggaatcccagaggccccgcccataaaTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATC                                                                                                                                                <  NB501016:230:HFVJLAFXY:2:11207:5297:17031/49‑1 (MQ=60)
         gtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCAT                                                                                                                                              >  NB501016:230:HFVJLAFXY:1:21307:6578:18554/97‑151 (MQ=60)
         gtcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCAT                                                                                                                                              >  NB501016:230:HFVJLAFXY:3:21603:14343:14356/97‑151 (MQ=60)
          tcagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATT                                                                                                                                             >  NB501016:230:HFVJLAFXY:3:11406:17036:8871/96‑151 (MQ=60)
           cagttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATT                                                                                                                                             >  NB501016:230:HFVJLAFXY:4:21507:14408:12622/95‑150 (MQ=60)
             gttcgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACC                                                                                                                                          >  NB501016:230:HFVJLAFXY:4:21605:22059:6390/93‑151 (MQ=60)
              AAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAAACAT                                                                                                                                              >  NB501016:230:HFVJLAFXY:3:21506:10832:15822/1‑146 (MQ=60)
                cgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTG                                                                                                                                       >  NB501016:230:HFVJLAFXY:4:11406:8339:17648/90‑151 (MQ=60)
                cgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTG                                                                                                                                       >  NB501016:230:HFVJLAFXY:4:11411:15694:13757/90‑151 (MQ=60)
                cgagcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGT                                                                                                                                        <  NB501016:230:HFVJLAFXY:3:21601:10620:5050/61‑1 (MQ=60)
                   gcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATA                                                                                                                                    <  NB501016:230:HFVJLAFXY:2:21309:17888:11885/65‑1 (MQ=60)
                   gcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGAT                                                                                                                                     >  NB501016:230:HFVJLAFXY:1:11304:20199:14119/87‑150 (MQ=60)
                   gcataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGAT                                                                                                                                     >  NB501016:230:HFVJLAFXY:3:21411:14109:19841/87‑150 (MQ=60)
                    cataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATA                                                                                                                                    >  NB501016:230:HFVJLAFXY:2:21312:25213:7558/86‑150 (MQ=60)
                    CGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAA                                                                                                                                   <  NB501016:230:HFVJLAFXY:4:21408:13651:13108/151‑1 (MQ=60)
                     ataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAAC                                                                                                                                  >  NB501016:230:HFVJLAFXY:1:11112:4771:10020/85‑151 (MQ=60)
                     ataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgagatTATCTCCTCTATGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAAC                                                                                                                                  >  NB501016:230:HFVJLAFXY:4:21402:6869:5043/85‑151 (MQ=60)
                     ataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtaagagagattgatgtatgagaagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTTCCTGTCATGGGCTGTTAATCATTACAGTGATAAC                                                                                                                                  >  NB501016:230:HFVJLAFXY:1:21204:2807:3570/85‑151 (MQ=60)
                     ataaggctgacagcctgattgcaaaattcaaagaagcgggcggaacgggcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAAC                                                                                                                                  <  NB501016:230:HFVJLAFXY:1:21309:3491:2986/67‑1 (MQ=60)
                        aggctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCC                                                                                                                               <  NB501016:230:HFVJLAFXY:1:21206:12312:2179/70‑1 (MQ=60)
                        aggctgacagcctgattgcaaaattcaaagaagcgggcgcaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCC                                                                                                                               >  NB501016:230:HFVJLAFXY:1:11202:10641:9049/82‑151 (MQ=60)
                         ggctgacagcctgattgcaaaattcaaagaagcgggcggcacggtcagagagagtgatgtatgagcagcgtccccgcgctTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCC                                                                                                                               >  NB501016:230:HFVJLAFXY:3:11607:4251:16643/81‑150 (MQ=60)
                          gctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCAT                                                                                                                             >  NB501016:230:HFVJLAFXY:3:21508:10035:9119/80‑151 (MQ=60)
                           ctgacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAAAGCCATT                                                                                                                            >  NB501016:230:HFVJLAFXY:2:21103:1497:17365/79‑151 (MQ=60)
                            ggacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtataagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATta                                                                                                                           <  NB501016:230:HFVJLAFXY:4:21601:7489:7949/74‑3 (MQ=60)
                             gacagcctgattgcaaaattcaaagaagcgtgcgtaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGTTTATCTGCATCATCGTCTGCCTGTCTTGGGATGTTAATCATTACCGTGATAACGCCATT                                                                                                                            <  NB501016:230:HFVJLAFXY:2:11107:21949:8212/73‑1 (MQ=60)
                             gacagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTTGTTATCTGCATCATCGTCTGCCTGTAATGGGCTGTTAATCATTACCGTGATAACGCCATta                                                                                                                           <  NB501016:230:HFVJLAFXY:3:11610:19288:3983/74‑3 (MQ=60)
                              acagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATtacc                                                                                                                         <  NB501016:230:HFVJLAFXY:2:21205:11179:2380/76‑5 (MQ=60)
                               cagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATtacct                                                                                                                        <  NB501016:230:HFVJLAFXY:4:11611:17939:17368/77‑6 (MQ=60)
                               cagcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATta                                                                                                                           >  NB501016:230:HFVJLAFXY:3:11508:4154:11519/75‑146 (MQ=60)
                                 gccttattgcaaaattcaaagaagcgggcggaacggtcagagagtttgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTAC                                                                                                                      <  NB501016:230:HFVJLAFXY:1:21304:2090:16561/79‑1 (MQ=60)
                                 gcctgattgcaaaattcaaagaagcgggcggatcggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCACTCTGGTTATCTGGATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTAC                                                                                                                      <  NB501016:230:HFVJLAFXY:3:21407:14298:13945/79‑1 (MQ=60)
                                 gcctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTAC                                                                                                                      <  NB501016:230:HFVJLAFXY:2:21107:11799:12561/79‑1 (MQ=60)
                                  cctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgttTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTAC                                                                                                                      >  NB501016:230:HFVJLAFXY:2:11111:15749:14325/72‑150 (MQ=21)
                                   ctgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgattatcgCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACA                                                                                                                     >  NB501016:230:HFVJLAFXY:4:11501:2758:1887/76‑150 (MQ=60)
                                    tgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAA                                                                                                                    >  NB501016:230:HFVJLAFXY:2:11106:16989:2630/70‑150 (MQ=60)
                                    tgattgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATA                                                                                                                                    >  NB501016:230:HFVJLAFXY:1:11308:14654:14023/70‑134 (MQ=60)
                                      attgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATta                                                                                                                           >  NB501016:230:HFVJLAFXY:1:21211:16285:13322/68‑139 (MQ=60)
                                      attgcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGC                                                                                                                 <  NB501016:230:HFVJLAFXY:4:21609:8927:10185/84‑1 (MQ=60)
                                         gcaaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGccca                                                                                                              <  NB501016:230:HFVJLAFXY:4:11506:5622:6032/87‑5 (MQ=60)
                                         gcaaaattcaaagaagcgggcggaacgggcagatagattgatgtatgagcagaggcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGA                                                                                                                                      <  NB501016:230:HFVJLAFXY:4:11601:25200:9289/63‑1 (MQ=60)
                                          caaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGcccag                                                                                                             <  NB501016:230:HFVJLAFXY:1:21212:18169:6361/88‑6 (MQ=60)
                                          caaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGccc                                                                                                               >  NB501016:230:HFVJLAFXY:2:21112:1678:12652/64‑146 (MQ=60)
                                          caaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGccc                                                                                                               >  NB501016:230:HFVJLAFXY:2:11201:6216:16375/64‑146 (MQ=60)
                                          caaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGCCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGccacg                                                                                                             >  NB501016:230:HFVJLAFXY:1:11210:6323:7104/64‑146 (MQ=60)
                                           aaaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGC                                                                                                            >  NB501016:230:HFVJLAFXY:1:21102:7943:1681/63‑151 (MQ=60)
                                            aaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCG                                                                                                           >  NB501016:230:HFVJLAFXY:1:21105:17199:5692/62‑151 (MQ=60)
                                            aaattcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCG                                                                                                           <  NB501016:230:HFVJLAFXY:4:21406:23783:16775/90‑1 (MQ=60)
                                              attcaaagaagcgggcggtacggtcagagagtttgatgtatgagcagagtcaccgcgatTATCTCCGCTCTTGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCG                                                                                                         <  NB501016:230:HFVJLAFXY:2:11311:14257:14921/92‑1 (MQ=60)
                                              attcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCG                                                                                                         <  NB501016:230:HFVJLAFXY:1:11106:12187:10490/92‑1 (MQ=60)
                                              attcaaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCG                                                                                                         >  NB501016:230:HFVJLAFXY:4:11408:16955:5939/60‑151 (MQ=60)
                                                  aaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgattatcgCCGCTCTGGTTATCTGCATCATCGTCTCCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACA                                                                                                      >  NB501016:230:HFVJLAFXY:3:11407:6308:4389/61‑150 (MQ=60)
                                                  aaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACA                                                                                                      >  NB501016:230:HFVJLAFXY:2:21304:21610:13803/56‑150 (MQ=60)
                                                  aaagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACA                                                                                                      >  NB501016:230:HFVJLAFXY:3:11510:8828:8056/56‑150 (MQ=60)
                                                   aagaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAA                                                                                                     >  NB501016:230:HFVJLAFXY:2:21306:21518:2976/55‑150 (MQ=60)
                                                    agaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAA                                                                                                    <  NB501016:230:HFVJLAFXY:4:11412:2088:17884/97‑1 (MQ=60)
                                                     gaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAA                                                                                                  <  NB501016:230:HFVJLAFXY:1:11205:19016:20086/99‑1 (MQ=60)
                                                     gaagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAA                                                                                                   >  NB501016:230:HFVJLAFXY:1:21209:26090:11987/53‑150 (MQ=60)
                                                      aagcgtgctgaacattcagatagattgatgtatgagcagagtcaccgcgatTATCTACGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAAT                                                                                                 <  NB501016:230:HFVJLAFXY:2:11204:26176:14047/100‑1 (MQ=60)
                                                      aagcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAAT                                                                                                 <  NB501016:230:HFVJLAFXY:1:21301:15971:3575/100‑1 (MQ=60)
                                                       agcgggcggaacgggcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCGTGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATG                                                                                                <  NB501016:230:HFVJLAFXY:3:11505:10694:4619/101‑1 (MQ=60)
                                                        gcgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAACCCCAGCGCGACAAAAATGC                                                                                               <  NB501016:230:HFVJLAFXY:2:21110:5160:16563/102‑1 (MQ=60)
                                                         cgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCC                                                                                              <  NB501016:230:HFVJLAFXY:4:11505:1119:7663/103‑1 (MQ=0)
                                                         cgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCC                                                                                              <  NB501016:230:HFVJLAFXY:2:21303:4559:9243/103‑1 (MQ=0)
                                                         cgggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCC                                                                                              <  NB501016:230:HFVJLAFXY:3:21511:6837:18959/103‑1 (MQ=0)
                                                           ggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAG                                                                                            <  NB501016:230:HFVJLAFXY:2:11310:16438:13063/105‑1 (MQ=0)
                                                           ggcggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAG                                                                                            <  NB501016:230:HFVJLAFXY:4:11405:22732:6350/105‑1 (MQ=0)
                                                              ggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAG                                                                                            >  NB501016:230:HFVJLAFXY:2:21307:3149:5276/44‑148 (MQ=0)
                                                              ggaacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACC                                                                                                                                          >  NB501016:230:HFVJLAFXY:4:21511:1633:14105/44‑102 (MQ=60)
                                                                aacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCTCTCTGGTTATCTGCATCATCTGATGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGACAAC                                                                                       <  NB501016:230:HFVJLAFXY:1:21107:22324:7767/110‑1 (MQ=0)
                                                                aacggtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATTCCAGAGAAC                                                                                       >  NB501016:230:HFVJLAFXY:2:21209:15108:2301/42‑151 (MQ=60)
                                                                 acggtcagatagcttgatgtatgagcagagtcaccgagatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACT                                                                                      <  NB501016:230:HFVJLAFXY:2:11101:3696:2795/111‑1 (MQ=9)
                                                                    gtcagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGA                                                                                    >  NB501016:230:HFVJLAFXY:3:11401:19217:7679/38‑150 (MQ=0)
                                                                      cagagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGC                                                                                 <  NB501016:230:HFVJLAFXY:1:21306:18709:17008/116‑1 (MQ=0)
                                                                       agagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCT                                                                                >  NB501016:230:HFVJLAFXY:1:11301:8505:16132/35‑151 (MQ=0)
                                                                        gagagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAG                                                                                            <  NB501016:230:HFVJLAFXY:2:21211:14913:4581/105‑1 (MQ=0)
                                                                          gagattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGC                                                                             <  NB501016:230:HFVJLAFXY:2:11307:24051:8700/120‑1 (MQ=5)
                                                                           agtttgttgtttgagcagagtcacctcgatTATCTCGGCTCTGGTTATCTGCATCTTGGTCTGCCTGTCTTGGGCTGTTAATCATTACCGTGATAAGGCCATTACCTACAATGCCCTGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCG                                                                            <  NB501016:230:HFVJLAFXY:3:21504:24238:15108/121‑1 (MQ=23)
                                                                           agattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGTG                                                                            <  NB501016:230:HFVJLAFXY:3:21612:17780:9544/121‑1 (MQ=3)
                                                                           agattgatgtatgagcagaggcaccgcgatTATCTCCGCTCTGGTTAGCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCG                                                                            <  NB501016:230:HFVJLAFXY:1:21201:15882:9551/121‑1 (MQ=19)
                                                                           ATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCG                                                                            <  NB501016:230:HFVJLAFXY:2:11304:11079:19849/151‑1 (MQ=60)
                                                                            gattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCG                                                                            >  NB501016:230:HFVJLAFXY:3:11404:11033:16863/30‑150 (MQ=7)
                                                                            gattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCG                                                                            >  NB501016:230:HFVJLAFXY:4:11508:19595:3866/30‑150 (MQ=7)
                                                                            gattgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCG                                                                            >  NB501016:230:HFVJLAFXY:4:21611:5862:4973/30‑150 (MQ=7)
                                                                             attgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAA                                                                          <  NB501016:230:HFVJLAFXY:2:21301:15159:17406/123‑1 (MQ=8)
                                                                              ttgatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTTCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGACCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCTAAC                                                                         >  NB501016:230:HFVJLAFXY:4:21401:19961:1307/28‑151 (MQ=17)
                                                                                gatgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGC                                                                       >  NB501016:230:HFVJLAFXY:3:21404:4893:8226/26‑151 (MQ=8)
                                                                                  tgtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGG                                                                     >  NB501016:230:HFVJLAFXY:2:11312:2531:15693/24‑151 (MQ=8)
                                                                                   gtatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGC                                                                    <  NB501016:230:HFVJLAFXY:2:11305:12081:3529/129‑1 (MQ=8)
                                                                                    tatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCCGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGC                                                                    >  NB501016:230:HFVJLAFXY:3:21408:20193:15335/22‑150 (MQ=8)
                                                                                    tatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGC                                                                    >  NB501016:230:HFVJLAFXY:4:11412:12117:7798/22‑150 (MQ=8)
                                                                                    tatgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGC                                                                    >  NB501016:230:HFVJLAFXY:3:11502:10057:7235/22‑150 (MQ=8)
                                                                                     atgagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGA                                                                                         >  NB501016:230:HFVJLAFXY:2:21302:2610:5884/21‑128 (MQ=0)
                                                                                       gagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATT                                                                >  NB501016:230:HFVJLAFXY:1:11107:22048:8295/19‑151 (MQ=10)
                                                                                       gagcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAAT                                                                 <  NB501016:230:HFVJLAFXY:1:21312:11031:17807/132‑1 (MQ=10)
                                                                                        agcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACGGAAGCTGGCGAACGCGGCAATT                                                                >  NB501016:230:HFVJLAFXY:3:21508:9234:18394/18‑150 (MQ=15)
                                                                                        agcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCAACAAAGCCCAGCGCGACAAAACTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTA                                                               <  NB501016:230:HFVJLAFXY:3:21605:4025:15486/134‑1 (MQ=60)
                                                                                        agcagagtcaccgcgatTATCTCCGCTCTGGGTATCTGCATCATCGTCTGCCTGTCATGGGCTGGTAATCAGTACCGTGCTAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGCACGCGGCAATT                                                                >  NB501016:230:HFVJLAFXY:3:21510:1842:18127/18‑150 (MQ=13)
                                                                                         gcagattcaccgcgataATCTCCGCTCTGGTTATCAGCATCATCGTCTGCCTGTCATCGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTAC                                                              <  NB501016:230:HFVJLAFXY:1:11209:16542:11501/134‑1 (MQ=28)
                                                                                         gcagagtcaccgcgatTATCTCCGCTCTGGTTTTCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAAGTGAAGCTGGCGAACGCGGCAATTAC                                                              >  NB501016:230:HFVJLAFXY:1:21203:26056:16936/17‑151 (MQ=16)
                                                                                         gcagagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTAC                                                              <  NB501016:230:HFVJLAFXY:1:11311:11665:3388/135‑1 (MQ=17)
                                                                                           agagtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTG                                                            <  NB501016:230:HFVJLAFXY:1:11203:26541:8204/137‑1 (MQ=21)
                                                                                             agtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAC                                                          <  NB501016:230:HFVJLAFXY:4:11602:26281:16243/139‑1 (MQ=20)
                                                                                              gtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGCCAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAC                                                          >  NB501016:230:HFVJLAFXY:4:11503:1761:16225/12‑150 (MQ=20)
                                                                                              gtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACA                                                         <  NB501016:230:HFVJLAFXY:4:21411:9975:8202/140‑1 (MQ=18)
                                                                                              gtcaccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCC                                                                                              >  NB501016:230:HFVJLAFXY:3:11612:4758:16986/12‑114 (MQ=11)
                                                                                                caccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGccatg                                                       >  NB501016:230:HFVJLAFXY:2:21203:10717:13710/10‑146 (MQ=22)
                                                                                                caccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATG                                                       >  NB501016:230:HFVJLAFXY:3:11602:20304:19527/10‑151 (MQ=15)
                                                                                                caccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTG                                                            >  NB501016:230:HFVJLAFXY:1:11312:20925:18099/10‑146 (MQ=22)
                                                                                                caccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATCCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATG                                                       >  NB501016:230:HFVJLAFXY:1:11109:20810:19664/10‑151 (MQ=14)
                                                                                                caccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTCTCATGGGCTGTTAATCATTACGGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATG                                                       <  NB501016:230:HFVJLAFXY:3:11506:15675:7074/142‑1 (MQ=14)
                                                                                                 accgcgattatctccggaCTGGTTATCTGCATCATCGTCTGCCTCTCAAGGGCTGTTAATCATTACCGTGATAAAGCCATTACCTACAAAGCCCAGCGAGACAAAAATGCCAGAGAACTGAAGCGGGCGATCGCGGCAATTACTGAcctgc                                                      >  NB501016:230:HFVJLAFXY:3:11402:18827:1102/19‑146 (MQ=17)
                                                                                                 accgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCTACAAAAATGCCAGAGAACT                                                                                      >  NB501016:230:HFVJLAFXY:2:11212:19086:18169/9‑119 (MQ=10)
                                                                                                  ccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGTGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGTGAACGCGGCAATTACTG                                                            <  NB501016:230:HFVJLAFXY:2:11102:5918:16724/137‑1 (MQ=23)
                                                                                                  ccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGATAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGC                                                      >  NB501016:230:HFVJLAFXY:1:11101:13282:18448/8‑150 (MQ=14)
                                                                                                  ccgcgatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTG                                                            <  NB501016:230:HFVJLAFXY:4:21606:3997:14842/137‑1 (MQ=24)
                                                                                                      gatTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATG                                                 <  NB501016:230:HFVJLAFXY:2:11104:13261:6136/148‑1 (MQ=21)
                                                                                                         TATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGT                                              >  NB501016:230:HFVJLAFXY:1:11103:20784:15863/1‑151 (MQ=25)
                                                                                                         TATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCG                                               <  NB501016:230:HFVJLAFXY:2:11209:25241:12170/150‑1 (MQ=25)
                                                                                                         TATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCG                                               <  NB501016:230:HFVJLAFXY:3:11604:17411:13562/150‑1 (MQ=25)
                                                                                                          ATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTC                                             >  NB501016:230:HFVJLAFXY:4:11603:9122:19088/1‑151 (MQ=23)
                                                                                                           TCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTC                                             >  NB501016:230:HFVJLAFXY:1:11102:19881:2593/1‑150 (MQ=21)
                                                                                                           TCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCA                                                     <  NB501016:230:HFVJLAFXY:4:21510:22637:14306/142‑1 (MQ=22)
                                                                                                            CTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAACGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTC                                             >  NB501016:230:HFVJLAFXY:1:21207:26770:12485/1‑149 (MQ=19)
                                                                                                            CTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTC                                             >  NB501016:230:HFVJLAFXY:2:11309:15801:6124/1‑149 (MQ=20)
                                                                                                            CTCCGCTCTGGTTATCAGCATCATCGTCTGCCTGACATGGGCTGTTAATCATTACCGTGATAACTCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAGCGAGGCAAATACTGACATGCAGATGCGTC                                             >  NB501016:230:HFVJLAFXY:1:21310:9303:13188/1‑149 (MQ=17)
                                                                                                             TCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGAC                                                                                                       <  NB501016:230:HFVJLAFXY:1:11305:20434:3872/90‑1 (MQ=60)
                                                                                                              CCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCA                                            <  NB501016:230:HFVJLAFXY:3:11601:6876:6518/148‑1 (MQ=18)
                                                                                                                   CTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATG                                    >  NB501016:230:HFVJLAFXY:1:11110:2340:11181/1‑151 (MQ=16)
                                                                                                                     GGTTATCTGCATCCTCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTT                                  >  NB501016:230:HFVJLAFXY:3:21604:1712:10619/1‑151 (MQ=21)
                                                                                                                     GGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGCCATGCAGATGCG                                               >  NB501016:230:HFVJLAFXY:1:21208:18499:13464/1‑138 (MQ=9)
                                                                                                                     GGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTT                                  >  NB501016:230:HFVJLAFXY:1:21303:14322:7582/1‑151 (MQ=9)
                                                                                                                     GGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTT                                  >  NB501016:230:HFVJLAFXY:2:11206:12168:3753/1‑151 (MQ=9)
                                                                                                                        TATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCT                               >  NB501016:230:HFVJLAFXY:2:21109:6733:5467/1‑151 (MQ=6)
                                                                                                                        AATCTGCATCATCGTCTGCCTGTCATGGGCTGTTGATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCAAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGATGCT                               >  NB501016:230:HFVJLAFXY:1:11207:3564:11771/1‑151 (MQ=5)
                                                                                                                         ATCTGCTTCATCGTCTGCCTGTCATGTGCTGTTAATCATTACCGGGATATCGCGATTACCTACATAGCGCTGGGCGACAAAAATGCGAGAGAACTGTATCTGGCGAACGCGGCAATTACTGAGATGCAGATTCGTCAGCGTGATGTTGCT                               <  NB501016:230:HFVJLAFXY:3:11605:11831:1066/150‑1 (MQ=17)
                                                                                                                          TCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGC                             >  NB501016:230:HFVJLAFXY:4:11608:3428:5698/1‑151 (MQ=6)
                                                                                                                           CTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCG                            >  NB501016:230:HFVJLAFXY:2:21208:10070:13535/1‑151 (MQ=6)
                                                                                                                           CTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCG                            >  NB501016:230:HFVJLAFXY:2:11301:4786:4928/1‑151 (MQ=6)
                                                                                                                            tgcctcatcggcttCCTGTCATGGGCTGTTTATCATTACCGTGCTAACGCCATTTCCTAACAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTAC                                                              <  NB501016:230:HFVJLAFXY:2:21102:24895:1804/102‑1 (MQ=18)
                                                                                                                            TGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGC                           >  NB501016:230:HFVJLAFXY:3:21506:6270:8580/1‑151 (MQ=6)
                                                                                                                            TGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAG                                                                                            >  NB501016:230:HFVJLAFXY:4:11403:23860:10709/1‑86 (MQ=7)
                                                                                                                             GCATCCTCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTCCAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTG                              <  NB501016:230:HFVJLAFXY:3:21610:19796:9761/147‑1 (MQ=17)
                                                                                                                             GCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCT                          <  NB501016:230:HFVJLAFXY:3:21503:13310:6440/151‑1 (MQ=6)
                                                                                                                             GCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTG                              <  NB501016:230:HFVJLAFXY:2:21311:7138:15294/147‑1 (MQ=6)
                                                                                                                             GCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTG                              <  NB501016:230:HFVJLAFXY:4:21409:20122:5342/147‑1 (MQ=6)
                                                                                                                             GCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTG                              <  NB501016:230:HFVJLAFXY:4:11402:4544:3873/147‑1 (MQ=6)
                                                                                                                              CATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTC                         >  NB501016:230:HFVJLAFXY:1:21103:9348:13205/1‑151 (MQ=6)
                                                                                                                               ATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCG                        >  NB501016:230:HFVJLAFXY:4:11604:7346:20194/1‑151 (MQ=8)
                                                                                                                                TCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCG                        >  NB501016:230:HFVJLAFXY:1:21111:14474:10739/1‑150 (MQ=10)
                                                                                                                                TCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCG                        >  NB501016:230:HFVJLAFXY:2:11211:20261:19389/1‑150 (MQ=10)
                                                                                                                                 CATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAT                      >  NB501016:230:HFVJLAFXY:4:21404:24692:15836/1‑151 (MQ=14)
                                                                                                                                  ATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATG                     <  NB501016:230:HFVJLAFXY:4:21411:13535:1115/151‑1 (MQ=16)
                                                                                                                                  ATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATG                     >  NB501016:230:HFVJLAFXY:3:21607:11808:17230/1‑151 (MQ=16)
                                                                                                                                    CGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATACCGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGC                    >  NB501016:230:HFVJLAFXY:3:21609:19230:18084/1‑150 (MQ=18)
                                                                                                                                         GCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAAT               <  NB501016:230:HFVJLAFXY:3:11606:25447:5913/150‑1 (MQ=6)
                                                                                                                                         GCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAAT               >  NB501016:230:HFVJLAFXY:2:11202:7327:4629/1‑150 (MQ=6)
                                                                                                                                         GCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAA                <  NB501016:230:HFVJLAFXY:3:21403:4999:3897/149‑1 (MQ=6)
                                                                                                                                         GCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAA                 <  NB501016:230:HFVJLAFXY:4:11512:1169:6860/148‑1 (MQ=6)
                                                                                                                                          CCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTGCCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATAC             <  NB501016:230:HFVJLAFXY:2:11306:2252:17121/151‑1 (MQ=6)
                                                                                                                                          CCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATAC             <  NB501016:230:HFVJLAFXY:4:21608:10506:9377/151‑1 (MQ=6)
                                                                                                                                            TGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACAC           >  NB501016:230:HFVJLAFXY:4:11401:18710:5611/1‑151 (MQ=6)
                                                                                                                                            TGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACAC           >  NB501016:230:HFVJLAFXY:3:21502:17735:10728/1‑151 (MQ=6)
                                                                                                                                             GTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAACATACACG          <  NB501016:230:HFVJLAFXY:4:21512:5511:18035/151‑1 (MQ=6)
                                                                                                                                             GTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACG          <  NB501016:230:HFVJLAFXY:1:21110:21647:2465/151‑1 (MQ=6)
                                                                                                                                             GTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACG          <  NB501016:230:HFVJLAFXY:3:11409:13397:10198/151‑1 (MQ=6)
                                                                                                                                             GTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACG          <  NB501016:230:HFVJLAFXY:4:11404:21998:14016/151‑1 (MQ=6)
                                                                                                                                             GTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCA                                                                   >  NB501016:230:HFVJLAFXY:4:21504:6330:16125/1‑94 (MQ=21)
                                                                                                                                               CATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCTAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGA         >  NB501016:230:HFVJLAFXY:2:21104:9119:10676/1‑150 (MQ=6)
                                                                                                                                               CATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGA         >  NB501016:230:HFVJLAFXY:1:11104:25599:16314/1‑150 (MQ=6)
                                                                                                                                               CATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGA         <  NB501016:230:HFVJLAFXY:4:11511:11407:16073/150‑1 (MQ=6)
                                                                                                                                               CATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGC                                                >  NB501016:230:HFVJLAFXY:4:21502:24939:15019/1‑111 (MQ=7)
                                                                                                                                               CATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGCCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGA         >  NB501016:230:HFVJLAFXY:2:21206:2782:8820/1‑150 (MQ=6)
                                                                                                                                                ATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGA         >  NB501016:230:HFVJLAFXY:4:11406:11460:18656/1‑149 (MQ=6)
                                                                                                                                                ATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGA         >  NB501016:230:HFVJLAFXY:1:21109:16023:17111/1‑149 (MQ=6)
                                                                                                                                                ATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGA         >  NB501016:230:HFVJLAFXY:1:21101:23128:1707/1‑149 (MQ=6)
                                                                                                                                                  GGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAagatcggaaagagcacacgtctgaactccagtcaccggc     >  NB501016:230:HFVJLAFXY:4:11510:11542:6667/1‑112 (MQ=5)
                                                                                                                                                  GGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGG      <  NB501016:230:HFVJLAFXY:4:11606:4490:19865/150‑1 (MQ=4)
                                                                                                                                                  GGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGG      <  NB501016:230:HFVJLAFXY:2:21310:11345:19307/150‑1 (MQ=4)
                                                                                                                                                  GGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCGCGATGCAAAATACACGAAGG      <  NB501016:230:HFVJLAFXY:1:21308:4424:8662/150‑1 (MQ=4)
                                                                                                                                                   GGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAG    <  NB501016:230:HFVJLAFXY:3:11608:7146:7353/151‑1 (MQ=4)
                                                                                                                                                   GGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAG    <  NB501016:230:HFVJLAFXY:4:11409:23288:11914/151‑1 (MQ=4)
                                                                                                                                                   GGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGGCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGCAGG      >  NB501016:230:HFVJLAFXY:3:21402:6206:17703/1‑149 (MQ=4)
                                                                                                                                                     CTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTT  >  NB501016:230:HFVJLAFXY:4:21509:17765:16696/1‑151 (MQ=4)
                                                                                                                                                     CTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTT  >  NB501016:230:HFVJLAFXY:1:11307:12285:17821/1‑151 (MQ=4)
                                                                                                                                                     |                                                                                                                                                      
CAATGGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTT  >  NC_000913/1423095‑1423394

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ 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: