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 F394 I0 R1 1342 215.0 7134094 96.2% 6862998 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,404,059 C→A 100% R141S (CGC→AGC)  dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,404,0590CA98.6% 695.5 / 0.7 216R141S (CGC→AGC) dusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (108/105);  total (109/107)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-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 216 total aligned reads displayed.

ATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG  >  CP009273/3403915‑3404203
                                                                                                                                                |                                                                                                                                                
attattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCa                                                                                                                                               <  2:3566161/148‑1 (MQ=255)
attattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCa                                                                                                                                               >  1:1855189/1‑148 (MQ=255)
 ttattGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCAc                                                                                                                                              <  1:3514823/148‑1 (MQ=255)
   attGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCg                                                                                                                                            >  1:731130/1‑148 (MQ=255)
    ttGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCgg                                                                                                                                           <  2:60781/148‑1 (MQ=255)
    ttGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCgg                                                                                                                                           >  1:2985431/1‑148 (MQ=255)
        tatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTgg                                                                                                                                       >  1:1556471/1‑148 (MQ=255)
        tatCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTgg                                                                                                                                       >  1:3165166/1‑148 (MQ=255)
          tCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGc                                                                                                                                     >  2:485178/1‑148 (MQ=255)
           cAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCa                                                                                                                                    <  1:3109939/148‑1 (MQ=255)
              tatGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCg                                                                                                                                 <  1:200671/148‑1 (MQ=255)
               atGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCgg                                                                                                                                >  2:2344994/1‑148 (MQ=255)
                 gggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGaa                                                                                                                              <  1:355392/148‑1 (MQ=255)
                  ggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAac                                                                                                                             <  1:1156865/148‑1 (MQ=255)
                  ggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAac                                                                                                                             >  2:2776188/1‑148 (MQ=255)
                  ggTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAac                                                                                                                             >  2:2174049/1‑148 (MQ=255)
                   gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAaca                                                                                                                            >  2:1819793/1‑148 (MQ=255)
                   gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAaca                                                                                                                            <  1:375477/148‑1 (MQ=255)
                    ttGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAacac                                                                                                                           >  1:791225/1‑148 (MQ=255)
                         cGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTaa                                                                                                                      <  2:622622/148‑1 (MQ=255)
                          ggCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAAc                                                                                                                     >  2:13008/1‑148 (MQ=255)
                           gCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACt                                                                                                                    >  1:3551935/1‑148 (MQ=255)
                             tAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGc                                                                                                                  >  2:3126636/1‑148 (MQ=255)
                              aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGAAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGAAACTGCg                                                                                                                 >  2:2155984/1‑148 (MQ=255)
                               aaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGa                                                                                                                <  2:899156/148‑1 (MQ=255)
                                 aaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAag                                                                                                              <  1:594747/148‑1 (MQ=255)
                                    gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAgaga                                                                                                           >  2:1826138/1‑148 (MQ=255)
                                      gAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAtt                                                                                                         >  2:711619/1‑148 (MQ=255)
                                        aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGc                                                                                                       <  1:1335198/148‑1 (MQ=255)
                                        aTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGc                                                                                                       >  1:779245/1‑148 (MQ=255)
                                         tCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGcc                                                                                                      <  2:2817638/148‑1 (MQ=255)
                                          cGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACTCCGTAACTGCGAAGAGATTGccc                                                                                                     <  1:641579/148‑1 (MQ=255)
                                          cGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGccc                                                                                                     >  1:252904/1‑148 (MQ=255)
                                            cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCaa                                                                                                   <  1:477317/148‑1 (MQ=255)
                                              aGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                 <  1:1900195/148‑1 (MQ=255)
                                                    cAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGa                                                                                           >  1:2448415/1‑148 (MQ=255)
                                                      ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAg                                                                                         >  2:693836/1‑148 (MQ=255)
                                                      ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAg                                                                                         <  1:2531940/148‑1 (MQ=255)
                                                      ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAg                                                                                         <  1:3496195/148‑1 (MQ=255)
                                                      ggTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAg                                                                                         >  2:3264464/1‑148 (MQ=255)
                                                             cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTgg                                                                                  >  1:2676031/1‑148 (MQ=255)
                                                             cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTgg                                                                                  >  2:2442006/1‑148 (MQ=255)
                                                             cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTgg                                                                                  >  2:2441296/1‑148 (MQ=255)
                                                             cccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTgg                                                                                  >  1:229175/1‑148 (MQ=255)
                                                              ccTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGc                                                                                 <  1:2184123/148‑1 (MQ=255)
                                                               ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCa                                                                                <  1:3206644/148‑1 (MQ=255)
                                                               ctctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCa                                                                                >  2:3044680/1‑148 (MQ=255)
                                                                tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAt                                                                               <  2:3436970/148‑1 (MQ=255)
                                                                tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAt                                                                               >  2:1203930/1‑148 (MQ=255)
                                                                 ctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                              >  1:171582/1‑148 (MQ=255)
                                                                 ctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                              <  1:991136/148‑1 (MQ=255)
                                                                   tGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCa                                                                            <  1:3493868/148‑1 (MQ=255)
                                                                     cAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAgg                                                                          >  2:2674926/1‑148 (MQ=255)
                                                                        tACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGctc                                                                       <  1:3391132/148‑1 (MQ=255)
                                                                             ggATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGAcc                                                                  <  2:1855189/148‑1 (MQ=255)
                                                                              gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCa                                                                 <  1:1769361/148‑1 (MQ=255)
                                                                              gATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCa                                                                 <  1:554630/148‑1 (MQ=255)
                                                                               aTGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAt                                                                >  1:3238555/1‑148 (MQ=255)
                                                                               aTGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAt                                                                <  1:1435775/148‑1 (MQ=255)
                                                                                tGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCCtt                                                               >  2:2774541/1‑148 (MQ=255)
                                                                                 gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTc                                                              >  2:3306834/1‑148 (MQ=255)
                                                                                 gTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTc                                                              >  2:279360/1‑148 (MQ=255)
                                                                                  tCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCa                                                             <  2:1093721/148‑1 (MQ=255)
                                                                                   cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAt                                                            >  1:1231345/1‑148 (MQ=255)
                                                                                   cGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAt                                                            >  1:3322502/1‑148 (MQ=255)
                                                                                        aaTCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGTACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCg                                                       >  1:2042299/1‑148 (MQ=255)
                                                                                          tCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTa                                                     >  2:2083711/1‑148 (MQ=255)
                                                                                                  tACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGcc                                             <  2:2066426/148‑1 (MQ=255)
                                                                                                    ccGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtg                                           >  1:2776699/1‑148 (MQ=255)
                                                                                                    ccGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtg                                           <  1:2798038/148‑1 (MQ=255)
                                                                                                      gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtt                                         <  1:3244216/148‑1 (MQ=255)
                                                                                                      gAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtt                                         <  1:2654096/148‑1 (MQ=255)
                                                                                                        ggTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgtttg                                       >  2:163348/1‑148 (MQ=255)
                                                                                                           cgtcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTc                                    <  2:748881/148‑1 (MQ=255)
                                                                                                             tcAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCaa                                  >  1:1934401/1‑148 (MQ=255)
                                                                                                              cAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAAt                                 >  2:2905028/1‑148 (MQ=255)
                                                                                                              cAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAAt                                 >  2:728865/1‑148 (MQ=255)
                                                                                                                aTGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATgg                               >  1:1856540/1‑148 (MQ=255)
                                                                                                                 tGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGga                              <  1:2377493/148‑1 (MQ=255)
                                                                                                                 tGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGga                              >  1:3261259/1‑148 (MQ=255)
                                                                                                                 tGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGga                              >  1:644141/1‑148 (MQ=255)
                                                                                                                  gCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGgag                             <  1:2530763/148‑1 (MQ=255)
                                                                                                                     gTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAg                          <  2:3509507/148‑1 (MQ=255)
                                                                                                                      tGGACGTTCCTGTTACCCAGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACAGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGATCTGACCATTAATGTCCGTACACGCGCCTGTTTGTTCAATGGAGAAGc                         >  2:2396146/1‑148 (MQ=255)
                                                                                                                        gACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTg                       <  2:3551945/148‑1 (MQ=255)
                                                                                                                        gACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTg                       <  2:3345421/148‑1 (MQ=255)
                                                                                                                         aCGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGa                      <  2:385521/148‑1 (MQ=255)
                                                                                                                            ttCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                   >  2:142064/1‑148 (MQ=255)
                                                                                                                             tCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                  >  2:2407424/1‑148 (MQ=255)
                                                                                                                                tGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGAc               <  2:1536882/148‑1 (MQ=255)
                                                                                                                                     cccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAt          <  1:485178/148‑1 (MQ=255)
                                                                                                                                      ccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt         <  2:2710251/148‑1 (MQ=255)
                                                                                                                                      ccTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAtt         <  2:226752/148‑1 (MQ=255)
                                                                                                                                       cTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTc        <  1:2747502/148‑1 (MQ=255)
                                                                                                                                        tGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCg       <  1:2772422/148‑1 (MQ=255)
                                                                                                                                         gAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCgg      >  1:1028736/1‑148 (MQ=255)
                                                                                                                                         gAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCgg      <  1:2794231/148‑1 (MQ=255)
                                                                                                                                           aGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGc    <  1:890877/148‑1 (MQ=255)
                                                                                                                                             aTTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAg  >  1:1310601/1‑148 (MQ=255)
                                                                                                                                             aTTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAg  >  2:2189034/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG  >  CP009273/3403915‑3404203

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAG  >  CP009273/3403910‑3404208
                                                                                                                                                     |                                                                                                                                                     
CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTA                                                                                                                                                       >  SRR19973066.441873/1‑150 (MQ=60)
CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTA                                                                                                                                                       >  SRR19973066.441771/1‑150 (MQ=60)
 CCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAG                                                                                                                                                      >  SRR19973066.2651409/1‑150 (MQ=60)
  CAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGC                                                                                                                                                     <  SRR19973066.3185798/150‑1 (MQ=60)
   AGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCA                                                                                                                                                    >  SRR19973066.1868463/1‑150 (MQ=60)
   AGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCA                                                                                                                                                    <  SRR19973066.895812/150‑1 (MQ=60)
   AGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCA                                                                                                                                                    <  SRR19973066.2693767/150‑1 (MQ=60)
     ATTATTGATATCAATATGGGTTGCCCGTCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACC                                                                                                                                                  <  SRR19973066.3307315/150‑1 (MQ=60)
      TTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCG                                                                                                                                                 >  SRR19973066.644956/1‑150 (MQ=60)
      TTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCG                                                                                                                                                 >  SRR19973066.736551/1‑150 (MQ=60)
      TTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCG                                                                                                                                                 <  SRR19973066.3538230/150‑1 (MQ=60)
      TTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCG                                                                                                                                                 <  SRR19973066.3538169/150‑1 (MQ=60)
       TATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGG                                                                                                                                                >  SRR19973066.3005163/1‑150 (MQ=60)
        ATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGC                                                                                                                                               >  SRR19973066.2105960/1‑150 (MQ=60)
          TGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTG                                                                                                                                             >  SRR19973066.1547844/1‑150 (MQ=60)
           GATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGG                                                                                                                                            >  SRR19973066.1567662/1‑150 (MQ=60)
           GATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGG                                                                                                                                            <  SRR19973066.2484250/150‑1 (MQ=60)
           GATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGG                                                                                                                                            >  SRR19973066.3186250/1‑150 (MQ=60)
             TATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGC                                                                                                                                          >  SRR19973066.2794148/1‑150 (MQ=60)
                CAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACC                                                                                                                                       <  SRR19973066.3130711/150‑1 (MQ=60)
                  ATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGG                                                                                                                                     >  SRR19973066.2025431/1‑150 (MQ=60)
                   TATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGA                                                                                                                                    <  SRR19973066.202374/150‑1 (MQ=60)
                   TATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGA                                                                                                                                    >  SRR19973066.1148800/1‑150 (MQ=60)
                      GGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACA                                                                                                                                 <  SRR19973066.714278/150‑1 (MQ=60)
                      GGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACA                                                                                                                                 <  SRR19973066.3207286/150‑1 (MQ=60)
                      GGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACA                                                                                                                                 <  SRR19973066.358154/150‑1 (MQ=60)
                       GGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACAC                                                                                                                                <  SRR19973066.1165090/150‑1 (MQ=60)
                       GGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACAC                                                                                                                                >  SRR19973066.797044/1‑150 (MQ=60)
                        GTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACC                                                                                                                               <  SRR19973066.348900/150‑1 (MQ=60)
                        GTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACC                                                                                                                               <  SRR19973066.378430/150‑1 (MQ=60)
                        GTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACC                                                                                                                               <  SRR19973066.378413/150‑1 (MQ=60)
                              CGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACT                                                                                                                         >  SRR19973066.3575504/1‑150 (MQ=60)
                              CGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACT                                                                                                                         >  SRR19973066.3575514/1‑150 (MQ=60)
                                      AAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAG                                                                                                                 <  SRR19973066.597342/150‑1 (MQ=60)
                                      AAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAG                                                                                                                 <  SRR19973066.598864/150‑1 (MQ=60)
                                       AAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGA                                                                                                                >  SRR19973066.3260002/1‑150 (MQ=60)
                                         GTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATT                                                                                                              >  SRR19973066.308215/1‑150 (MQ=60)
                                          TGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTG                                                                                                             <  SRR19973066.1323705/150‑1 (MQ=60)
                                           GAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGC                                                                                                            >  SRR19973066.784999/1‑150 (MQ=60)
                                             ATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCC                                                                                                          <  SRR19973066.1344889/150‑1 (MQ=60)
                                             ATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCC                                                                                                          >  SRR19973066.254833/1‑150 (MQ=60)
                                               CGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACTCCGTAACTGCGAAGAGATTGCCCAA                                                                                                        <  SRR19973066.646219/150‑1 (MQ=60)
                                                 CAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACT                                                                                                      <  SRR19973066.480917/150‑1 (MQ=60)
                                                 CAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACT                                                                                                      <  SRR19973066.479244/150‑1 (MQ=60)
                                                 CAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTCGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACT                                                                                                      <  SRR19973066.479250/150‑1 (MQ=60)
                                                   AGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGG                                                                                                    <  SRR19973066.98944/150‑1 (MQ=60)
                                                   AGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGG                                                                                                    <  SRR19973066.1913752/150‑1 (MQ=60)
                                                    GCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGC                                                                                                   <  SRR19973066.2658505/150‑1 (MQ=60)
                                                     CTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCT                                                                                                  <  SRR19973066.2773702/150‑1 (MQ=60)
                                                       CGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGA                                                                                                >  SRR19973066.2465099/1‑150 (MQ=60)
                                                       CGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGA                                                                                                >  SRR19973066.228574/1‑150 (MQ=60)
                                                         CAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAG                                                                                              >  SRR19973066.758750/1‑150 (MQ=60)
                                                           GGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGAC                                                                                            <  SRR19973066.2548648/150‑1 (MQ=60)
                                                           GGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGAC                                                                                            <  SRR19973066.2549153/150‑1 (MQ=60)
                                                           GGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGAC                                                                                            <  SRR19973066.3519503/150‑1 (MQ=60)
                                                           GGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGAC                                                                                            <  SRR19973066.3519420/150‑1 (MQ=60)
                                                                AGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGG                                                                                       >  SRR19973066.231004/1‑150 (MQ=60)
                                                                AGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGG                                                                                       >  SRR19973066.3195739/1‑150 (MQ=60)
                                                                AGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGG                                                                                       >  SRR19973066.1803335/1‑150 (MQ=60)
                                                                AGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGG                                                                                       >  SRR19973066.2693903/1‑150 (MQ=60)
                                                                AGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGG                                                                                       >  SRR19973066.1605973/1‑150 (MQ=60)
                                                                   CCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAT                                                                                    <  SRR19973066.2199277/150‑1 (MQ=60)
                                                                    CTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATT                                                                                   >  SRR19973066.173113/1‑150 (MQ=60)
                                                                    CTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATT                                                                                   <  SRR19973066.3227916/150‑1 (MQ=60)
                                                                      CTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCA                                                                                 >  SRR19973066.484192/1‑150 (MQ=60)
                                                                      CTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCA                                                                                 <  SRR19973066.998444/150‑1 (MQ=60)
                                                                        TGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGG                                                                               <  SRR19973066.3517168/150‑1 (MQ=60)
                                                                             TACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTG                                                                          <  SRR19973066.3413765/150‑1 (MQ=60)
                                                                                CCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACC                                                                       >  SRR19973066.3160232/1‑150 (MQ=60)
                                                                                  GGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAT                                                                     >  SRR19973066.3259931/1‑150 (MQ=60)
                                                                                   GATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATT                                                                    <  SRR19973066.3485846/150‑1 (MQ=60)
                                                                                   GATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATT                                                                    <  SRR19973066.1782038/150‑1 (MQ=60)
                                                                                   GATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATT                                                                    <  SRR19973066.558609/150‑1 (MQ=60)
                                                                                    ATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTC                                                                   <  SRR19973066.1446008/150‑1 (MQ=60)
                                                                                    ATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTC                                                                   <  SRR19973066.3228640/150‑1 (MQ=60)
                                                                                      GTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAT                                                                 >  SRR19973066.3344572/1‑150 (MQ=60)
                                                                                      GTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAT                                                                 >  SRR19973066.1240096/1‑150 (MQ=60)
                                                                                      GTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAT                                                                 >  SRR19973066.1959045/1‑150 (MQ=60)
                                                                                           TAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGTACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCG                                                            >  SRR19973066.2056516/1‑150 (MQ=60)
                                                                                                    CCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGC                                                   >  SRR19973066.727451/1‑150 (MQ=60)
                                                                                                       TACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTG                                                >  SRR19973066.2795391/1‑150 (MQ=60)
                                                                                                       TACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTG                                                >  SRR19973066.2798434/1‑150 (MQ=60)
                                                                                                       TACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTG                                                >  SRR19973066.3543739/1‑150 (MQ=60)
                                                                                                         CCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTT                                              <  SRR19973066.2818010/150‑1 (MQ=60)
                                                                                                         CCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTT                                              <  SRR19973066.2816873/150‑1 (MQ=60)
                                                                                                           GAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTG                                            <  SRR19973066.2671988/150‑1 (MQ=60)
                                                                                                           GAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTG                                            <  SRR19973066.2671901/150‑1 (MQ=60)
                                                                                                           GAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTG                                            <  SRR19973066.3265610/150‑1 (MQ=60)
                                                                                                           GAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTG                                            <  SRR19973066.792246/150‑1 (MQ=60)
                                                                                                             GGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTT                                          >  SRR19973066.1909938/1‑150 (MQ=60)
                                                                                                                CGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAA                                       >  SRR19973066.1948130/1‑150 (MQ=60)
                                                                                                                CGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAA                                       >  SRR19973066.1948950/1‑150 (MQ=60)
                                                                                                                   CAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGG                                    <  SRR19973066.3347866/150‑1 (MQ=60)
                                                                                                                   CAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGG                                    >  SRR19973066.1869695/1‑150 (MQ=60)
                                                                                                                   CAATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGG                                    >  SRR19973066.1869820/1‑150 (MQ=60)
                                                                                                                    AATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGA                                   >  SRR19973066.648808/1‑150 (MQ=60)
                                                                                                                    AATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGA                                   >  SRR19973066.648803/1‑150 (MQ=60)
                                                                                                                    AATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGA                                   >  SRR19973066.1669512/1‑150 (MQ=60)
                                                                                                                    AATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGA                                   >  SRR19973066.3282707/1‑150 (MQ=60)
                                                                                                                     ATGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAG                                  >  SRR19973066.3277761/1‑150 (MQ=60)
                                                                                                                      TGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGA                                 <  SRR19973066.2393544/150‑1 (MQ=60)
                                                                                                                      TGCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGA                                 >  SRR19973066.2314769/1‑150 (MQ=60)
                                                                                                                       GCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAA                                <  SRR19973066.2547674/150‑1 (MQ=60)
                                                                                                                       GCAGTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAA                                <  SRR19973066.2547972/150‑1 (MQ=60)
                                                                                                                          GTGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACAGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCT                             >  SRR19973066.1741557/1‑150 (MQ=60)
                                                                                                                           TGGACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTG                            >  SRR19973066.1409756/1‑150 (MQ=60)
                                                                                                                             GACGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAG                          <  SRR19973066.2737735/150‑1 (MQ=60)
                                                                                                                               CGTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTA                        >  SRR19973066.2505741/1‑150 (MQ=60)
                                                                                                                                GTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAC                       >  SRR19973066.1932233/1‑150 (MQ=60)
                                                                                                                                GTTCCTGTTACCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAC                       <  SRR19973066.1212304/150‑1 (MQ=60)
                                                                                                                                          CCCTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTC             <  SRR19973066.488821/150‑1 (MQ=60)
                                                                                                                                            CTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGG           >  SRR19973066.1036310/1‑150 (MQ=60)
                                                                                                                                            CTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGG           <  SRR19973066.2765990/150‑1 (MQ=60)
                                                                                                                                            CTGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGG           <  SRR19973066.2765988/150‑1 (MQ=60)
                                                                                                                                             TGAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGG          <  SRR19973066.2791089/150‑1 (MQ=60)
                                                                                                                                              GAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGC         <  SRR19973066.2813034/150‑1 (MQ=60)
                                                                                                                                              GAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGC         <  SRR19973066.2813042/150‑1 (MQ=60)
                                                                                                                                              GAAGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGC         >  SRR19973066.2441263/1‑150 (MQ=60)
                                                                                                                                                AGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG       >  SRR19973066.1320113/1‑150 (MQ=60)
                                                                                                                                                AGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG       >  SRR19973066.1320106/1‑150 (MQ=60)
                                                                                                                                                AGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG       <  SRR19973066.897184/150‑1 (MQ=60)
                                                                                                                                                AGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG       <  SRR19973066.897187/150‑1 (MQ=60)
                                                                                                                                                AGATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAG       >  SRR19973066.552396/1‑150 (MQ=60)
                                                                                                                                                  ATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTT     >  SRR19973066.2447097/1‑150 (MQ=60)
                                                                                                                                                  ATTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTT     >  SRR19973066.877656/1‑150 (MQ=60)
                                                                                                                                                   TTAGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTA    >  SRR19973066.579400/1‑150 (MQ=60)
                                                                                                                                                     AGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAG  >  SRR19973066.3576594/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
CCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAG  >  CP009273/3403910‑3404208

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: