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,002,903 T→C S261P (TCA→CCA)  ygeY → YgeY family selenium metabolism‑linked hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,002,9030TC92.7% 190.0 / 5.9 68S261P (TCA→CCA) ygeYYgeY family selenium metabolism‑linked hydrolase
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (30/33);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

GAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAAC  >  NZ_CP009273/3002816‑3002989
                                                                                       |                                                                                      
gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa                                                                                      >  2:212234/1‑90 (MQ=255)
gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa                                                                                      >  1:264588/1‑90 (MQ=255)
gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa                                                                                      >  1:1265666/1‑90 (MQ=255)
 aaTTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCaa                                                                                     >  2:97743/1‑90 (MQ=255)
 aaTTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCaa                                                                                     >  1:1268865/1‑90 (MQ=255)
        cGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACcg                                                                              <  2:697830/90‑1 (MQ=255)
        cGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACcg                                                                              <  1:541778/90‑1 (MQ=255)
        cGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcg                                                                              <  2:295379/90‑1 (MQ=255)
         ggCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcgt                                                                             <  2:59317/90‑1 (MQ=255)
             aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                         <  2:390943/90‑1 (MQ=255)
             aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                         <  2:1159967/90‑1 (MQ=255)
             aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                         <  1:126165/90‑1 (MQ=255)
             aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                         <  1:440912/90‑1 (MQ=255)
                  aCCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGAcc                                                                    >  1:1435309/1‑90 (MQ=255)
                    ccTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTg                                                                  <  2:1214834/90‑1 (MQ=255)
                    ccTCACCGTTTCTGAAATCTTCTTCACACCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTg                                                                  >  2:1001104/1‑90 (MQ=255)
                       cACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTgggg                                                               >  2:881365/1‑90 (MQ=255)
                          cGTTTCTGAAATCTTCTTCACTTCCCCAAGCCGTTGCGCTGCAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                            <  1:1350391/90‑1 (MQ=255)
                          cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                            <  1:212234/90‑1 (MQ=255)
                          cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                            <  1:1318857/90‑1 (MQ=255)
                            tttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGaaa                                                          <  2:1331399/90‑1 (MQ=255)
                            tttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGaaa                                                          <  1:595090/90‑1 (MQ=255)
                             ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAAc                                                         >  1:559011/1‑90 (MQ=255)
                             ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAc                                                         >  1:480611/1‑90 (MQ=255)
                             ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAc                                                         >  2:1275997/1‑90 (MQ=255)
                             ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAc                                                         >  2:184341/1‑90 (MQ=255)
                             ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGACGAAAc                                                         >  1:1526765/1‑90 (MQ=255)
                                           tcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcgcgc                                           >  1:989735/1‑90 (MQ=255)
                                           tcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcgcgc                                           >  1:975098/1‑90 (MQ=255)
                                            cacaTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCt                                          <  2:159438/90‑1 (MQ=255)
                                            cacaTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCt                                          <  2:264588/90‑1 (MQ=255)
                                                tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGAc                                      >  1:847192/1‑90 (MQ=255)
                                                   ccAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGaa                                   <  1:913059/90‑1 (MQ=255)
                                                        ccGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcg                              >  1:853051/1‑90 (MQ=255)
                                                          gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                            >  2:757533/1‑90 (MQ=255)
                                                          gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                            >  1:1483002/1‑90 (MQ=255)
                                                          gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                            >  2:1203887/1‑90 (MQ=255)
                                                           ttGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCgcgc                           >  1:1452184/1‑90 (MQ=255)
                                                           ttGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCgcgc                           >  1:629347/1‑90 (MQ=255)
                                                            tGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGcc                          <  1:869123/90‑1 (MQ=255)
                                                            tGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGcc                          <  1:206737/90‑1 (MQ=255)
                                                            tGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGcc                          <  1:574991/90‑1 (MQ=255)
                                                             gcgcTGTAGCAGACAGCTGCGCAGTTTCAAGCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGccc                         <  1:1046468/90‑1 (MQ=255)
                                                             gcgcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGccc                         <  2:393938/90‑1 (MQ=255)
                                                               gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                       <  2:853051/90‑1 (MQ=255)
                                                               gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                       <  2:1570242/90‑1 (MQ=255)
                                                                 tGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcc                     <  1:893438/90‑1 (MQ=255)
                                                                 tGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcc                     <  1:1004685/90‑1 (MQ=255)
                                                                  gTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcct                    <  2:680590/90‑1 (MQ=255)
                                                                    aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc                  <  1:1203887/90‑1 (MQ=255)
                                                                    aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc                  <  2:173774/90‑1 (MQ=255)
                                                                    aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc                  <  2:1270660/90‑1 (MQ=255)
                                                                    aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc                  <  2:477724/90‑1 (MQ=255)
                                                                      cAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAg                <  2:1175837/90‑1 (MQ=255)
                                                                        gACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTa              >  1:713240/1‑90 (MQ=255)
                                                                          cAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACa            >  2:1338884/1‑90 (MQ=255)
                                                                          cAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACa            >  1:1321741/1‑90 (MQ=255)
                                                                           aGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAg           <  2:576068/90‑1 (MQ=255)
                                                                            gCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          >  1:789871/1‑90 (MQ=255)
                                                                            gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          >  1:1318398/1‑90 (MQ=255)
                                                                            gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          <  2:1265666/90‑1 (MQ=255)
                                                                            gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          >  1:374494/1‑90 (MQ=255)
                                                                            gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          >  1:569860/1‑90 (MQ=255)
                                                                            gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          >  2:1213476/1‑90 (MQ=255)
                                                                            gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          >  1:726088/1‑90 (MQ=255)
                                                                            gCTGCGCAGTTCCAAACGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa          <  1:757163/90‑1 (MQ=255)
                                                                               gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAg       >  2:584261/1‑90 (MQ=255)
                                                                                    gTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAAc  <  2:1067407/90‑1 (MQ=255)
                                                                                       |                                                                                      
GAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAAC  >  NZ_CP009273/3002816‑3002989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ 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
GGGTTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTT  >  NZ_CP009273/3002806‑3002999
                                                                                                 |                                                                                                
GGGTTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCA                                                                                                >  SRR3722078.267821/1‑100 (MQ=60)
GGGTTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCA                                                                                                >  SRR3722078.680745/1‑100 (MQ=60)
GGGTTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCA                                                                                                >  SRR3722078.1285751/1‑100 (MQ=60)
 GGTTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAA                                                                                               >  SRR3722078.1288996/1‑100 (MQ=60)
    TATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCG                                                                                            <  SRR3722078.1118927/100‑1 (MQ=60)
      TGACGAATTCCTCGGCAAAGGCCCCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGAC                                                                                          <  SRR3722078.1387325/100‑1 (MQ=60)
             gaccaccgagatctacacgaccgccatcgtcggcagcgtcagatgtgtataagagacaGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTC                                                                                   <  SRR3722078.703101/42‑1 (MQ=60)
                  CGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACC                                                                              <  SRR3722078.548178/100‑1 (MQ=60)
                  CGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACC                                                                              >  SRR3722078.1458659/1‑100 (MQ=60)
                   ccaccgatatctacacgaccgccctcgtcggcagcgtcagatgtgtataagagacagGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCT                                                                             <  SRR3722078.580963/43‑1 (MQ=60)
                       AAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGG                                                                         <  SRR3722078.446201/100‑1 (MQ=60)
                       AAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGG                                                                         <  SRR3722078.127626/100‑1 (MQ=60)
                             CCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAAC                                                                   >  SRR3722078.565590/1‑100 (MQ=60)
                             CCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAC                                                                   >  SRR3722078.486308/1‑100 (MQ=60)
                             CCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGACGAAAC                                                                   >  SRR3722078.1551641/1‑100 (MQ=60)
                                    CGTTTCTGAAATCTTCTTCACTTCCCCAAGCCGTTGCGCTGCAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTCGGAA                                                            <  SRR3722078.1372073/100‑1 (MQ=60)
                                    CGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAA                                                            <  SRR3722078.1339941/100‑1 (MQ=60)
                                    CGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAA                                                            <  SRR3722078.214802/100‑1 (MQ=60)
                                      TTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGG                                                          <  SRR3722078.601990/100‑1 (MQ=60)
                                           GAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGC                                                     >  SRR3722078.989101/1‑100 (MQ=60)
                                           GAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGC                                                     >  SRR3722078.1004102/1‑100 (MQ=60)
                                                CTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGAC                                                >  SRR3722078.858502/1‑100 (MQ=60)
                                                        CATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCG                                        >  SRR3722078.864443/1‑100 (MQ=60)
                                                          TCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCG                                      >  SRR3722078.1507186/1‑100 (MQ=60)
                                                           CCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGC                                     >  SRR3722078.1475863/1‑100 (MQ=60)
                                                           CCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGC                                     >  SRR3722078.636589/1‑100 (MQ=60)
                                                             CCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCC                                   <  SRR3722078.925419/100‑1 (MQ=60)
                                                                      TGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAG                          <  SRR3722078.209232/100‑1 (MQ=60)
                                                                      TGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAG                          <  SRR3722078.880762/100‑1 (MQ=60)
                                                                      TGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAG                          <  SRR3722078.581702/100‑1 (MQ=60)
                                                                       GCGCTGTAGCAGACAGCTGCGCAGTTTCAAGCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGT                         <  SRR3722078.1062120/100‑1 (MQ=60)
                                                                        CGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTA                        >  SRR3722078.721987/1‑100 (MQ=60)
                                                                          CTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACA                      >  SRR3722078.1342863/1‑100 (MQ=60)
                                                                           TGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAG                     <  SRR3722078.1019399/100‑1 (MQ=60)
                                                                           TGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAG                     <  SRR3722078.905445/100‑1 (MQ=60)
                                                                            GTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGA                    >  SRR3722078.800202/1‑100 (MQ=60)
                                                                            GTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGA                    >  SRR3722078.1339475/1‑100 (MQ=60)
                                                                            GTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGA                    >  SRR3722078.379070/1‑100 (MQ=60)
                                                                            GTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGA                    >  SRR3722078.735137/1‑100 (MQ=60)
                                                                            GTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGA                    >  SRR3722078.576534/1‑100 (MQ=60)
                                                                              AGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAA                  <  SRR3722078.1222766/100‑1 (MQ=60)
                                                                                      GCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAActgtctcttatacacatctgacgctgccgacgatggcggtcgtgtagatc          >  SRR3722078.465569/1‑50 (MQ=60)
                                                                                      GCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGC          >  SRR3722078.539326/1‑100 (MQ=60)
                                                                                      GCTGCGCAGTTCCAAACGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGC          <  SRR3722078.766860/100‑1 (MQ=60)
                                                                                              GTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTT  >  SRR3722078.791840/1‑100 (MQ=60)
                                                                                              GTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTT  >  SRR3722078.825478/1‑100 (MQ=60)
                                                                                              GTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTT  >  SRR3722078.1056820/1‑100 (MQ=60)
                                                                                              GTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTT  >  SRR3722078.622659/1‑100 (MQ=60)
                                                                                                 |                                                                                                
GGGTTATGACGAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCTAACGCGGTTGTTT  >  NZ_CP009273/3002806‑3002999

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