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
A2 F5 I175 R1 328 64.7 2582366 96.2% 2484236 120.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,043,563 A→C G195G (GGT→GGG gfcC ← capsule biosynthesis GfcC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,043,5630AC100.0% 177.8 / NA 49G195G (GGT→GGGgfcCcapsule biosynthesis GfcC family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/27);  total (22/27)

CAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCA  >  NZ_CP009273/1043427‑1043693
                                                                                                                                        |                                                                                                                                  
cAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTcccc                                                                                                                                  <  2:1215684/139‑1 (MQ=255)
  gACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCt                                                                                                                                <  1:581152/139‑1 (MQ=255)
  gACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCt                                                                                                                                <  2:419843/139‑1 (MQ=255)
  gACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCt                                                                                                                                <  2:612650/139‑1 (MQ=255)
  gACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCt                                                                                                                                >  2:752537/1‑139 (MQ=255)
   acacAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCt                                                                                                                                >  1:572010/1‑138 (MQ=255)
    cacaATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCa                                                                                                                              <  1:714723/139‑1 (MQ=255)
             gTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTa                                                                                                                        <  2:69163/136‑1 (MQ=255)
                 tttAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCa                                                                                                                 <  2:347743/139‑1 (MQ=255)
                 tttAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCa                                                                                                                 >  2:1247306/1‑139 (MQ=255)
                 tttAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCa                                                                                                                 >  2:538702/1‑139 (MQ=255)
                      aTCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACg                                                                                                            <  2:515523/139‑1 (MQ=255)
                                         aGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTg                                                                                         <  2:572010/139‑1 (MQ=255)
                                                  gcAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCggg                                                                                >  2:463688/1‑139 (MQ=255)
                                                     gAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCgg                                                                                 >  1:850756/1‑135 (MQ=255)
                                                     gAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCgg                                                                                 <  2:850756/135‑1 (MQ=255)
                                                     gAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGt                                                                             >  1:1042306/1‑139 (MQ=255)
                                                     gAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGt                                                                             >  2:930629/1‑139 (MQ=255)
                                                      aaaaGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTg                                                                            <  2:811434/139‑1 (MQ=255)
                                                                              cAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCt                                                    <  1:201123/139‑1 (MQ=255)
                                                                                       gtTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGgcctg                                            <  2:283328/138‑1 (MQ=255)
                                                                                                  cTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAAt                                                              >  1:802169/1‑109 (MQ=255)
                                                                                                  cTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAAt                                                              <  2:802169/109‑1 (MQ=255)
                                                                                                   ttttATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATc                                                                         <  1:975998/97‑1 (MQ=255)
                                                                                                   ttttATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATc                                                                         >  2:975998/1‑97 (MQ=255)
                                                                                                        ttCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGcc                          <  1:1102125/139‑1 (MQ=255)
                                                                                                            acaGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGgcctg                                            >  2:583179/1‑117 (MQ=255)
                                                                                                            acaGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGgcctg                                            <  1:583179/117‑1 (MQ=255)
                                                                                                            acaGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGc                      >  1:598152/1‑139 (MQ=255)
                                                                                                             caGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTa                                                         >  1:298680/1‑103 (MQ=255)
                                                                                                             caGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTa                                                         <  2:298680/103‑1 (MQ=255)
                                                                                                             caGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGcc                     >  1:5479/1‑139 (MQ=255)
                                                                                                             caGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGcc                     >  2:1246851/1‑139 (MQ=255)
                                                                                                                ggCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGccc                    >  1:1275532/1‑137 (MQ=255)
                                                                                                                     cAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGt                                                          >  1:623995/1‑94 (MQ=255)
                                                                                                                     cAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGt                                                          <  2:623995/94‑1 (MQ=255)
                                                                                                                         ggCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCa         <  1:463688/139‑1 (MQ=255)
                                                                                                                          gCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAc        >  2:1052237/1‑139 (MQ=255)
                                                                                                                          gCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAc        >  2:251284/1‑139 (MQ=255)
                                                                                                                           cgacgaCAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTg                                                                                         <  2:490042/57‑1 (MQ=255)
                                                                                                                           cgacgaCAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTg                                                                                         >  1:490042/1‑57 (MQ=255)
                                                                                                                           cgacgaCAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAcc       <  1:251284/139‑1 (MQ=255)
                                                                                                                            gacgacAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCAccc      <  1:845395/139‑1 (MQ=255)
                                                                                                                              cgacAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCcag    <  2:773581/139‑1 (MQ=255)
                                                                                                                              cgacAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCcag    <  2:278437/139‑1 (MQ=255)
                                                                                                                                acAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGcc                          <  1:256666/115‑1 (MQ=255)
                                                                                                                                acAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGcc                          >  2:256666/1‑115 (MQ=255)
                                                                                                                                acAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCcagca  <  2:598152/139‑1 (MQ=255)
                                                                                                                                    tctcCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGaaa              <  2:348170/123‑1 (MQ=255)
                                                                                                                                    tctcCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGaaa              >  1:348170/1‑123 (MQ=255)
                                                                                                                                        |                                                                                                                                  
CAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCA  >  NZ_CP009273/1043427‑1043693

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAAT  >  NZ_CP009273/1043419‑1043700
                                                                                                                                                |                                                                                                                                         
CGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTC                                                                                                                                      >  SRR3722129.580473/1‑150 (MQ=60)
         AGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATC                                                                                                                            <  SRR3722129.589730/151‑1 (MQ=60)
           ACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCAC                                                                                                                          <  SRR3722129.725215/151‑1 (MQ=60)
              CAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCctgtctcttatacacat                                                                                                                       >  SRR3722129.766708/1‑134 (MQ=60)
                                                CAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtat                                                                                     >  SRR3722129.1028515/1‑100 (MQ=60)
                                                  GTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTG                                                                                   >  SRR3722129.1056947/1‑151 (MQ=60)
                                                  GTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTG                                                                                   >  SRR3722129.862916/1‑151 (MQ=60)
                                                             tttttcgtttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGA                                                                         <  SRR3722129.496860/72‑1 (MQ=60)
                                                                  cgtgggctcggagatgtgtataagagacaGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAAT                                                                     <  SRR3722129.989812/120‑1 (MQ=60)
                                                                                     CCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCT                                                <  SRR3722129.203935/151‑1 (MQ=60)
                                                                                              ataagagacaGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGC                                        <  SRR3722129.591788/140‑1 (MQ=60)
                                                                                               GTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACtgtctcttatacacatctc                                      >  SRR3722129.813749/1‑132 (MQ=60)
                                                                                                         TCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCC                            >  SRR3722129.607079/1‑151 (MQ=60)
                                                                                                          CTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCctgtctcttatacacatctccgagcc                           >  SRR3722129.302952/1‑125 (MQ=60)
                                                                                                          CTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCC                           >  SRR3722129.5549/1‑151 (MQ=60)
                                                                                                             TTATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCC                          >  SRR3722129.1293243/1‑149 (MQ=60)
                                                                                                               ATTCCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAA                      <  SRR3722129.1117491/151‑1 (MQ=60)
                                                                                                                gtataagagacagCAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAA                      <  SRR3722129.260289/137‑1 (MQ=60)
                                                                                                                  CCACAGGGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGtctcttatacacatctccgagcccacgagact                   >  SRR3722129.633269/1‑119 (MQ=60)
                                                                                                                        GGCAACAGGGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGctgtctcttatacacatctccgagcccacgagactaaggcgaatatcgtatgccgtcttctgcttgaaaaa              >  SRR3722129.497120/1‑79 (MQ=60)
                                                                                                                                GGGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGT     <  SRR3722129.470431/151‑1 (MQ=60)
                                                                                                                                 GGCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCctgtct    >  SRR3722129.353183/1‑145 (MQ=60)
                                                                                                                                  GCGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAA   <  SRR3722129.254825/151‑1 (MQ=60)
                                                                                                                                   CGACGACAGTCTCCCCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAAT  <  SRR3722129.857501/151‑1 (MQ=60)
                                                                                                                                                |                                                                                                                                         
CGTTAAAACAGACACAATCTGGTCATTTAAATCAGCATATTTTTCGGGCAGTACATGCGCAGAAAAGCCCAACCAAAGCTGGCTGCCAGGCGGTGGTTCCACATGTCTTTTATTCCACAGGGCAACAGGGGCGACGACAGTCTCACCCTCAGGCGTAATCACCATAACGTTATTCTTGTCCGCGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAAT  >  NZ_CP009273/1043419‑1043700

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