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
A2 F5 I179 R1 182 40.3 1547118 96.8% 1497610 120.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 583,595 A→G V939A (GTC→GCC)  nfrA ← bacteriophage adsorption protein NfrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273583,5950AG95.5% 140.7 / ‑0.8 44V939A (GTC→GCC) nfrAbacteriophage adsorption protein NfrA
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (23/19);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTG  >  NZ_CP009273/583460‑583728
                                                                                                                                       |                                                                                                                                     
agaAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTc                                                                                                                                    >  1:563521/1‑139 (MQ=255)
        gttgttGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGc                                                                                                                            >  2:684568/1‑139 (MQ=255)
                         tACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGcc                                                                                                           >  2:510429/1‑139 (MQ=255)
                            gTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCttt                                                                                                        <  1:337962/139‑1 (MQ=255)
                                aTTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAg                                                                                                               <  1:530049/128‑1 (MQ=255)
                                aTTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAg                                                                                                               >  2:530049/1‑128 (MQ=255)
                                aTTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCa                                                                                                    >  2:527160/1‑139 (MQ=255)
                                         cTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           <  2:649839/139‑1 (MQ=255)
                                          ttAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTAt                                                                                                      >  1:341089/1‑127 (MQ=255)
                                          ttAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTAt                                                                                                      <  2:341089/127‑1 (MQ=255)
                                                   tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGa                                                                                                                <  1:50026/108‑1 (MQ=255)
                                                   tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGa                                                                                                                >  2:50026/1‑108 (MQ=255)
                                                   tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGAt                                                                                                >  2:591188/1‑124 (MQ=255)
                                                   tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGAt                                                                                                <  1:591188/124‑1 (MQ=255)
                                                   tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGaa                                                                                 >  2:576455/1‑139 (MQ=255)
                                                   tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGaa                                                                                 >  2:230271/1‑139 (MQ=255)
                                                   tGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGaa                                                                                 >  2:11754/1‑139 (MQ=255)
                                                              gACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTa                                                                       <  1:510429/138‑1 (MQ=255)
                                                               aCGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAg                                                                     >  1:544982/1‑139 (MQ=255)
                                                                cGCCGAGACTGACTTTGGGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCCCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGCAAgg                                                                    >  2:497742/1‑139 (MQ=255)
                                                                            ctTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATACTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTg                                                        <  2:544982/139‑1 (MQ=255)
                                                                                gtgCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCg                                                    >  1:634659/1‑139 (MQ=255)
                                                                                     gCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGc                                               >  2:269625/1‑139 (MQ=255)
                                                                                          gCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCt                                          >  2:156055/1‑139 (MQ=255)
                                                                                               gtAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGa                                     >  1:339804/1‑139 (MQ=255)
                                                                                                    gCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           <  1:79767/80‑1 (MQ=255)
                                                                                                    gCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           >  2:79767/1‑80 (MQ=255)
                                                                                                     cGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAg                               >  1:217072/1‑139 (MQ=255)
                                                                                                       tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           >  2:139910/1‑77 (MQ=255)
                                                                                                       tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCg                                                                                           <  1:139910/77‑1 (MQ=255)
                                                                                                       tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCt                             >  2:493381/1‑139 (MQ=255)
                                                                                                       tctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCt                             <  1:527160/139‑1 (MQ=255)
                                                                                                        ctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGt                                                                              >  1:765761/1‑89 (MQ=255)
                                                                                                        ctcGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGt                                                                              <  2:765761/89‑1 (MQ=255)
                                                                                                           gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGAcg                         <  1:684568/139‑1 (MQ=255)
                                                                                                           gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGAcg                         <  2:560305/139‑1 (MQ=255)
                                                                                                           gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGAcg                         <  1:450996/139‑1 (MQ=255)
                                                                                                           gCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGACAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGa                           <  1:657835/137‑1 (MQ=255)
                                                                                                                           cAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTc          <  2:353134/138‑1 (MQ=255)
                                                                                                                           cAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTc          <  1:558349/138‑1 (MQ=255)
                                                                                                                              cGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTa                                                                       <  2:511092/74‑1 (MQ=255)
                                                                                                                              cGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTa                                                                       >  1:511092/1‑74 (MQ=255)
                                                                                                                                   cccGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTg  >  2:350058/1‑138 (MQ=255)
                                                                                                                                   cccGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTg  >  1:555906/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                     
AGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTG  >  NZ_CP009273/583460‑583728

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAAT  >  NZ_CP009273/583449‑583731
                                                                                                                                                  |                                                                                                                                        
CTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCC                                                                                                                                      >  SRR3722133.571004/1‑151 (MQ=60)
                              ggagatgtgtataagagacaGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATC                                                                                                        <  SRR3722133.50654/131‑1 (MQ=60)
                               GCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCA                                                                                                       <  SRR3722133.537145/151‑1 (MQ=60)
                                      CGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGC                                                                                                <  SRR3722133.342717/151‑1 (MQ=60)
                                          GATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGCT                                                                                            >  SRR3722133.345875/1‑151 (MQ=60)
                                              gacaGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCC                                                                                        <  SRR3722133.598997/147‑1 (MQ=60)
                                                   catacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAAC                                                                                   <  SRR3722133.141760/99‑1 (MQ=60)
                                                   acgagatctagtacggtctcgtgggctcggagatgtgtataagagacagGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAAC                                                                                   <  SRR3722133.80776/102‑1 (MQ=60)
                                                               GTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGG                                                                       >  SRR3722133.552239/1‑151 (MQ=60)
                                                                        CGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAG                                                               <  SRR3722133.517301/150‑1 (MQ=60)
                                                                                AGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGT                                                      >  SRR3722133.642978/1‑151 (MQ=60)
                                                                                               GGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGAT                                       >  SRR3722133.344574/1‑151 (MQ=60)
                                                                                                     GCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGC                                 >  SRR3722133.219999/1‑151 (MQ=60)
                                                                                                        TCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCctgtctcttatacacatctccgagcccacgagaccgtact                              >  SRR3722133.775657/1‑111 (MQ=60)
                                                                                                                 GTCTCGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAAT                     <  SRR3722133.534237/151‑1 (MQ=60)
                                                                                                                     CGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGC                 <  SRR3722133.457125/151‑1 (MQ=60)
                                                                                                                     CGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGC                 <  SRR3722133.693426/151‑1 (MQ=60)
                                                                                                                     CGCCGGTCCAGATATTCCAGCGGACCCCGGCTCCGACAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCT                   <  SRR3722133.666399/149‑1 (MQ=60)
                                                                                                                              AGATATTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccg        >  SRR3722133.517976/1‑96 (MQ=60)
                                                                                                                                   TTCCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGA    >  SRR3722133.563322/1‑150 (MQ=60)
                                                                                                                                     CCAGCGGACCCCGGCTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAAT  <  SRR3722133.565783/150‑1 (MQ=60)
                                                                                                                                                  |                                                                                                                                        
CTCCAATGGTGAGAAACGCGTTGTTGCGCTCTCCGTTACGTTGATTAATCGCCTTAAAGGTATGTTGATACTCGACGCCGAGACTGACTTTGTGCGGCCAGGCGTCGTAGTGCGTCTCGCCGGTCCAGATATTCCAGCGGACCCCGACTCCGCCAAGCTGCGCGCCCTGAGTGCCTTTATCACGATAGCCGTTGTCCTGAACGTGAGCGTAAGGCTCAATAGTCTGTCCGTTAGCTACCTTCTGATGCCAGCTGACGCGATAATCTGCCGTCCACGCCTGAAT  >  NZ_CP009273/583449‑583731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: