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 F361 I0 R1 1353 250.7 8230166 96.0% 7900959 148.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,144,158 C→A 100% intergenic (+253/+19) yghA → / ← exbD putative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,144,1580CA99.5% 684.2 / ‑4.1 209intergenic (+253/+19)yghA/exbDputative oxidoreductase/membrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (114/94);  total (115/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 209 total aligned reads displayed.

TGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTG  >  CP009273/3144013‑3144303
                                                                                                                                                 |                                                                                                                                                 
tGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCa                                                                                                                                                <  2:2386163/149‑1 (MQ=255)
 gCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAt                                                                                                                                               >  1:3661893/1‑149 (MQ=255)
 gCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAt                                                                                                                                               >  2:10434/1‑149 (MQ=255)
 gCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAt                                                                                                                                               >  1:1614407/1‑149 (MQ=255)
   aaTAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCa                                                                                                                                             >  1:3804390/1‑149 (MQ=255)
       gATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCCTCAGGCa                                                                                                                                         >  2:3192968/1‑149 (MQ=255)
       gATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCa                                                                                                                                         >  2:887455/1‑149 (MQ=255)
       gATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCa                                                                                                                                         >  2:3193068/1‑149 (MQ=255)
        aTTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCaa                                                                                                                                        >  2:1210267/1‑149 (MQ=255)
         ttttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAt                                                                                                                                       <  1:3500364/149‑1 (MQ=255)
         ttttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAt                                                                                                                                       <  1:1879302/149‑1 (MQ=255)
         ttttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAt                                                                                                                                       <  1:3289454/149‑1 (MQ=255)
          tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAtt                                                                                                                                      <  2:8374/149‑1 (MQ=255)
          tttGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAAtt                                                                                                                                      <  1:2600263/149‑1 (MQ=255)
             gTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGACGATTCAACATGTTGCAATAGCTTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCta                                                                                                                                   >  2:3041129/1‑149 (MQ=255)
             gTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCta                                                                                                                                   <  1:143292/149‑1 (MQ=255)
               aGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtact                                                                                                                                 >  2:3769549/1‑149 (MQ=255)
                  ccTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtacttta                                                                                                                              <  1:632517/149‑1 (MQ=255)
                  ccTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtacttta                                                                                                                              <  2:1340357/149‑1 (MQ=255)
                   cTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGACCGCCTTATCCGGCGTACAAACGCACGAAGCTTCAACATGTTGCACTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCtactttac                                                                                                                             >  2:1839961/1‑149 (MQ=38)
                       tAAGCGTAGCGCATCAGGCATTGAGCCCCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGTGCGTCGAGCATCAGGCAATTCTCCTGTACGTCg                                                                                                                         >  2:2497363/1‑149 (MQ=25)
                       tAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCg                                                                                                                         >  1:3339852/1‑149 (MQ=255)
                         aGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACATTACTTCGCt                                                                                                                       >  2:220899/1‑149 (MQ=255)
                          gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                      >  1:1901883/1‑149 (MQ=255)
                          gCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                      >  1:1303955/1‑149 (MQ=255)
                          gCGGAGCGCAACAGGCATTGAGCACCTACTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGAGGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCtt                                                                                                                      <  1:4013193/149‑1 (MQ=25)
                           cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCttt                                                                                                                     <  1:2254481/149‑1 (MQ=255)
                           cGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCttt                                                                                                                     <  2:3591514/149‑1 (MQ=255)
                            gTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTg                                                                                                                    <  2:615812/149‑1 (MQ=255)
                               gcgcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCg                                                                                                                 >  1:3740271/1‑149 (MQ=255)
                               gcgcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCg                                                                                                                 >  1:1104323/1‑149 (MQ=255)
                                cgcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTtggtgg                                                                                                                <  2:2983413/149‑1 (MQ=255)
                                 gcctaaGGCATTGAGCACATTATGACGGGTGCTTCGTTAACGGCTTTTCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGt                                                                                                               <  2:380487/144‑1 (MQ=21)
                                 gcATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGt                                                                                                               >  1:1187324/1‑149 (MQ=255)
                                   aTCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGttt                                                                                                             >  1:2848698/1‑149 (MQ=255)
                                    tCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTAACTTCGCTTTGGCGGTttc                                                                                                            >  2:1455102/1‑149 (MQ=255)
                                     cAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttct                                                                                                           <  1:3102094/149‑1 (MQ=255)
                                       ggCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttcttc                                                                                                         <  2:97074/149‑1 (MQ=255)
                                       ggCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGAGTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTttcttc                                                                                                         <  1:3050945/149‑1 (MQ=255)
                                          aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGcc                                                                                                      >  2:2327418/1‑149 (MQ=255)
                                          aTTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGcc                                                                                                      >  1:2435175/1‑149 (MQ=255)
                                           ttGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCg                                                                                                     >  1:2948537/1‑149 (MQ=255)
                                               gCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCa                                                                                                 >  1:1989307/1‑149 (MQ=255)
                                                 aCCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGa                                                                                               <  1:2553344/149‑1 (MQ=255)
                                                 aCCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGa                                                                                               <  1:2552275/149‑1 (MQ=255)
                                                    tAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCt                                                                                            >  2:4014710/1‑149 (MQ=255)
                                                     aaTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTa                                                                                           <  1:212133/149‑1 (MQ=255)
                                                     aaTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTa                                                                                           <  2:68292/149‑1 (MQ=255)
                                                       tGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATc                                                                                         <  2:2577765/149‑1 (MQ=255)
                                                        gCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCt                                                                                        <  1:3072800/149‑1 (MQ=255)
                                                        gCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCt                                                                                        <  1:1125789/149‑1 (MQ=255)
                                                            gATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAg                                                                                    >  1:1752075/1‑149 (MQ=255)
                                                                 gtcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAc                                                                               >  2:3069929/1‑149 (MQ=255)
                                                                  tcgtGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAAcc                                                                              >  1:571403/1‑149 (MQ=255)
                                                                       aaCGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCt                                                                         <  1:3218026/149‑1 (MQ=255)
                                                                          gCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGAt                                                                      >  2:3290884/1‑149 (MQ=255)
                                                                            cTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGc                                                                    <  1:1210267/149‑1 (MQ=255)
                                                                            cTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACGTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGc                                                                    <  2:3248759/149‑1 (MQ=255)
                                                                             ttATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCa                                                                   <  2:344710/149‑1 (MQ=255)
                                                                             ttATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCa                                                                   <  2:344820/149‑1 (MQ=255)
                                                                                tCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTCCTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCg                                                                >  2:1996497/1‑149 (MQ=255)
                                                                                     ccTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATcc                                                           <  2:288657/149‑1 (MQ=255)
                                                                                        aCAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAtt                                                        <  1:2405501/149‑1 (MQ=255)
                                                                                        aCAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCAtt                                                        <  1:1469965/149‑1 (MQ=255)
                                                                                          aaaCGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTAc                                                      <  2:454981/149‑1 (MQ=255)
                                                                                            aCGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCt                                                    <  2:1240131/149‑1 (MQ=255)
                                                                                              gCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtc                                                  <  1:2297532/149‑1 (MQ=255)
                                                                                               cACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTtca                                                 >  2:3692351/1‑149 (MQ=255)
                                                                                                   aaGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCaa                                             >  1:1364904/1‑149 (MQ=255)
                                                                                                    aGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAAc                                            >  2:1504997/1‑149 (MQ=255)
                                                                                                    aGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAAc                                            <  1:2988120/149‑1 (MQ=255)
                                                                                                      aTTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATGAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATCACCTTCATCAACGt                                          <  1:884064/149‑1 (MQ=255)
                                                                                                      aTTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGt                                          >  1:2966222/1‑149 (MQ=255)
                                                                                                      aTTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGt                                          <  1:928865/149‑1 (MQ=255)
                                                                                                       ttCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtc                                         >  2:2225648/1‑149 (MQ=255)
                                                                                                       ttCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtc                                         >  1:511107/1‑149 (MQ=255)
                                                                                                       ttCAACATGTTGCAATAGATGTTGTAGGCCTGATGAGCGTAGAGCATCAGGCACGTCTACGTTACGTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTCTCTTCAGGGACCCCGCCTGATGCAGCGTATCCATTACCTGCAGCAACGtc                                         >  2:1961012/1‑149 (MQ=37)
                                                                                                         cAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtctc                                       >  1:1041630/1‑149 (MQ=255)
                                                                                                          aaCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGtctct                                      >  2:3534347/1‑148 (MQ=255)
                                                                                                           aCATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGt                                     >  2:168735/1‑149 (MQ=255)
                                                                                                            cATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTa                                    >  1:3707798/1‑149 (MQ=255)
                                                                                                              tGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAAt                                  <  1:3079253/149‑1 (MQ=255)
                                                                                                               gTTGCCAGAGATTTTGGAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATc                                 <  1:3376036/149‑1 (MQ=255)
                                                                                                                   cAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACg                             <  2:1901883/149‑1 (MQ=255)
                                                                                                                    aaTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGTCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACgg                            >  2:1535040/1‑149 (MQ=255)
                                                                                                                    aaTAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACgg                            >  1:810834/1‑149 (MQ=255)
                                                                                                                      tAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtt                          <  1:3311354/149‑1 (MQ=255)
                                                                                                                        gATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                        <  2:1429128/149‑1 (MQ=255)
                                                                                                                        gATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGtttt                        <  2:3930220/149‑1 (MQ=255)
                                                                                                                         aTTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTg                       <  1:4086979/149‑1 (MQ=255)
                                                                                                                           tttGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTc                     <  1:619553/149‑1 (MQ=255)
                                                                                                                                 ggCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCg               >  1:3278939/1‑149 (MQ=255)
                                                                                                                                  gCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCGTCAACGTCTCGTAATCGACGGTTTTGTCCGCTCgg              >  2:3994506/1‑149 (MQ=255)
                                                                                                                                      gatgaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaaga          >  1:866370/1‑149 (MQ=255)
                                                                                                                                       atgaGCGTAGAGCATCAGGAAATTCTAGTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGgaagaa         <  1:361093/149‑1 (MQ=255)
                                                                                                                                         gaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAagaaga       <  1:3192971/149‑1 (MQ=255)
                                                                                                                                         gaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAagaaga       <  1:3193075/149‑1 (MQ=255)
                                                                                                                                         gaGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAagaaga       >  1:3562650/1‑149 (MQ=255)
                                                                                                                                           gCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtg     >  2:2995601/1‑149 (MQ=255)
                                                                                                                                             gTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAtggt   <  1:2680344/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTG  >  CP009273/3144013‑3144303

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGTGT  >  CP009273/3144011‑3144307
                                                                                                                                                   |                                                                                                                                                     
GTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGC                                                                                                                                                     <  SRR19973068.4030020/150‑1 (MQ=60)
  TGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAT                                                                                                                                                   >  SRR19973068.3684040/1‑150 (MQ=60)
  TGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCAT                                                                                                                                                   >  SRR19973068.1624971/1‑150 (MQ=60)
    CAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCA                                                                                                                                                 >  SRR19973068.3827212/1‑150 (MQ=60)
     AATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAG                                                                                                                                                <  SRR19973068.2357584/150‑1 (MQ=60)
     AATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAG                                                                                                                                                <  SRR19973068.2357201/150‑1 (MQ=60)
     AATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAG                                                                                                                                                <  SRR19973068.2357190/150‑1 (MQ=60)
     AATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAG                                                                                                                                                <  SRR19973068.2357316/150‑1 (MQ=60)
     AATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAG                                                                                                                                                <  SRR19973068.2357583/150‑1 (MQ=60)
           TTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATT                                                                                                                                          <  SRR19973068.1891764/150‑1 (MQ=60)
           TTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATT                                                                                                                                          <  SRR19973068.3521637/150‑1 (MQ=60)
           TTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATT                                                                                                                                          <  SRR19973068.3309622/150‑1 (MQ=60)
            TTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTC                                                                                                                                         <  SRR19973068.2616673/150‑1 (MQ=60)
               GTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTAC                                                                                                                                      <  SRR19973068.144552/150‑1 (MQ=60)
                   GCCTGATAAGCGTGGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTA                                                                                                                                  <  SRR19973068.3195322/150‑1 (MQ=60)
                   GCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTA                                                                                                                                  <  SRR19973068.3213207/150‑1 (MQ=60)
                    CCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTAC                                                                                                                                 >  SRR19973068.3249526/1‑150 (MQ=60)
                    CCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTAC                                                                                                                                 <  SRR19973068.636688/150‑1 (MQ=60)
                      TGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTT                                                                                                                               <  SRR19973068.845896/150‑1 (MQ=60)
                        ATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCG                                                                                                                             >  SRR19973068.3360214/1‑150 (MQ=60)
                          AAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCT                                                                                                                           <  SRR19973068.2989014/150‑1 (MQ=60)
                           AGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTT                                                                                                                          >  SRR19973068.1313147/1‑150 (MQ=60)
                           AGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTT                                                                                                                          >  SRR19973068.1914515/1‑150 (MQ=60)
                           AGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTT                                                                                                                          >  SRR19973068.1312496/1‑150 (MQ=60)
                           AGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTT                                                                                                                          >  SRR19973068.1914503/1‑150 (MQ=60)
                            GCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTT                                                                                                                         <  SRR19973068.4037707/150‑1 (MQ=60)
                            GCGGAGCGCAACAGGCATTGAGCACCTACTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGAGGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTC                                                                                                                         <  SRR19973068.4037429/150‑1 (MQ=60)
                             CGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTG                                                                                                                        <  SRR19973068.2269150/150‑1 (MQ=60)
                                AGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCG                                                                                                                     >  SRR19973068.1111678/1‑150 (MQ=60)
                                AGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCG                                                                                                                     >  SRR19973068.3762682/1‑150 (MQ=60)
                                  CGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGT                                                                                                                   >  SRR19973068.1195257/1‑150 (MQ=60)
                                  CGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGT                                                                                                                   >  SRR19973068.1195323/1‑150 (MQ=60)
                                   GCATCAGGCATTGAGCACCTACTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTT                                                                                                                  <  SRR19973068.1101354/150‑1 (MQ=60)
                                    CATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTT                                                                                                                 >  SRR19973068.2866445/1‑150 (MQ=60)
                                     ATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTC                                                                                                                <  SRR19973068.1124311/150‑1 (MQ=60)
                                     ATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTC                                                                                                                <  SRR19973068.1614084/150‑1 (MQ=60)
                                     ATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTC                                                                                                                >  SRR19973068.2982423/1‑150 (MQ=60)
                                       CAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTT                                                                                                              <  SRR19973068.3120983/150‑1 (MQ=60)
                                       CAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTT                                                                                                              <  SRR19973068.3120563/150‑1 (MQ=60)
                                       CAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTT                                                                                                              <  SRR19973068.2947844/150‑1 (MQ=60)
                                         GGCATTGAGCACCTACTGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCG                                                                                                            <  SRR19973068.3132737/150‑1 (MQ=60)
                                         GGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCG                                                                                                            <  SRR19973068.1065226/150‑1 (MQ=60)
                                         GGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGAGTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCG                                                                                                            <  SRR19973068.3069639/150‑1 (MQ=60)
                                         GGCATTGAGCACCTAAGGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCG                                                                                                            <  SRR19973068.3069652/150‑1 (MQ=60)
                                           CATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCC                                                                                                          >  SRR19973068.2450657/1‑150 (MQ=60)
                                           CATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCC                                                                                                          >  SRR19973068.2448527/1‑150 (MQ=60)
                                            ATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCG                                                                                                         >  SRR19973068.2966831/1‑150 (MQ=60)
                                                AGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCA                                                                                                     >  SRR19973068.2002440/1‑150 (MQ=60)
                                                AGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCA                                                                                                     >  SRR19973068.1375406/1‑150 (MQ=60)
                                                 GCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAG                                                                                                    <  SRR19973068.649720/150‑1 (MQ=60)
                                                   ACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGAC                                                                                                  <  SRR19973068.2568461/150‑1 (MQ=60)
                                                   ACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGAC                                                                                                  <  SRR19973068.2569478/150‑1 (MQ=60)
                                                   ACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGAC                                                                                                  <  SRR19973068.2569532/150‑1 (MQ=60)
                                                      TAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTA                                                                                               >  SRR19973068.2292593/1‑150 (MQ=60)
                                                       AATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTAT                                                                                              <  SRR19973068.213739/150‑1 (MQ=60)
                                                          GCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTT                                                                                           <  SRR19973068.1133307/150‑1 (MQ=60)
                                                          GCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTT                                                                                           <  SRR19973068.3091560/150‑1 (MQ=60)
                                                          GCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTT                                                                                           <  SRR19973068.3944882/150‑1 (MQ=60)
                                                             GGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAG                                                                                        >  SRR19973068.1896536/1‑150 (MQ=60)
                                                             GGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAG                                                                                        >  SRR19973068.1763304/1‑150 (MQ=60)
                                                             GGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAG                                                                                        >  SRR19973068.3816072/1‑150 (MQ=60)
                                                              GATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGG                                                                                       <  SRR19973068.10596/150‑1 (MQ=60)
                                                                TGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTA                                                                                     >  SRR19973068.2987749/1‑150 (MQ=60)
                                                                   GTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACC                                                                                  >  SRR19973068.575321/1‑150 (MQ=60)
                                                                         AACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTG                                                                            >  SRR19973068.2027462/1‑150 (MQ=60)
                                                                         AACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTG                                                                            <  SRR19973068.3235784/150‑1 (MQ=60)
                                                                         AACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTG                                                                            <  SRR19973068.3237806/150‑1 (MQ=60)
                                                                              CTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCA                                                                       <  SRR19973068.977539/150‑1 (MQ=60)
                                                                              CTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCA                                                                       <  SRR19973068.1218333/150‑1 (MQ=60)
                                                                                     GGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATC                                                                >  SRR19973068.1163768/1‑150 (MQ=60)
                                                                                     GGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATC                                                                >  SRR19973068.123202/1‑150 (MQ=60)
                                                                                      GCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCC                                                               >  SRR19973068.916318/1‑150 (MQ=60)
                                                                                         TACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATT                                                            <  SRR19973068.3092076/150‑1 (MQ=60)
                                                                                          ACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTA                                                           <  SRR19973068.1479623/150‑1 (MQ=60)
                                                                                          ACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTA                                                           <  SRR19973068.2420773/150‑1 (MQ=60)
                                                                                             AACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCT                                                        <  SRR19973068.3486889/150‑1 (MQ=60)
                                                                                                GCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCA                                                     <  SRR19973068.2312392/150‑1 (MQ=60)
                                                                                                   CGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCA                                                  <  SRR19973068.2538515/150‑1 (MQ=60)
                                                                                                    GAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAA                                                 >  SRR19973068.1373644/1‑150 (MQ=60)
                                                                                                     AAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAAC                                                >  SRR19973068.2609101/1‑150 (MQ=60)
                                                                                                      ANATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACG                                               >  SRR19973068.3732078/1‑150 (MQ=60)
                                                                                                      AGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACG                                               <  SRR19973068.3006600/150‑1 (MQ=60)
                                                                                                      AGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACG                                               >  SRR19973068.3824932/1‑150 (MQ=60)
                                                                                                      AGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACG                                               <  SRR19973068.112457/150‑1 (MQ=60)
                                                                                                       GATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGT                                              >  SRR19973068.2984609/1‑150 (MQ=60)
                                                                                                        ATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATGAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATCACCTTCATCAACGTC                                             <  SRR19973068.890048/150‑1 (MQ=60)
                                                                                                        ATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTC                                             >  SRR19973068.514660/1‑150 (MQ=60)
                                                                                                        ATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTC                                             >  SRR19973068.514670/1‑150 (MQ=60)
                                                                                                        ATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTC                                             <  SRR19973068.935097/150‑1 (MQ=60)
                                                                                                          TCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTC                                           >  SRR19973068.1048303/1‑150 (MQ=60)
                                                                                                           CAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCG                                          <  SRR19973068.3270661/150‑1 (MQ=60)
                                                                                                             ACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTA                                        >  SRR19973068.3730101/1‑150 (MQ=60)
                                                                                                             ACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTA                                        >  SRR19973068.3730095/1‑150 (MQ=60)
                                                                                                             ACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTA                                        >  SRR19973068.1802956/1‑150 (MQ=60)
                                                                                                               ATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAAT                                      >  SRR19973068.1978812/1‑150 (MQ=60)
                                                                                                                TGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATC                                     <  SRR19973068.3098032/150‑1 (MQ=60)
                                                                                                                TGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATC                                     >  SRR19973068.1748274/1‑150 (MQ=60)
                                                                                                                 GTTGCCAGAGATTTTGGAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCG                                    <  SRR19973068.3396523/150‑1 (MQ=60)
                                                                                                                 GTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCG                                    <  SRR19973068.3396199/150‑1 (MQ=60)
                                                                                                                 GTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCG                                    <  SRR19973068.2251222/150‑1 (MQ=60)
                                                                                                                   TGCAATAGATTTTGTAGGCCTGATGAGCTTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGAC                                  <  SRR19973068.1473112/150‑1 (MQ=60)
                                                                                                                     CAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGG                                >  SRR19973068.3026762/1‑150 (MQ=60)
                                                                                                                     CAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGG                                >  SRR19973068.816357/1‑150 (MQ=60)
                                                                                                                     CAATAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGG                                >  SRR19973068.367534/1‑150 (MQ=60)
                                                                                                                        TAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTT                             <  SRR19973068.1564084/150‑1 (MQ=60)
                                                                                                                        TAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTT                             <  SRR19973068.3332582/150‑1 (MQ=60)
                                                                                                                        TAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTT                             <  SRR19973068.3331593/150‑1 (MQ=60)
                                                                                                                        TAGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTT                             <  SRR19973068.3331485/150‑1 (MQ=60)
                                                                                                                         AGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTT                            <  SRR19973068.2703999/150‑1 (MQ=60)
                                                                                                                         AGATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTT                            <  SRR19973068.2703986/150‑1 (MQ=60)
                                                                                                                           ATTTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGT                          <  SRR19973068.4111487/150‑1 (MQ=60)
                                                                                                                             TTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCC                        <  SRR19973068.623687/150‑1 (MQ=60)
                                                                                                                             TTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCC                        <  SRR19973068.623676/150‑1 (MQ=60)
                                                                                                                             TTTGTAGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCC                        <  SRR19973068.623330/150‑1 (MQ=60)
                                                                                                                                  AGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCG                   >  SRR19973068.3299076/1‑150 (MQ=60)
                                                                                                                                  AGGCCTGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCG                   >  SRR19973068.3299051/1‑150 (MQ=60)
                                                                                                                                       TGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGA              >  SRR19973068.3572148/1‑150 (MQ=60)
                                                                                                                                       TGATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGA              >  SRR19973068.872256/1‑150 (MQ=60)
                                                                                                                                         ATGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAG            <  SRR19973068.363273/150‑1 (MQ=60)
                                                                                                                                         ATGAGCGTAGAGCATCAGGAAATTCTAGTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAG            <  SRR19973068.363486/150‑1 (MQ=60)
                                                                                                                                          TGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGA           >  SRR19973068.438253/1‑150 (MQ=60)
                                                                                                                                          TGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGA           >  SRR19973068.3072967/1‑150 (MQ=60)
                                                                                                                                          TGAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGA           >  SRR19973068.3584289/1‑150 (MQ=60)
                                                                                                                                           GAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAT          <  SRR19973068.3212617/150‑1 (MQ=60)
                                                                                                                                           GAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAT          <  SRR19973068.3212620/150‑1 (MQ=60)
                                                                                                                                           GAGCGTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAT          <  SRR19973068.3212724/150‑1 (MQ=60)
                                                                                                                                           GAGCGTAGAGCATCAGGCAATCCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGAT          <  SRR19973068.3212717/150‑1 (MQ=60)
                                                                                                                                               GTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTG      <  SRR19973068.2697045/150‑1 (MQ=60)
                                                                                                                                               GTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTG      <  SRR19973068.2697046/150‑1 (MQ=60)
                                                                                                                                               GTAGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTG      >  SRR19973068.873259/1‑150 (MQ=60)
                                                                                                                                                 AGAGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGT    <  SRR19973068.1852693/150‑1 (MQ=60)
                                                                                                                                                   AGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGTGT  <  SRR19973068.570138/150‑1 (MQ=60)
                                                                                                                                                   AGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGTGT  >  SRR19973068.1660655/1‑150 (MQ=60)
                                                                                                                                                   AGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGTGT  >  SRR19973068.1660756/1‑150 (MQ=60)
                                                                                                                                                   |                                                                                                                                                     
GTTGCAATAGATTTTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTGAGCACCTAATGCCGGATGCGTCGTGAACGCCTTATCCGGCCTACAAACGCACGAAGATTCAACATGTTGCAATAGATTTTGTAGGCCTGATGAGCGTAGCGCATCAGGCAATTCTACTTTACTTCGCTTTGGCGGTTTCTTCGCCGACCAGACCTATCTTCAGGTAACCCGCCTGATGCAGCGTATCCATTACCTTCATCAACGTCTCGTAATCGACGGTTTTGTCCGCTCGGAAGAAGATGGTGGTGT  >  CP009273/3144011‑3144307

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