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 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,618,693 G→T 89.1% E74* (GAA→TAA)  ppx → exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,618,6930GT89.1% 389.3 / 45.3 156E74* (GAA→TAA) ppxexopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base G (7/10);  new base T (72/67);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 156 total aligned reads displayed.

GGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAG  >  CP009273/2618548‑2618837
                                                                                                                                                 |                                                                                                                                                
ggTCATAGCCCGTGGGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTaa                                                                                                                                                >  2:2603386/1‑148 (MQ=255)
  tCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACg                                                                                                                                              >  1:600453/1‑148 (MQ=255)
  tCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACg                                                                                                                                              >  2:311323/1‑148 (MQ=255)
   cATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACgg                                                                                                                                             >  2:102799/1‑148 (MQ=255)
   cATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACgg                                                                                                                                             <  1:1887841/148‑1 (MQ=255)
   cATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACgg                                                                                                                                             <  1:1520004/148‑1 (MQ=255)
   cATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACgg                                                                                                                                             <  1:1911772/148‑1 (MQ=255)
    aTAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGc                                                                                                                                            >  1:79509/1‑148 (MQ=255)
    aTAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGc                                                                                                                                            >  1:1754391/1‑148 (MQ=255)
       gCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTAc                                                                                                                                         >  1:6325/1‑148 (MQ=255)
           gtgtGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAgg                                                                                                                                     <  1:2165332/148‑1 (MQ=255)
             gtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGt                                                                                                                                   <  2:3154881/148‑1 (MQ=255)
             gtgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGt                                                                                                                                   <  2:1789144/148‑1 (MQ=255)
              tgGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGtt                                                                                                                                  >  1:96059/1‑148 (MQ=255)
                   gATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTtctc                                                                                                                             <  2:729601/148‑1 (MQ=255)
                     tGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCt                                                                                                                           >  1:1728664/1‑148 (MQ=255)
                       gTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGc                                                                                                                         >  1:1498474/1‑148 (MQ=255)
                        tGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGcc                                                                                                                        <  2:2130704/148‑1 (MQ=255)
                         gCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCa                                                                                                                       >  2:363377/1‑148 (MQ=255)
                          ccATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAg                                                                                                                      >  1:1167887/1‑148 (MQ=255)
                              gCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTc                                                                                                                  >  2:559861/1‑148 (MQ=255)
                                 gATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGt                                                                                                               >  2:2708032/1‑148 (MQ=255)
                                  attattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGta                                                                                                              >  1:2536119/1‑148 (MQ=255)
                                   ttattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtat                                                                                                             >  2:1565336/1‑148 (MQ=255)
                                    tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                                                                                            >  1:2185467/1‑148 (MQ=255)
                                    tattGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGtata                                                                                                            >  1:3105152/1‑148 (MQ=255)
                                     attGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAg                                                                                                           >  2:1330248/1‑148 (MQ=255)
                                      ttGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGGGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGGTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGt                                                                                                          >  1:1282540/1‑148 (MQ=255)
                                         gccgccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAAGATGTTGAGTGAAGAGTCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCAAGCGTCTGTATAGTTgg                                                                                                       >  1:829423/1‑148 (MQ=255)
                                             ccTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTAcc                                                                                                   >  1:1288052/1‑148 (MQ=255)
                                                  aaCAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATAc                                                                                              >  1:3075809/1‑148 (MQ=255)
                                                        gggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCg                                                                                        <  1:2227037/148‑1 (MQ=255)
                                                         ggTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGt                                                                                       <  1:2177522/148‑1 (MQ=255)
                                                           tGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCa                                                                                     >  2:1895722/1‑148 (MQ=255)
                                                                 tGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCt                                                                               <  1:1595771/148‑1 (MQ=255)
                                                                   gCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGa                                                                             <  1:2142981/148‑1 (MQ=255)
                                                                      gACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACg                                                                          >  1:3212241/1‑148 (MQ=255)
                                                                       aCGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTTCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGc                                                                         >  2:475677/1‑148 (MQ=255)
                                                                         ggCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCa                                                                       <  2:956646/148‑1 (MQ=255)
                                                                          gCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCAc                                                                      >  1:2301422/1‑148 (MQ=255)
                                                                           ccTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTTCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACt                                                                     >  1:2036497/1‑148 (MQ=255)
                                                                           ccTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACt                                                                     <  1:380657/148‑1 (MQ=255)
                                                                            cTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTg                                                                    <  2:2623729/148‑1 (MQ=255)
                                                                             tGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGa                                                                   >  2:1515635/1‑148 (MQ=255)
                                                                               gggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACt                                                                 >  2:1801836/1‑148 (MQ=255)
                                                                               gggCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACt                                                                 >  2:563567/1‑148 (MQ=255)
                                                                                 gCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACttt                                                               <  2:1381815/148‑1 (MQ=255)
                                                                                   cAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCt                                                             <  1:2473546/148‑1 (MQ=255)
                                                                                      ataataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGaa                                                          >  2:2927280/1‑148 (MQ=255)
                                                                                       taataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGaaa                                                         >  1:2747420/1‑148 (MQ=255)
                                                                                       taataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGaaa                                                         >  2:3189788/1‑148 (MQ=255)
                                                                                       taataTGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGaaa                                                         >  1:2295340/1‑148 (MQ=255)
                                                                                         atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                                                       <  1:1329026/148‑1 (MQ=255)
                                                                                         atatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcg                                                       >  1:3156878/1‑148 (MQ=255)
                                                                                          tatGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAAcgc                                                      >  2:2051580/1‑148 (MQ=255)
                                                                                                 aGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGaaa                                               <  2:1884405/148‑1 (MQ=255)
                                                                                                   tGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAg                                             <  2:1589230/148‑1 (MQ=255)
                                                                                                      agagGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTc                                          <  2:759578/148‑1 (MQ=255)
                                                                                                       gagGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGACAAGGTCa                                         >  2:24804/1‑148 (MQ=255)
                                                                                                        agGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCAt                                        <  1:528251/148‑1 (MQ=255)
                                                                                                        agGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCAt                                        <  2:888502/148‑1 (MQ=255)
                                                                                                          gCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCCTACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTc                                      >  2:2335776/1‑148 (MQ=255)
                                                                                                          gCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTc                                      >  2:2825737/1‑148 (MQ=255)
                                                                                                          gCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTc                                      <  2:683484/148‑1 (MQ=255)
                                                                                                           cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                     >  2:1488425/1‑148 (MQ=255)
                                                                                                           cAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTcc                                     <  2:2129851/148‑1 (MQ=255)
                                                                                                            aaTGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTccc                                    >  1:3094545/1‑148 (MQ=255)
                                                                                                              tGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTa                                  >  1:2049400/1‑148 (MQ=255)
                                                                                                               gACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAc                                 <  1:1079845/148‑1 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGACAAGGTCATTCCCTAcc                                >  2:2318419/1‑148 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                                >  1:2691705/1‑148 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                                <  2:857337/148‑1 (MQ=255)
                                                                                                                aCGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTAcc                                >  2:1295652/1‑148 (MQ=255)
                                                                                                                   cgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGa                             <  2:6325/148‑1 (MQ=255)
                                                                                                                   cgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGa                             <  2:2074560/148‑1 (MQ=255)
                                                                                                                   cgcgGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGa                             <  2:2074476/148‑1 (MQ=255)
                                                                                                                        tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaa                        <  1:102752/148‑1 (MQ=255)
                                                                                                                        tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaa                        >  1:2666719/1‑148 (MQ=255)
                                                                                                                        tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaa                        <  2:1498474/148‑1 (MQ=255)
                                                                                                                        tttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTCCGTCAGGCGCAGAACGCCACTGACTTTCGGAAACGCGCGGAAAAGGTCATTCCCTACACGATTGaa                        >  1:1014629/1‑148 (MQ=255)
                                                                                                                         ttAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaaa                       >  1:1853709/1‑148 (MQ=255)
                                                                                                                         ttAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaaa                       <  2:351233/148‑1 (MQ=255)
                                                                                                                         ttAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGaaa                       <  2:441721/148‑1 (MQ=255)
                                                                                                                          tAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAat                      <  1:1072913/148‑1 (MQ=255)
                                                                                                                             aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattat                   <  1:3278661/148‑1 (MQ=255)
                                                                                                                             aCTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattat                   >  1:1635602/1‑148 (MQ=255)
                                                                                                                              ctgtctgtCGCTGTTTTCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCAAACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                  >  1:192253/1‑148 (MQ=255)
                                                                                                                              ctgtctgtCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAattatt                  <  2:1919358/148‑1 (MQ=255)
                                                                                                                               tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt                 >  1:1186320/1‑148 (MQ=255)
                                                                                                                               tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt                 <  2:813386/148‑1 (MQ=255)
                                                                                                                               tgtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAttt                 >  2:1070504/1‑148 (MQ=255)
                                                                                                                                gtctgtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTc                >  1:3148384/1‑148 (MQ=255)
                                                                                                                                 tctgtcGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTcc               <  2:2908825/148‑1 (MQ=255)
                                                                                                                                    gtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGt            <  2:1122952/148‑1 (MQ=255)
                                                                                                                                    gtcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGt            <  1:192401/148‑1 (MQ=255)
                                                                                                                                     tcGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTa           <  1:2946625/148‑1 (MQ=255)
                                                                                                                                       gCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAAt         <  2:302871/148‑1 (MQ=255)
                                                                                                                                        cTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATg        >  1:3132801/1‑148 (MQ=255)
                                                                                                                                          gTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaa      <  2:1710018/148‑1 (MQ=255)
                                                                                                                                             tGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATgaagaa   <  2:312174/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAG  >  CP009273/2618548‑2618837

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
CATGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAGCC  >  CP009273/2618545‑2618839
                                                                                                                                                    |                                                                                                                                                  
CATGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTA                                                                                                                                                   >  SRR19973070.305004/1‑150 (MQ=60)
  TGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAAC                                                                                                                                                 <  SRR19973070.3045432/150‑1 (MQ=60)
   GGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACG                                                                                                                                                >  SRR19973070.604636/1‑150 (MQ=60)
     TCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGC                                                                                                                                              >  SRR19973070.80342/1‑150 (MQ=60)
     TCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGC                                                                                                                                              >  SRR19973070.1766549/1‑150 (MQ=60)
      CATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCT                                                                                                                                             >  SRR19973070.2248744/1‑150 (MQ=60)
      CATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCT                                                                                                                                             <  SRR19973070.1900652/150‑1 (MQ=60)
      CATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCT                                                                                                                                             <  SRR19973070.2726323/150‑1 (MQ=60)
      CATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCT                                                                                                                                             <  SRR19973070.1846101/150‑1 (MQ=60)
      CATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCT                                                                                                                                             <  SRR19973070.1530670/150‑1 (MQ=60)
      CATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCT                                                                                                                                             <  SRR19973070.1924689/150‑1 (MQ=60)
        TAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTAC                                                                                                                                           >  SRR19973070.6419/1‑150 (MQ=60)
              GTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGT                                                                                                                                     <  SRR19973070.2179810/150‑1 (MQ=60)
              GTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGT                                                                                                                                     <  SRR19973070.2179701/150‑1 (MQ=60)
               TGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTT                                                                                                                                    >  SRR19973070.96999/1‑150 (MQ=60)
                    TAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTC                                                                                                                               <  SRR19973070.1277700/150‑1 (MQ=60)
                      GATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCT                                                                                                                             >  SRR19973070.1740644/1‑150 (MQ=60)
                      GATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACGGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCT                                                                                                                             <  SRR19973070.2196094/150‑1 (MQ=60)
                       ATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTG                                                                                                                            <  SRR19973070.1941391/150‑1 (MQ=60)
                        TGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGC                                                                                                                           >  SRR19973070.1508967/1‑150 (MQ=60)
                           TGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAG                                                                                                                        >  SRR19973070.1176063/1‑150 (MQ=60)
                           TGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAG                                                                                                                        >  SRR19973070.1175723/1‑150 (MQ=60)
                               ATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTC                                                                                                                    >  SRR19973070.2665183/1‑150 (MQ=60)
                                 GCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTG                                                                                                                  <  SRR19973070.3259209/150‑1 (MQ=60)
                                   AGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTA                                                                                                                >  SRR19973070.2552696/1‑150 (MQ=60)
                                   AGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTA                                                                                                                >  SRR19973070.2149684/1‑150 (MQ=60)
                                     ATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATA                                                                                                              >  SRR19973070.3125179/1‑150 (MQ=60)
                                     ATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATA                                                                                                              >  SRR19973070.3125718/1‑150 (MQ=60)
                                     ATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATA                                                                                                              >  SRR19973070.2199916/1‑150 (MQ=60)
                                       gcTTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGGGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGGTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGT                                                                                                            >  SRR19973070.1291496/3‑150 (MQ=60)
                                          TGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAAGATGTTGAGTGAAGAGTCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCAAGCGTCTGTATAGTTGG                                                                                                         >  SRR19973070.835053/1‑150 (MQ=60)
                                              CGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACC                                                                                                     >  SRR19973070.1297035/1‑150 (MQ=60)
                                                   GAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATAC                                                                                                >  SRR19973070.3095563/1‑150 (MQ=60)
                                                    AAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACG                                                                                               <  SRR19973070.1373957/150‑1 (MQ=60)
                                                           GGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTC                                                                                        <  SRR19973070.2241679/150‑1 (MQ=60)
                                                            GGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCA                                                                                       <  SRR19973070.2191941/150‑1 (MQ=60)
                                                                 ATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGC                                                                                  >  SRR19973070.863268/1‑150 (MQ=60)
                                                                    TGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGA                                                                               <  SRR19973070.1606923/150‑1 (MQ=60)
                                                                      GCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAAC                                                                             <  SRR19973070.2157271/150‑1 (MQ=60)
                                                                       CGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACG                                                                            >  SRR19973070.3233012/1‑150 (MQ=60)
                                                                           CGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCAC                                                                        >  SRR19973070.2316729/1‑150 (MQ=60)
                                                                            GGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTTCCTAACGGCTACAAGGGTTTTCTACTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACA                                                                       >  SRR19973070.2050275/1‑150 (MQ=60)
                                                                              CCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGA                                                                     <  SRR19973070.1924744/150‑1 (MQ=60)
                                                                              CCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGA                                                                     <  SRR19973070.383448/150‑1 (MQ=60)
                                                                                 GGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTT                                                                  <  SRR19973070.640741/150‑1 (MQ=60)
                                                                                     CCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTG                                                              <  SRR19973070.1576348/150‑1 (MQ=60)
                                                                                      CAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGA                                                             <  SRR19973070.2489662/150‑1 (MQ=60)
                                                                                        GATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAA                                                           >  SRR19973070.2765078/1‑150 (MQ=60)
                                                                                        GATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAA                                                           >  SRR19973070.2310606/1‑150 (MQ=60)
                                                                                          TAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACG                                                         >  SRR19973070.1721827/1‑150 (MQ=60)
                                                                                          TAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACG                                                         >  SRR19973070.3176393/1‑150 (MQ=60)
                                                                                          TAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACG                                                         >  SRR19973070.3177287/1‑150 (MQ=60)
                                                                                            ATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCG                                                       <  SRR19973070.1338194/150‑1 (MQ=60)
                                                                                              ATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCG                                                     <  SRR19973070.54846/150‑1 (MQ=60)
                                                                                                       GAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTC                                            <  SRR19973070.1908570/150‑1 (MQ=60)
                                                                                                           AGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTC                                        <  SRR19973070.531820/150‑1 (MQ=60)
                                                                                                           AGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTC                                        <  SRR19973070.3101466/150‑1 (MQ=60)
                                                                                                             GCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCC                                      >  SRR19973070.3114481/1‑150 (MQ=60)
                                                                                                               AATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTA                                    <  SRR19973070.349885/150‑1 (MQ=60)
                                                                                                               AATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTA                                    >  SRR19973070.2063260/1‑150 (MQ=60)
                                                                                                                 TGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACC                                  >  SRR19973070.2709141/1‑150 (MQ=60)
                                                                                                                 TGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACC                                  >  SRR19973070.2434245/1‑150 (MQ=60)
                                                                                                                  GACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCC                                 <  SRR19973070.1087352/150‑1 (MQ=60)
                                                                                                                         GGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAA                          >  SRR19973070.2684040/1‑150 (MQ=60)
                                                                                                                         GGTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTCCGTCAGGCGCAGAACGCCACTGACTTTCGGAAACGCGCGGAAAAGGTCATTCCCTACACGATTGAA                          >  SRR19973070.1021835/1‑150 (MQ=60)
                                                                                                                          GTTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAA                         >  SRR19973070.1866407/1‑150 (MQ=60)
                                                                                                                           TTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAAT                        <  SRR19973070.103774/150‑1 (MQ=60)
                                                                                                                           TTTAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAAT                        <  SRR19973070.103727/150‑1 (MQ=60)
                                                                                                                             TAAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTA                      <  SRR19973070.1080401/150‑1 (MQ=60)
                                                                                                                              AAACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTAT                     >  SRR19973070.1646942/1‑150 (MQ=60)
                                                                                                                               AACTGTCTGTCGCTGTTTTCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCAAACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATT                    >  SRR19973070.193734/1‑150 (MQ=60)
                                                                                                                                ACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTT                   <  SRR19973070.3299728/150‑1 (MQ=60)
                                                                                                                                ACTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTT                   >  SRR19973070.1194659/1‑150 (MQ=60)
                                                                                                                                 CTGTCTGTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTC                  >  SRR19973070.3168730/1‑150 (MQ=60)
                                                                                                                                  TGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCC                 >  SRR19973070.1302410/1‑150 (MQ=60)
                                                                                                                                       GTCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAA            <  SRR19973070.193886/150‑1 (MQ=60)
                                                                                                                                        TCGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAAT           <  SRR19973070.2965597/150‑1 (MQ=60)
                                                                                                                                         CGCTGTTTGCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATG          >  SRR19973070.3153036/1‑150 (MQ=60)
                                                                                                                                                 GCCTAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAGCC  <  SRR19973070.2946174/150‑1 (MQ=60)
                                                                                                                                                    |                                                                                                                                                  
CATGGTCATAGCCCGTGTGGTAGATGGTGCCATGCAGATTATTGGCCGCCTGAAACAGCGGGTGCATCTGGCGGACGGCCTGGGGCCAGATAATATGTTGAGTGAAGAGGCAATGACGCGCGGTTTAAACTGTCTGTCGCTGTTTGCCGAACGGCTACAAGGGTTTTCTCCTGCCAGCGTCTGTATAGTTGGTACCCATACGCTGCGTCAGGCGCTGAACGCCACTGACTTTCTGAAACGCGCGGAAAAGGTCATTCCCTACCCGATTGAAATTATTTCCGGTAATGAAGAAGCC  >  CP009273/2618545‑2618839

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: