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 3,785,687 G→C G75G (GGC→GGG kbl ← glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,785,6870GC100.0% 180.6 / NA 54G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (27/27);  total (27/27)

TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785603‑3785771
                                                                                    |                                                                                    
ttCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGaa                                                                                 >  2:480002/1‑90 (MQ=255)
         cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                        <  2:474918/90‑1 (MQ=255)
         cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                        >  1:471027/1‑90 (MQ=255)
          aGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGaa                                                                       >  1:706952/1‑90 (MQ=255)
              aGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCa                                                                   >  1:935784/1‑90 (MQ=255)
               ggCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAt                                                                  >  1:580915/1‑90 (MQ=255)
                gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                 >  1:441064/1‑90 (MQ=255)
                 ccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTc                                                                <  2:992878/90‑1 (MQ=255)
                 ccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTc                                                                >  2:1040509/1‑90 (MQ=255)
                   gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                              >  1:1215319/1‑90 (MQ=255)
                                cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                 <  2:649172/90‑1 (MQ=255)
                                cAAGCTCTTTATGGCTGTCCTGAGTCCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                 <  1:815044/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTTCCGCAAATAAAACGCACCGAAGCCATCCCGAACCCGTGGGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  1:194841/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  1:1040509/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:975007/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:852073/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:706952/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:441064/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:1215319/90‑1 (MQ=255)
                                     tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  1:194842/90‑1 (MQ=255)
                                          aTGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAatcagatc                                       <  1:728542/90‑1 (MQ=255)
                                            ggCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatcag                                     >  2:850908/1‑90 (MQ=255)
                                             gCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAgg                                    <  1:397857/90‑1 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgccgcc                                                  >  1:262784/1‑73 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgccgcc                                                  <  2:262784/73‑1 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            >  1:38/1‑79 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  1:1115961/79‑1 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:38/79‑1 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            >  2:1115961/1‑79 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:76989/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:813020/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:45602/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:287221/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:778603/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:932823/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:360492/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  2:1062820/1‑90 (MQ=255)
                                                gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                 >  1:981088/1‑90 (MQ=255)
                                                     gAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCg                            >  2:349339/1‑90 (MQ=255)
                                                       gtcccGCAAATAAAACGCGCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                          <  2:1034339/87‑1 (MQ=255)
                                                          ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                       >  1:543106/1‑90 (MQ=255)
                                                          ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                       >  2:724560/1‑90 (MQ=255)
                                                          ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                       >  2:547135/1‑90 (MQ=255)
                                                            gCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     <  1:104571/90‑1 (MQ=255)
                                                            gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     >  2:934137/1‑90 (MQ=255)
                                                            gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                     <  1:349339/90‑1 (MQ=255)
                                                                 tAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                <  2:574257/90‑1 (MQ=255)
                                                                        ccccGAAGCCACCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt         <  1:611821/88‑1 (MQ=255)
                                                                        cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                <  2:77147/83‑1 (MQ=255)
                                                                        cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                >  1:77147/1‑83 (MQ=255)
                                                                               gCCATCCCGAACCCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  2:580915/90‑1 (MQ=255)
                                                                               gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  1:891590/90‑1 (MQ=255)
                                                                               gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  2:883556/90‑1 (MQ=255)
                                                                               cccATCCCGAACCCGGGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  1:549436/89‑1 (MQ=255)
                                                                                    |                                                                                    
TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785603‑3785771

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  NZ_CP009273/3785588‑3785785
                                                                                                   |                                                                                                  
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                    >  SRR3722089.1171551/1‑100 (MQ=60)
       gagacagCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                             <  SRR3722089.978275/93‑1 (MQ=60)
       ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                             <  SRR3722089.406350/100‑1 (MQ=60)
              CTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGA                                                                                      >  SRR3722089.476656/1‑100 (MQ=60)
               TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAA                                                                                     >  SRR3722089.716159/1‑100 (MQ=60)
                   ATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCA                                                                                 >  SRR3722089.949481/1‑100 (MQ=60)
                    TCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAT                                                                                >  SRR3722089.587874/1‑100 (MQ=60)
                     CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                               >  SRR3722089.446358/1‑100 (MQ=60)
                        CAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCC                                                                            >  SRR3722089.1233807/1‑100 (MQ=60)
                                               CAAGCTCTTTATGGCTGTCCTGAGTCCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                     <  SRR3722089.826449/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTTCCGCAAATAAAACGCACCGAAGCCATCCCGAACCCGTGGGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                <  SRR3722089.197292/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                <  SRR3722089.1056089/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                <  SRR3722089.197293/100‑1 (MQ=60)
                                                    GCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAACCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                <  SRR3722089.1132865/100‑1 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCtgtctc                                               >  SRR3722089.266149/1‑94 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722089.824387/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722089.789196/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722089.77929/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722089.946462/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722089.46176/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722089.290898/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722089.995614/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATC                                               >  SRR3722089.38/1‑100 (MQ=60)
                                                         ATGGCTGTCCTGAGTGCCGCAAATAAAACGCACGGAAGCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTC                                           <  SRR3722089.738017/100‑1 (MQ=60)
                                                            GCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCC                                        <  SRR3722089.402698/100‑1 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                     >  SRR3722089.549432/1‑100 (MQ=60)
                                                                           GCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                         <  SRR3722089.105864/100‑1 (MQ=60)
                                                                           GCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAAAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                         <  SRR3722089.353676/100‑1 (MQ=60)
                                                                             AAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTT                       >  SRR3722089.78089/1‑100 (MQ=60)
                                                                                       CCCCGAAGCCACCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG             <  SRR3722089.619447/100‑1 (MQ=60)
                                                                                         CCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCctgtctcttatacacatctgacgctgccgacg           >  SRR3722089.709800/1‑68 (MQ=60)
                                                                                              GCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA      <  SRR3722089.904477/100‑1 (MQ=60)
                                                                                              CCCATCCCGAACCCGGGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA      <  SRR3722089.555819/100‑1 (MQ=60)
                                                                                               CCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGAC     >  SRR3722089.955832/1‑100 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  SRR3722089.481236/1‑100 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  SRR3722089.1073305/1‑100 (MQ=60)
                                                                                                   |                                                                                                  
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  NZ_CP009273/3785588‑3785785

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