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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,785,687 G→C G75G (GGC→GGG kbl ← glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,785,6870GC96.5% 184.0 / ‑1.1 58G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (34/22);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

ATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785600‑3785772
                                                                                       |                                                                                     
aTCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATccc                                                                                     <  1:237216/90‑1 (MQ=255)
    tCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGaaa                                                                                 <  2:271528/90‑1 (MQ=255)
      cATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAcc                                                                               <  2:312031/90‑1 (MQ=255)
            cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  1:104790/1‑90 (MQ=255)
            cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  2:351774/1‑90 (MQ=255)
             aGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                          <  2:1514676/88‑1 (MQ=255)
             aGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                          >  1:1514676/1‑88 (MQ=255)
              ggAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAAt                                                                       >  1:610263/1‑90 (MQ=255)
               gAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      <  1:887667/90‑1 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:503734/1‑90 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  2:1413590/1‑90 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  2:1559861/1‑90 (MQ=255)
                   gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:1259543/1‑90 (MQ=255)
                    ccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTc                                                                 >  2:582333/1‑90 (MQ=255)
                    ccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTc                                                                 >  2:104111/1‑90 (MQ=255)
                      gccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTccc                                                               >  1:300620/1‑90 (MQ=255)
                              gtGTTAAAGCTCTTTATGGCGGCCCTGGGTCCCGCAAAAAAACCGCCCCGAACCCACCCCGAACCCGTGAGAACCCATTCCCGCCTTTgc                                                       <  1:940904/89‑1 (MQ=255)
                               ttttCAAAGTCTTTATGGCTGTCCTGAGTGGCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgcc                                                      <  2:1519358/88‑1 (MQ=255)
                                   cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  1:1413590/90‑1 (MQ=255)
                                       ctctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAat                                              <  2:22435/90‑1 (MQ=255)
                                        tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:542739/90‑1 (MQ=255)
                                        tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:522062/90‑1 (MQ=255)
                                        tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:331784/90‑1 (MQ=255)
                                         ctTTATGGCTGTCCTGAGGGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcc                                            >  2:940904/1‑89 (MQ=255)
                                           ttATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcaga                                          <  1:578263/90‑1 (MQ=255)
                                             aTGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                        <  1:900241/90‑1 (MQ=255)
                                             aTGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                        <  1:72816/90‑1 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:693726/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGGGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:771640/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:1443894/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:1344437/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:8975/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:183336/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1225138/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:125159/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1457007/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1470409/1‑90 (MQ=255)
                                                   gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:201369/1‑90 (MQ=255)
                                                   gTCCTGAGCGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:100873/1‑90 (MQ=255)
                                                     ccTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatgat                                >  1:975922/1‑90 (MQ=255)
                                                          gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           <  2:300620/90‑1 (MQ=255)
                                                          gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           >  1:466212/1‑90 (MQ=255)
                                                          gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           >  2:641283/1‑90 (MQ=255)
                                                             ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  2:1147856/1‑90 (MQ=255)
                                                             ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:883425/1‑90 (MQ=255)
                                                             ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:395842/1‑90 (MQ=255)
                                                               gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:784740/1‑90 (MQ=255)
                                                               gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  1:486104/90‑1 (MQ=255)
                                                                  aaTAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCgag                   >  2:182418/1‑90 (MQ=255)
                                                                    tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:404673/1‑85 (MQ=255)
                                                                    tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:404673/85‑1 (MQ=255)
                                                                    tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                 <  2:975922/90‑1 (MQ=255)
                                                                      aaccGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAg               <  2:1099546/87‑1 (MQ=255)
                                                                       aaCGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgt              >  2:386060/1‑90 (MQ=255)
                                                                           cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          <  2:784740/90‑1 (MQ=255)
                                                                           cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          >  1:1369669/1‑90 (MQ=255)
                                                                           cACCGAAGCCATCCCGAAACCGTGAAATTCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          <  1:1487604/90‑1 (MQ=255)
                                                                                  gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  2:435266/90‑1 (MQ=255)
                                                                                   ccATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaaa  <  2:1259543/90‑1 (MQ=255)
                                                                                       |                                                                                     
ATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785600‑3785772

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ 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
ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  NZ_CP009273/3785595‑3785785
                                                                                            |                                                                                                  
ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                             <  SRR3722078.1201463/100‑1 (MQ=60)
ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                             <  SRR3722078.1006150/100‑1 (MQ=60)
     ATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGA                                                                                        <  SRR3722078.240109/100‑1 (MQ=60)
      gatctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacagGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAG                                                                                       <  SRR3722078.176951/50‑1 (MQ=60)
       CTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGA                                                                                      >  SRR3722078.106033/1‑100 (MQ=60)
        TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAA                                                                                     >  SRR3722078.1539381/1‑100 (MQ=60)
         TCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAAT                                                                                    >  SRR3722078.617307/1‑100 (MQ=60)
              CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                               >  SRR3722078.1279508/1‑100 (MQ=60)
              CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                               >  SRR3722078.509688/1‑100 (MQ=60)
                 CAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCC                                                                            >  SRR3722078.304299/1‑100 (MQ=60)
                    GAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCC                                                                         <  SRR3722078.899595/100‑1 (MQ=60)
                                   GTGTTAAAGCTCTTTATGGCGGCCCTGGGTCCCGCAAAAAAACCGCCCCGAACCCACCCCGAACCCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATC                                                          <  SRR3722078.953978/100‑1 (MQ=22)
                                        agacagCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCACCCCAAACCCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                     <  SRR3722078.1203433/94‑1 (MQ=60)
                                        CAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                     <  SRR3722078.1436538/100‑1 (MQ=60)
                                             TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                <  SRR3722078.335841/100‑1 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGGGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722078.781599/1‑100 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722078.203793/1‑100 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722078.1244430/1‑100 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722078.1494391/1‑100 (MQ=60)
                                              CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                               >  SRR3722078.1480763/1‑100 (MQ=60)
                                                TTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGAT                                             >  SRR3722078.989943/1‑100 (MQ=60)
                                                TTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGAT                                             <  SRR3722078.585003/100‑1 (MQ=60)
                                                  ATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTC                                           <  SRR3722078.73681/100‑1 (MQ=60)
                                                  ATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTC                                           <  SRR3722078.912325/100‑1 (MQ=60)
                                                     GCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCC                                        >  SRR3722078.471772/1‑100 (MQ=60)
                                                        GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                     >  SRR3722078.400647/1‑100 (MQ=60)
                                                        GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                     >  SRR3722078.895299/1‑100 (MQ=60)
                                                          CCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCC                                   >  SRR3722078.794969/1‑100 (MQ=60)
                                                               gcccccCAAATAAAACGCACCGACCCCATGCCAAACCCGTGAAAATCCATCCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAT                              <  SRR3722078.975491/94‑1 (MQ=60)
                                                               GTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAT                              >  SRR3722078.409580/1‑100 (MQ=60)
                                                                    caaaaaAAAACGCACCGAAGCCATCCCAAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                         <  SRR3722078.1472338/94‑1 (MQ=60)
                                                                    GCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                         <  SRR3722078.491874/100‑1 (MQ=60)
                                                                      AAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTT                       >  SRR3722078.1391719/1‑100 (MQ=60)
                                                                                CACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATctgtctcttatacacatctgacgctgccgacgatg             >  SRR3722078.1383757/1‑65 (MQ=60)
                                                                                CACCGAAGCCATCCCGAAACCGTGAAATTCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG             <  SRR3722078.1511849/100‑1 (MQ=60)
                                                                                 ACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGT            >  SRR3722078.709714/1‑100 (MQ=60)
                                                                                           TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  SRR3722078.184129/1‑100 (MQ=60)
                                                                                           TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  SRR3722078.418178/1‑100 (MQ=60)
                                                                                           TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  SRR3722078.748688/1‑100 (MQ=60)
                                                                                            |                                                                                                  
ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG  >  NZ_CP009273/3785595‑3785785

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