Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F368 I0 R1 1472 214.4 7116750 95.9% 6824963 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,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,7030CA96.7% 572.5 / 9.1 184A900S (GCG→TCG) gcvPglycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base A (82/96);  total (84/100)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 189 total aligned reads displayed.

CAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC  >  CP009273/3039558‑3039850
                                                                                                                                                 |                                                                                                                                                   
   aaCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCa                                                                                                                                                >  2:18936/1‑148 (MQ=255)
   aaCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCa                                                                                                                                                <  2:846391/148‑1 (MQ=255)
   aaCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCa                                                                                                                                                >  2:18584/1‑148 (MQ=255)
    aCAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAc                                                                                                                                               >  1:431147/1‑148 (MQ=255)
           aCGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCggg                                                                                                                                        <  1:2201167/148‑1 (MQ=255)
               tCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTAt                                                                                                                                    >  1:2949252/1‑148 (MQ=255)
                   cGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTc                                                                                                                                <  1:2255109/148‑1 (MQ=255)
                      aaaCATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAg                                                                                                                             >  2:2392914/1‑148 (MQ=255)
                          atcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGc                                                                                                                         <  1:2011301/148‑1 (MQ=255)
                           tcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGcc                                                                                                                        >  2:3212681/1‑148 (MQ=255)
                           tcatcCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGcc                                                                                                                        >  2:2632917/1‑148 (MQ=255)
                                cAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGaca                                                                                                                   <  1:2387416/148‑1 (MQ=255)
                                 aGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGacac                                                                                                                  <  1:234000/148‑1 (MQ=255)
                                  gACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACAcc                                                                                                                 >  1:3066211/1‑148 (MQ=255)
                                      tttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCt                                                                                                             <  1:318207/148‑1 (MQ=255)
                                       ttCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCtt                                                                                                            <  2:1765520/148‑1 (MQ=255)
                                        tCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCttt                                                                                                           >  2:3554129/1‑148 (MQ=255)
                                        tCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCttt                                                                                                           >  2:2434421/1‑148 (MQ=255)
                                          accgTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACTGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGTTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCtttaa                                                                                                         <  1:1055991/145‑3 (MQ=255)
                                          aCTGTCGGCCAGTATTTGTCTGCCTCACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGAGGGTTATCTTACAGCGGCCAGACACCGGCTTTCa                                                                                                         <  1:2088082/148‑1 (MQ=255)
                                           cTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAc                                                                                                        <  2:2565505/148‑1 (MQ=255)
                                            tGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                                       <  1:3544631/148‑1 (MQ=255)
                                            tGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                                       <  1:413882/148‑1 (MQ=255)
                                            tGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                                       <  1:414622/148‑1 (MQ=255)
                                            tGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCAcc                                                                                                       >  1:2536283/1‑148 (MQ=255)
                                             gTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCt                                                                                                      >  2:2561645/1‑148 (MQ=255)
                                             gTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCt                                                                                                      >  2:779660/1‑148 (MQ=255)
                                                ggCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGt                                                                                                   <  2:678039/148‑1 (MQ=255)
                                                 gCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTc                                                                                                  >  2:1684273/1‑148 (MQ=255)
                                                 gCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTc                                                                                                  >  2:1684272/1‑148 (MQ=255)
                                                  ccAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGcc                                                                                                                        <  2:2402146/125‑1 (MQ=255)
                                                  ccAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                                 >  2:500230/1‑148 (MQ=255)
                                                  ccAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                                 >  1:1371949/1‑148 (MQ=255)
                                                  ccAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                                 >  1:1845525/1‑148 (MQ=255)
                                                  ccAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCa                                                                                                 >  2:2285739/1‑148 (MQ=255)
                                                   cAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCaa                                                                                                >  2:504177/1‑148 (MQ=255)
                                                    aGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAt                                                                                               >  1:1075701/1‑148 (MQ=255)
                                                    aGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAt                                                                                               <  1:508539/148‑1 (MQ=255)
                                                    aGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAt                                                                                               <  1:145370/148‑1 (MQ=255)
                                                          tGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGGACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTgc                                                                                         <  2:2480210/148‑1 (MQ=255)
                                                           gTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTgcg                                                                                        <  1:1420227/148‑1 (MQ=255)
                                                           gTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTgcg                                                                                        <  1:1967223/148‑1 (MQ=255)
                                                                ccACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAt                                                                                   >  2:1022783/1‑148 (MQ=255)
                                                                ccACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGAt                                                                                   <  2:2285034/148‑1 (MQ=255)
                                                                  acacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAg                                                                                 <  1:543615/148‑1 (MQ=255)
                                                                   cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                                >  2:2706441/1‑148 (MQ=255)
                                                                   cacCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGc                                                                                >  2:2135759/1‑148 (MQ=255)
                                                                    acCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGcc                                                                               <  1:746560/148‑1 (MQ=255)
                                                                    acCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGcc                                                                               <  1:2530036/148‑1 (MQ=255)
                                                                     ccTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCa                                                                              >  1:2180647/1‑148 (MQ=255)
                                                                      cTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAg                                                                             >  1:391325/1‑148 (MQ=255)
                                                                        gCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCa                                                                           >  1:2375039/1‑148 (MQ=255)
                                                                        gCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCa                                                                           <  2:611073/148‑1 (MQ=255)
                                                                            ggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgc                                                                       >  1:658112/1‑148 (MQ=255)
                                                                            ggAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATcgc                                                                       >  1:1790329/1‑148 (MQ=255)
                                                                                        cTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCg                                                           >  2:1241899/1‑148 (MQ=255)
                                                                                            aCGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                                       <  2:559719/148‑1 (MQ=255)
                                                                                            aCGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATcc                                                       >  2:59378/1‑148 (MQ=255)
                                                                                             cGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCa                                                      >  1:127157/1‑148 (MQ=255)
                                                                                              ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAg                                                     <  2:472412/148‑1 (MQ=255)
                                                                                              ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAg                                                     <  2:1048560/148‑1 (MQ=255)
                                                                                              ggCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAg                                                     <  2:215859/148‑1 (MQ=255)
                                                                                                cTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGtt                                                   <  2:2776072/148‑1 (MQ=255)
                                                                                                  gTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTcc                                                 <  2:1494663/148‑1 (MQ=255)
                                                                                                   tACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCa                                                >  2:1923725/1‑148 (MQ=255)
                                                                                                   tACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCa                                                >  2:1682914/1‑148 (MQ=255)
                                                                                                    aCGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCAc                                               >  2:1415955/1‑148 (MQ=255)
                                                                                                     cGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACt                                              <  2:2949252/148‑1 (MQ=255)
                                                                                                      ggATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACtt                                             <  2:415967/148‑1 (MQ=255)
                                                                                                      ggATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACtt                                             <  1:2881813/148‑1 (MQ=255)
                                                                                                        aTGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                           >  2:2866277/1‑148 (MQ=255)
                                                                                                        aTGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTg                                           <  2:1330069/148‑1 (MQ=255)
                                                                                                           cgcCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCtt                                        <  2:2874882/148‑1 (MQ=255)
                                                                                                             cccACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTc                                      <  2:585591/148‑1 (MQ=255)
                                                                                                                 cTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAt                                  <  1:766611/148‑1 (MQ=255)
                                                                                                                  tCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGAtt                                 <  2:2223057/148‑1 (MQ=255)
                                                                                                                   cGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTc                                >  2:3411971/1‑148 (MQ=255)
                                                                                                                     gCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAg                              <  2:2536283/148‑1 (MQ=255)
                                                                                                                      cGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGt                             >  1:939243/1‑148 (MQ=255)
                                                                                                                       gACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTc                            <  2:2617373/148‑1 (MQ=255)
                                                                                                                        aCCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCg                           >  2:2978174/1‑148 (MQ=255)
                                                                                                                        aCCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCg                           <  2:1594637/148‑1 (MQ=255)
                                                                                                                          cAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGt                         >  1:3286220/1‑148 (MQ=255)
                                                                                                                          cAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGt                         >  1:2752859/1‑148 (MQ=255)
                                                                                                                             ttCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCa                      <  2:2335577/148‑1 (MQ=255)
                                                                                                                              tCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCaa                     <  1:18936/148‑1 (MQ=255)
                                                                                                                              tCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCaa                     <  1:18584/148‑1 (MQ=255)
                                                                                                                               cGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAAc                    <  1:2251431/148‑1 (MQ=255)
                                                                                                                                 ctctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCAGTTTGCTTTCAGATTCAGTCGGTTCAACCa                  >  2:1790046/1‑148 (MQ=255)
                                                                                                                                 ctctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCa                  >  2:1124418/1‑148 (MQ=255)
                                                                                                                                 ctctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCa                  >  2:465323/1‑148 (MQ=255)
                                                                                                                                   ctGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATc                <  1:2915961/148‑1 (MQ=255)
                                                                                                                                     gAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAg              >  2:3028808/1‑148 (MQ=255)
                                                                                                                                       aTGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCg            >  2:2457878/1‑148 (MQ=255)
                                                                                                                                         gtgCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTa          <  2:3469172/148‑1 (MQ=255)
                                                                                                                                           gcggcgAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTAcc        >  1:2429802/1‑148 (MQ=255)
                                                                                                                                           gcggcgAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTNTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTAcn        <  1:2281801/148‑2 (MQ=255)
                                                                                                                                              gcgaGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGc     >  1:1380128/1‑148 (MQ=255)
                                                                                                                                                 cgTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAc  >  1:1633438/1‑148 (MQ=255)
                                                                                                                                                 cgTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAc  >  1:3546231/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC  >  CP009273/3039558‑3039850

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
AGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCG  >  CP009273/3039556‑3039852
                                                                                                                                                   |                                                                                                                                                     
AGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGT                                                                                                                                                     >  SRR19973067.764973/1‑150 (MQ=60)
  CAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTC                                                                                                                                                   <  SRR19973067.2296936/150‑1 (MQ=60)
    GAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAC                                                                                                                                                 >  SRR19973067.418854/1‑150 (MQ=60)
    GAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCAC                                                                                                                                                 >  SRR19973067.434144/1‑150 (MQ=60)
      ACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAGCTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCA                                                                                                                                               <  SRR19973067.1905583/150‑1 (MQ=60)
       CAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAG                                                                                                                                              <  SRR19973067.2835405/150‑1 (MQ=60)
          GTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGG                                                                                                                                           >  SRR19973067.2950221/1‑150 (MQ=60)
           TTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGG                                                                                                                                          >  SRR19973067.541022/1‑150 (MQ=60)
             ACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTT                                                                                                                                        <  SRR19973067.2215850/150‑1 (MQ=60)
               GGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTAT                                                                                                                                      >  SRR19973067.2967915/1‑150 (MQ=60)
                  CGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTT                                                                                                                                   <  SRR19973067.2501584/150‑1 (MQ=60)
                     CGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCA                                                                                                                                <  SRR19973067.2270079/150‑1 (MQ=60)
                            ATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCA                                                                                                                         <  SRR19973067.2025577/150‑1 (MQ=60)
                            ATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCA                                                                                                                         <  SRR19973067.2024835/150‑1 (MQ=60)
                             TCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAG                                                                                                                        >  SRR19973067.835613/1‑150 (MQ=60)
                                  CAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACC                                                                                                                   <  SRR19973067.2402754/150‑1 (MQ=60)
                                  CAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACC                                                                                                                   >  SRR19973067.3085729/1‑150 (MQ=60)
                                  CAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACC                                                                                                                   <  SRR19973067.2402922/150‑1 (MQ=60)
                                   AGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCG                                                                                                                  <  SRR19973067.235662/150‑1 (MQ=60)
                                       GTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTT                                                                                                              <  SRR19973067.1397872/150‑1 (MQ=60)
                                        TTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTT                                                                                                             <  SRR19973067.320354/150‑1 (MQ=60)
                                            ACTGTCGGCCAGTATTTGTCTGCCTCACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGAGGGTTATCTTACAGCGGCCAGACACCGGCTTTCACC                                                                                                         <  SRR19973067.2101901/150‑1 (MQ=60)
                                            ACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACC                                                                                                         >  SRR19973067.2552893/1‑150 (MQ=60)
                                            ACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACC                                                                                                         <  SRR19973067.2101914/150‑1 (MQ=60)
                                            ACCGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACTGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGTTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTAACC                                                                                                         <  SRR19973067.1063231/150‑1 (MQ=60)
                                             CTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCT                                                                                                        >  SRR19973067.2591859/1‑150 (MQ=60)
                                             CTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGACCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCT                                                                                                        <  SRR19973067.1432462/150‑1 (MQ=60)
                                              TGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTG                                                                                                       <  SRR19973067.3566999/150‑1 (MQ=60)
                                              TGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTG                                                                                                       <  SRR19973067.417500/150‑1 (MQ=60)
                                              TGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCNG                                                                                                       <  SRR19973067.416758/150‑1 (MQ=60)
                                                  GGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGAgatcggaagagcacacgtct                                                                                                   >  SRR19973067.2417824/1‑130 (MQ=60)
                                                  GGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCA                                                                                                   >  SRR19973067.1380904/1‑150 (MQ=60)
                                                  GGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCA                                                                                                   >  SRR19973067.1381316/1‑150 (MQ=60)
                                                  GGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCA                                                                                                   >  SRR19973067.1381612/1‑150 (MQ=60)
                                                  GGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCA                                                                                                   >  SRR19973067.1858105/1‑150 (MQ=60)
                                                    CCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAAT                                                                                                 >  SRR19973067.1083077/1‑150 (MQ=60)
                                                      AGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTT                                                                                               <  SRR19973067.146636/150‑1 (MQ=60)
                                                      AGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTT                                                                                               <  SRR19973067.512006/150‑1 (MQ=60)
                                                      AGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTT                                                                                               <  SRR19973067.511867/150‑1 (MQ=60)
                                                      AGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTT                                                                                               <  SRR19973067.1142123/150‑1 (MQ=60)
                                                            TGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGC                                                                                         <  SRR19973067.39478/150‑1 (MQ=60)
                                                             GTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCG                                                                                        <  SRR19973067.1980577/150‑1 (MQ=60)
                                                             GTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCG                                                                                        <  SRR19973067.1429825/150‑1 (MQ=60)
                                                             GTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCG                                                                                        >  SRR19973067.1696861/1‑150 (MQ=60)
                                                              TCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGG                                                                                       <  SRR19973067.2370767/150‑1 (MQ=60)
                                                               CTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGA                                                                                      >  SRR19973067.202093/1‑150 (MQ=60)
                                                                    ACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCC                                                                                 <  SRR19973067.547220/150‑1 (MQ=60)
                                                                     CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                                >  SRR19973067.2195227/1‑150 (MQ=60)
                                                                     CACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCA                                                                                >  SRR19973067.2195381/1‑150 (MQ=60)
                                                                      ACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAG                                                                               <  SRR19973067.2546612/150‑1 (MQ=60)
                                                                      ACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAG                                                                               <  SRR19973067.751743/150‑1 (MQ=60)
                                                                      ACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAG                                                                               >  SRR19973067.394064/1‑150 (MQ=60)
                                                                      ACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAG                                                                               >  SRR19973067.394071/1‑150 (MQ=60)
                                                                        CTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCA                                                                             >  SRR19973067.2390477/1‑150 (MQ=60)
                                                                          GCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATC                                                                           >  SRR19973067.3189384/1‑150 (MQ=60)
                                                                            CGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGC                                                                         >  SRR19973067.662585/1‑150 (MQ=60)
                                                                            CGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGC                                                                         >  SRR19973067.1802394/1‑150 (MQ=60)
                                                                            CGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGC                                                                         >  SRR19973067.663010/1‑150 (MQ=60)
                                                                               GAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTC                                                                      <  SRR19973067.16485/150‑1 (MQ=60)
                                                                                AATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCG                                                                     <  SRR19973067.659125/150‑1 (MQ=60)
                                                                                 ATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGA                                                                    >  SRR19973067.2547842/1‑150 (MQ=60)
                                                                                             CACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCA                                                        >  SRR19973067.128325/1‑150 (MQ=60)
                                                                                                  CTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCC                                                   >  SRR19973067.3553698/1‑150 (MQ=60)
                                                                                                        GGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTG                                             <  SRR19973067.2900242/150‑1 (MQ=60)
                                                                                                              GCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCA                                       <  SRR19973067.3250486/150‑1 (MQ=60)
                                                                                                                   CTCGGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTC                                  <  SRR19973067.771934/150‑1 (MQ=60)
                                                                                                                      GGCGACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGT                               >  SRR19973067.945631/1‑150 (MQ=60)
                                                                                                                         GACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGG                            >  SRR19973067.1911500/1‑150 (MQ=60)
                                                                                                                          ACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGT                           >  SRR19973067.3306852/1‑150 (MQ=60)
                                                                                                                          ACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGT                           >  SRR19973067.1831176/1‑150 (MQ=60)
                                                                                                                          ACCAGTTCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGT                           >  SRR19973067.2770644/1‑150 (MQ=60)
                                                                                                                                TCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACC                     <  SRR19973067.18862/150‑1 (MQ=60)
                                                                                                                                TCGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACC                     <  SRR19973067.19221/150‑1 (MQ=60)
                                                                                                                                 CGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCA                    <  SRR19973067.2266379/150‑1 (MQ=60)
                                                                                                                                 CGCTCTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCA                    <  SRR19973067.2266438/150‑1 (MQ=60)
                                                                                                                                     CTGAATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAG                <  SRR19973067.2934519/150‑1 (MQ=60)
                                                                                                                                        AATGTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGT             <  SRR19973067.221866/150‑1 (MQ=60)
                                                                                                                                           GTGCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACC          >  SRR19973067.2445748/1‑150 (MQ=60)
                                                                                                                                             GCGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTNTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACNCG        <  SRR19973067.2296888/150‑1 (MQ=60)
                                                                                                                                              CGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGC       >  SRR19973067.1389535/1‑150 (MQ=60)
                                                                                                                                              CGGCGAGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGC       >  SRR19973067.3077229/1‑150 (MQ=60)
                                                                                                                                                 CGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC    >  SRR19973067.3568606/1‑150 (MQ=60)
                                                                                                                                                 CGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC    >  SRR19973067.1644674/1‑150 (MQ=60)
                                                                                                                                                 CGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCAC    >  SRR19973067.3568585/1‑150 (MQ=60)
                                                                                                                                                   AGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCG  <  SRR19973067.3417023/150‑1 (MQ=60)
                                                                                                                                                   |                                                                                                                                                     
AGCAGAACAGGTTACGGTCGCCGTAAACATCATCCAGACGTTTCACTGTCGGCCAGTATTTGTCTGCCACACCTGCCGGGAATACCGCAACTTCACGGCTGTACGGATGCGCCCACTCGGCGACCAGTTCGCTCTGAATGTGCGGCGCGTTCACCAGCGGGTTATCTTCCAGCGGCCAGACACCGGCTTTCACCTGGTCAATTTCTGCGCGGATAGCCAGCATCGCGTCGATAAAGCGATCCAGTTCCACTTTGCTTTCAGATTCAGTCGGTTCAACCATCAGCGTACCCGCCACCG  >  CP009273/3039556‑3039852

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