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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,057,664 T→C G95G (GGA→GGG cobT ← nicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,057,6640TC100.0% 167.7 / NA 49G95G (GGA→GGGcobTnicotinate‑nucleotide‑‑dimethylbenzimidazole phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/23);  total (26/23)

GGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCAC  >  NZ_CP009273/2057577‑2057749
                                                                                       |                                                                                     
ggggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTccc                                                                                     >  1:411300/1‑90 (MQ=255)
ggggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTccc                                                                                     >  1:160289/1‑90 (MQ=255)
  ggATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCAc                                                                                   <  2:896414/90‑1 (MQ=255)
      aGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGt                                                                               >  1:972151/1‑90 (MQ=255)
      aGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGt                                                                               >  1:1128992/1‑90 (MQ=255)
         cTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCat                                                                            <  1:369526/90‑1 (MQ=255)
           cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                          >  2:1115914/1‑90 (MQ=255)
           cagcagTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCatat                                                                          >  1:71640/1‑90 (MQ=255)
               agTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTc                                                                      <  2:1177033/90‑1 (MQ=255)
               agTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTc                                                                      <  2:1045885/90‑1 (MQ=255)
               agTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTc                                                                      >  1:54581/1‑90 (MQ=255)
                gTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCa                                                                     >  1:536820/1‑90 (MQ=255)
                gTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCa                                                                     >  1:769190/1‑90 (MQ=255)
                 tATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                    <  1:401541/90‑1 (MQ=255)
                 tATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAg                                                                    <  1:58047/90‑1 (MQ=255)
                  aTCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGc                                                                   >  1:1215979/1‑90 (MQ=255)
                  aTCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGc                                                                   >  1:846961/1‑90 (MQ=255)
                   tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                  <  2:1108845/90‑1 (MQ=255)
                   tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                  <  1:1062007/90‑1 (MQ=255)
                   tCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGcc                                                                  <  2:507067/90‑1 (MQ=255)
                       tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                              <  1:1036549/90‑1 (MQ=255)
                       tACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTa                                                              >  1:295487/1‑90 (MQ=255)
                          cAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATgg                                                           >  2:1194937/1‑90 (MQ=255)
                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                         >  1:831536/1‑90 (MQ=255)
                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                         >  2:686864/1‑90 (MQ=255)
                            aCATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCt                                                         >  1:645500/1‑90 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                              >  1:1043505/1‑90 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                              >  2:449014/1‑90 (MQ=255)
                                       gTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCttt                                              >  2:433732/1‑90 (MQ=255)
                                        tGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                             <  2:769190/90‑1 (MQ=255)
                                        tGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCtttt                                             <  2:54581/90‑1 (MQ=255)
                                          gACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTgg                                           >  1:930601/1‑90 (MQ=255)
                                          gACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTgg                                           >  1:1007817/1‑90 (MQ=255)
                                                       gCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                              >  2:1004306/1‑90 (MQ=255)
                                                       gCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                              <  2:824673/90‑1 (MQ=255)
                                                       gCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGAcc                              <  1:68160/90‑1 (MQ=255)
                                                           gTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCtt                          <  2:1043505/90‑1 (MQ=255)
                                                                gCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTAt                                                             >  1:873410/1‑50 (MQ=255)
                                                                gCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTAt                                                             <  2:873410/50‑1 (MQ=255)
                                                                         cacacGCCGGTTGTCCCCCGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCg            <  2:823107/90‑1 (MQ=255)
                                                                          acacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGt           >  1:1138021/1‑90 (MQ=255)
                                                                          acacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGt           >  1:818054/1‑90 (MQ=255)
                                                                          acacGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGt           >  2:1103614/1‑90 (MQ=255)
                                                                             cGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGAt        >  2:1156826/1‑90 (MQ=255)
                                                                             cGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGAt        >  2:480359/1‑90 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  2:273396/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  2:1168030/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  1:232904/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  2:579577/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  2:972151/90‑1 (MQ=255)
                                                                                 ggTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGc    <  2:995600/90‑1 (MQ=255)
                                                                                   ttgtcCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGcac  <  1:433732/90‑1 (MQ=255)
                                                                                       |                                                                                     
GGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCAC  >  NZ_CP009273/2057577‑2057749

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTA  >  NZ_CP009273/2057567‑2057761
                                                                                                 |                                                                                                 
TGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCC                                                                                                 >  SRR3722089.416273/1‑100 (MQ=60)
TGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCC                                                                                                 >  SRR3722089.162284/1‑100 (MQ=60)
    accgagatctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGG                                                                                             <  SRR3722089.1163360/45‑1 (MQ=60)
     AGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGG                                                                                            <  SRR3722089.125754/100‑1 (MQ=60)
      GCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGT                                                                                           >  SRR3722089.986512/1‑100 (MQ=60)
      GCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGT                                                                                           >  SRR3722089.1146148/1‑100 (MQ=60)
           GGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATAT                                                                                      >  SRR3722089.72512/1‑100 (MQ=60)
               TAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTC                                                                                  >  SRR3722089.55257/1‑100 (MQ=60)
                AGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCA                                                                                 >  SRR3722089.779559/1‑100 (MQ=60)
                AGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCA                                                                                 >  SRR3722089.543074/1‑100 (MQ=60)
                  GCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGC                                                                               >  SRR3722089.1234479/1‑100 (MQ=60)
                  GCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGC                                                                               >  SRR3722089.859008/1‑100 (MQ=60)
                   CTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCC                                                                              <  SRR3722089.374067/100‑1 (MQ=60)
                       GCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTA                                                                          >  SRR3722089.299260/1‑100 (MQ=60)
                           agatctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGC                                                                      <  SRR3722089.482184/49‑1 (MQ=60)
                           TATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGC                                                                      <  SRR3722089.406418/100‑1 (MQ=60)
                           TATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGC                                                                      <  SRR3722089.58760/100‑1 (MQ=60)
                            ATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCT                                                                     >  SRR3722089.653753/1‑100 (MQ=60)
                            ATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCT                                                                     >  SRR3722089.843285/1‑100 (MQ=60)
                             TCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTG                                                                    <  SRR3722089.1077944/100‑1 (MQ=60)
                                 TACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCAC                                                                <  SRR3722089.1052055/100‑1 (MQ=60)
                                       CATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTT                                                          >  SRR3722089.1059131/1‑100 (MQ=60)
                                          CAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGG                                                       >  SRR3722089.944193/1‑100 (MQ=60)
                                          CAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGG                                                       >  SRR3722089.1022819/1‑100 (MQ=60)
                                                                CGCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTctgtctcttatacacatctgacgccgccga                                 >  SRR3722089.885953/1‑70 (MQ=60)
                                                                 GCTTGTTCTGCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCC                                <  SRR3722089.68982/100‑1 (MQ=60)
                                                                          GCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGT                       >  SRR3722089.829531/1‑100 (MQ=60)
                                                                          GCCAGCACACACACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGT                       >  SRR3722089.1155336/1‑100 (MQ=60)
                                                                                     CACGCCGGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACAC            >  SRR3722089.90174/1‑100 (MQ=60)
                                                                                           GGTTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACC      <  SRR3722089.235862/100‑1 (MQ=60)
                                                                                             TTGTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAG    <  SRR3722089.438951/100‑1 (MQ=60)
                                                                                               GTCCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTA  >  SRR3722089.1013213/1‑100 (MQ=60)
                                                                                                 |                                                                                                 
TGATAAGCCCGGGGATAGGCTCAGCAGTATCAATACCAACATCAATTACGTGGACGTTAGCGCCCGCTTGTTCTGCCAGCACACACACGCCGGTTGTTCCACGGGTCATATTTTCAGCCTGTATGGCTGTCACTTCTTTTGGAGAAATAGCGACCCCTTCCTCCCAGACGCCGTGATCGGCACACATAACCAGTA  >  NZ_CP009273/2057567‑2057761

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