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 F3 I96 R1 172 34.7 2005794 96.2% 1929573 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,527,151 T→C Q87Q (CAA→CAG nanS ← 9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,527,1510TC96.7% 196.0 / ‑1.1 60Q87Q (CAA→CAGnanS9‑O‑acetyl‑N‑acetylneuraminic acid deacetylase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (28/30);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527072‑4527236
                                                                               |                                                                                     
gAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCg                                                                             <  1:548316/90‑1 (MQ=255)
      cTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCa                                                                       >  2:193971/1‑90 (MQ=255)
        gCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCaga                                                                     >  2:155938/1‑90 (MQ=255)
           ttATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  <  2:34355/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATg                                                                <  1:193971/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                >  2:94394/1‑90 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  2:920005/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  2:699903/90‑1 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                >  2:62702/1‑90 (MQ=255)
             aTCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATg                                                                <  1:210414/90‑1 (MQ=255)
              tCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATgg                                                               >  2:406409/1‑90 (MQ=255)
               cAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGt                                                              >  1:46655/1‑90 (MQ=255)
                  gAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAt                                                           >  1:383458/1‑90 (MQ=255)
                  gAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAt                                                           <  2:352135/90‑1 (MQ=255)
                      gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  1:433686/90‑1 (MQ=255)
                      gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       >  1:296877/1‑90 (MQ=255)
                      gAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  2:383458/90‑1 (MQ=255)
                       aaGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAccc                                                      >  1:424718/1‑90 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:179257/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:27758/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:46655/90‑1 (MQ=255)
                           gCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  2:664979/90‑1 (MQ=255)
                             aGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  2:604468/90‑1 (MQ=255)
                             aGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  1:754535/90‑1 (MQ=255)
                             aGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  1:505019/90‑1 (MQ=255)
                                 aTTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCt                                            <  1:444099/90‑1 (MQ=255)
                                   ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      <  1:617121/62‑1 (MQ=255)
                                   ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCag                                                                      >  2:617121/1‑62 (MQ=255)
                                   ttCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGa                                          >  1:607922/1‑90 (MQ=255)
                                      cGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       <  1:406409/90‑1 (MQ=255)
                                        tGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATc                                     >  1:453368/1‑90 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       <  1:255420/84‑1 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACa                                       >  2:255420/1‑84 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCgtgt                                 <  1:8312/90‑1 (MQ=255)
                                            atatGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCgtgt                                 >  2:206312/1‑90 (MQ=255)
                                                   gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          >  1:783411/1‑90 (MQ=255)
                                                   gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          >  1:244826/1‑90 (MQ=255)
                                                   gTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGt                          <  1:301123/90‑1 (MQ=255)
                                                    tGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCNAGAGGATGGTGATAACNCTTNATANCCTTAACATCGTGTGGGCAGTg                         >  2:224631/1‑90 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  >  2:700424/1‑45 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGa                                                                  <  1:700424/45‑1 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       >  1:430781/1‑56 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAcc                                                       <  2:430781/56‑1 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     >  1:283441/1‑90 (MQ=255)
                                                        tGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTc                     <  1:323105/90‑1 (MQ=255)
                                                                 gTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  <  1:665345/52‑1 (MQ=255)
                                                                 gTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCa                                                  >  2:665345/1‑52 (MQ=255)
                                                                   gCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaataat          >  2:468991/1‑90 (MQ=255)
                                                                   gCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGaataat          >  2:756420/1‑90 (MQ=255)
                                                                     cGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGt        <  2:424718/90‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  2:854943/44‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                >  2:864935/1‑44 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                >  1:854943/1‑44 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCata                                                <  1:864935/44‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:902688/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:117313/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  2:847137/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  <  1:468991/90‑1 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  >  1:778736/1‑90 (MQ=255)
                                                                           gtgtCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTaa  <  2:117313/90‑1 (MQ=255)
                                                                               |                                                                                     
GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA  >  NZ_CP009273/4527072‑4527236

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  >  NZ_CP009273/4527057‑4527246
                                                                                              |                                                                                               
AACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATG                                                                                            <  SRR3722236.619517/100‑1 (MQ=60)
               GAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGG                                                                             <  SRR3722236.556768/100‑1 (MQ=60)
                    CCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGT                                                                        >  SRR3722236.47330/1‑100 (MQ=60)
                      gatgtgtataagagacagGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGA                                                                      <  SRR3722236.627059/82‑1 (MQ=60)
                       GCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGAT                                                                     >  SRR3722236.388839/1‑100 (MQ=60)
                          cttcgtcggcagcgtcagatgtgtataagagacagATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAAC                                                                  <  SRR3722236.712255/65‑1 (MQ=60)
                           TATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACC                                                                 >  SRR3722236.301174/1‑100 (MQ=60)
                            ATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCC                                                                <  SRR3722236.196845/100‑1 (MQ=60)
                            ATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCC                                                                <  SRR3722236.213572/100‑1 (MQ=60)
                            ATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCC                                                                >  SRR3722236.430646/1‑100 (MQ=60)
                                 ctacacctaagccttcgtcggcagcgtcagatgtgtataagagacaGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCAT                                                           <  SRR3722236.967386/54‑1 (MQ=60)
                                 ctacacctaagccttcgtcggcagcgtcagatgtgtataagagacaGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCAT                                                           <  SRR3722236.726005/54‑1 (MQ=60)
                                     GAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCC                                                       <  SRR3722236.439770/100‑1 (MQ=60)
                                        GGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGA                                                    >  SRR3722236.617631/1‑100 (MQ=60)
                                          ggcagcgtcagatgtgtataagagacagCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAAC                                                  <  SRR3722236.676398/72‑1 (MQ=60)
                                            ccttcgtcggcagcgtcagatgtgtataagagacaGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACAT                                                <  SRR3722236.880108/65‑1 (MQ=60)
                                            AGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACAT                                                <  SRR3722236.512667/100‑1 (MQ=60)
                                            AGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACAT                                                <  SRR3722236.767478/100‑1 (MQ=60)
                                             GTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATC                                               >  SRR3722236.459934/1‑100 (MQ=60)
                                                ATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTG                                            <  SRR3722236.450466/100‑1 (MQ=60)
                                                     CGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGC                                       <  SRR3722236.412093/100‑1 (MQ=60)
                                                     CGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGC                                       <  SRR3722236.259186/100‑1 (MQ=60)
                                                        GCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGT                                    >  SRR3722236.796939/1‑100 (MQ=60)
                                                        GCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGT                                    >  SRR3722236.248470/1‑100 (MQ=60)
                                                           ATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAG                                 <  SRR3722236.8423/100‑1 (MQ=60)
                                                             ATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCctgtctcttatacacatctgacgc                               >  SRR3722236.436804/1‑76 (MQ=60)
                                                             ATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTC                               >  SRR3722236.287541/1‑100 (MQ=60)
                                                                  GTGCCTGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGG                          <  SRR3722236.305477/100‑1 (MQ=60)
                                                                       TGGCCAACGGTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAA                     <  SRR3722236.327719/100‑1 (MQ=60)
                                                                                GTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCtgtctcttatacacatctgacgctgccgacgaagg            >  SRR3722236.869916/1‑65 (MQ=60)
                                                                                GTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA            >  SRR3722236.792160/1‑100 (MQ=60)
                                                                                GTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA            >  SRR3722236.918609/1‑100 (MQ=60)
                                                                                GTGCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAA            >  SRR3722236.118955/1‑100 (MQ=60)
                                                                                  GCCGTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAG          >  SRR3722236.957316/1‑100 (MQ=60)
                                                                                     GTACTGTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCctgtctcttatacacatctgacgctgccgacgaaggcttaggtgtagatctcggtgg       >  SRR3722236.499733/1‑43 (MQ=60)
                                                                                          GTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  >  SRR3722236.547097/1‑100 (MQ=60)
                                                                                          GTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  <  SRR3722236.475921/100‑1 (MQ=60)
                                                                                          GTGTCTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGGGGAATAATGTCGTTAAAGTGACATGG  >  SRR3722236.456493/1‑100 (MQ=60)
                                                                                              |                                                                                               
AACACGGAACGATAAGAATCCCTGCATTATCAGGAATGAAGGGCAGTAATTTCCGTGCAATATGCAGTGCCTGGCCAACGGTGCCGTACTGTGTTTGATGATTCGTTGCCAGAGGATGGTGATAACCCTGCATATCCTGAACATCGTGTGGGCAGTGAGTCAGTGGAATAATGTCGTTAAAGTGACATGG  >  NZ_CP009273/4527057‑4527246

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