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 F224 I0 R1 1457 284.0 9290224 95.7% 8890744 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,266,646 C→A 100% G91C (GGT→TGT)  garL ← alpha‑dehydro‑beta‑deoxy‑D‑glucarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,266,6460CA98.6% 966.6 / 4.0 282G91C (GGT→TGT) garLalpha‑dehydro‑beta‑deoxy‑D‑glucarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base A (158/120);  total (161/121)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
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 282 total aligned reads displayed.

CGAAATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACCGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAG  >  CP009273/3266499‑3266791
                                                                                                                                                   |                                                                                                                                                 
    ataatCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATc                                                                                                                                             <  2:2193962/150‑1 (MQ=255)
     taatCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGACGTTATAGAAACAGATATcc                                                                                                                                            >  2:553945/1‑150 (MQ=255)
     taatCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATcc                                                                                                                                            >  1:585055/1‑150 (MQ=255)
      aatCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCa                                                                                                                                           >  2:3804136/1‑150 (MQ=255)
      aatCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCa                                                                                                                                           >  1:2798648/1‑150 (MQ=255)
       atCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCag                                                                                                                                          >  1:401881/1‑150 (MQ=255)
            ccACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGAc                                                                                                                                     <  2:2574443/150‑1 (MQ=255)
            ccACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGAc                                                                                                                                     >  2:1847452/1‑150 (MQ=255)
            ccACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCATCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGAc                                                                                                                                     <  1:565842/150‑1 (MQ=255)
             cACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACg                                                                                                                                    >  2:1660067/1‑150 (MQ=255)
              aCGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGc                                                                                                                                   >  2:586160/1‑150 (MQ=255)
               cGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCt                                                                                                                                  >  2:2238980/1‑150 (MQ=255)
                ggTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCtt                                                                                                                                 <  2:3002239/150‑1 (MQ=255)
                  tGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTtaa                                                                                                                               >  2:121409/1‑150 (MQ=255)
                    ccAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTtaata                                                                                                                             <  1:3999737/150‑1 (MQ=255)
                      aaaCATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTaataat                                                                                                                           >  2:4543444/1‑150 (MQ=255)
                      aaaCATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTaataat                                                                                                                           >  2:739269/1‑150 (MQ=255)
                        aCATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTa                                                                                                                         >  2:252092/1‑150 (MQ=255)
                          atatTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTAcc                                                                                                                       >  2:2334923/1‑150 (MQ=255)
                          atatTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTAcc                                                                                                                       >  1:200965/1‑150 (MQ=255)
                           tatTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCg                                                                                                                      >  2:397400/1‑150 (MQ=255)
                           tatTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCg                                                                                                                      >  1:942704/1‑150 (MQ=255)
                           tatTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCg                                                                                                                      <  1:4294020/150‑1 (MQ=255)
                               ggCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTc                                                                                                                  >  1:2387399/1‑150 (MQ=255)
                               ggCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTc                                                                                                                  >  1:1835211/1‑150 (MQ=255)
                               ggCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTAACGGcaa                                                                                                                  >  2:4257746/1‑148 (MQ=255)
                                  gcgGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGtt                                                                                                               >  2:3437005/1‑150 (MQ=255)
                                   cgGTGAGACACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTg                                                                                                              >  2:2663231/1‑150 (MQ=255)
                                       gagaAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCg                                                                                                          <  1:2703595/150‑1 (MQ=255)
                                         gaAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGc                                                                                                        >  1:1318158/1‑150 (MQ=255)
                                            aCGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACt                                                                                                     <  1:877718/150‑1 (MQ=255)
                                             cGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTc                                                                                                    >  1:4445077/1‑150 (MQ=255)
                                              ggAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCg                                                                                                   <  2:1275164/150‑1 (MQ=255)
                                                agaCGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCa                                                                                                 >  2:3769801/1‑150 (MQ=255)
                                                 gaCGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCac                                                                                                >  2:1116566/1‑150 (MQ=255)
                                                   cgccgcGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCacta                                                                                              <  1:4155118/150‑1 (MQ=255)
                                                        cgAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGgc                                                                                         <  1:2056266/150‑1 (MQ=255)
                                                        cgAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGgc                                                                                         <  1:460874/150‑1 (MQ=255)
                                                         gAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGgcg                                                                                        >  2:2109992/1‑150 (MQ=255)
                                                           aTGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGgcgcg                                                                                      <  2:875330/150‑1 (MQ=255)
                                                            tGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGgcgcgc                                                                                     <  1:3485893/150‑1 (MQ=255)
                                                             gCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgct                                                                                    >  2:2831655/1‑150 (MQ=255)
                                                             gCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgct                                                                                    <  1:3706512/150‑1 (MQ=255)
                                                             gCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgct                                                                                    <  2:81702/150‑1 (MQ=255)
                                                              ccTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctg                                                                                   <  1:1815282/150‑1 (MQ=255)
                                                              ccTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctg                                                                                   <  2:4621558/150‑1 (MQ=255)
                                                                ttCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctggc                                                                                 <  2:2934813/150‑1 (MQ=255)
                                                                 tCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctggcg                                                                                >  2:3633202/1‑150 (MQ=255)
                                                                  ccGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctggcgc                                                                               >  1:3868084/1‑150 (MQ=255)
                                                                   cGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctggcgct                                                                              <  1:4303370/150‑1 (MQ=255)
                                                                   cGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctggcgct                                                                              <  2:3674072/150‑1 (MQ=255)
                                                                   cGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctggcgct                                                                              <  1:1424727/150‑1 (MQ=255)
                                                                    ggTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCgcgctggcgctg                                                                             <  1:2887736/150‑1 (MQ=255)
                                                                      tGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGcc                                                                           >  2:4038403/1‑150 (MQ=255)
                                                                        ggTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCtt                                                                         >  2:693726/1‑150 (MQ=255)
                                                                         gTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCttt                                                                        <  2:2435215/150‑1 (MQ=255)
                                                                           aaCGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCa                                                                      >  2:3682524/1‑150 (MQ=255)
                                                                            aCGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCaa                                                                     <  2:2721581/150‑1 (MQ=255)
                                                                              gggaagATTCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCTTTGGTCTGAACTCGCACTACTGGCGCGCTGGCGCTGACTTTCAAgg                                                                   >  2:503552/6‑150 (MQ=255)
                                                                                  tGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCAt                                                               <  2:3294985/150‑1 (MQ=255)
                                                                                   gATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCAtt                                                              >  2:2756274/1‑150 (MQ=255)
                                                                                    aTGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTa                                                             <  2:1442801/150‑1 (MQ=255)
                                                                                      gCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAAc                                                           <  1:1721578/150‑1 (MQ=255)
                                                                                        cACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTg                                                         >  1:1712413/1‑150 (MQ=255)
                                                                                         aCCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGc                                                        <  1:1169705/150‑1 (MQ=255)
                                                                                         aCCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGc                                                        >  1:2042886/1‑150 (MQ=255)
                                                                                           cGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCgg                                                      >  2:4460020/1‑150 (MQ=255)
                                                                                            gCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGa                                                     >  1:4065227/1‑150 (MQ=255)
                                                                                              cAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAAt                                                   <  2:430985/150‑1 (MQ=255)
                                                                                              cAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAAt                                                   >  1:4186629/1‑150 (MQ=255)
                                                                                               aGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATa                                                  >  1:3138217/1‑150 (MQ=255)
                                                                                               aGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATa                                                  >  1:2168743/1‑150 (MQ=255)
                                                                                                   ctGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACg                                              >  2:3376351/1‑150 (MQ=255)
                                                                                                    tGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGt                                             >  1:3216812/1‑150 (MQ=255)
                                                                                                    tGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGt                                             <  1:2628587/150‑1 (MQ=255)
                                                                                                       ttCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGga                                          <  1:1237902/150‑1 (MQ=255)
                                                                                                       ttCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGga                                          <  1:3172284/150‑1 (MQ=255)
                                                                                                         ccTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGgaga                                        >  1:3726928/1‑150 (MQ=255)
                                                                                                            ctTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGatat                                     <  1:1764299/150‑1 (MQ=255)
                                                                                                                 gTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGttt                                <  2:3916461/150‑1 (MQ=255)
                                                                                                                    tCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGgc                             >  2:252964/1‑150 (MQ=255)
                                                                                                                    tCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGgc                             >  1:2520902/1‑150 (MQ=255)
                                                                                                                       aCAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCa                          <  2:2044419/150‑1 (MQ=255)
                                                                                                                        cAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCAt                         <  2:1975321/150‑1 (MQ=255)
                                                                                                                         aaaaGGAATCAGGAAGTTATAGAAACCGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATg                        >  2:1774763/1‑150 (MQ=255)
                                                                                                                         aaaaGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATg                        >  2:2313484/1‑150 (MQ=255)
                                                                                                                         aaaaGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATg                        >  1:541285/1‑150 (MQ=255)
                                                                                                                           aaGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGtt                      >  2:1554392/1‑150 (MQ=255)
                                                                                                                            aGGAATCAGGAAGTTAGAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTc                     <  2:1599823/150‑1 (MQ=255)
                                                                                                                             ggAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCg                    <  2:3603598/150‑1 (MQ=255)
                                                                                                                              gAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGc                   >  2:1904756/1‑150 (MQ=255)
                                                                                                                                 tCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCAt                >  2:1865009/1‑150 (MQ=255)
                                                                                                                                 tCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCAt                >  2:2024970/1‑150 (MQ=255)
                                                                                                                                   aGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATcc              >  1:1356847/1‑150 (MQ=255)
                                                                                                                                    ggAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCa             <  1:4543444/150‑1 (MQ=255)
                                                                                                                                     gAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAg            >  2:3946219/1‑150 (MQ=255)
                                                                                                                                          tataGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCAcca       <  1:2606081/150‑1 (MQ=255)
                                                                                                                                           ataGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCAccag      >  2:3740848/1‑150 (MQ=255)
                                                                                                                                           ataGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATAAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCAccag      >  1:3961945/1‑150 (MQ=255)
                                                                                                                                              gAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCAccagcca   >  1:3525839/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGAAATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACCGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAG  >  CP009273/3266499‑3266791

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
AGCGAAATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACCGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAGTCAA  >  CP009273/3266497‑3266795
                                                                                                                                                     |                                                                                                                                                     
AGCGAAATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACA                                                                                                                                                       >  SRR19973072.3396078/1‑150 (MQ=60)
   GAAATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGAT                                                                                                                                                    <  SRR19973072.1831360/150‑1 (MQ=60)
   GAAATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGAT                                                                                                                                                    <  SRR19973072.1831405/150‑1 (MQ=60)
      ATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATC                                                                                                                                                 <  SRR19973072.1793926/150‑1 (MQ=60)
       TAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCC                                                                                                                                                >  SRR19973072.586416/1‑150 (MQ=60)
       TAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCC                                                                                                                                                <  SRR19973072.4109655/150‑1 (MQ=60)
        AATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCA                                                                                                                                               <  SRR19973072.1747880/150‑1 (MQ=60)
        AATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCA                                                                                                                                               >  SRR19973072.2804634/1‑150 (MQ=60)
         ATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAG                                                                                                                                              >  SRR19973072.402736/1‑150 (MQ=60)
         ATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAG                                                                                                                                              >  SRR19973072.402749/1‑150 (MQ=60)
         ATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAG                                                                                                                                              >  SRR19973072.402767/1‑150 (MQ=60)
            CGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAG                                                                                                                                           <  SRR19973072.3694224/150‑1 (MQ=60)
            CGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAG                                                                                                                                           <  SRR19973072.1499359/150‑1 (MQ=60)
              CCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCATCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGAC                                                                                                                                         <  SRR19973072.567164/150‑1 (MQ=60)
               CACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACG                                                                                                                                        >  SRR19973072.4205201/1‑150 (MQ=60)
                   GTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGTAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTA                                                                                                                                    <  SRR19973072.4488515/150‑1 (MQ=60)
                      CCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATA                                                                                                                                 <  SRR19973072.1723140/150‑1 (MQ=60)
                      CCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATA                                                                                                                                 <  SRR19973072.4008115/150‑1 (MQ=60)
                      CCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATA                                                                                                                                 <  SRR19973072.4007592/150‑1 (MQ=60)
                            ATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACC                                                                                                                           >  SRR19973072.201468/1‑150 (MQ=60)
                            ATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACC                                                                                                                           >  SRR19973072.201445/1‑150 (MQ=60)
                             TATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCG                                                                                                                          >  SRR19973072.944763/1‑150 (MQ=60)
                             TATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCG                                                                                                                          <  SRR19973072.4302972/150‑1 (MQ=60)
                             TATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCG                                                                                                                          <  SRR19973072.1968515/150‑1 (MQ=60)
                               TTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGC                                                                                                                        <  SRR19973072.580511/150‑1 (MQ=60)
                                 GGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTC                                                                                                                      >  SRR19973072.1839199/1‑150 (MQ=60)
                                 GGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTC                                                                                                                      >  SRR19973072.4313422/1‑150 (MQ=60)
                                 GGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTC                                                                                                                      >  SRR19973072.2392562/1‑150 (MQ=60)
                                 GGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTC                                                                                                                      >  SRR19973072.2392332/1‑150 (MQ=60)
                                  GCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCG                                                                                                                     >  SRR19973072.1388778/1‑150 (MQ=60)
                                   CGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGT                                                                                                                    >  SRR19973072.3522907/1‑150 (MQ=60)
                                     CGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTG                                                                                                                  <  SRR19973072.388807/150‑1 (MQ=60)
                                      GGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGG                                                                                                                 >  SRR19973072.4311573/1‑150 (MQ=60)
                                      GGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGG                                                                                                                 >  SRR19973072.1352599/1‑150 (MQ=60)
                                       GTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGT                                                                                                                <  SRR19973072.4164651/150‑1 (MQ=60)
                                         GAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCG                                                                                                              <  SRR19973072.2709347/150‑1 (MQ=60)
                                         GAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTATTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCG                                                                                                              <  SRR19973072.2709248/150‑1 (MQ=60)
                                           GAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGC                                                                                                            >  SRR19973072.1321037/1‑150 (MQ=60)
                                             AACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCAC                                                                                                          <  SRR19973072.3586915/150‑1 (MQ=60)
                                              ACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACT                                                                                                         <  SRR19973072.879602/150‑1 (MQ=60)
                                               CGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTC                                                                                                        <  SRR19973072.2363066/150‑1 (MQ=60)
                                               CGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTC                                                                                                        >  SRR19973072.4454452/1‑150 (MQ=60)
                                                     CGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTA                                                                                                  <  SRR19973072.4163829/150‑1 (MQ=60)
                                                          CGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGC                                                                                             <  SRR19973072.461941/150‑1 (MQ=60)
                                                          CGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGC                                                                                             <  SRR19973072.2060741/150‑1 (MQ=60)
                                                              TGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGC                                                                                         <  SRR19973072.3493180/150‑1 (MQ=60)
                                                              TGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGC                                                                                         <  SRR19973072.2782433/150‑1 (MQ=60)
                                                               GCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCT                                                                                        <  SRR19973072.3714321/150‑1 (MQ=60)
                                                                CCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTG                                                                                       >  SRR19973072.2026264/1‑150 (MQ=60)
                                                                CCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTG                                                                                       >  SRR19973072.2616505/1‑150 (MQ=60)
                                                                CCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTG                                                                                       >  SRR19973072.2026284/1‑150 (MQ=60)
                                                                CCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTG                                                                                       <  SRR19973072.1819235/150‑1 (MQ=60)
                                                                 CTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTGAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGG                                                                                      >  SRR19973072.4548011/1‑150 (MQ=60)
                                                                    CCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGC                                                                                   >  SRR19973072.3876287/1‑150 (MQ=60)
                                                                    CCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGC                                                                                   >  SRR19973072.3876128/1‑150 (MQ=60)
                                                                     CGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCT                                                                                  <  SRR19973072.4312288/150‑1 (MQ=60)
                                                                     CGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCT                                                                                  <  SRR19973072.1427815/150‑1 (MQ=60)
                                                                     CGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCT                                                                                  <  SRR19973072.4312360/150‑1 (MQ=60)
                                                                      GGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTG                                                                                 <  SRR19973072.3229745/150‑1 (MQ=60)
                                                                      GGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTG                                                                                 <  SRR19973072.3229772/150‑1 (MQ=60)
                                                                      GGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTG                                                                                 <  SRR19973072.2893875/150‑1 (MQ=60)
                                                                      GGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTG                                                                                 <  SRR19973072.2893926/150‑1 (MQ=60)
                                                                           GTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTT                                                                            >  SRR19973072.3863220/1‑150 (MQ=60)
                                                                              ACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGAGCTGGCGCTGCCTTTCAA                                                                         <  SRR19973072.252007/150‑1 (MQ=60)
                                                                              ACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGAGCTGGCGCTGCCTTTCAA                                                                         <  SRR19973072.403264/150‑1 (MQ=60)
                                                                                GGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGG                                                                       <  SRR19973072.176330/150‑1 (MQ=60)
                                                                                 GGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGC                                                                      >  SRR19973072.3157037/1‑150 (MQ=60)
                                                                                  GTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCC                                                                     <  SRR19973072.71190/150‑1 (MQ=60)
                                                                                     GATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATT                                                                  >  SRR19973072.2674436/1‑150 (MQ=60)
                                                                                      ATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTA                                                                 >  SRR19973072.794106/1‑150 (MQ=60)
                                                                                        GCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAAC                                                               >  SRR19973072.569762/1‑150 (MQ=60)
                                                                                        GCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAAC                                                               <  SRR19973072.1725349/150‑1 (MQ=60)
                                                                                         CCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACT                                                              <  SRR19973072.2780218/150‑1 (MQ=60)
                                                                                         CCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACT                                                              <  SRR19973072.371540/150‑1 (MQ=60)
                                                                                         CCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACT                                                              <  SRR19973072.371622/150‑1 (MQ=60)
                                                                                          CACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTG                                                             >  SRR19973072.1716164/1‑150 (MQ=60)
                                                                                          CACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTG                                                             >  SRR19973072.2158043/1‑150 (MQ=60)
                                                                                           ACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGC                                                            <  SRR19973072.1172161/150‑1 (MQ=60)
                                                                                           ACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGC                                                            >  SRR19973072.2047324/1‑150 (MQ=60)
                                                                                           ACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGC                                                            >  SRR19973072.2047325/1‑150 (MQ=60)
                                                                                             CGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGG                                                          >  SRR19973072.552551/1‑150 (MQ=60)
                                                                                             CGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGG                                                          >  SRR19973072.552495/1‑150 (MQ=60)
                                                                                              GCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGA                                                         >  SRR19973072.4073757/1‑150 (MQ=60)
                                                                                               CCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGACACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAA                                                        >  SRR19973072.4441068/1‑150 (MQ=60)
                                                                                                CAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAAT                                                       >  SRR19973072.4195391/1‑150 (MQ=60)
                                                                                                 AGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATA                                                      >  SRR19973072.2173470/1‑150 (MQ=60)
                                                                                                 AGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATA                                                      >  SRR19973072.3144837/1‑150 (MQ=60)
                                                                                                   CTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAA                                                    >  SRR19973072.4002110/1‑150 (MQ=60)
                                                                                                      TGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGT                                                 >  SRR19973072.3223622/1‑150 (MQ=60)
                                                                                                      TGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGT                                                 <  SRR19973072.2634235/150‑1 (MQ=60)
                                                                                                       GCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTG                                                >  SRR19973072.929136/1‑150 (MQ=60)
                                                                                                         TTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACCGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGA                                              <  SRR19973072.3515619/150‑1 (MQ=60)
                                                                                                         TTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGA                                              <  SRR19973072.1240542/150‑1 (MQ=60)
                                                                                                         TTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGA                                              <  SRR19973072.1240726/150‑1 (MQ=60)
                                                                                                         TTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGA                                              <  SRR19973072.3178977/150‑1 (MQ=60)
                                                                                                         TTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGA                                              <  SRR19973072.3178972/150‑1 (MQ=60)
                                                                                                         TTCCTCTTTTGTTACTAAATAAGGAATCATGTAGTTATAGAAAAATATATCCAGAAGAATCTTAATAATTACCTGCTCGTTGGTCGTCACTCTCACTACTGGCTCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGA                                              <  SRR19973072.3179077/150‑1 (MQ=60)
                                                                                                           CCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGGGGAGA                                            >  SRR19973072.3734770/1‑150 (MQ=60)
                                                                                                              CTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATAT                                         <  SRR19973072.4047512/150‑1 (MQ=60)
                                                                                                              CTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATAT                                         <  SRR19973072.1768155/150‑1 (MQ=60)
                                                                                                               TTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATC                                        >  SRR19973072.1173948/1‑150 (MQ=60)
                                                                                                                TTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCG                                       >  SRR19973072.2596103/1‑150 (MQ=60)
                                                                                                                   GTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTT                                    >  SRR19973072.2931223/1‑150 (MQ=60)
                                                                                                                   GTTTCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTT                                    >  SRR19973072.2490414/1‑150 (MQ=60)
                                                                                                                      TCTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGC                                 >  SRR19973072.2526380/1‑150 (MQ=60)
                                                                                                                       CTACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCG                                >  SRR19973072.17954/1‑150 (MQ=60)
                                                                                                                         ACAAAAGGAATCAGGAAGTTATAGAAACCGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCA                              >  SRR19973072.2820272/1‑150 (MQ=60)
                                                                                                                         ACAAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCA                              >  SRR19973072.4111488/1‑150 (MQ=60)
                                                                                                                           AAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATG                            >  SRR19973072.542568/1‑150 (MQ=60)
                                                                                                                           AAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATG                            >  SRR19973072.542131/1‑150 (MQ=60)
                                                                                                                           AAAAGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATG                            >  SRR19973072.542125/1‑150 (MQ=60)
                                                                                                                              AGGAATCAGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTC                         <  SRR19973072.252699/150‑1 (MQ=60)
                                                                                                                                     AGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCC                  <  SRR19973072.2121488/150‑1 (MQ=60)
                                                                                                                                     AGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCC                  >  SRR19973072.1359817/1‑150 (MQ=60)
                                                                                                                                     AGGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCC                  >  SRR19973072.1359923/1‑150 (MQ=60)
                                                                                                                                      GGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCA                 <  SRR19973072.4553027/150‑1 (MQ=60)
                                                                                                                                      GGAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCA                 <  SRR19973072.4553405/150‑1 (MQ=60)
                                                                                                                                       GAAGTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAAAGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAG                >  SRR19973072.2969776/1‑150 (MQ=60)
                                                                                                                                         AGTTATAGAAACAGATATCCAGAAGACGCTTAATAGTTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGGGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCA              >  SRR19973072.1602472/1‑150 (MQ=60)
                                                                                                                                          GTTATAGAAACAGATATCCAGAAGACGCTTNATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCAC             >  SRR19973072.4292686/1‑150 (MQ=60)
                                                                                                                                          GTTATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCAC             >  SRR19973072.1333516/1‑150 (MQ=60)
                                                                                                                                            TATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCA           <  SRR19973072.2611700/150‑1 (MQ=60)
                                                                                                                                             ATAGAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATAAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAG          >  SRR19973072.3970228/1‑150 (MQ=60)
                                                                                                                                                GAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCA       <  SRR19973072.1260978/150‑1 (MQ=60)
                                                                                                                                                GAAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCA       >  SRR19973072.3533240/1‑150 (MQ=60)
                                                                                                                                                 AAACAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAG      >  SRR19973072.933324/1‑150 (MQ=60)
                                                                                                                                                    CAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAGTCA   <  SRR19973072.1457480/150‑1 (MQ=60)
                                                                                                                                                    CAGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAGTCA   <  SRR19973072.3227568/150‑1 (MQ=60)
                                                                                                                                                     AGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAGTCAA  >  SRR19973072.3255397/1‑150 (MQ=60)
                                                                                                                                                     AGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAGTCAA  >  SRR19973072.3437081/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
AGCGAAATAATCCGCCACGGTGCCAAACATATTGGCGCGGTGAGAAACGGAGACGCCGCGAATGCCTTCCGGTGGGTAACGGGTTGATGCCACCGCCAGCTCTGCTTCCTCTTTTGTTTCTACAAAAGGAATCAGGAAGTTATAGAAACCGATATCCAGAAGACGCTTAATAATTACCGGCTCGTTGGTCGGCACTCGCACTACTGGCGCGCTGGCGCTGCCTTTCAAGGCCATTAACTGCGGAATAAACGTGGAGATATCGTTTGGCGCATGTTCGCCATCCAGCACCAGCCAGTCAA  >  CP009273/3266497‑3266795

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: