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 F368 I0 R1 1472 214.4 7116750 95.9% 6824963 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,727,580 C→A 100% A167A (GCC→GCA xylH → D‑xylose ABC transporter permease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,727,5800CA97.1% 536.1 / 6.7 173A167A (GCC→GCAxylHD‑xylose ABC transporter permease subunit
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (81/87);  total (84/89)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 178 total aligned reads displayed.

GGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTA  >  CP009273/3727436‑3727722
                                                                                                                                                |                                                                                                                                              
gggACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACcagcgcagcg                                                                                                                                             >  2:2982471/1‑148 (MQ=255)
gggACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACcagcgcagcg                                                                                                                                             >  2:2982445/1‑148 (MQ=255)
gggACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACcagcgcagcg                                                                                                                                             <  1:1420151/148‑1 (MQ=255)
  gACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGAt                                                                                                                                           >  1:3142438/1‑148 (MQ=255)
   aCTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATg                                                                                                                                          >  2:37514/1‑148 (MQ=255)
   aCTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATg                                                                                                                                          <  2:3441/148‑1 (MQ=255)
     tGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTc                                                                                                                                        >  2:149008/1‑148 (MQ=255)
     tGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTc                                                                                                                                        >  1:867702/1‑148 (MQ=255)
         tctcGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACaa                                                                                                                                    >  1:3519645/1‑148 (MQ=255)
         tctcGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACaa                                                                                                                                    >  1:1762793/1‑148 (MQ=255)
          ctcGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACaaa                                                                                                                                   >  1:2992228/1‑148 (MQ=255)
           tcGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAAt                                                                                                                                  <  2:1450369/148‑1 (MQ=255)
            cGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAAtt                                                                                                                                 >  1:3098789/1‑148 (MQ=255)
              gTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCAGGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTgg                                                                                                                               >  1:2949283/1‑148 (MQ=255)
                  cTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaa                                                                                                                           <  2:2251981/148‑1 (MQ=255)
                  cTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaa                                                                                                                           <  2:2245285/148‑1 (MQ=255)
                   tGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaaa                                                                                                                          <  2:2869943/148‑1 (MQ=255)
                   tGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaaa                                                                                                                          >  2:1719882/1‑148 (MQ=255)
                   tGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaaa                                                                                                                          <  1:333832/148‑1 (MQ=255)
                   tGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaaa                                                                                                                          >  2:2266499/1‑148 (MQ=255)
                   tGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaaa                                                                                                                          >  1:209868/1‑148 (MQ=255)
                   tGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaaa                                                                                                                          >  2:737645/1‑148 (MQ=255)
                   tGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCaaa                                                                                                                          >  2:737717/1‑148 (MQ=255)
                    ggAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAg                                                                                                                         <  2:3370448/148‑1 (MQ=255)
                      aaCGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCt                                                                                                                       <  2:2522752/148‑1 (MQ=255)
                      aaCGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTAGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCt                                                                                                                       <  1:1128181/148‑1 (MQ=255)
                       aCGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCCCCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCACCGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTa                                                                                                                      <  1:2110064/148‑1 (MQ=255)
                        cGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTAt                                                                                                                     <  2:1035425/148‑1 (MQ=255)
                           aTGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTAtctc                                                                                                                  >  2:1340382/1‑148 (MQ=255)
                             ggtggGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTccc                                                                                                                <  2:537884/148‑1 (MQ=255)
                               tggGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCg                                                                                                              >  2:3254081/1‑148 (MQ=255)
                                  gTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCa                                                                                                           <  1:1700654/148‑1 (MQ=255)
                                          ccGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGc                                                                                                   <  1:2337772/148‑1 (MQ=255)
                                            gTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCtt                                                                                                 <  1:371897/148‑1 (MQ=255)
                                             tAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTtc                                                                                                >  1:531877/1‑148 (MQ=255)
                                                 gTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTtcatca                                                                                            >  1:1769508/1‑148 (MQ=255)
                                                  tCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTcatcat                                                                                           <  2:2695558/148‑1 (MQ=255)
                                                  tCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTcatcat                                                                                           >  2:3320550/1‑148 (MQ=255)
                                                  tCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTcatcat                                                                                           >  2:3320143/1‑148 (MQ=255)
                                                    ccTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTg                                                                                         <  1:2857342/148‑1 (MQ=255)
                                                      ttCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGgc                                                                                       >  2:2206283/1‑148 (MQ=255)
                                                        cATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGgcgc                                                                                     >  2:3532509/1‑148 (MQ=255)
                                                            tATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCtt                                                                                 <  2:1430920/148‑1 (MQ=255)
                                                             aTTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTg                                                                                <  2:1506256/148‑1 (MQ=255)
                                                              ttGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTgg                                                                               >  1:877973/1‑148 (MQ=255)
                                                                 tCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCtt                                                                            >  1:3468161/1‑148 (MQ=255)
                                                                 tCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCtt                                                                            >  2:1960590/1‑148 (MQ=255)
                                                                 tCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCtt                                                                            >  1:2392958/1‑148 (MQ=255)
                                                                  cACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTa                                                                           <  2:2586742/148‑1 (MQ=255)
                                                                  cACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTa                                                                           >  1:2146183/1‑148 (MQ=255)
                                                                  cACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTa                                                                           <  1:776096/148‑1 (MQ=255)
                                                                   aCCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTaa                                                                          <  1:2237323/148‑1 (MQ=255)
                                                                      cTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATgg                                                                       <  2:3163402/148‑1 (MQ=255)
                                                                      cTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATgg                                                                       <  2:867702/148‑1 (MQ=255)
                                                                        cgcgGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCt                                                                     <  2:3061622/148‑1 (MQ=255)
                                                                        cgcgGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCt                                                                     >  1:1021375/1‑148 (MQ=255)
                                                                        cgcgGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCt                                                                     >  1:1023842/1‑148 (MQ=255)
                                                                            ggCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTttg                                                                 <  2:901888/148‑1 (MQ=255)
                                                                            ggCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTttg                                                                 >  1:3480215/1‑148 (MQ=255)
                                                                             gCATGTTGGCATTTCGCGGCATACTCATTGGCATCTCCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTttgt                                                                >  2:1859784/1‑148 (MQ=255)
                                                                              cATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTttgtt                                                               <  2:2059304/148‑1 (MQ=255)
                                                                                  ttGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTgttggtt                                                           >  1:3006541/1‑148 (MQ=255)
                                                                                   tGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTgttggttg                                                          <  2:2257592/148‑1 (MQ=255)
                                                                                     gCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGc                                                        <  1:2467430/148‑1 (MQ=255)
                                                                                       aTTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCaa                                                      <  1:3320143/148‑1 (MQ=255)
                                                                                       aTTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCaa                                                      <  1:3320550/148‑1 (MQ=255)
                                                                                        tttCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAAt                                                     <  1:1076024/148‑1 (MQ=255)
                                                                                        tttCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAAt                                                     <  2:1370377/148‑1 (MQ=255)
                                                                                         ttCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATg                                                    <  1:2665623/148‑1 (MQ=255)
                                                                                              ggCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGcgcg                                               <  1:1282451/148‑1 (MQ=255)
                                                                                              ggCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGcgcg                                               <  1:1817916/148‑1 (MQ=255)
                                                                                               gCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCgg                                              <  2:3284794/148‑1 (MQ=255)
                                                                                                 aTACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGgaa                                            >  2:852110/1‑148 (MQ=255)
                                                                                                    cTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGgaagaa                                         >  2:123978/1‑148 (MQ=255)
                                                                                                     tCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAAt                                        >  2:1891145/1‑148 (MQ=255)
                                                                                                     tCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAAt                                        >  1:167592/1‑148 (MQ=255)
                                                                                                        ttGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATgcg                                     >  2:2929851/1‑148 (MQ=255)
                                                                                                            cATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGt                                 >  2:3044386/1‑148 (MQ=255)
                                                                                                             aTCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTc                                >  1:1591755/1‑148 (MQ=255)
                                                                                                             aTCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTc                                >  1:1591795/1‑148 (MQ=255)
                                                                                                             aTCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTc                                <  2:2489201/148‑1 (MQ=255)
                                                                                                              tCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCa                               <  2:2174577/148‑1 (MQ=255)
                                                                                                               cACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAg                              >  1:1984337/1‑148 (MQ=255)
                                                                                                                 ccAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGc                            <  2:3519645/148‑1 (MQ=255)
                                                                                                                 ccAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGc                            <  2:1762793/148‑1 (MQ=255)
                                                                                                                   aaCGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCtt                          >  1:137093/1‑148 (MQ=255)
                                                                                                                   aaCGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCtt                          <  2:1176297/148‑1 (MQ=255)
                                                                                                                    aCGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCttt                         <  2:391100/148‑1 (MQ=255)
                                                                                                                        cACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGt                     <  1:14864/148‑1 (MQ=255)
                                                                                                                         aCGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGtt                    <  1:2781330/148‑1 (MQ=255)
                                                                                                                         aCGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTGGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGtt                    <  1:2450370/148‑1 (MQ=255)
                                                                                                                              tGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAg               <  1:888333/148‑1 (MQ=255)
                                                                                                                               gTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGt              <  2:2949283/148‑1 (MQ=255)
                                                                                                                                tATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGtc             <  2:1646313/148‑1 (MQ=255)
                                                                                                                                 aTCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGtct            >  1:2834138/1‑148 (MQ=255)
                                                                                                                                 aTCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGtct            >  1:2834146/1‑148 (MQ=255)
                                                                                                                                 aTCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTGGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGtct            <  2:1782722/148‑1 (MQ=255)
                                                                                                                                  tCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGtctc           <  2:2504666/148‑1 (MQ=255)
                                                                                                                                   cccccACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTcc          >  1:3237310/1‑148 (MQ=255)
                                                                                                                                           aGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTa  >  1:1229890/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTA  >  CP009273/3727436‑3727722

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
CTGGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTACCGCAGT  >  CP009273/3727434‑3727729
                                                                                                                                                  |                                                                                                                                                     
CTGGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCG                                                                                                                                                    <  SRR19973067.951262/150‑1 (MQ=60)
  GGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGAT                                                                                                                                                  <  SRR19973067.1429749/150‑1 (MQ=60)
  GGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGAT                                                                                                                                                  <  SRR19973067.3348651/150‑1 (MQ=60)
  GGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGAT                                                                                                                                                  >  SRR19973067.3162386/1‑150 (MQ=60)
     ACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTC                                                                                                                                               >  SRR19973067.873508/1‑150 (MQ=60)
         CTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAA                                                                                                                                           >  SRR19973067.1774786/1‑150 (MQ=60)
         CTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAA                                                                                                                                           >  SRR19973067.3541912/1‑150 (MQ=60)
          TTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAA                                                                                                                                          >  SRR19973067.3011185/1‑150 (MQ=60)
            CTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATT                                                                                                                                        >  SRR19973067.3118173/1‑150 (MQ=60)
            CTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATT                                                                                                                                        >  SRR19973067.3118506/1‑150 (MQ=60)
              CGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGG                                                                                                                                      >  SRR19973067.2105929/1‑150 (MQ=60)
              CGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCAGGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATACCCCACCAGCGCAGCGATGTCACAAATTGG                                                                                                                                      >  SRR19973067.2967946/1‑150 (MQ=60)
                  GCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAA                                                                                                                                  >  SRR19973067.1184078/1‑150 (MQ=60)
                   CCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAA                                                                                                                                 >  SRR19973067.211437/1‑150 (MQ=60)
                   CCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAA                                                                                                                                 >  SRR19973067.211453/1‑150 (MQ=60)
                     TGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGC                                                                                                                               <  SRR19973067.336118/150‑1 (MQ=60)
                        AACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTAGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTAT                                                                                                                            <  SRR19973067.1135796/150‑1 (MQ=60)
                         ACGGATGGGGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCCCCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCACCGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATC                                                                                                                           <  SRR19973067.2124090/150‑1 (MQ=60)
                          CGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCT                                                                                                                          <  SRR19973067.1498717/150‑1 (MQ=60)
                          CGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCT                                                                                                                          <  SRR19973067.1961948/150‑1 (MQ=60)
                             ATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCC                                                                                                                       <  SRR19973067.827296/150‑1 (MQ=60)
                                    GTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGT                                                                                                                <  SRR19973067.1712344/150‑1 (MQ=60)
                                            CCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTT                                                                                                        <  SRR19973067.2353079/150‑1 (MQ=60)
                                             CGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTC                                                                                                       >  SRR19973067.535442/1‑150 (MQ=60)
                                              GTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCA                                                                                                      <  SRR19973067.2811880/150‑1 (MQ=60)
                                              GTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCA                                                                                                      <  SRR19973067.374485/150‑1 (MQ=60)
                                                AAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATC                                                                                                    <  SRR19973067.1635410/150‑1 (MQ=60)
                                                 AAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCA                                                                                                   <  SRR19973067.1746265/150‑1 (MQ=60)
                                                 AAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCA                                                                                                   >  SRR19973067.1781524/1‑150 (MQ=60)
                                                      CCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGC                                                                                              <  SRR19973067.2875666/150‑1 (MQ=60)
                                                       CTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCG                                                                                             >  SRR19973067.907834/1‑150 (MQ=60)
                                                        TTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGC                                                                                            >  SRR19973067.3548250/1‑150 (MQ=60)
                                                          CATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGC                                                                                          >  SRR19973067.2189120/1‑150 (MQ=60)
                                                           ATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCT                                                                                         >  SRR19973067.591661/1‑150 (MQ=60)
                                                           ATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCT                                                                                         >  SRR19973067.3299059/1‑150 (MQ=60)
                                                             TTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCCCCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCGTG                                                                                       <  SRR19973067.1325716/150‑1 (MQ=60)
                                                              TATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGG                                                                                      >  SRR19973067.883828/1‑150 (MQ=60)
                                                                 TGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTT                                                                                   >  SRR19973067.3490636/1‑150 (MQ=60)
                                                                 TGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTT                                                                                   >  SRR19973067.3490201/1‑150 (MQ=60)
                                                                 TGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTT                                                                                   >  SRR19973067.2408530/1‑150 (MQ=60)
                                                                  GTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTA                                                                                  >  SRR19973067.2160530/1‑150 (MQ=60)
                                                                  GTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTA                                                                                  >  SRR19973067.1701410/1‑150 (MQ=60)
                                                                    CACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAAT                                                                                <  SRR19973067.781475/150‑1 (MQ=60)
                                                                     ACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATG                                                                               <  SRR19973067.2252204/150‑1 (MQ=60)
                                                                        CTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCT                                                                            >  SRR19973067.1030882/1‑150 (MQ=60)
                                                                        CTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCT                                                                            >  SRR19973067.1028392/1‑150 (MQ=60)
                                                                          CGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTT                                                                          <  SRR19973067.2553750/150‑1 (MQ=60)
                                                                            CGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTG                                                                        >  SRR19973067.3502302/1‑150 (MQ=60)
                                                                              GGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTT                                                                      >  SRR19973067.89164/1‑150 (MQ=60)
                                                                                  TGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTT                                                                  >  SRR19973067.3025646/1‑150 (MQ=60)
                                                                                   GTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTG                                                                 >  SRR19973067.457988/1‑150 (MQ=60)
                                                                                   GTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTG                                                                 <  SRR19973067.1469525/150‑1 (MQ=60)
                                                                                    TTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGG                                                                >  SRR19973067.3220063/1‑150 (MQ=60)
                                                                                     TGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGC                                                               >  SRR19973067.2689363/1‑150 (MQ=60)
                                                                                       GCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAA                                                             <  SRR19973067.2483653/150‑1 (MQ=60)
                                                                                         ATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATG                                                           >  SRR19973067.3460416/1‑150 (MQ=60)
                                                                                         ATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATG                                                           <  SRR19973067.3341550/150‑1 (MQ=60)
                                                                                         ATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATG                                                           <  SRR19973067.3341138/150‑1 (MQ=60)
                                                                                          TTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGG                                                          <  SRR19973067.1083400/150‑1 (MQ=60)
                                                                                           TTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGC                                                         <  SRR19973067.2682751/150‑1 (MQ=60)
                                                                                                GGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGA                                                    <  SRR19973067.1291191/150‑1 (MQ=60)
                                                                                                GGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGA                                                    <  SRR19973067.1830221/150‑1 (MQ=60)
                                                                                                    TACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAA                                                >  SRR19973067.2399884/1‑150 (MQ=60)
                                                                                                     ACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAAT                                               >  SRR19973067.168953/1‑150 (MQ=60)
                                                                                                        CATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCG                                            >  SRR19973067.1114561/1‑150 (MQ=60)
                                                                                                             GCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTC                                       >  SRR19973067.1602705/1‑150 (MQ=60)
                                                                                                             GCATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTC                                       >  SRR19973067.1602665/1‑150 (MQ=60)
                                                                                                               ATCACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAG                                     >  SRR19973067.1997768/1‑150 (MQ=60)
                                                                                                                 CACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGC                                   >  SRR19973067.1051818/1‑150 (MQ=60)
                                                                                                                 CACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGC                                   >  SRR19973067.2947868/1‑150 (MQ=60)
                                                                                                                  ACCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCT                                  <  SRR19973067.2290759/150‑1 (MQ=60)
                                                                                                                   CCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTT                                 >  SRR19973067.182974/1‑150 (MQ=60)
                                                                                                                   CCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTT                                 >  SRR19973067.138836/1‑150 (MQ=60)
                                                                                                                   CCAACGGCACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTT                                 >  SRR19973067.138314/1‑150 (MQ=60)
                                                                                                                          CACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTT                          <  SRR19973067.15097/150‑1 (MQ=60)
                                                                                                                           ACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTA                         <  SRR19973067.2799283/150‑1 (MQ=60)
                                                                                                                           ACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTGGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTA                         <  SRR19973067.2466489/150‑1 (MQ=60)
                                                                                                                           ACGACTGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTGGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTA                         <  SRR19973067.2466474/150‑1 (MQ=60)
                                                                                                                                TGTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTC                    <  SRR19973067.894226/150‑1 (MQ=60)
                                                                                                                                 GTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCT                   >  SRR19973067.2852358/1‑150 (MQ=60)
                                                                                                                                 GTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCT                   >  SRR19973067.2852350/1‑150 (MQ=60)
                                                                                                                                 GTATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCT                   >  SRR19973067.2852444/1‑150 (MQ=60)
                                                                                                                                   ATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCC                 >  SRR19973067.3257611/1‑150 (MQ=60)
                                                                                                                                   ATCCCCCACCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCC                 <  SRR19973067.3157505/150‑1 (MQ=60)
                                                                                                                                           CCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTA         >  SRR19973067.1238117/1‑150 (MQ=60)
                                                                                                                                           CCAGCGCAGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTA         >  SRR19973067.2987582/1‑150 (MQ=60)
                                                                                                                                                  AGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTACCGCAGT  >  SRR19973067.1759763/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
CTGGGACTGCTTCTCGGTGCCTGGAACGGATGGTGGGTCGCGTACCGTAAAGTCCCTTCATTTATTGTCACCCTCGCGGGCATGTTGGCATTTCGCGGCATACTCATTGGCATCACCAACGGCACGACTGTATCCCCCACCAGCGCCGCGATGTCACAAATTGGGCAAAGCTATCTCCCCGCCAGTACCGGCTTCATCATTGGCGCGCTTGGCTTAATGGCTTTTGTTGGTTGGCAATGGCGCGGAAGAATGCGCCGTCAGGCTTTGGGTTTACAGTCTCCGGCCTCTACCGCAGT  >  CP009273/3727434‑3727729

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: