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 I176 R1 180 42.9 1710034 96.5% 1650182 119.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,616,195 A→G E156G (GAA→GGA)  radA → DNA repair protein RadA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,616,1950AG100.0% 148.7 / NA 40E156G (GAA→GGA) radADNA repair protein RadA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/22);  total (18/22)

CTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACG  >  NZ_CP009273/4616062‑4616328
                                                                                                                                     |                                                                                                                                     
ctACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAAc                                                                                                                                    >  1:313499/1‑137 (MQ=255)
 tACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAg                                                                                                                                 >  2:270194/1‑139 (MQ=255)
   cTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCa                                                                                                                               <  2:417225/139‑1 (MQ=255)
   cTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCa                                                                                                                               <  2:169642/139‑1 (MQ=255)
       aaaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGa                                                                                                                           <  2:313499/139‑1 (MQ=255)
       aaaCGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGa                                                                                                                           <  2:560892/139‑1 (MQ=255)
               tgCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCt                                                                                                                   >  2:704000/1‑139 (MQ=255)
                     cTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGa                                                                                                             >  2:820682/1‑139 (MQ=255)
                        gCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGAt                                                                                                            <  2:361774/137‑1 (MQ=255)
                                  tGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAagaaga                                                                                                <  2:606632/139‑1 (MQ=255)
                                      nnngCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAagaaga                                                                                                <  2:64907/132‑1 (MQ=255)
                                                       gCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAAt                                                                           >  2:448707/1‑139 (MQ=255)
                                                          aaGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGa                                                                        <  2:359788/139‑1 (MQ=255)
                                                            gagTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACt                                                                      <  1:492084/139‑1 (MQ=255)
                                                            gagTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACt                                                                      <  1:429341/139‑1 (MQ=255)
                                                                  cTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATc                                                                 >  2:637569/1‑138 (MQ=255)
                                                                         aGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGAt                                                         >  2:181899/1‑139 (MQ=255)
                                                                          ggTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATg                                                        <  1:448707/139‑1 (MQ=255)
                                                                               cAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCatat                                                   <  1:780649/139‑1 (MQ=255)
                                                                                  tGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGc                                                <  2:578061/139‑1 (MQ=255)
                                                                                  tGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGc                                                <  1:40527/139‑1 (MQ=255)
                                                                                           aTCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAc                                        <  1:181899/138‑1 (MQ=255)
                                                                                               ccTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCa                                   <  2:11696/139‑1 (MQ=255)
                                                                                                   ggCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACt                                                                      >  1:814243/1‑100 (MQ=255)
                                                                                                   ggCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACt                                                                      <  2:814243/100‑1 (MQ=255)
                                                                                                   ggCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGc                               >  2:261353/1‑139 (MQ=255)
                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGa                                                                                                                           >  2:142879/1‑41 (MQ=255)
                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGa                                                                                                                           <  1:142879/41‑1 (MQ=255)
                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGAt                                                         >  2:360941/1‑107 (MQ=255)
                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGAt                                                         <  1:360941/107‑1 (MQ=255)
                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTgg                           <  2:736364/137‑1 (MQ=255)
                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTgg                           >  2:592393/1‑137 (MQ=255)
                                                                                                         ccGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTgg                           >  1:736364/1‑137 (MQ=255)
                                                                                                            aCTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCg                      <  1:94314/139‑1 (MQ=255)
                                                                                                                      tcAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGt            >  1:734245/1‑139 (MQ=255)
                                                                                                                          tatGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGcc                              >  2:822786/1‑117 (MQ=255)
                                                                                                                          tatGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGcc                              <  1:822786/117‑1 (MQ=255)
                                                                                                                              ttgtCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGaaa    >  1:788330/1‑139 (MQ=255)
                                                                                                                               tgtCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAAc   >  1:536911/1‑139 (MQ=255)
                                                                                                                                gtCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAAcg  >  1:178426/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACG  >  NZ_CP009273/4616062‑4616328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ 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
GAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTG  >  NZ_CP009273/4616049‑4616340
                                                                                                                                                  |                                                                                                                                                 
GAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACC                                                                                                                                               <  SRR3722130.596895/151‑1 (MQ=60)
  AATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACC                                                                                                                                               >  SRR3722130.318383/1‑149 (MQ=60)
                   ggattgtgcgttttttttttttcaagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCC                                                                                                                            <  SRR3722130.144993/63‑1 (MQ=60)
                                                                        agcagaagacggcacacgatattcgccttagtctcgtgggctcggagatgtgtataagagacagATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTTATTGACTCGATCCAGGTG                                                                       <  SRR3722130.205190/87‑1 (MQ=60)
                                                                        AGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTG                                                                       <  SRR3722130.499701/151‑1 (MQ=60)
                                                                        AGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTG                                                                       <  SRR3722130.436117/151‑1 (MQ=60)
                                                                            cagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGC                                                                   <  SRR3722130.709512/89‑1 (MQ=60)
                                                                                   gcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAATCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGC                                                            <  SRR3722130.766988/98‑1 (MQ=60)
                                                                                     cggagatgtgtagaagagacagGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGG                                                          <  SRR3722130.366581/129‑1 (MQ=60)
                                                                                      AGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGA                                                         <  SRR3722130.455729/151‑1 (MQ=60)
                                                                                           GCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTAC                                                    <  SRR3722130.792046/151‑1 (MQ=60)
                                                                                              ATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGT                                                 <  SRR3722130.41088/151‑1 (MQ=60)
                                                                                                     CTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGctgtctcttatacacatctccgaacccac                                          >  SRR3722130.826105/1‑122 (MQ=60)
                                                                                                       CATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCT                                         <  SRR3722130.184630/150‑1 (MQ=60)
                                                                                                           GCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGT                                      >  SRR3722130.747173/1‑149 (MQ=60)
                                                                                                                 ataagagacagGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGG                               <  SRR3722130.834745/139‑1 (MQ=60)
                                                                                                                        GACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTG                       >  SRR3722130.745021/1‑151 (MQ=60)
                                                                                                                        GACTGACAATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTG                       <  SRR3722130.95629/151‑1 (MQ=60)
                                                                                                                                ATCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAAC               >  SRR3722130.799823/1‑151 (MQ=60)
                                                                                                                                 TCTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACG              >  SRR3722130.545087/1‑151 (MQ=60)
                                                                                                                                  CTCAATATGTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGG             >  SRR3722130.181112/1‑151 (MQ=60)
                                                                                                                                          GTTGTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTctgtctcttatacacatctccgagcccacg        >  SRR3722130.102382/1‑118 (MQ=60)
                                                                                                                                             GTCGGGAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGtctcttatacacatctccgagcccacgagactaaggcg  >  SRR3722130.241205/1‑113 (MQ=60)
                                                                                                                                                  |                                                                                                                                                 
GAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGAAGAGTCGCTGCAACAGGTGGCAATGCGCGCTCATCGCCTTGGCCTGCCGACTGACAATCTCAATATGTTGTCGGAAACCAGCATCGAACAGATCTGCCTGATTGCCGAAGAAGAGCAACCGAAGCTGATGGTAATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTG  >  NZ_CP009273/4616049‑4616340

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