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 I102 R1 278 39.1 2251254 95.8% 2156701 87.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 552,665 T→C E178G (GAA→GGA)  folD ← bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase FolD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273552,6650TC100.0% 154.9 / NA 42E178G (GAA→GGA) folDbifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase FolD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/22);  total (20/22)

TCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGC  >  NZ_CP009273/552579‑552754
                                                                                      |                                                                                         
tCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCAt                                                                                       <  1:789244/91‑1 (MQ=255)
      gATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCAt                                                                                 >  2:44218/1‑91 (MQ=255)
      gATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCAt                                                                                 >  1:996428/1‑91 (MQ=255)
        tCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCg                                                                               >  2:877276/1‑91 (MQ=255)
        tCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCg                                                                               >  1:1090486/1‑91 (MQ=255)
         cGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCgg                                                                              >  2:1033045/1‑91 (MQ=255)
         cGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCgg                                                                              <  1:361930/91‑1 (MQ=255)
         cGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCgg                                                                              >  1:595151/1‑91 (MQ=255)
           gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGggc                                                                            >  1:100861/1‑91 (MQ=255)
           gCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGggc                                                                            <  1:877276/91‑1 (MQ=255)
              ttttCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGggcggc                                                                         <  1:139066/91‑1 (MQ=255)
               tttCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGcc                                                                        <  1:848457/91‑1 (MQ=255)
               tttCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGcc                                                                        <  1:232309/91‑1 (MQ=255)
                       tgatgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGatat                                                                >  1:459730/1‑91 (MQ=255)
                                aGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATgcgc                                                       <  1:253138/91‑1 (MQ=255)
                                aGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATgcgc                                                       <  2:182545/91‑1 (MQ=255)
                                aGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATgcgc                                                       <  2:68159/91‑1 (MQ=255)
                                aGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATgcgc                                                       <  2:996428/91‑1 (MQ=255)
                                           gAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACg                                            <  1:1058481/91‑1 (MQ=255)
                                             aGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGc                                          <  1:871874/91‑1 (MQ=255)
                                                ggTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGtt                                       >  2:915409/1‑91 (MQ=255)
                                                 gTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTg                                      <  1:1033045/91‑1 (MQ=255)
                                                   tcgTCACTGTAGTGGTGCAACGTGCCAGCAGGACGCCCGTGCTCCTCGCGCGGCCAACGATATTCGATGCGCAAACCACCACGGCGTTGac                                    >  2:855559/3‑90 (MQ=255)
                                                     gTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGc                                  >  2:130927/1‑91 (MQ=255)
                                                      tCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGcc                                 <  1:130927/91‑1 (MQ=255)
                                                         cTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGaa                              <  2:339317/91‑1 (MQ=255)
                                                              gtggtgCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTAt                         >  2:1118246/1‑91 (MQ=255)
                                                               tggtgCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATc                        <  1:18738/91‑1 (MQ=255)
                                                               tggtgCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATc                        <  1:73158/91‑1 (MQ=255)
                                                               tggtgCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATc                        <  2:29904/91‑1 (MQ=255)
                                                                 gtgCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCaa                      >  2:831432/1‑91 (MQ=255)
                                                                  tgCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAAt                     <  1:285094/91‑1 (MQ=255)
                                                                   gCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAg                             <  2:162939/82‑1 (MQ=255)
                                                                   gCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAg                             >  1:162939/1‑82 (MQ=255)
                                                                   gCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAAtg                    >  1:668858/1‑91 (MQ=255)
                                                                   gCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAAtg                    >  2:898488/1‑91 (MQ=255)
                                                                   gCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAAtg                    >  1:672916/1‑91 (MQ=255)
                                                                     aaCCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAAtgtt                  >  2:225725/1‑91 (MQ=255)
                                                                     aaCCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAAtgtt                  >  1:690660/1‑91 (MQ=255)
                                                                     aaCCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAAtgtt                  >  1:868178/1‑91 (MQ=255)
                                                                        cTGCCAGCAGCATGCCCATGCTCATCGGGCGCCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTa               <  2:322676/91‑1 (MQ=255)
                                                                                 gcagTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTa                          <  2:1045269/71‑1 (MQ=255)
                                                                                 gcagTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTa                          >  1:1045269/1‑71 (MQ=255)
                                                                                  cagTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAagc     <  2:100861/91‑1 (MQ=255)
                                                                                     ttCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAagcagc  <  2:267250/91‑1 (MQ=255)
                                                                                      |                                                                                         
TCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGC  >  NZ_CP009273/552579‑552754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ 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
GCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACC  >  NZ_CP009273/552565‑552765
                                                                                                    |                                                                                                    
GCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTC                                                                                                      >  SRR3722125.69413/1‑101 (MQ=60)
 CCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCC                                                                                                     <  SRR3722125.415974/101‑1 (MQ=60)
      cagaagacggcatacgagatcatgcctagtctcgtgggctcggagatgtgtataagagacAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGC                                                                                                <  SRR3722125.500987/41‑1 (MQ=60)
      GGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGC                                                                                                <  SRR3722125.1112166/101‑1 (MQ=60)
      GGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGC                                                                                                <  SRR3722125.413757/101‑1 (MQ=60)
      GGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGC                                                                                                <  SRR3722125.277785/101‑1 (MQ=60)
          ACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCAT                                                                                            >  SRR3722125.1018674/1‑101 (MQ=60)
          ACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCAT                                                                                            <  SRR3722125.118188/101‑1 (MQ=60)
           CGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATC                                                                                           <  SRR3722125.475313/101‑1 (MQ=60)
            GATCAAAAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCG                                                                                          >  SRR3722125.1115017/1‑101 (MQ=60)
             ATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGG                                                                                         >  SRR3722125.607951/1‑101 (MQ=60)
              TCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGG                                                                                        <  SRR3722125.806463/101‑1 (MQ=60)
               CAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGC                                                                                       >  SRR3722125.102752/1‑101 (MQ=60)
                       CGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACG                                                                               <  SRR3722125.369515/101‑1 (MQ=60)
                         GCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGAT                                                                             <  SRR3722125.896684/101‑1 (MQ=60)
                           ATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATAT                                                                           >  SRR3722125.469430/1‑101 (MQ=60)
                            TTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATT                                                                          <  SRR3722125.141711/101‑1 (MQ=60)
                             TTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTC                                                                         <  SRR3722125.867129/101‑1 (MQ=60)
                             TTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTC                                                                         <  SRR3722125.236954/101‑1 (MQ=60)
                                              AGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCAC                                                        <  SRR3722125.258240/101‑1 (MQ=60)
                                                         GAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGC                                             <  SRR3722125.1082298/101‑1 (MQ=60)
                                                           AGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCG                                           <  SRR3722125.891140/101‑1 (MQ=60)
                                                               GTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGG                                       <  SRR3722125.1056219/101‑1 (MQ=60)
                                                                    TCACTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCA                                  <  SRR3722125.133398/101‑1 (MQ=60)
                                                                       CTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATG                               >  SRR3722125.166059/1‑101 (MQ=60)
                                                                       CTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATG                               >  SRR3722125.683255/1‑101 (MQ=60)
                                                                       CTGTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATG                               >  SRR3722125.687406/1‑101 (MQ=60)
                                                                         GTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTT                             >  SRR3722125.887350/1‑101 (MQ=60)
                                                                         GTAGTGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTT                             >  SRR3722125.705544/1‑101 (MQ=60)
                                                                             TGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAA                         <  SRR3722125.19079/101‑1 (MQ=60)
                                                                             TGGTGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAA                         <  SRR3722125.74518/101‑1 (MQ=60)
                                                                                TGCAACCTGCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGC                      <  SRR3722125.290887/101‑1 (MQ=60)
                                                                                     CCTGCCAGCAGCAGTCCCAAGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTctgtctctta                 >  SRR3722125.1068722/1‑91 (MQ=60)
                                                                                        GCCAGCAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACtgtctcttat              >  SRR3722125.918993/1‑91 (MQ=60)
                                                                                             CAGCAGTCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGA         >  SRR3722125.94009/1‑101 (MQ=60)
                                                                                                   TCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATAC   <  SRR3722125.1068837/101‑1 (MQ=60)
                                                                                                   TCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATAC   <  SRR3722125.511133/101‑1 (MQ=60)
                                                                                                   TCCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATAC   <  SRR3722125.230214/101‑1 (MQ=60)
                                                                                                    CCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACC  >  SRR3722125.849890/1‑101 (MQ=60)
                                                                                                    CCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACC  >  SRR3722125.183748/1‑101 (MQ=60)
                                                                                                    |                                                                                                    
GCCAACGGCAACGATCAATAGATCGGCATTTTCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACC  >  NZ_CP009273/552565‑552765

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: