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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 661,892 A→C G594G (GGT→GGG mrdA ← peptidoglycan DD‑transpeptidase MrdA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273661,8920AC100.0% 191.7 / NA 55G594G (GGT→GGGmrdApeptidoglycan DD‑transpeptidase MrdA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (30/25);  total (30/25)

TTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGC  >  NZ_CP009273/661805‑661976
                                                                                       |                                                                                    
ttCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCaccccc                                                                                    >  2:1054256/1‑90 (MQ=255)
  ccgccgGCAGATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                  <  1:682221/90‑1 (MQ=255)
  ccGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCGTCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                  <  2:368584/90‑1 (MQ=255)
     caggcagATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTtct                                                                               >  1:537225/1‑90 (MQ=255)
       ggcagATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTcc                                                                             <  1:874480/90‑1 (MQ=255)
         cagATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAg                                                                           <  1:761940/90‑1 (MQ=255)
           gATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGaa                                                                         >  1:864906/1‑90 (MQ=255)
                gtgtTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATg                                                                    <  1:54221/90‑1 (MQ=255)
                     gttATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                  >  1:1173930/1‑71 (MQ=255)
                     gttATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                  <  2:1173930/71‑1 (MQ=255)
                     gttATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGAc                                                               >  1:407651/1‑90 (MQ=255)
                        aTCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGc                                                            <  2:743656/90‑1 (MQ=255)
                          cACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCa                                                          >  1:824298/1‑90 (MQ=255)
                              cAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTg                                                      >  2:1232416/1‑90 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                     >  1:388548/1‑90 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                     >  1:886266/1‑90 (MQ=255)
                               aGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                     >  1:1168602/1‑90 (MQ=255)
                                 cATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                   <  1:43032/90‑1 (MQ=255)
                                 cATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                   <  1:126025/90‑1 (MQ=255)
                                 cATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                   <  2:480619/90‑1 (MQ=255)
                                     aTGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAttg                                               >  1:321171/1‑90 (MQ=255)
                                     aTGTGGTCGAGGATCTGGCACATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAttg                                               >  1:1139263/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACATCCACTTGCGTATTGTTGTa                                          >  1:146396/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                     <  2:1010199/79‑1 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGc                                                     >  1:1010199/1‑79 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                   <  1:964211/81‑1 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                   <  2:653339/81‑1 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                   >  2:964211/1‑81 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCgg                                                   >  1:653339/1‑81 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                          >  2:86979/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                          >  1:1031725/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                          >  2:100705/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                          >  2:734241/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                          >  1:394542/1‑90 (MQ=255)
                                          gTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTa                                          >  1:196759/1‑90 (MQ=255)
                                           tCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                                         <  1:929776/90‑1 (MQ=255)
                                           tCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                                         <  1:1011803/90‑1 (MQ=255)
                                           tCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAt                                         <  2:1131933/90‑1 (MQ=255)
                                                 aTCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCa                                   >  1:44701/1‑90 (MQ=255)
                                                    tGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                  <  1:449281/40‑1 (MQ=255)
                                                    tGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGt                                                                                  >  2:449281/1‑40 (MQ=255)
                                                        gcATCAGTGTACCAACCGCCGGACCCGCCCCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCgg                            <  1:1155864/90‑1 (MQ=255)
                                                             aGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc                       <  2:1182132/90‑1 (MQ=255)
                                                             aGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGtcatc                       <  2:715340/90‑1 (MQ=255)
                                                                     aaCCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTtgtg               >  1:1140349/1‑90 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGcc                                                           >  2:1141874/1‑43 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGcc                                                           <  1:1141874/43‑1 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc            <  2:192362/90‑1 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc            <  1:446312/90‑1 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc            >  2:1033681/1‑90 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc            >  1:494855/1‑90 (MQ=255)
                                                                        cgccgGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTc            >  1:549017/1‑90 (MQ=255)
                                                                                ccGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAAc    >  2:683418/1‑90 (MQ=255)
                                                                                  gCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                 >  1:962870/1‑43 (MQ=255)
                                                                                  gCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAt                                                 <  2:962870/43‑1 (MQ=255)
                                                                                  gCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGc  <  2:1168602/90‑1 (MQ=255)
                                                                                       |                                                                                    
TTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGC  >  NZ_CP009273/661805‑661976

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCTGGATTTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAA  >  NZ_CP009273/661797‑661983
                                                                                               |                                                                                           
GCTGGATTTTCCGCAGGCAGATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGT                                                                                         <  SRR3722089.1090946/100‑1 (MQ=60)
   GGATTTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCT                                                                                      >  SRR3722089.543484/1‑100 (MQ=60)
       TTTCCGCAGGCAGATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAG                                                                                  <  SRR3722089.329434/100‑1 (MQ=60)
         TCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAA                                                                                >  SRR3722089.877288/1‑100 (MQ=60)
          aagatggtcgtcggcagcgtcagatgtgtataagagacaGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAAT                                                                               <  SRR3722089.454671/61‑1 (MQ=60)
          CCGCCGGCAGATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAAT                                                                               <  SRR3722089.691085/100‑1 (MQ=60)
               GGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGG                                                                          <  SRR3722089.887037/100‑1 (MQ=60)
                 CAGATCGGTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCG                                                                        <  SRR3722089.772138/100‑1 (MQ=60)
                   GATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCtgtctctt                                                                      >  SRR3722089.1191810/1‑92 (MQ=60)
                   GATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGAC                                                                      >  SRR3722089.412585/1‑100 (MQ=60)
                        GTGTTGTTATCCCCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCA                                                                 <  SRR3722089.54893/100‑1 (MQ=60)
                        GTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCA                                                                 >  SRR3722089.835915/1‑100 (MQ=60)
                             GTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGC                                                            >  SRR3722089.1186383/1‑100 (MQ=60)
                             GTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGC                                                            >  SRR3722089.393279/1‑100 (MQ=60)
                             GTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGC                                                            >  SRR3722089.899053/1‑100 (MQ=60)
                                 atggtcgtcggcagcgtcagatgtgtataagagacagGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGAT                                                        <  SRR3722089.1159241/63‑1 (MQ=60)
                                   ACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTG                                                      >  SRR3722089.325217/1‑100 (MQ=60)
                                   ACCCAGCATAATGTGGTCGAGGATCTGGCACATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTG                                                      >  SRR3722089.1156588/1‑100 (MQ=60)
                                        GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACATCCACTTGCGTATTGTTGTA                                                 >  SRR3722089.148211/1‑100 (MQ=60)
                                        GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTC                                                 >  SRR3722089.1025249/1‑100 (MQ=60)
                                        GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                 >  SRR3722089.661733/1‑100 (MQ=60)
                                        GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                 >  SRR3722089.199244/1‑100 (MQ=60)
                                        GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                 >  SRR3722089.399350/1‑100 (MQ=60)
                                        GCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTA                                                 >  SRR3722089.1047135/1‑100 (MQ=60)
                                         CATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT                                                <  SRR3722089.43570/100‑1 (MQ=60)
                                         CATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT                                                <  SRR3722089.978429/100‑1 (MQ=60)
                                         CATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTAT                                                <  SRR3722089.127564/100‑1 (MQ=60)
                                               GTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCA                                          >  SRR3722089.45260/1‑100 (MQ=60)
                                                   TCGAGGATCTGGCGCATCAGTGTCCCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGG                                      <  SRR3722089.943355/100‑1 (MQ=60)
                                                   TCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGG                                      <  SRR3722089.1026874/100‑1 (MQ=60)
                                                                GCATCAGTGTACCAACCGCCGGACCCGCCCCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTT                         <  SRR3722089.1173455/100‑1 (MQ=60)
                                                                   TCAGTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTG                      >  SRR3722089.1157692/1‑100 (MQ=60)
                                                                      GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC                   >  SRR3722089.555395/1‑100 (MQ=60)
                                                                      GTGTACCAACCGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTC                   >  SRR3722089.500720/1‑100 (MQ=60)
                                                                                CGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCtgtctcttatacacatctgacgctgccgacgaccat         >  SRR3722089.977066/1‑64 (MQ=60)
                                                                                CGCCGGACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGC         <  SRR3722089.451663/100‑1 (MQ=60)
                                                                                      ACCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGGCATCAGTTTGTGGTCACGTAAACGCTCGGCA   >  SRR3722089.1103502/1‑100 (MQ=60)
                                                                                       CCCGCACCCCCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAA  >  SRR3722089.639826/1‑100 (MQ=60)
                                                                                               |                                                                                           
GCTGGATTTTCCGCAGGCAGATCGGTGTTGTTATCACCCAGCATAATGTGGTCGAGGATCTGGCGCATCAGTGTACCAACCGCCGGACCCGCACCACCGTTCTCCAGAATCATGGCGACAGCCACTTGCGGATTGTTGTATGGCGCAAAGGCGGTCATCAGTTTGTGGTCACGTAAACGCTCGGCAA  >  NZ_CP009273/661797‑661983

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: