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 F394 I0 R1 1342 215.0 7134094 96.2% 6862998 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,717,788 C→A 100% P297P (CCG→CCT)  glyQ ← glycine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,717,7880CA99.4% 566.3 / ‑4.2 176P297P (CCG→CCT) glyQglycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (102/73);  total (102/74)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 176 total aligned reads displayed.

CAGCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATCGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTT  >  CP009273/3717643‑3717931
                                                                                                                                                 |                                                                                                                                               
 aGCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAATAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAggg                                                                                                                                              <  2:372326/148‑1 (MQ=255)
  gCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGa                                                                                                                                             <  2:1118819/148‑1 (MQ=255)
   cGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGaa                                                                                                                                            >  2:425839/1‑148 (MQ=255)
      tATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGcc                                                                                                                                         <  2:2606637/148‑1 (MQ=255)
         ccAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAg                                                                                                                                      >  2:2616066/1‑148 (MQ=255)
            gCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGc                                                                                                                                   <  2:3169759/148‑1 (MQ=255)
                 gCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCAc                                                                                                                              <  1:2530508/148‑1 (MQ=255)
                  cAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACg                                                                                                                             <  1:922292/148‑1 (MQ=255)
                     tAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGa                                                                                                                          >  2:1205001/1‑148 (MQ=255)
                     tAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGa                                                                                                                          >  2:1625687/1‑148 (MQ=255)
                     tAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGa                                                                                                                          >  2:1625576/1‑148 (MQ=255)
                     tAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGa                                                                                                                          >  2:43756/1‑148 (MQ=255)
                          ttCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCgt                                                                                                                     >  1:978491/1‑148 (MQ=255)
                             gcagcaAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCgtagt                                                                                                                  <  2:6495/148‑1 (MQ=255)
                             gcagcaAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCgtagt                                                                                                                  <  2:6455/148‑1 (MQ=255)
                              cagcaAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCgtagta                                                                                                                 <  1:1218182/148‑1 (MQ=255)
                                   aaGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTTCCGATTTACACCAGAAAAGTTTTCTAAGACATAGCCGCCTCTTACTTATCTTTGTATCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGGAGTATGCtt                                                                                                            >  2:3068797/1‑148 (MQ=255)
                                   aaGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCtt                                                                                                            >  2:3068884/1‑148 (MQ=255)
                                   aaGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCtt                                                                                                            >  2:1831502/1‑148 (MQ=255)
                                    aGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTc                                                                                                           >  2:808772/1‑148 (MQ=255)
                                    aGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTc                                                                                                           >  1:2847091/1‑148 (MQ=255)
                                    aGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCGTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGACGCGTAGTATGCTTc                                                                                                           >  2:808724/1‑148 (MQ=255)
                                         tCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCa                                                                                                      <  1:2560828/148‑1 (MQ=255)
                                         tCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCa                                                                                                      <  2:32485/148‑1 (MQ=255)
                                          cAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCAc                                                                                                     >  1:2758895/1‑148 (MQ=255)
                                          cAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCAc                                                                                                     >  1:2758847/1‑148 (MQ=255)
                                            gCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTg                                                                                                   <  1:1462587/148‑1 (MQ=255)
                                            gCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTg                                                                                                   <  2:1583313/148‑1 (MQ=255)
                                             ccAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGc                                                                                                  <  2:1167856/148‑1 (MQ=255)
                                             ccAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGc                                                                                                  >  1:674408/1‑148 (MQ=255)
                                             ccAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGc                                                                                                  >  2:388063/1‑148 (MQ=255)
                                             ccAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGc                                                                                                  <  2:1167897/148‑1 (MQ=255)
                                              cAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCt                                                                                                 <  1:1063839/148‑1 (MQ=255)
                                               aGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCtt                                                                                                <  2:1525492/148‑1 (MQ=255)
                                                ggCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCttt                                                                                               >  2:921133/1‑148 (MQ=255)
                                                  cTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTgg                                                                                             <  1:458179/148‑1 (MQ=255)
                                                   tGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGt                                                                                            <  1:3493570/148‑1 (MQ=255)
                                                       cAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAgg                                                                                        >  2:1762853/1‑148 (MQ=255)
                                                          tGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTg                                                                                     >  2:3357737/1‑148 (MQ=255)
                                                           gCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGc                                                                                    >  2:2173191/1‑148 (MQ=255)
                                                           gCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGc                                                                                    >  1:3303358/1‑148 (MQ=255)
                                                             ttttGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGa                                                                                  <  1:352050/148‑1 (MQ=255)
                                                              tttGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGaa                                                                                 >  2:2384684/1‑148 (MQ=255)
                                                              tttGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGaa                                                                                 <  1:2813503/148‑1 (MQ=255)
                                                               ttGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAAt                                                                                <  1:3377486/148‑1 (MQ=255)
                                                                  gtggCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATgcg                                                                             <  1:1181039/148‑1 (MQ=255)
                                                                   tggCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATgcgc                                                                            <  1:1584628/148‑1 (MQ=255)
                                                                      cAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGa                                                                         <  1:3440239/148‑1 (MQ=255)
                                                                        gCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAAt                                                                       >  2:3546396/1‑148 (MQ=255)
                                                                        gCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAAt                                                                       <  1:1541019/148‑1 (MQ=255)
                                                                           ctTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTa                                                                    >  2:912742/1‑148 (MQ=255)
                                                                           ctTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTa                                                                    >  2:2943095/1‑148 (MQ=255)
                                                                            ttCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAg                                                                   <  1:1756540/148‑1 (MQ=255)
                                                                             tCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAgc                                                                  >  2:3525717/1‑148 (MQ=255)
                                                                              cAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAgcg                                                                 <  2:978491/148‑1 (MQ=255)
                                                                                gTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGGAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGACTGCGCAGAATGTAGCGCt                                                               >  2:905415/1‑148 (MQ=255)
                                                                                gTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGAACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCt                                                               >  1:3069786/1‑148 (MQ=255)
                                                                                  gCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGa                                                             <  1:2447437/148‑1 (MQ=255)
                                                                                     gATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGc                                                          <  2:101546/148‑1 (MQ=255)
                                                                                      aTTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCt                                                         >  2:1539207/1‑148 (MQ=255)
                                                                                      aTTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCt                                                         >  2:1872720/1‑148 (MQ=255)
                                                                                       tttCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTc                                                        >  2:291610/1‑148 (MQ=255)
                                                                                       tttCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTAGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTc                                                        >  1:1365973/1‑148 (MQ=255)
                                                                                        ttCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCg                                                       <  2:3324639/148‑1 (MQ=255)
                                                                                            accaGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGa                                                   >  2:806443/1‑148 (MQ=255)
                                                                                             ccaGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGAc                                                  >  1:2721588/1‑148 (MQ=255)
                                                                                                 aaaaGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGgag                                              >  2:2457879/1‑148 (MQ=255)
                                                                                                   aaGTTTTCTCAGACATAGCCGCCTCTTACTGATCTTTGTTGCACATAGGGAAGCCGAGGGATTCACGGGAAGCGTAGTATGCTTCTGCAACTGCTTTGGTCAGGGTGCGAATGCGAAGAATGTAGCGCTGACGCTCGGTGACGGAGAt                                            >  2:1108720/1‑148 (MQ=255)
                                                                                                    aGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATg                                           >  1:1079329/1‑148 (MQ=255)
                                                                                                      ttttCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGc                                         >  2:2612194/1‑148 (MQ=255)
                                                                                                      ttttCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGc                                         >  2:1238326/1‑148 (MQ=255)
                                                                                                      ttttCTCAGACATAGCCGCCTATTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGc                                         >  2:2727533/1‑148 (MQ=255)
                                                                                                       tttCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCt                                        >  1:3327310/1‑148 (MQ=255)
                                                                                                         tctcAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCttt                                      >  1:1894633/1‑148 (MQ=255)
                                                                                                           tcAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTAc                                    >  1:3308715/1‑148 (MQ=255)
                                                                                                               aCATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACgcgc                                >  2:3475274/1‑148 (MQ=255)
                                                                                                               aCATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACgcgc                                >  2:3475263/1‑148 (MQ=255)
                                                                                                               aCATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACgcgc                                >  2:1553310/1‑148 (MQ=255)
                                                                                                                cATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCa                               >  2:1490585/1‑148 (MQ=255)
                                                                                                                cATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCa                               >  2:1490390/1‑148 (MQ=255)
                                                                                                                   aGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATcc                            <  1:2173191/148‑1 (MQ=255)
                                                                                                                    gccgccTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCca                           >  2:1265716/1‑148 (MQ=255)
                                                                                                                     ccgccTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCcag                          >  2:1859656/1‑148 (MQ=255)
                                                                                                                       gccTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCcagca                        >  2:1742861/1‑148 (MQ=255)
                                                                                                                        ccTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCcagcag                       >  1:2532305/1‑148 (MQ=255)
                                                                                                                              acttaTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGaa                 >  2:1604418/1‑148 (MQ=255)
                                                                                                                              acttaTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGaa                 <  2:2847091/148‑1 (MQ=255)
                                                                                                                              acttaTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGaa                 <  1:3172454/148‑1 (MQ=255)
                                                                                                                               cttaTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAg                <  2:2407402/148‑1 (MQ=255)
                                                                                                                               cttaTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAg                <  1:921133/148‑1 (MQ=255)
                                                                                                                                ttaTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGc               >  1:2989660/1‑148 (MQ=255)
                                                                                                                                ttaTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGc               >  1:266013/1‑148 (MQ=255)
                                                                                                                                 taTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCt              >  2:157714/1‑148 (MQ=255)
                                                                                                                                 taTCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCt              >  2:3253183/1‑148 (MQ=255)
                                                                                                                                   tCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCtgt            <  2:1097224/148‑1 (MQ=255)
                                                                                                                                    cTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCtgtg           >  1:2176626/1‑148 (MQ=255)
                                                                                                                                     tttGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTgg          <  1:3389222/148‑1 (MQ=255)
                                                                                                                                      ttgttgCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTggg         <  1:1625687/148‑1 (MQ=255)
                                                                                                                                      ttgttgCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTggg         <  1:1625576/148‑1 (MQ=255)
                                                                                                                                         ttgCACAGAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGggcgg      <  2:1084432/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAGCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATCGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTT  >  CP009273/3717643‑3717931

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
AGGCCAGCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATCGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAAT  >  CP009273/3717639‑3717937
                                                                                                                                                     |                                                                                                                                                     
AGGCCAGCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATA                                                                                                                                                       >  SRR19973066.1092351/1‑150 (MQ=60)
   CCAGCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGG                                                                                                                                                    <  SRR19973066.3365226/150‑1 (MQ=60)
        GTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCC                                                                                                                                               >  SRR19973066.1337414/1‑150 (MQ=60)
             CCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATCGGGAAGCCGAGGG                                                                                                                                          <  SRR19973066.3456530/150‑1 (MQ=60)
                GCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGACGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTT                                                                                                                                       >  SRR19973066.3498018/1‑150 (MQ=60)
                    CGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACG                                                                                                                                   >  SRR19973066.2788894/1‑150 (MQ=60)
                     GCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGG                                                                                                                                  <  SRR19973066.2547716/150‑1 (MQ=60)
                      CAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGG                                                                                                                                 <  SRR19973066.928940/150‑1 (MQ=60)
                            AGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGT                                                                                                                           >  SRR19973066.985782/1‑150 (MQ=60)
                            AGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGT                                                                                                                           >  SRR19973066.985705/1‑150 (MQ=60)
                                  CAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATG                                                                                                                     <  SRR19973066.1226764/150‑1 (MQ=60)
                                    GCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGGAGTATGCT                                                                                                                   >  SRR19973066.2469350/1‑150 (MQ=60)
                                     CAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTT                                                                                                                  >  SRR19973066.3027432/1‑150 (MQ=60)
                                     CAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTT                                                                                                                  >  SRR19973066.2725691/1‑150 (MQ=60)
                                      AAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTC                                                                                                                 >  SRR19973066.2866220/1‑150 (MQ=60)
                                            CTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCAC                                                                                                           >  SRR19973066.2777420/1‑150 (MQ=60)
                                            CTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCAC                                                                                                           >  SRR19973066.2777468/1‑150 (MQ=60)
                                             TCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACT                                                                                                          <  SRR19973066.2578209/150‑1 (MQ=60)
                                               AGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGC                                                                                                        >  SRR19973066.679345/1‑150 (MQ=60)
                                                GCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCT                                                                                                       <  SRR19973066.1472927/150‑1 (MQ=60)
                                                 CCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTT                                                                                                      >  SRR19973066.2665690/1‑150 (MQ=60)
                                                  CAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTT                                                                                                     <  SRR19973066.1071636/150‑1 (MQ=60)
                                                      CTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTC                                                                                                 >  SRR19973066.1124795/1‑150 (MQ=60)
                                                      CTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTC                                                                                                 <  SRR19973066.461693/150‑1 (MQ=60)
                                                       TGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCA                                                                                                <  SRR19973066.3516869/150‑1 (MQ=60)
                                                             GTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGC                                                                                          >  SRR19973066.3325386/1‑150 (MQ=60)
                                                             GTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGC                                                                                          >  SRR19973066.3325210/1‑150 (MQ=60)
                                                                 TTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAAT                                                                                      <  SRR19973066.354793/150‑1 (MQ=60)
                                                                  TTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATG                                                                                     <  SRR19973066.2833075/150‑1 (MQ=60)
                                                                  TTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATG                                                                                     <  SRR19973066.2832425/150‑1 (MQ=60)
                                                                   TTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGC                                                                                    <  SRR19973066.3400010/150‑1 (MQ=60)
                                                                      GTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCA                                                                                 <  SRR19973066.1189341/150‑1 (MQ=60)
                                                                       TGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAG                                                                                <  SRR19973066.1596101/150‑1 (MQ=60)
                                                                       TGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAG                                                                                <  SRR19973066.1596086/150‑1 (MQ=60)
                                                                       TGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAG                                                                                <  SRR19973066.1596064/150‑1 (MQ=60)
                                                                          CAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAAT                                                                             <  SRR19973066.3463252/150‑1 (MQ=60)
                                                                            GCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGT                                                                           <  SRR19973066.1552032/150‑1 (MQ=60)
                                                                               CTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGC                                                                        >  SRR19973066.930780/1‑150 (MQ=60)
                                                                               CTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGC                                                                        >  SRR19973066.2549199/1‑150 (MQ=60)
                                                                               CTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGC                                                                        >  SRR19973066.3404639/1‑150 (MQ=60)
                                                                                TTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCG                                                                       <  SRR19973066.1769161/150‑1 (MQ=60)
                                                                                  CAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGAACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCT                                                                     >  SRR19973066.3090311/1‑150 (MQ=60)
                                                                                      GCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACG                                                                 <  SRR19973066.2464109/150‑1 (MQ=60)
                                                                                         GATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTAGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTC                                                              >  SRR19973066.1375851/1‑150 (MQ=60)
                                                                                               CACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGAC                                                        >  SRR19973066.2739896/1‑150 (MQ=60)
                                                                                                     AAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGAT                                                  <  SRR19973066.291741/150‑1 (MQ=60)
                                                                                                      AAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATG                                                 >  SRR19973066.1087218/1‑150 (MQ=60)
                                                                                                       AAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGG                                                <  SRR19973066.791203/150‑1 (MQ=60)
                                                                                                         GTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCT                                              >  SRR19973066.3349395/1‑150 (MQ=60)
                                                                                                         GTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCT                                              >  SRR19973066.3349413/1‑150 (MQ=60)
                                                                                                          TTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTT                                             <  SRR19973066.2751207/150‑1 (MQ=60)
                                                                                                           TTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTT                                            >  SRR19973066.1908143/1‑150 (MQ=60)
                                                                                                            TTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTA                                           <  SRR19973066.209978/150‑1 (MQ=60)
                                                                                                             TCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTAC                                          >  SRR19973066.3330639/1‑150 (MQ=60)
                                                                                                                       AGCCGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAG                                <  SRR19973066.2188283/150‑1 (MQ=60)
                                                                                                                          CGCCTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAG                             >  SRR19973066.2549520/1‑150 (MQ=60)
                                                                                                                             CTCTTACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTT                          <  SRR19973066.1746198/150‑1 (MQ=60)
                                                                                                                                  ACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGC                     >  SRR19973066.268002/1‑150 (MQ=60)
                                                                                                                                  ACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGC                     <  SRR19973066.3193868/150‑1 (MQ=60)
                                                                                                                                  ACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGC                     <  SRR19973066.3193579/150‑1 (MQ=60)
                                                                                                                                  ACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGC                     >  SRR19973066.3009448/1‑150 (MQ=60)
                                                                                                                                  ACTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGC                     <  SRR19973066.1846897/150‑1 (MQ=60)
                                                                                                                                   CTTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCT                    <  SRR19973066.927767/150‑1 (MQ=60)
                                                                                                                                    TTATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTG                   <  SRR19973066.575826/150‑1 (MQ=60)
                                                                                                                                      ATCTTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTG                 >  SRR19973066.2191739/1‑150 (MQ=60)
                                                                                                                                         TTTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGC              <  SRR19973066.3411843/150‑1 (MQ=60)
                                                                                                                                          TTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCG             <  SRR19973066.1637361/150‑1 (MQ=60)
                                                                                                                                          TTGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCG             <  SRR19973066.1637473/150‑1 (MQ=60)
                                                                                                                                           TGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGG            >  SRR19973066.1203178/1‑150 (MQ=60)
                                                                                                                                           TGTTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGG            >  SRR19973066.381119/1‑150 (MQ=60)
                                                                                                                                             TTGCACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCT          >  SRR19973066.2241432/1‑150 (MQ=60)
                                                                                                                                                 ACATAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCA      <  SRR19973066.1544049/150‑1 (MQ=60)
                                                                                                                                                    TAGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAA   >  SRR19973066.1653681/1‑150 (MQ=60)
                                                                                                                                                     AGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAAT  <  SRR19973066.1327039/150‑1 (MQ=60)
                                                                                                                                                     AGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAAT  >  SRR19973066.3164345/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                     
AGGCCAGCGTTATCCAGCTCCGCAGTAAAGTTCGCAGCAAAGGACTCAGCCAGGCTGCGCAGTGCTTTTGGTGGCAGCTCTTCAGTGCCGATTTCCACCAGAAAAGTTTTCTCAGACATAGCCGCCTCTTACTTATCTTTGTTGCACATCGGGAAGCCGAGGGCTTCACGGGAAGCGTAGTATGCTTCTGCCACTGCTTTGGTCAGGGTGCGAATGCGCAGAATGTAGCGCTGACGCTCGGTGACGGAGATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAAT  >  CP009273/3717639‑3717937

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