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 F342 I0 R1 1331 225.9 7641038 95.8% 7320114 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,039,703 C→A 100% A900S (GCG→TCG)  gcvP ← glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,039,7030CA98.1% 668.0 / 3.2 210A900S (GCG→TCG) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base A (86/120);  total (86/124)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-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 210 total aligned reads displayed.

GAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGT  >  CP009273/3039560‑3039841
                                                                                                                                               |                                                                                                                                          
  aCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCAc                                                                                                                                      <  1:3642178/148‑1 (MQ=255)
  aCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAc                                                                                                                                      >  1:1860761/1‑148 (MQ=255)
  aCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAc                                                                                                                                      >  1:3311217/1‑148 (MQ=255)
   cAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAcc                                                                                                                                     >  2:405034/1‑148 (MQ=255)
   cAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAcc                                                                                                                                     >  2:1864992/1‑148 (MQ=255)
       ttACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCCGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCg                                                                                                                                 >  2:2021909/1‑148 (MQ=255)
       ttACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGACTGTGCGGCGAGTTCACCAGCg                                                                                                                                 >  2:89822/1‑148 (MQ=255)
                ccGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCtt                                                                                                                        <  2:3445274/148‑1 (MQ=255)
                 cGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTc                                                                                                                       <  2:3397666/148‑1 (MQ=255)
                      aCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCg                                                                                                                  >  1:3394499/1‑148 (MQ=255)
                       catcatCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCgg                                                                                                                 >  1:1200937/1‑148 (MQ=255)
                        atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGc                                                                                                                <  1:904064/148‑1 (MQ=255)
                         tcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCTGTTCGCTCTGAATGTGCGGCGATTTCACCAGCGGTTTATCTTCCAGCGGcc                                                                                                               <  1:755012/148‑1 (MQ=255)
                         tcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGcc                                                                                                               >  1:2698613/1‑148 (MQ=255)
                            tcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGa                                                                                                            <  2:3599753/148‑1 (MQ=255)
                            tcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGa                                                                                                            <  2:41158/148‑1 (MQ=255)
                             ccAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGac                                                                                                           <  1:85867/148‑1 (MQ=255)
                             ccAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGac                                                                                                           <  1:2852820/148‑1 (MQ=255)
                                gACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAcc                                                                                                        >  1:2155829/1‑148 (MQ=255)
                                gACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAcc                                                                                                        >  1:2592430/1‑148 (MQ=255)
                                 aCGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCg                                                                                                       >  2:1857701/1‑148 (MQ=255)
                                  cGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCgg                                                                                                      <  1:437308/148‑1 (MQ=255)
                                  cGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCgg                                                                                                      <  1:2815187/148‑1 (MQ=255)
                                   gTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGc                                                                                                     <  2:3089587/148‑1 (MQ=255)
                                     ttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCtt                                                                                                   <  1:2721789/148‑1 (MQ=255)
                                      tCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCttt                                                                                                  <  1:285830/148‑1 (MQ=255)
                                      tCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCttt                                                                                                  <  1:284341/148‑1 (MQ=255)
                                       cACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTc                                                                                                 <  2:2899724/148‑1 (MQ=255)
                                         cTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAc                                                                                               <  1:2451491/148‑1 (MQ=255)
                                          tGTCGGCCAGTATGTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                              <  1:1301713/148‑1 (MQ=255)
                                             cGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTgg                                                                                           <  1:2357148/148‑1 (MQ=255)
                                                ccAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                        >  2:736140/1‑148 (MQ=255)
                                                ccAGTATTTGTCTGCCACACCTGCCGGGAAAACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCCGACACCGGCTTTCACCTGGTCa                                                                                        >  2:2005548/1‑148 (MQ=255)
                                                  aGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAt                                                                                      <  2:1914803/148‑1 (MQ=255)
                                                    tagttGTCTGCCACACCAGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                    <  2:705541/145‑1 (MQ=255)
                                                    tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                    >  1:1445699/1‑148 (MQ=255)
                                                    tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                    <  2:336923/148‑1 (MQ=255)
                                                    tATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAttt                                                                                    >  2:3353362/1‑148 (MQ=255)
                                                            tGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCgg                                                                            <  1:545633/148‑1 (MQ=255)
                                                            tGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCgg                                                                            <  1:545534/148‑1 (MQ=255)
                                                             gCCACACATGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGa                                                                           >  2:1389003/1‑148 (MQ=255)
                                                              ccACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAt                                                                          <  1:758702/148‑1 (MQ=255)
                                                               cacaCCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGTGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATa                                                                         <  1:224969/148‑1 (MQ=255)
                                                                 cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                       <  1:2516431/148‑1 (MQ=255)
                                                                 cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                       <  1:415023/148‑1 (MQ=255)
                                                                 cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                       <  1:580210/148‑1 (MQ=255)
                                                                 cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                       <  1:3363747/148‑1 (MQ=255)
                                                                 cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                       >  2:3736435/1‑148 (MQ=255)
                                                                 cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                       <  1:3137561/148‑1 (MQ=255)
                                                                  acCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGcc                                                                      <  2:366430/148‑1 (MQ=255)
                                                                     tGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGc                                                                   >  2:3592282/1‑148 (MQ=255)
                                                                      gCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCa                                                                  <  1:1790470/148‑1 (MQ=255)
                                                                      gCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCa                                                                  >  2:618608/1‑148 (MQ=255)
                                                                       ccGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCAt                                                                 <  2:2031573/148‑1 (MQ=255)
                                                                              tACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCg                                                          <  2:1052431/148‑1 (MQ=255)
                                                                              tACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCg                                                          <  1:2241600/148‑1 (MQ=255)
                                                                               aCCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGa                                                         >  2:3110573/1‑148 (MQ=255)
                                                                               aCCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGa                                                         <  2:3573105/148‑1 (MQ=255)
                                                                               aCCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGa                                                         >  2:3554462/1‑148 (MQ=255)
                                                                                   cAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATaaa                                                     >  2:1018369/1‑148 (MQ=255)
                                                                                      cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                  <  2:3285645/148‑1 (MQ=255)
                                                                                      cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                  >  2:164261/1‑148 (MQ=255)
                                                                                      cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                  <  1:2476519/148‑1 (MQ=255)
                                                                                      cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                  >  1:781813/1‑148 (MQ=255)
                                                                                         cACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATc                                               <  1:3812305/148‑1 (MQ=255)
                                                                                          aCGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                              >  2:914861/1‑148 (MQ=255)
                                                                                          aCGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                              <  1:511054/148‑1 (MQ=255)
                                                                                          aCGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                              <  2:1189184/148‑1 (MQ=255)
                                                                                           cGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCa                                             >  2:2981598/1‑148 (MQ=255)
                                                                                             gCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGt                                           <  1:1860102/148‑1 (MQ=255)
                                                                                                 tACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCa                                       <  1:3346035/148‑1 (MQ=255)
                                                                                                   cGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACt                                     <  1:2333737/148‑1 (MQ=255)
                                                                                                      aTGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                  <  2:761111/148‑1 (MQ=255)
                                                                                                      aTGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                  <  1:1146138/148‑1 (MQ=255)
                                                                                                      aTGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                  <  2:2675525/148‑1 (MQ=255)
                                                                                                        gcgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCt                                <  2:3299234/148‑1 (MQ=255)
                                                                                                         cgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATATCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCtt                               >  1:504904/1‑148 (MQ=255)
                                                                                                              aCTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGa                          <  1:1558409/148‑1 (MQ=255)
                                                                                                              aCTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGa                          >  1:2650300/1‑148 (MQ=255)
                                                                                                              aCTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGa                          <  1:3099150/148‑1 (MQ=255)
                                                                                                               cTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCACCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGATTTCAGAt                         >  2:3292477/1‑148 (MQ=255)
                                                                                                                tCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAtt                        <  1:3174195/148‑1 (MQ=255)
                                                                                                                tCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAtt                        >  2:3198713/1‑148 (MQ=255)
                                                                                                                 cGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTc                       <  2:1750857/148‑1 (MQ=255)
                                                                                                                  ggCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCa                      >  1:2768705/1‑148 (MQ=255)
                                                                                                                    cGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGt                    <  1:2945589/148‑1 (MQ=255)
                                                                                                                      aCCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCg                  >  1:2072370/1‑148 (MQ=255)
                                                                                                                      aCCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCg                  >  2:1152521/1‑148 (MQ=255)
                                                                                                                      aCCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCg                  >  1:721853/1‑148 (MQ=255)
                                                                                                                       ccAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCgg                 <  1:3173393/148‑1 (MQ=255)
                                                                                                                           ttCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCa             >  1:848908/1‑148 (MQ=255)
                                                                                                                             cGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAAc           >  2:3492055/1‑148 (MQ=255)
                                                                                                                                tctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCAt        <  2:2155829/148‑1 (MQ=255)
                                                                                                                                tctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCAt        >  1:1879089/1‑148 (MQ=255)
                                                                                                                                 ctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATc       <  1:2129460/148‑1 (MQ=255)
                                                                                                                                 ctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATc       >  1:3561877/1‑148 (MQ=255)
                                                                                                                                 ctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATc       >  1:1526485/1‑148 (MQ=255)
                                                                                                                                    aaTGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGc    >  1:3490036/1‑148 (MQ=255)
                                                                                                                                      tgtgCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAGAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGt  >  2:3482714/1‑148 (MQ=255)
                                                                                                                                      tgtgCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGt  >  2:16607/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                          
GAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGT  >  CP009273/3039560‑3039841

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
GAGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC  >  CP009273/3039555‑3039850
                                                                                                                                                    |                                                                                                                                                   
GAGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAG                                                                                                                                                    <  SRR19973069.1913514/150‑1 (MQ=60)
  GCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCCGAATGTGCGGCGAGTT                                                                                                                                                  <  SRR19973069.2861464/150‑1 (MQ=60)
   CAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTC                                                                                                                                                 >  SRR19973069.751413/1‑150 (MQ=60)
     GAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTTTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAC                                                                                                                                               >  SRR19973069.1873999/1‑150 (MQ=60)
     GAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAC                                                                                                                                               >  SRR19973069.3332918/1‑150 (MQ=60)
     GAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAC                                                                                                                                               >  SRR19973069.1873750/1‑150 (MQ=60)
       ACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCA                                                                                                                                             <  SRR19973069.3666040/150‑1 (MQ=60)
           GTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGG                                                                                                                                         >  SRR19973069.2163019/1‑150 (MQ=60)
           GTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGG                                                                                                                                         >  SRR19973069.1763180/1‑150 (MQ=60)
            TTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGG                                                                                                                                        <  SRR19973069.502210/150‑1 (MQ=60)
            TTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGG                                                                                                                                        <  SRR19973069.404157/150‑1 (MQ=60)
                      CGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCA                                                                                                                              <  SRR19973069.806594/150‑1 (MQ=60)
                         AAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCG                                                                                                                           >  SRR19973069.3416658/1‑150 (MQ=60)
                         AAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCG                                                                                                                           >  SRR19973069.3416924/1‑150 (MQ=60)
                          AACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGG                                                                                                                          >  SRR19973069.1207540/1‑150 (MQ=60)
                          AACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGG                                                                                                                          >  SRR19973069.1209196/1‑150 (MQ=60)
                            CATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCC                                                                                                                        >  SRR19973069.2716695/1‑150 (MQ=60)
                             ATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCA                                                                                                                       <  SRR19973069.87157/150‑1 (MQ=60)
                             ATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCA                                                                                                                       <  SRR19973069.910268/150‑1 (MQ=60)
                              TCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCTGTTCGCTCTGAATGTGCGGCGATTTCACCAGCGGTTTATCTTCCAGCGGCCAG                                                                                                                      <  SRR19973069.760352/150‑1 (MQ=60)
                              TCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAG                                                                                                                      <  SRR19973069.760353/150‑1 (MQ=60)
                                  CCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAC                                                                                                                  <  SRR19973069.2871593/150‑1 (MQ=60)
                                  CCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAC                                                                                                                  <  SRR19973069.86774/150‑1 (MQ=60)
                                  CCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAC                                                                                                                  <  SRR19973069.2871613/150‑1 (MQ=60)
                                  CCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAC                                                                                                                  <  SRR19973069.2871601/150‑1 (MQ=60)
                                   CAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACC                                                                                                                 >  SRR19973069.2609716/1‑150 (MQ=60)
                                   CAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACC                                                                                                                 >  SRR19973069.2170591/1‑150 (MQ=60)
                                       CGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCT                                                                                                             <  SRR19973069.2833830/150‑1 (MQ=60)
                                       CGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCT                                                                                                             <  SRR19973069.2833828/150‑1 (MQ=60)
                                       CGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCT                                                                                                             <  SRR19973069.440511/150‑1 (MQ=60)
                                       CGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCT                                                                                                             <  SRR19973069.440497/150‑1 (MQ=60)
                                          TTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTC                                                                                                          <  SRR19973069.2740037/150‑1 (MQ=60)
                                          TTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTC                                                                                                          <  SRR19973069.2740026/150‑1 (MQ=60)
                                           TCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCA                                                                                                         <  SRR19973069.64932/150‑1 (MQ=60)
                                           TCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCA                                                                                                         <  SRR19973069.286307/150‑1 (MQ=60)
                                           TCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCA                                                                                                         <  SRR19973069.287804/150‑1 (MQ=60)
                                           TCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCA                                                                                                         <  SRR19973069.2112689/150‑1 (MQ=60)
                                              CTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCT                                                                                                      <  SRR19973069.2467979/150‑1 (MQ=60)
                                               TGTCGGCCTGTATGTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTG                                                                                                     <  SRR19973069.1310528/150‑1 (MQ=60)
                                               TGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTG                                                                                                     <  SRR19973069.1310460/150‑1 (MQ=60)
                                                  CGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTC                                                                                                  >  SRR19973069.2794676/1‑150 (MQ=60)
                                                  CGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTC                                                                                                  <  SRR19973069.2373169/150‑1 (MQ=60)
                                                  CGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTC                                                                                                  <  SRR19973069.2373132/150‑1 (MQ=60)
                                                     CCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAT                                                                                               >  SRR19973069.2301768/1‑150 (MQ=60)
                                                     CCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAT                                                                                               >  SRR19973069.516679/1‑150 (MQ=60)
                                                       AGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTT                                                                                             >  SRR19973069.3371381/1‑150 (MQ=60)
                                                       AGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTT                                                                                             >  SRR19973069.1455785/1‑150 (MQ=60)
                                                        GTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTC                                                                                            >  SRR19973069.2884286/1‑150 (MQ=60)
                                                        GTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTACAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTC                                                                                            <  SRR19973069.1354002/150‑1 (MQ=60)
                                                                 TGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAT                                                                                   <  SRR19973069.2481976/150‑1 (MQ=60)
                                                                 TGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAT                                                                                   <  SRR19973069.549360/150‑1 (MQ=60)
                                                                 TGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAT                                                                                   <  SRR19973069.549369/150‑1 (MQ=60)
                                                                 TGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAT                                                                                   <  SRR19973069.549459/150‑1 (MQ=60)
                                                                   CCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAG                                                                                 <  SRR19973069.764057/150‑1 (MQ=60)
                                                                    CACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGTGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGC                                                                                <  SRR19973069.226706/150‑1 (MQ=60)
                                                                    CACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGTGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGC                                                                                <  SRR19973069.2563801/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.2533178/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.584192/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.418068/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.3385751/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.3158163/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.3385830/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.3158115/150‑1 (MQ=60)
                                                                      CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                              <  SRR19973069.2599572/150‑1 (MQ=60)
                                                                       ACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCACTTTCTGCGCGGATAGCCAG                                                                             <  SRR19973069.365678/150‑1 (MQ=60)
                                                                           GCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATC                                                                         <  SRR19973069.1803087/150‑1 (MQ=60)
                                                                              GGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTTCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCG                                                                      >  SRR19973069.2415207/1‑150 (MQ=60)
                                                                              GGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCG                                                                      <  SRR19973069.185150/150‑1 (MQ=60)
                                                                              GGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCG                                                                      <  SRR19973069.556715/150‑1 (MQ=60)
                                                                                   TACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGAT                                                                 <  SRR19973069.2256705/150‑1 (MQ=60)
                                                                                   TACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGAT                                                                 <  SRR19973069.315498/150‑1 (MQ=60)
                                                                                   GACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGAT                                                                 <  SRR19973069.2256711/150‑1 (MQ=60)
                                                                                        CAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGC                                                            <  SRR19973069.2630738/150‑1 (MQ=60)
                                                                                         AACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCG                                                           >  SRR19973069.787358/1‑150 (MQ=60)
                                                                                         AACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCG                                                           <  SRR19973069.2532960/150‑1 (MQ=60)
                                                                                           CTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGAT                                                         <  SRR19973069.2493114/150‑1 (MQ=60)
                                                                                              CACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCA                                                      <  SRR19973069.3836979/150‑1 (MQ=60)
                                                                                               ACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAG                                                     <  SRR19973069.2039603/150‑1 (MQ=60)
                                                                                               ACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAG                                                     <  SRR19973069.514689/150‑1 (MQ=60)
                                                                                                  GCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTC                                                  <  SRR19973069.1873089/150‑1 (MQ=60)
                                                                                                  GCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTC                                                  <  SRR19973069.1873099/150‑1 (MQ=60)
                                                                                                      TACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACT                                              <  SRR19973069.3367934/150‑1 (MQ=60)
                                                                                                        CGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTT                                            <  SRR19973069.2349618/150‑1 (MQ=60)
                                                                                                        CGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTT                                            >  SRR19973069.3320842/1‑150 (MQ=60)
                                                                                                           ATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCT                                         <  SRR19973069.1154139/150‑1 (MQ=60)
                                                                                                            TGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATATCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTT                                        >  SRR19973069.508515/1‑150 (MQ=60)
                                                                                                            TGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTT                                        >  SRR19973069.508813/1‑150 (MQ=60)
                                                                                                            TGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTT                                        <  SRR19973069.779044/150‑1 (MQ=60)
                                                                                                             GCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTT                                       <  SRR19973069.3615718/150‑1 (MQ=60)
                                                                                                                 CCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGA                                   >  SRR19973069.2669020/1‑150 (MQ=60)
                                                                                                                 CCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGA                                   >  SRR19973069.2668053/1‑150 (MQ=60)
                                                                                                                 CCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCGTTCAGA                                   >  SRR19973069.3015813/1‑150 (MQ=60)
                                                                                                                  CACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAT                                  >  SRR19973069.984075/1‑150 (MQ=60)
                                                                                                                  CACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAT                                  >  SRR19973069.745299/1‑150 (MQ=60)
                                                                                                                   ACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATT                                 <  SRR19973069.3765891/150‑1 (MQ=60)
                                                                                                                   ACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATT                                 <  SRR19973069.3561219/150‑1 (MQ=60)
                                                                                                                   ACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATT                                 <  SRR19973069.3119568/150‑1 (MQ=60)
                                                                                                                   ACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATT                                 <  SRR19973069.2199314/150‑1 (MQ=60)
                                                                                                                   ACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATT                                 <  SRR19973069.1569052/150‑1 (MQ=60)
                                                                                                                   ACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATT                                 <  SRR19973069.1569051/150‑1 (MQ=60)
                                                                                                                     TCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCA                               <  SRR19973069.3194418/150‑1 (MQ=60)
                                                                                                                     TCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCA                               >  SRR19973069.2787162/1‑150 (MQ=60)
                                                                                                                     TCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCA                               <  SRR19973069.3194857/150‑1 (MQ=60)
                                                                                                                     TCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCA                               <  SRR19973069.3336252/150‑1 (MQ=60)
                                                                                                                        GCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTC                            <  SRR19973069.272246/150‑1 (MQ=60)
                                                                                                                         CGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCG                           >  SRR19973069.726944/1‑150 (MQ=60)
                                                                                                                         CGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCG                           >  SRR19973069.2086727/1‑150 (MQ=60)
                                                                                                                         CGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCG                           <  SRR19973069.2965134/150‑1 (MQ=60)
                                                                                                                            CCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTT                        <  SRR19973069.3194055/150‑1 (MQ=60)
                                                                                                                              AGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCA                      >  SRR19973069.1492365/1‑150 (MQ=60)
                                                                                                                              AGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCA                      >  SRR19973069.854859/1‑150 (MQ=60)
                                                                                                                               GTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAA                     >  SRR19973069.860345/1‑150 (MQ=60)
                                                                                                                               GTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAA                     >  SRR19973069.311277/1‑150 (MQ=60)
                                                                                                                                   GCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCAT                 >  SRR19973069.1892155/1‑150 (MQ=60)
                                                                                                                                    CTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATC                >  SRR19973069.1536994/1‑150 (MQ=60)
                                                                                                                                    CTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATC                >  SRR19973069.3585065/1‑150 (MQ=60)
                                                                                                                                      CTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAG              <  SRR19973069.2144106/150‑1 (MQ=60)
                                                                                                                                      CTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAG              >  SRR19973069.207715/1‑150 (MQ=60)
                                                                                                                                       TGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGC             >  SRR19973069.3512548/1‑150 (MQ=60)
                                                                                                                                                  CGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC  >  SRR19973069.916936/1‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                                                   
GAGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC  >  CP009273/3039555‑3039850

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