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 2,639,116 C→A 100% A523A (GCG→GCT)  pbpC ← penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,639,1160CA97.2% 503.0 / 6.6 181A523A (GCG→GCT) pbpCpenicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (94/82);  total (98/83)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-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 181 total aligned reads displayed.

GCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC  >  CP009273/2638971‑2639260
                                                                                                                                                 |                                                                                                                                                
gCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGc                                                                                                                                                <  2:1508672/148‑1 (MQ=255)
gCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGc                                                                                                                                                >  2:2321595/1‑148 (MQ=255)
 cAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGcc                                                                                                                                               <  2:86820/148‑1 (MQ=255)
 cAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGcc                                                                                                                                               <  1:1652106/148‑1 (MQ=255)
    aTTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAg                                                                                                                                            >  2:3398058/1‑148 (MQ=255)
     ttCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAgg                                                                                                                                           >  1:3350596/1‑148 (MQ=255)
      tCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGc                                                                                                                                          >  1:2673728/1‑148 (MQ=255)
         cACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTc                                                                                                                                       <  1:2117392/148‑1 (MQ=255)
           cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                     >  2:3352647/1‑148 (MQ=255)
           cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                     >  1:2651762/1‑148 (MQ=255)
           cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                     >  1:3548241/1‑148 (MQ=255)
           cTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAc                                                                                                                                     >  1:3239176/1‑148 (MQ=255)
            tGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACg                                                                                                                                    <  2:3041276/148‑1 (MQ=255)
                 gACAGTAAGATATTATTGACCTGTTTGAACAATGGTTCGGCACTGGCAATGCCAAACTGAACAACAACTGGCGTGCCGTGCGGTCTGCCTGTCCAGTTGTCAATGTCATAGCGTGCTTTAACCCCAATATGCCAGTCGTATCGATAGc                                                                                                                               <  2:3243598/148‑1 (MQ=255)
                  aCAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGcc                                                                                                                              >  2:1250229/1‑148 (MQ=255)
                       aaGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGc                                                                                                                         >  1:695109/1‑148 (MQ=255)
                        aGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCt                                                                                                                        >  1:987997/1‑148 (MQ=255)
                          atatTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTgg                                                                                                                      >  1:1974936/1‑148 (MQ=255)
                           tattatTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGt                                                                                                                     >  2:2508742/1‑148 (MQ=255)
                             ttattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGc                                                                                                                   >  1:2710673/1‑148 (MQ=255)
                              tattGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGATAACCCCAATAGCCCAGGCGACACGATAGCCATAGCTGGTGcc                                                                                                                  >  1:2348664/1‑148 (MQ=255)
                                ttGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCg                                                                                                                <  2:1166689/148‑1 (MQ=255)
                                ttGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCg                                                                                                                <  2:2817454/148‑1 (MQ=255)
                                ttGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCg                                                                                                                <  2:3392894/148‑1 (MQ=255)
                                ttGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCg                                                                                                                >  1:3184156/1‑148 (MQ=255)
                                ttGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCg                                                                                                                >  2:1148567/1‑148 (MQ=255)
                                  gACCTGATTCAACACTGGTACGGCACTGGCAAAGCCAAACTGACCATCAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCttt                                                                                                              <  2:2396339/148‑1 (MQ=255)
                                   aCCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGttt                                                                                                             >  1:2964391/1‑148 (MQ=255)
                                   aCCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGttt                                                                                                             <  2:1476453/148‑1 (MQ=255)
                                    ccTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGtttt                                                                                                            >  2:2604326/1‑148 (MQ=255)
                                     cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTc                                                                                                           >  2:2790316/1‑148 (MQ=255)
                                     cTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTc                                                                                                           >  1:851494/1‑148 (MQ=255)
                                      tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                          >  1:3496127/1‑148 (MQ=255)
                                      tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                          >  1:1390025/1‑148 (MQ=255)
                                      tGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTcc                                                                                                          <  2:2859531/148‑1 (MQ=255)
                                        aTTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAt                                                                                                        >  2:2153788/1‑148 (MQ=255)
                                        aTTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCAt                                                                                                        >  2:2153526/1‑148 (MQ=255)
                                          tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCTTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      <  2:703383/148‑1 (MQ=255)
                                          tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  2:2657047/1‑148 (MQ=255)
                                          tCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGc                                                                                                      >  2:829514/1‑148 (MQ=255)
                                           caacaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGcc                                                                                                     >  1:2023652/1‑148 (MQ=255)
                                           caacaaTGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGcc                                                                                                     >  1:2996087/1‑148 (MQ=255)
                                                 tGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTggg                                                                                               >  1:48613/1‑148 (MQ=255)
                                                 tGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTggg                                                                                               >  1:1725625/1‑148 (MQ=255)
                                                        gCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgc                                                                                        >  1:2469855/1‑148 (MQ=255)
                                                         cACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGAcgcg                                                                                       >  2:2534068/1‑148 (MQ=255)
                                                              gCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCCGTGGGGCGACGCGCGGCa                                                                                  >  2:2237620/1‑148 (MQ=255)
                                                              gCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCa                                                                                  >  2:2785658/1‑148 (MQ=255)
                                                               cAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCaa                                                                                 >  1:1824773/1‑148 (MQ=255)
                                                                   gCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCa                                                                             <  2:279392/148‑1 (MQ=255)
                                                                   gCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCa                                                                             <  2:3505455/148‑1 (MQ=255)
                                                                    ccAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                                                                            >  1:1245288/1‑148 (MQ=255)
                                                                    ccAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAc                                                                            >  1:2544127/1‑148 (MQ=255)
                                                                     cAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACt                                                                           >  1:147395/1‑148 (MQ=255)
                                                                      aaccTGACCAACAAAGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCGCGTTTTCCATGCCTGTGGGGCGACGCGCGGCAAGGCACTa                                                                          <  2:1607336/145‑1 (MQ=255)
                                                                      aaaCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTa                                                                          <  2:1602306/148‑1 (MQ=255)
                                                                        aCTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCACGGCAAGGCACTATc                                                                        >  2:2674867/1‑148 (MQ=255)
                                                                         cTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATcc                                                                       >  2:1130352/1‑148 (MQ=255)
                                                                         cTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATcc                                                                       <  2:703191/148‑1 (MQ=255)
                                                                          tGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCg                                                                      >  1:850473/1‑148 (MQ=255)
                                                                              caacaaCGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCaa                                                                  <  1:2598604/148‑1 (MQ=255)
                                                                              caacaaCGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCaa                                                                  >  2:3471113/1‑148 (MQ=255)
                                                                                      ggCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTgcg                                                          >  1:6291/1‑148 (MQ=255)
                                                                                         gTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCtt                                                       <  1:2282022/148‑1 (MQ=255)
                                                                                          tGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtc                                                      <  2:905526/148‑1 (MQ=255)
                                                                                           gCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTtca                                                     <  1:3508513/148‑1 (MQ=255)
                                                                                                ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                >  2:2453141/1‑148 (MQ=255)
                                                                                                ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                <  1:590934/148‑1 (MQ=255)
                                                                                                ccGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGc                                                >  1:3040520/1‑148 (MQ=255)
                                                                                                 cGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGcc                                               <  2:3249852/148‑1 (MQ=255)
                                                                                                 cGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGcc                                               >  2:1018869/1‑148 (MQ=255)
                                                                                                    tCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATa                                            <  2:2964391/148‑1 (MQ=255)
                                                                                                       gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         >  1:3491659/1‑148 (MQ=255)
                                                                                                       gCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATc                                         <  1:2185528/148‑1 (MQ=255)
                                                                                                         cAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCg                                       >  2:3538311/1‑148 (MQ=255)
                                                                                                          aGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCgg                                      <  1:1116549/148‑1 (MQ=255)
                                                                                                          aGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCgg                                      <  1:3527901/148‑1 (MQ=255)
                                                                                                             ccAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   <  2:1737136/148‑1 (MQ=255)
                                                                                                             ccAGATCCCAATGACATAGCGAGCGTTAACCCCAAGAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   <  2:2983970/148‑1 (MQ=255)
                                                                                                             ccAGATCCCAATGACATAGCGAGAGTTAACCCCAATAGCCCAGTCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGa                                   >  2:693687/1‑148 (MQ=255)
                                                                                                              cAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGAGAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGaa                                  <  2:1716095/148‑1 (MQ=255)
                                                                                                                gATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                <  2:533707/148‑1 (MQ=255)
                                                                                                                gATCCCAATGACATAGCGAGCGGTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATg                                <  2:181045/148‑1 (MQ=255)
                                                                                                                 aTCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGACTGa                               >  2:2165628/1‑148 (MQ=255)
                                                                                                                   cccAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATc                             >  1:1101736/1‑148 (MQ=255)
                                                                                                                    ccAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATcc                            <  1:2913841/148‑1 (MQ=255)
                                                                                                                         gACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgcc                       >  2:1937876/1‑148 (MQ=255)
                                                                                                                          aCATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAcgccg                      >  2:36869/1‑148 (MQ=255)
                                                                                                                             tAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGccc                   <  1:2790316/148‑1 (MQ=255)
                                                                                                                              agcgagcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTTTCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGcccc                  >  2:3510445/1‑148 (MQ=255)
                                                                                                                                  agcgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgac              >  1:738328/1‑148 (MQ=255)
                                                                                                                                    cgTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCcgacga            >  1:1934208/1‑148 (MQ=255)
                                                                                                                                      ttAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACa          <  1:1454093/148‑1 (MQ=255)
                                                                                                                                      ttAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACa          <  2:2079277/148‑1 (MQ=255)
                                                                                                                                       tAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAt         >  1:1462546/1‑148 (MQ=255)
                                                                                                                                         aCCCCACTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTa       <  1:515481/148‑1 (MQ=255)
                                                                                                                                          ccccAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTaa      >  2:1860301/1‑148 (MQ=255)
                                                                                                                                             cAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAgg   <  2:664658/148‑1 (MQ=255)
                                                                                                                                              aaTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGc  >  2:1413984/1‑148 (MQ=255)
                                                                                                                                              aaTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGc  <  2:2995903/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGC  >  CP009273/2638971‑2639260

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
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGT  >  CP009273/2638972‑2639263
                                                                                                                                                |                                                                                                                                                   
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCA                                                                                                                                                <  SRR19973067.1663483/150‑1 (MQ=60)
  GATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGG                                                                                                                                              >  SRR19973067.3371841/1‑150 (MQ=60)
   ATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGC                                                                                                                                             >  SRR19973067.2690882/1‑150 (MQ=60)
       GCAATGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCA                                                                                                                                         >  SRR19973067.417456/1‑150 (MQ=60)
        CACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAC                                                                                                                                        >  SRR19973067.3570621/1‑150 (MQ=60)
        CACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAC                                                                                                                                        <  SRR19973067.2131505/150‑1 (MQ=60)
        CACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAC                                                                                                                                        >  SRR19973067.1723908/1‑150 (MQ=60)
        CACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAC                                                                                                                                        >  SRR19973067.2668856/1‑150 (MQ=60)
        CACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAC                                                                                                                                        >  SRR19973067.2668907/1‑150 (MQ=60)
        CACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCAC                                                                                                                                        >  SRR19973067.3259488/1‑150 (MQ=60)
             CGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAG                                                                                                                                   <  SRR19973067.2611725/150‑1 (MQ=60)
             CGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAG                                                                                                                                   <  SRR19973067.199147/150‑1 (MQ=60)
                    GTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGC                                                                                                                            >  SRR19973067.699907/1‑150 (MQ=60)
                    GTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGC                                                                                                                            >  SRR19973067.699909/1‑150 (MQ=60)
                     TAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCT                                                                                                                           >  SRR19973067.2654064/1‑150 (MQ=60)
                     TAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCT                                                                                                                           <  SRR19973067.3215121/150‑1 (MQ=60)
                     TAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCT                                                                                                                           >  SRR19973067.994794/1‑150 (MQ=60)
                     TAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCT                                                                                                                           <  SRR19973067.1902844/150‑1 (MQ=60)
                     TAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCT                                                                                                                           >  SRR19973067.994793/1‑150 (MQ=60)
                       AGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGG                                                                                                                         >  SRR19973067.1988319/1‑150 (MQ=60)
                          TATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGC                                                                                                                      >  SRR19973067.2728134/1‑150 (MQ=60)
                           ATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGATAACCCCAATAGCCCAGGCGACACGATAGCCATAGCTGGTGCC                                                                                                                     >  SRR19973067.2364001/1‑150 (MQ=60)
                            TTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCC                                                                                                                    <  SRR19973067.2973880/150‑1 (MQ=60)
                             TATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCG                                                                                                                   >  SRR19973067.3204295/1‑150 (MQ=60)
                                TGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTT                                                                                                                >  SRR19973067.2983113/1‑150 (MQ=60)
                                  ACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTC                                                                                                              >  SRR19973067.857246/1‑150 (MQ=60)
                                   CCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCC                                                                                                             >  SRR19973067.3518307/1‑150 (MQ=60)
                                   CCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCC                                                                                                             >  SRR19973067.3518313/1‑150 (MQ=60)
                                   CCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCC                                                                                                             >  SRR19973067.1399472/1‑150 (MQ=60)
                                        TTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCC                                                                                                        >  SRR19973067.2037233/1‑150 (MQ=60)
                                        TTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCC                                                                                                        >  SRR19973067.3014894/1‑150 (MQ=60)
                                        TTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCC                                                                                                        >  SRR19973067.3015079/1‑150 (MQ=60)
                                           AACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGT                                                                                                     >  SRR19973067.3283028/1‑150 (MQ=60)
                                             CAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGG                                                                                                   >  SRR19973067.611523/1‑150 (MQ=60)
                                              AATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGG                                                                                                  >  SRR19973067.49258/1‑150 (MQ=60)
                                              AATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGG                                                                                                  >  SRR19973067.1737452/1‑150 (MQ=60)
                                                     CGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGC                                                                                           >  SRR19973067.2486092/1‑150 (MQ=60)
                                                      GGCACTGGCAAAGACAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCG                                                                                          >  SRR19973067.2093063/1‑150 (MQ=60)
                                                           TGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCA                                                                                     >  SRR19973067.3385563/1‑150 (MQ=60)
                                                           TGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCA                                                                                     >  SRR19973067.669196/1‑150 (MQ=60)
                                                            GGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAA                                                                                    >  SRR19973067.1836647/1‑150 (MQ=60)
                                                            GGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAA                                                                                    >  SRR19973067.1837159/1‑150 (MQ=60)
                                                              CAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCACTCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGTGGCGACGCGCGGCAAGG                                                                                  <  SRR19973067.2525200/150‑1 (MQ=60)
                                                                 AGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAC                                                                               >  SRR19973067.1253684/1‑150 (MQ=60)
                                                                 AGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAC                                                                               >  SRR19973067.2560777/1‑150 (MQ=60)
                                                                 AGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCAC                                                                               >  SRR19973067.767251/1‑150 (MQ=60)
                                                                  GCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACT                                                                              >  SRR19973067.148670/1‑150 (MQ=60)
                                                                     AAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATC                                                                           <  SRR19973067.1258672/150‑1 (MQ=60)
                                                                       ACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCG                                                                         >  SRR19973067.972181/1‑150 (MQ=60)
                                                                       ACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCG                                                                         >  SRR19973067.856223/1‑150 (MQ=60)
                                                                             CAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGG                                                                   <  SRR19973067.2615503/150‑1 (MQ=60)
                                                                             CAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGG                                                                   >  SRR19973067.2143756/1‑150 (MQ=60)
                                                                                 AACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTG                                                               <  SRR19973067.1188807/150‑1 (MQ=60)
                                                                                   CGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCG                                                             >  SRR19973067.6384/1‑150 (MQ=60)
                                                                                        GTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCA                                                        >  SRR19973067.2741033/1‑150 (MQ=60)
                                                                                        GTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCA                                                        <  SRR19973067.2297112/150‑1 (MQ=60)
                                                                                          GCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATC                                                      <  SRR19973067.1501139/150‑1 (MQ=60)
                                                                                          GCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATC                                                      <  SRR19973067.3530733/150‑1 (MQ=60)
                                                                                             GTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCACGGGTTGCGCTTCATCAGC                                                   >  SRR19973067.3059887/1‑150 (MQ=60)
                                                                                               CCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCA                                                 <  SRR19973067.594696/150‑1 (MQ=60)
                                                                                                    CTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATC                                            >  SRR19973067.3513819/1‑150 (MQ=60)
                                                                                                      GCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCG                                          <  SRR19973067.2200129/150‑1 (MQ=60)
                                                                                                         AGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCG                                       <  SRR19973067.1124119/150‑1 (MQ=60)
                                                                                                         AGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCG                                       <  SRR19973067.3550208/150‑1 (MQ=60)
                                                                                                         AGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCG                                       <  SRR19973067.3550109/150‑1 (MQ=60)
                                                                                                             CAGATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATG                                   <  SRR19973067.3145187/150‑1 (MQ=60)
                                                                                                                ATCCCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATC                                >  SRR19973067.1109230/1‑150 (MQ=60)
                                                                                                                  CCCAATGACATCGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCA                              <  SRR19973067.2427069/150‑1 (MQ=60)
                                                                                                                   CCAATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCAC                             <  SRR19973067.2932392/150‑1 (MQ=60)
                                                                                                                      ATGACATAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCC                          >  SRR19973067.2219260/1‑150 (MQ=60)
                                                                                                                            TAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCC                    <  SRR19973067.2808321/150‑1 (MQ=60)
                                                                                                                            TAGCGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCC                    <  SRR19973067.2808549/150‑1 (MQ=60)
                                                                                                                               CGAGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGAC                 >  SRR19973067.743451/1‑150 (MQ=60)
                                                                                                                                 AGCGTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGA               >  SRR19973067.1947403/1‑150 (MQ=60)
                                                                                                                                    GTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAT            >  SRR19973067.3551307/1‑150 (MQ=60)
                                                                                                                                    GTTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACAT            >  SRR19973067.1472314/1‑150 (MQ=60)
                                                                                                                                     TTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATT           <  SRR19973067.1465037/150‑1 (MQ=60)
                                                                                                                                     TTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATT           <  SRR19973067.1463825/150‑1 (MQ=60)
                                                                                                                                     TTAACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATT           <  SRR19973067.1463775/150‑1 (MQ=60)
                                                                                                                                        ACCCCACTAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAA        <  SRR19973067.518979/150‑1 (MQ=60)
                                                                                                                                        ACCCCAATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAA        <  SRR19973067.518980/150‑1 (MQ=60)
                                                                                                                                              ATAGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGT  <  SRR19973067.361305/150‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                                   
CAGATTCGCACTGCGCGACAGTAAGATATTATTGACCTGATTCAACAATGGTACGGCACTGGCAAAGCCAAACTGACCAACAACGGGCGTGCCGTCCGGTCTGCCAGTCCAGATCCCAATGACATAGCGAGCGTTAACCCCAATCGCCCAGGCGTCACGATAGCCATAGCTGGTGCCCGTTTTCCATGCCAGTGGGGCGACGCGCGGCAAGGCACTATCCGGCAAGGGTTGCGCTTCATCAGCCATAATCCGGCGAATGATCCACGCCGCCCCCGACGACATTAAAGGCCGT  >  CP009273/2638972‑2639263

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: ‑