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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 832,557 T→C S251P (TCC→CCC)  hcxB → hydroxycarboxylate dehydrogenase HcXB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273832,5570TC97.8% 148.2 / NA 46S251P (TCC→CCC) hcxBhydroxycarboxylate dehydrogenase HcXB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/22);  minor base A (1/0);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGACGTTAACGGCGTGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGA  >  NZ_CP009273/832425‑832681
                                                                                                                                    |                                                                                                                            
aTTGACGTTAACGGCGTGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCcgg                                                                                                                          <  2:281774/137‑1 (MQ=255)
                tGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCAt                                                                                                          >  2:754932/1‑137 (MQ=255)
                tGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCa                                                                                                        >  2:382557/1‑139 (MQ=255)
                tGCCGACGACCAATCCGGCGGGAATGCAGGAGGCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTCCGCCCTTGCAGCGATGCGTGTAATTTTTTGCGGGGCGCTTACCGGCGGGACAACGACGcct                                                                                                          >  2:328773/1‑135 (MQ=255)
                     acgacCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   <  1:382557/139‑1 (MQ=255)
                          cAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTa                                                                                              <  2:843890/139‑1 (MQ=255)
                                 gcggTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAAc                                                                                         >  1:145932/1‑137 (MQ=255)
                                          cAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATg                                                                              <  1:620119/139‑1 (MQ=255)
                                               gTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCAtt                                                                         >  2:102566/1‑139 (MQ=255)
                                                  gCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   <  1:109066/110‑1 (MQ=255)
                                                  gCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   >  2:109066/1‑110 (MQ=255)
                                                  gCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCt                                                                       >  1:756563/1‑138 (MQ=255)
                                                          gTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              <  1:328773/139‑1 (MQ=255)
                                                          gTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              <  2:613044/139‑1 (MQ=255)
                                                             cGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAACCCAGTCCCGATGCCATTCTTCACTGCATgcc                                                           >  2:836052/1‑137 (MQ=255)
                                                              gCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGAc                                                           >  1:5324/1‑138 (MQ=255)
                                                              gCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGAc                                                           <  2:5324/138‑1 (MQ=255)
                                                                     aCCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatca                                                   >  1:68435/1‑139 (MQ=255)
                                                                          tGCTGATCATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatcat                                               <  2:243308/138‑1 (MQ=255)
                                                                          tGCCGTACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCCTCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatcat                                               >  1:137068/1‑138 (MQ=255)
                                                                          tGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacgac                                                                                                              >  1:542545/1‑75 (MQ=255)
                                                                          tGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacgac                                                                                                              <  2:542545/75‑1 (MQ=255)
                                                                          tGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTatcatcat                                               >  1:822552/1‑138 (MQ=255)
                                                                                           aCGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCgg                             >  1:215032/1‑139 (MQ=255)
                                                                                            cGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCa                                                         <  2:422221/110‑1 (MQ=255)
                                                                                            cGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCa                                                         >  1:422221/1‑110 (MQ=255)
                                                                                             gCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTgcggcg                           <  2:137068/139‑1 (MQ=255)
                                                                                               ccTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAACCCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGgcgcg                         <  1:102566/139‑1 (MQ=255)
                                                                                                  tGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCa                                                               <  1:468587/98‑1 (MQ=255)
                                                                                                  tGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCa                                                               >  2:468587/1‑98 (MQ=255)
                                                                                                         aTGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGt                >  1:23081/1‑138 (MQ=255)
                                                                                                         aTGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGt                >  1:533328/1‑138 (MQ=255)
                                                                                                         aTGTGTGAAATTCTTGGCGGGGCGCCTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGt                >  2:124308/1‑138 (MQ=255)
                                                                                                            tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacga                                                                                                               >  1:29267/1‑40 (MQ=255)
                                                                                                            tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAacga                                                                                                               <  2:29267/40‑1 (MQ=255)
                                                                                                            tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              >  1:465322/1‑89 (MQ=255)
                                                                                                            tgtgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCAt                                                              <  2:465322/89‑1 (MQ=255)
                                                                                                              tgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg          <  2:68435/139‑1 (MQ=255)
                                                                                                              tgAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAAcgcg          >  2:833300/1‑139 (MQ=255)
                                                                                                                aaaTTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACgcgc         <  1:833300/138‑1 (MQ=255)
                                                                                                                        ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   <  2:482826/40‑1 (MQ=255)
                                                                                                                        ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGaaa                                                                                                   >  1:482826/1‑40 (MQ=255)
                                                                                                                        ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGa                                                                                <  2:4779/59‑1 (MQ=255)
                                                                                                                        ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGa                                                                                >  1:4779/1‑59 (MQ=255)
                                                                                                                        ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCg   <  1:434211/136‑1 (MQ=255)
                                                                                                                        ggCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGa  <  1:754932/137‑1 (MQ=255)
                                                                                                                                    |                                                                                                                            
ATTGACGTTAACGGCGTGCCGACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGA  >  NZ_CP009273/832425‑832681

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCG  >  NZ_CP009273/832445‑832692
                                                                                                                |                                                                                                                                       
agatgtgtataagagacagATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACC                                                                                                   <  SRR3722120.110514/132‑1 (MQ=60)
GACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACC                                                                                                   <  SRR3722120.388179/151‑1 (MQ=60)
  CGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACC                                                                                                   >  SRR3722120.147902/1‑149 (MQ=60)
                   ATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTT                                                                                 >  SRR3722120.766781/1‑150 (MQ=60)
                     GCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAAC                                                                              <  SRR3722120.628744/151‑1 (MQ=60)
                               CCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACC                                                                     >  SRR3722120.5387/1‑150 (MQ=60)
                                     cgtgggctcggagatgtgtataagagacagGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATC                                                               <  SRR3722120.475382/120‑1 (MQ=60)
                                     GGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCA                                                              <  SRR3722120.333493/151‑1 (MQ=60)
                                      GTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCAT                                                             >  SRR3722120.69334/1‑151 (MQ=60)
                                           CTGTTGACCTTTGCCGTACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCCTCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATC                                                         >  SRR3722120.138921/1‑150 (MQ=60)
                                           CTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgcc                                                        >  SRR3722120.550198/1‑98 (MQ=60)
                                           CTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATC                                                         >  SRR3722120.833617/1‑150 (MQ=60)
                                                            ACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGC                                       >  SRR3722120.218007/1‑151 (MQ=60)
                                                             CATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCctgtctcttatacacatct                                      >  SRR3722120.428400/1‑132 (MQ=60)
                                                                          CCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAACCCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAA                         <  SRR3722120.103937/151‑1 (MQ=60)
                                                                          CCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTA                          >  SRR3722120.23342/1‑150 (MQ=60)
                                                                          CCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTA                          >  SRR3722120.540895/1‑150 (MQ=60)
                                                                             TTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAActgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcctctgcttgaaaaaaaaaaaagcacgatcggg                      >  SRR3722120.29597/1‑62 (MQ=60)
                                                                             TTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCActgtctcttatacacatctccgagcccacgagactaaggc                      >  SRR3722120.472074/1‑111 (MQ=60)
                                                                              TGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAtgcggcgaggctttcacccactcggcaaaggct                       <  SRR3722120.24804/149‑34 (MQ=60)
                                                                                         GTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCctgtctcttatacacatctccgagcccacgagactaaggcgaatatcgtatgccgtcttctgcttgaaaaaaaaaaaacagaacgtgg           >  SRR3722120.489779/1‑62 (MQ=60)
                                                                                         GTGAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTActgtctcttatacacatctccgagcccacgggactaaggcgaaaatcgtatgccgtcttctgcttgaaa           >  SRR3722120.4835/1‑81 (MQ=60)
                                                                                           GAAATTCTTGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCC         <  SRR3722120.844491/150‑1 (MQ=60)
                                                                                                   TGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCG  <  SRR3722120.765135/149‑1 (MQ=60)
                                                                                                   TGGCGGGGCGCTTCCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCC   <  SRR3722120.440571/148‑1 (MQ=60)
                                                                                                                |                                                                                                                                       
GACGACCAATCCGGCGGTAATGCAGGAGTCGCCGTTGGGTTCGCTGTTGACCTTTGCCGAACATAAAGGCTACGCCCTTGCAGCGATGTGTGAAATTCTTGGCGGGGCGCTTTCCGGCGGTAAAACGACGCATCAGGAAACGTTACAAACCAGTCCCGATGCCATTCTTAACTGCATGACCACTATCATCATCAACCCGGAACTCTTCGGCGCGCCGGATTGTAACGCGCAGACCGAAGCCTTTGCCG  >  NZ_CP009273/832445‑832692

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