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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,084,653 A→G G358G (GGT→GGC frvB ← PTS fructose‑like transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,084,6530AG97.8% 125.6 / ‑3.5 46G358G (GGT→GGCfrvBPTS fructose‑like transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (21/24);  minor base C (1/0);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGA  >  NZ_CP009273/4084575‑4084738
                                                                              |                                                                                     
ggCCCCGACCACAATGGCACTATTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGAc                                                                            <  1:846849/90‑1 (MQ=255)
ggCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGAc                                                                            <  1:416765/90‑1 (MQ=255)
      gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                      >  1:926968/1‑90 (MQ=255)
      gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGggcgg                                                                      >  2:1005420/1‑90 (MQ=255)
      gACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGCCCCGCCGCGACGggcgg                                                                      >  2:1514051/1‑90 (MQ=255)
          acaATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATc                                                                  <  1:531762/90‑1 (MQ=255)
          acaATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATc                                                                  <  1:701457/90‑1 (MQ=255)
             aTGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAAGTCGCAAGGCCCGCCGCGACGGGCGGCATCAgcg                                                              >  1:855380/1‑90 (MQ=255)
                gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                            >  1:77924/1‑90 (MQ=255)
                gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                            >  1:1052603/1‑90 (MQ=255)
                gCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGcgcg                                                            >  2:1250871/1‑90 (MQ=255)
                     actGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAAc                                                       >  1:538089/1‑90 (MQ=255)
                        ggCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGcc                                                    >  2:1255634/1‑90 (MQ=255)
                            gCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcgacg                                                <  2:302930/90‑1 (MQ=255)
                             ccTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCcgacga                                               >  1:657429/1‑90 (MQ=255)
                              ctctTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGAt                                              <  1:152908/90‑1 (MQ=255)
                              ctctTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGAt                                              <  2:235654/90‑1 (MQ=255)
                                ctTTTTCCTGTGGGGTAAACAGTTTGGGCCCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATgg                                            <  1:1268004/90‑1 (MQ=255)
                                 tttttCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGc                                           >  1:1234072/1‑90 (MQ=255)
                                     tCCGGTCGGTTAAACAGTTTGGGCGCAACGACCGCCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTAc                                       <  1:1221284/90‑1 (MQ=255)
                                     cccTGTGGGGTAAACAGTTTGGGCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTAc                                       <  1:927243/89‑1 (MQ=255)
                                       ccgccGGGTAACAAGTTAGGCCGCAATGATCGGCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTCCCa                                     <  1:810563/86‑1 (MQ=255)
                                       cTGTCGGGTAAACAGTTTGGGGGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCa                                     <  1:1437279/90‑1 (MQ=255)
                                        tGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAg                                    >  2:1069378/1‑90 (MQ=255)
                                        tGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAg                                    >  1:1306823/1‑90 (MQ=255)
                                         gTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGt                                   >  2:320742/1‑90 (MQ=255)
                                          tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGAc                                                          <  2:1349382/66‑1 (MQ=255)
                                          tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGAc                                                          >  1:1349382/1‑66 (MQ=255)
                                          tCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTc                                  >  1:1220537/1‑90 (MQ=255)
                                                     gTTTGGGCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGtt                       <  2:885950/90‑1 (MQ=255)
                                                     gTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGtt                       >  1:103474/1‑90 (MQ=255)
                                                      tttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGttt                      >  1:554166/1‑90 (MQ=255)
                                                      tttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGttt                      >  1:1052076/1‑90 (MQ=255)
                                                       ttGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCa                                     >  1:695000/1‑74 (MQ=255)
                                                       ttGGGCGCAAGGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCa                                     <  2:695000/74‑1 (MQ=255)
                                                          ggCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTCCCAGTCATAAATGTGTTTTTcc                  <  1:1224763/90‑1 (MQ=255)
                                                               aaTGAAGGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag             <  1:543247/90‑1 (MQ=255)
                                                               aaTGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag             >  1:116932/1‑90 (MQ=255)
                                                               aaTGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCcagcag             <  2:80821/90‑1 (MQ=255)
                                                                 tGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTg           <  2:538089/90‑1 (MQ=255)
                                                                     cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa       <  1:23827/90‑1 (MQ=255)
                                                                     cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa       <  1:1209637/90‑1 (MQ=255)
                                                                     cGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGaa       <  2:659180/90‑1 (MQ=255)
                                                                        cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGTGTACCAGTCATAAATGTGTTTTTCcagcag             <  1:1471550/81‑1 (MQ=255)
                                                                        cGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGTGTACCAGTCATAAATGTGTTTTTCcagcag             >  2:1471550/1‑81 (MQ=255)
                                                                          cAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCaga  <  2:116932/90‑1 (MQ=255)
                                                                              |                                                                                     
GGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGA  >  NZ_CP009273/4084575‑4084738

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTT  >  NZ_CP009273/4084554‑4084752
                                                                                                   |                                                                                                   
AGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGG                                                                                                     <  SRR3722078.753190/100‑1 (MQ=60)
AGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGG                                                                                                     <  SRR3722078.1551421/100‑1 (MQ=60)
 GCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGC                                                                                                    <  SRR3722078.21189/100‑1 (MQ=60)
                 CTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGG                                                                                    >  SRR3722078.939705/1‑100 (MQ=60)
                     GGCCCCGACCACAATGGCACTATTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATC                                                                                <  SRR3722078.858153/100‑1 (MQ=60)
                     GGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATC                                                                                <  SRR3722078.421804/100‑1 (MQ=60)
                        ACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAAGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCG                                                                            >  SRR3722078.866811/1‑100 (MQ=60)
                           GACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCG                                                                          >  SRR3722078.1068400/1‑100 (MQ=60)
                           GACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCG                                                                          >  SRR3722078.78844/1‑100 (MQ=60)
                               ACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAA                                                                      <  SRR3722078.709909/100‑1 (MQ=60)
                               ACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAA                                                                      <  SRR3722078.538046/100‑1 (MQ=60)
                                CAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAAC                                                                     >  SRR3722078.544445/1‑100 (MQ=60)
                                        CTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGA                                                             >  SRR3722078.664959/1‑100 (MQ=60)
                                            TGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGC                                                         >  SRR3722078.1253538/1‑100 (MQ=60)
                                                   CTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAG                                                  >  SRR3722078.1327673/1‑100 (MQ=60)
                                                   CTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAG                                                  <  SRR3722078.154682/100‑1 (MQ=60)
                                                     CTTTTTCCTGTGGGGTAAACAGTTTGGGCCCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTC                                                <  SRR3722078.1288119/100‑1 (MQ=60)
                                                     CTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGctgtctcttataca                                                >  SRR3722078.1371036/1‑86 (MQ=60)
                                                     CTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTC                                                >  SRR3722078.1239742/1‑100 (MQ=60)
                                                          TCCGGTCGGTTAAACAGTTTGGGCGCAACGACCGCCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAA                                           <  SRR3722078.1240502/100‑1 (MQ=60)
                                                          CCCTGTGGGGTAAACAGTTTGGGCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAA                                           <  SRR3722078.939988/100‑1 (MQ=60)
                                                            ccgccgggtaacaagttaggcCGCAATGATCGGCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTCCCAGTCATAAATG                                         <  SRR3722078.821281/79‑1 (MQ=60)
                                                            CTGTCGGGTAAACAGTTTGGGGGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATG                                         <  SRR3722078.1460672/100‑1 (MQ=60)
                                                                CGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTT                                     >  SRR3722078.104703/1‑100 (MQ=60)
                                                                CGGGTAAAAAGGTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGAGGCATCAGCGCGACAACGCCGACGATGGCGgccccgtaggacacaagga                                     >  SRR3722078.1079248/1‑81 (MQ=60)
                                                                 GGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTT                                    >  SRR3722078.1067858/1‑100 (MQ=60)
                                                                 GGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTT                                    >  SRR3722078.560702/1‑100 (MQ=60)
                                                                  GGTAAACAGTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGctgtct                                   >  SRR3722078.703293/1‑94 (MQ=60)
                                                                          GTTTGGGCGCAATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAG                           >  SRR3722078.118300/1‑100 (MQ=60)
                                                                               GGCGCAATGACCGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTCCCAGTCATAAATGTGTTTTTCCAGCAGTGAGA                      <  SRR3722078.1244046/100‑1 (MQ=60)
                                                                                    AATGAAGGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAG                 <  SRR3722078.549667/100‑1 (MQ=60)
                                                                                    AATGAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGTGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAG                 <  SRR3722078.1495559/100‑1 (MQ=60)
                                                                                       GAACGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAG              >  SRR3722078.1259139/1‑100 (MQ=60)
                                                                                          CGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAAC           <  SRR3722078.24122/100‑1 (MQ=60)
                                                                                          CGTCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAAC           <  SRR3722078.1228644/100‑1 (MQ=60)
                                                                                            TCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAActgtctcttatacacatctgacgctgccgacgat         >  SRR3722078.971445/1‑66 (MQ=60)
                                                                                            TCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATctgtctcttatacacatctgacgctgccgacga         >  SRR3722078.872279/1‑67 (MQ=60)
                                                                                            TCGCAAGGCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCA         >  SRR3722078.134870/1‑100 (MQ=60)
                                                                                                   GCCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTT  >  SRR3722078.907684/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
AGCCGGTTCGGTTGCGACTGTGGCACCGACCACAATGGCACTACTGGCAGCCTCTTTTTCCTGTCGGGTAAACAGTTTGGGCGCAATGAACGTCGCAAGACCCGCCGCGACGGGCGGCATCAGCGCGACAACGCCGACGATGGCGTACCAGTCATAAATGTGTTTTTCCAGCAGTGAGAAGCAGAAGAACCAGGCGGTT  >  NZ_CP009273/4084554‑4084752

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