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,536,106 T→C E261G (GAA→GGA)  cysP ← thiosulfate/sulfate ABC transporter substrate‑binding protein CysP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,536,1060TC100.0% 232.1 / NA 68E261G (GAA→GGA) cysPthiosulfate/sulfate ABC transporter substrate‑binding protein CysP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (43/25);  total (43/25)

GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA  >  NZ_CP009273/2536021‑2536190
                                                                                     |                                                                                    
gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  <  1:1115791/90‑1 (MQ=255)
gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  <  2:84430/90‑1 (MQ=255)
gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  <  1:941805/90‑1 (MQ=255)
gTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGt                                                                                  <  1:790689/90‑1 (MQ=255)
  tATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAc                                                                                >  1:362152/1‑90 (MQ=255)
   aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAcc                                                                               >  2:261987/1‑90 (MQ=255)
   aTTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAcc                                                                               >  1:848161/1‑90 (MQ=255)
        cgcgGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTgg                                                                          >  2:131464/1‑90 (MQ=255)
             tagtaATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGc                                                                     >  2:79751/1‑90 (MQ=255)
                  aTAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGtt                                                                <  2:568672/90‑1 (MQ=255)
                     gTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGttttt                                                             >  1:102775/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:299364/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:928286/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:1028841/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:1053450/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  2:229551/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:1179446/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:419386/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:796500/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:48508/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  1:794776/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  2:1188165/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTa                                                            >  2:1124470/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGGTTTTa                                                            >  2:873780/1‑90 (MQ=255)
                      tCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGGACCGTTGGCCTGCACGTTTTTa                                                            >  1:806864/1‑90 (MQ=255)
                        ggTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATc                                                          >  2:880140/1‑90 (MQ=255)
                        ggTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATc                                                          >  2:955056/1‑90 (MQ=255)
                        ggTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATc                                                          <  2:51841/90‑1 (MQ=255)
                         gTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCa                                                         >  1:458715/1‑90 (MQ=255)
                         gTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCa                                                         >  1:41581/1‑90 (MQ=255)
                            atgatgGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAAcc                                                      >  1:975824/1‑90 (MQ=255)
                              gatgGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCa                                                    >  1:54566/1‑90 (MQ=255)
                               atgGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCAc                                                   <  2:287853/90‑1 (MQ=255)
                                 ggTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGc                                                 <  1:114604/90‑1 (MQ=255)
                                   tttGCGCCTGCGTGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCACCCCACGCCa                                               <  1:955056/90‑1 (MQ=255)
                                       cgcCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCgg                                           <  2:1179446/90‑1 (MQ=255)
                                           tGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAt                                       >  2:44799/1‑90 (MQ=255)
                                           tGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAt                                       >  1:340607/1‑90 (MQ=255)
                                            gCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAtt                                      >  2:117294/1‑90 (MQ=255)
                                            gCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAtt                                      <  1:724032/90‑1 (MQ=255)
                                            gCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAtt                                      <  1:390599/90‑1 (MQ=255)
                                                   ataGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAg                               <  2:848161/90‑1 (MQ=255)
                                                    taGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTCCCGTTGGACGGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGa                              >  1:138809/1‑90 (MQ=255)
                                                      gagCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAAt                            >  2:833541/1‑90 (MQ=255)
                                                        gCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGt                          >  2:836031/1‑90 (MQ=255)
                                                           aGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTgg                       >  1:180109/1‑90 (MQ=255)
                                                            gTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGt                      <  2:245679/90‑1 (MQ=255)
                                                             ttCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGtt                     >  2:1161608/1‑90 (MQ=255)
                                                             ttCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGtt                     >  1:316946/1‑90 (MQ=255)
                                                             ttCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGtt                     >  2:793293/1‑90 (MQ=255)
                                                                 gATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCa                                                                    >  2:195609/1‑39 (MQ=39)
                                                                 gATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCa                                                                    <  1:195609/39‑1 (MQ=39)
                                                                 gATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGc                                                 <  2:263757/58‑1 (MQ=255)
                                                                 gATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGc                                                 >  1:263757/1‑58 (MQ=255)
                                                                 gATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCg                 >  2:361408/1‑90 (MQ=255)
                                                                    aGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGaa              <  2:975824/90‑1 (MQ=255)
                                                                    aGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGaa              <  2:1078708/90‑1 (MQ=255)
                                                                      gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATc            >  2:731339/1‑90 (MQ=255)
                                                                      gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATc            >  2:1228825/1‑90 (MQ=255)
                                                                      gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATc            >  1:1139833/1‑90 (MQ=255)
                                                                      gCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCAGGAATTCCGCCAGAATGTTGGTTTTCGGAATc            >  1:935532/1‑90 (MQ=255)
                                                                         ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  2:706567/90‑1 (MQ=255)
                                                                         ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  2:796500/90‑1 (MQ=255)
                                                                         ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  2:1403/90‑1 (MQ=255)
                                                                         ttGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcacc         <  1:1010514/90‑1 (MQ=255)
                                                                           ggcggcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcaccac       <  2:41581/90‑1 (MQ=255)
                                                                           ggcggcTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATcaccac       <  1:793293/90‑1 (MQ=255)
                                                                                cTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCaa  <  2:180109/90‑1 (MQ=255)
                                                                                     |                                                                                    
GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAA  >  NZ_CP009273/2536021‑2536190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ 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
CCATCACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCC  >  NZ_CP009273/2536008‑2536195
                                                                                                  |                                                                                         
CCATCACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCC                                                                                          <  SRR3722089.853125/100‑1 (MQ=60)
CCATCACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCC                                                                                          <  SRR3722089.129311/100‑1 (MQ=60)
     cagatgtgtataagagacaGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAC                                                                                     <  SRR3722089.787509/81‑1 (MQ=60)
     ACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAC                                                                                     <  SRR3722089.963697/100‑1 (MQ=60)
     ACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTAC                                                                                     >  SRR3722089.366613/1‑100 (MQ=60)
      CCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACC                                                                                    >  SRR3722089.860225/1‑100 (MQ=60)
             GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCC                                                                             <  SRR3722089.955600/100‑1 (MQ=60)
             GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCC                                                                             <  SRR3722089.1132692/100‑1 (MQ=60)
             GTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCC                                                                             <  SRR3722089.801559/100‑1 (MQ=60)
                        GGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTT                                                                  >  SRR3722089.104052/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACnnnnnnnnnnnnnnnnnnnn                                                                 >  SRR3722089.759089/1‑80 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.424450/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.941844/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.807485/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.805728/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.49115/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.1044202/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.303172/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.1069257/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.1197415/1‑100 (MQ=60)
                         GTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGGACCGTTGGCCTGCACGTTTTTA                                                                 >  SRR3722089.818093/1‑100 (MQ=60)
                           taagatggtcgtcggcagcgtcagatgtgtataagagacAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATC                                                               <  SRR3722089.198072/61‑1 (MQ=60)
                            GTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCA                                                              >  SRR3722089.42101/1‑100 (MQ=60)
                            GTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCA                                                              >  SRR3722089.464197/1‑100 (MQ=60)
                               ATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACC                                                           >  SRR3722089.990246/1‑100 (MQ=60)
                                 AGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCA                                                         >  SRR3722089.55242/1‑100 (MQ=60)
                                              GGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAT                                            >  SRR3722089.344853/1‑100 (MQ=60)
                                              GGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAAT                                            <  SRR3722089.116006/100‑1 (MQ=60)
                                                TTTTCGCCTGCGTGCTATAGAGCCATTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCACCCCACGCCACCGGGAATTC                                          <  SRR3722089.969095/100‑1 (MQ=60)
                                                       CTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTCCCGTTGGACGGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGA                                   >  SRR3722089.140501/1‑100 (MQ=60)
                                                         GCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAAT                                 <  SRR3722089.395357/100‑1 (MQ=60)
                                                         GCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAAT                                 <  SRR3722089.733453/100‑1 (MQ=60)
                                                              CTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGG                            >  SRR3722089.182362/1‑100 (MQ=60)
                                                                ATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTT                          >  SRR3722089.320964/1‑100 (MQ=60)
                                                                    AGCCAGTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATctgtctcttatacacatctgac                      >  SRR3722089.267136/1‑78 (MQ=60)
                                                                         GTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATC                 >  SRR3722089.1157172/1‑100 (MQ=60)
                                                                         GTTCAGATAGGCTTTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCAGGAATTCCGCCAGAATGTTGGTTTTCGGAATC                 >  SRR3722089.949227/1‑100 (MQ=60)
                                                                                      TTGGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGC    <  SRR3722089.1025566/100‑1 (MQ=60)
                                                                                        GGCGGCTTTTCCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCC  <  SRR3722089.804220/100‑1 (MQ=60)
                                                                                                  |                                                                                         
CCATCACCTCCGGGTTATTCACGCGGTAGTAATAGTCGGTGATGATGGTTTGCGCCTGCGGGCTATAGAGCCAGTTCAGATAGGCTTTGGCGGCTTTTTCCGTACCGTTGGCCTGCACGTTTTTATCAACCCACGCCACCGGGAATTCCGCCAGAATGTTGGTTTTCGGAATCACCACTTCAAAGCCC  >  NZ_CP009273/2536008‑2536195

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