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 728,467 T→C I1143I (ATT→ATC rhsC → rhs element protein RhsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273728,4670TC100.0% 142.9 / NA 42I1143I (ATT→ATCrhsCrhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/23);  total (19/23)

CTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTG  >  NZ_CP009273/728333‑728593
                                                                                                                                      |                                                                                                                               
cTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCAc                                                                                                                              <  2:144802/138‑1 (MQ=21)
cTCGACCGGCTGGAAAGTGAAACCCTGGCTGACCGGGTTAGGTAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCAc                                                                                                                              <  2:124041/138‑1 (MQ=11)
 tCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCt                                                                                                                            >  2:475493/1‑139 (MQ=21)
        gCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCAc                                                                                                                     <  2:694861/139‑1 (MQ=21)
                        cTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGcc                                                                                                     <  1:651837/139‑1 (MQ=21)
                                ccGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTg                                                                                             >  1:852225/1‑139 (MQ=21)
                                       agtgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCtt                                                                                       <  1:213418/138‑1 (MQ=37)
                                        gtgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTAt                                                                                     <  2:162136/139‑1 (MQ=37)
                                        gtgagGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTAt                                                                                     <  2:13722/139‑1 (MQ=37)
                                           agGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAg                                                                                  >  1:526766/1‑139 (MQ=37)
                                                 gccgccGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGa                                                                            <  1:594319/139‑1 (MQ=255)
                                                           gctggcATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGacaaca                                                                  <  2:336598/139‑1 (MQ=255)
                                                                cATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTg                                                             <  1:629892/139‑1 (MQ=255)
                                                                cATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTg                                                             <  2:528714/139‑1 (MQ=255)
                                                                     tGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACg                                                        <  2:204231/139‑1 (MQ=255)
                                                                     tGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACg                                                        <  1:475493/139‑1 (MQ=255)
                                                                     tGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACg                                                        <  1:21646/139‑1 (MQ=255)
                                                                     tGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACg                                                        >  1:85732/1‑139 (MQ=255)
                                                                                  tGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGa                                           <  2:208656/139‑1 (MQ=255)
                                                                                  tGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGa                                           <  2:85732/139‑1 (MQ=255)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCt                                                                                                    >  1:67031/1‑80 (MQ=255)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCt                                                                                                    <  2:67031/80‑1 (MQ=255)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGcc                                                                                          >  1:812111/1‑90 (MQ=37)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGcc                                                                                          <  2:812111/90‑1 (MQ=37)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTg                                                             >  1:822899/1‑119 (MQ=255)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTg                                                             <  2:822899/119‑1 (MQ=255)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAAt                                         >  1:189376/1‑139 (MQ=255)
                                                                                    gcgcAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAAt                                         >  2:603615/1‑139 (MQ=255)
                                                                                       cAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATggg                                      >  1:93530/1‑139 (MQ=255)
                                                                                           tGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTaa                                  <  1:568033/139‑1 (MQ=255)
                                                                                            gCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGCAAc                                 >  2:630430/1‑139 (MQ=255)
                                                                                                aGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACctgc                             <  1:603615/139‑1 (MQ=255)
                                                                                                          aCCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCtgaatgaa                    >  1:561429/1‑137 (MQ=255)
                                                                                                                tataCACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAAccc            >  1:449631/1‑139 (MQ=255)
                                                                                                                     acGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATc       <  2:346646/139‑1 (MQ=255)
                                                                                                                         cGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGt                                    >  2:40335/1‑107 (MQ=255)
                                                                                                                         cGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGt                                    <  1:40335/107‑1 (MQ=255)
                                                                                                                         cGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCt   >  1:607821/1‑139 (MQ=255)
                                                                                                                          ggCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAAATGAAGAGAACCCGCATCAGCt   >  1:64503/1‑139 (MQ=255)
                                                                                                                          ggCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTg  >  2:244629/1‑139 (MQ=255)
                                                                                                                          ggCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTg  >  2:204326/1‑139 (MQ=255)
                                                                                                                          ggCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTg  >  1:360819/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                               
CTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTG  >  NZ_CP009273/728333‑728593

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ 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
ATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGC  >  NZ_CP009273/728354‑728594
                                                                                                                 |                                                                                                                                
ATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGG                                                                                             >  SRR3722120.863642/1‑151 (MQ=21)
  CCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCC                                                                                           <  SRR3722120.660813/151‑1 (MQ=23)
           CCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGC                                                                                  >  SRR3722120.534260/1‑151 (MQ=31)
                 GAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGA                                                                             <  SRR3722120.216360/150‑1 (MQ=40)
                           AGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAG                                                                  <  SRR3722120.602608/151‑1 (MQ=42)
                                     GGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGC                                                        >  SRR3722120.86829/1‑151 (MQ=60)
                                          GCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAAT                                                   <  SRR3722120.638612/151‑1 (MQ=60)
                                               GTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGAT                                              <  SRR3722120.21888/151‑1 (MQ=60)
                                               GTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGAT                                              <  SRR3722120.482368/151‑1 (MQ=60)
                                                    GCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgt                                         >  SRR3722120.67911/1‑102 (MQ=29)
                                                    GCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACctgtctcttatacacatctccgagcccacgagactaagg                                         >  SRR3722120.823048/1‑112 (MQ=38)
                                                    GCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATctgtctctt                                          >  SRR3722120.833968/1‑141 (MQ=60)
                                                    GCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATG                                         >  SRR3722120.192004/1‑151 (MQ=60)
                                                       TGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGG                                      >  SRR3722120.94781/1‑151 (MQ=60)
                                                                   cggagatgtgtataagagacagGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAA                         <  SRR3722120.40826/129‑1 (MQ=53)
                                                                     ATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATG                       <  SRR3722120.575970/151‑1 (MQ=60)
                                                                          AAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAG                  <  SRR3722120.611997/151‑1 (MQ=59)
                                                                          AAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAG                    >  SRR3722120.569289/1‑149 (MQ=59)
                                                                                GATGGACCCGGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCG            >  SRR3722120.456194/1‑151 (MQ=52)
                                                                                         GGTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTG   >  SRR3722120.616251/1‑151 (MQ=50)
                                                                                          GTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAAATGAAGAGAACCCGCATCAGCTG   >  SRR3722120.65341/1‑151 (MQ=46)
                                                                                          GTATACACGCCGGCGCGAAAAATCCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGC  >  SRR3722120.366059/1‑151 (MQ=48)
                                                                                                                 |                                                                                                                                
ATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACGGTGGCGCAGATGCAAAGCCAGATGGACCCGGTATACACGCCGGCGCGAAAAATTCACCTGTACCACTGCGACCATCGCGGCCTGCCGCTGGCCCTTATCAGCACGGAAGGGACAACAGCGTGGTACGCAGAATACGATGAATGGGGTAACCTGCTGAATGAAGAGAACCCGCATCAGCTGC  >  NZ_CP009273/728354‑728594

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