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 F3 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,782,276 T→G V192G (GTG→GGG)  ppsR → posphoenolpyruvate synthetase regulatory kinase/phosphorylase PpsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,782,2760TG100.0% 241.3 / NA 64V192G (GTG→GGG) ppsRposphoenolpyruvate synthetase regulatory kinase/phosphorylase PpsR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (46/18);  total (46/18)

CCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAA  >  NZ_CP009273/1782193‑1782361
                                                                                   |                                                                                     
ccccACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTAcc                                                                                 >  1:37332/1‑90 (MQ=255)
 cccACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTAccc                                                                                >  2:1188990/1‑90 (MQ=255)
  ccaccaGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCcg                                                                               >  2:2790318/1‑90 (MQ=255)
  ccaccaGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCcg                                                                               >  1:571037/1‑90 (MQ=255)
   caccaGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTAcc                                                                                 >  2:2469305/1‑87 (MQ=255)
   caccaGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTAcc                                                                                 <  1:2469305/87‑1 (MQ=255)
   caccaGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCcgc                                                                              >  2:453749/1‑90 (MQ=255)
   caccaGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCcgc                                                                              >  1:1995323/1‑90 (MQ=255)
   caccaGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCcgc                                                                              >  2:2668553/1‑90 (MQ=255)
        gTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGc                                                                         >  1:55623/1‑90 (MQ=255)
        gTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGc                                                                         >  1:220436/1‑90 (MQ=255)
             tATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCaaa                                                                    >  2:2707517/1‑90 (MQ=255)
               tCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAAcc                                                                  >  1:2305485/1‑90 (MQ=255)
               tCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAAcc                                                                  >  1:412142/1‑90 (MQ=255)
               tCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAAcc                                                                  <  2:1062774/90‑1 (MQ=255)
                 tGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGc                                                                >  1:2028249/1‑90 (MQ=255)
                 tGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGc                                                                >  1:2014695/1‑90 (MQ=255)
                 tGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGc                                                                >  1:1772643/1‑90 (MQ=255)
                  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTc                                                                           <  1:2601816/78‑1 (MQ=255)
                  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTc                                                                           >  2:2601816/1‑78 (MQ=255)
                  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCt                                                               >  1:600218/1‑90 (MQ=255)
                  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCt                                                               >  1:1755281/1‑90 (MQ=255)
                  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCt                                                               >  1:351113/1‑90 (MQ=255)
                  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCt                                                               >  2:86246/1‑90 (MQ=255)
                  ggCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCt                                                               >  1:1313572/1‑90 (MQ=255)
                   gCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCtt                                                              <  2:55623/90‑1 (MQ=255)
                   gCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCtt                                                              <  1:2461950/90‑1 (MQ=255)
                   gCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCtt                                                              <  1:1188089/90‑1 (MQ=255)
                   gCAATGCAGTTTGGTATCCCCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCtt                                                              <  1:2979606/90‑1 (MQ=255)
                     aaTGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCa                                                            <  1:2870675/90‑1 (MQ=255)
                      aTGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAg                                                           >  1:419905/1‑90 (MQ=255)
                      aTGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAg                                                           >  1:1995793/1‑90 (MQ=255)
                      aTGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAg                                                           >  1:267554/1‑90 (MQ=255)
                      aTGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTAGGGCTACCCGCGTCGCTCAAACCGCTTCAg                                                           >  1:832824/1‑90 (MQ=255)
                            tttGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATaaa                                                     >  2:607400/1‑90 (MQ=255)
                                  aTCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTc                                               <  1:625135/90‑1 (MQ=255)
                                     cgcgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCtt                                                              <  2:2350036/72‑1 (MQ=255)
                                     cgcgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCtt                                                              >  1:2350036/1‑72 (MQ=255)
                                     cgcgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGc                                            >  1:2719818/1‑90 (MQ=255)
                                     cgcgcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGc                                            >  1:2202864/1‑90 (MQ=255)
                                        gcgGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTg                                         <  1:62323/90‑1 (MQ=255)
                                                   ccccTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGAccc                              >  1:2985664/1‑90 (MQ=255)
                                                     ccTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCaa                                                                     <  1:2718365/49‑1 (MQ=255)
                                                     ccTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCaa                                                                     >  2:2718365/1‑49 (MQ=255)
                                                     ccTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCgg                            >  1:655935/1‑90 (MQ=255)
                                                       tttATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGaa                          >  1:45013/1‑90 (MQ=255)
                                                       tttATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGaa                          <  2:655935/90‑1 (MQ=255)
                                                        ttATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCa                                                         <  2:1598023/58‑1 (MQ=255)
                                                        ttATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCa                                                         >  1:1598023/1‑58 (MQ=255)
                                                          aTTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGt                       >  1:1890363/1‑90 (MQ=255)
                                                            tGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCt                     >  1:425567/1‑90 (MQ=255)
                                                            tGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCt                     >  1:2050858/1‑90 (MQ=255)
                                                               cgacgaTATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTggc                  >  2:903406/1‑90 (MQ=255)
                                                                     tatGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGAt            <  2:1995793/90‑1 (MQ=255)
                                                                       tGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTc          >  2:692636/1‑90 (MQ=255)
                                                                            aaTCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGaa                          >  2:367959/1‑69 (MQ=255)
                                                                            aaTCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGaa                          <  1:367959/69‑1 (MQ=255)
                                                                             aTCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAgg    >  2:296995/1‑90 (MQ=255)
                                                                             aTCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAgg    >  1:1328033/1‑90 (MQ=255)
                                                                             aTCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAgg    >  1:1125208/1‑90 (MQ=255)
                                                                             aTCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAgg    >  2:2081028/1‑90 (MQ=255)
                                                                              tCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa   <  2:2305485/90‑1 (MQ=255)
                                                                              tCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGa   >  2:2030789/1‑90 (MQ=255)
                                                                               cTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGaa  <  2:2050858/90‑1 (MQ=255)
                                                                                   |                                                                                     
CCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAA  >  NZ_CP009273/1782193‑1782361

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ 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
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAG  >  NZ_CP009273/1782179‑1782375
                                                                                                 |                                                                                                   
gtcagatgtgtataagagacagGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC                                                                                                   <  SRR3722030.2481986/78‑1 (MQ=60)
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC                                                                                                   <  SRR3722030.332291/100‑1 (MQ=60)
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC                                                                                                   <  SRR3722030.2455691/100‑1 (MQ=60)
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC                                                                                                   <  SRR3722030.2466803/100‑1 (MQ=60)
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC                                                                                                   <  SRR3722030.2387539/100‑1 (MQ=60)
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC                                                                                                   <  SRR3722030.575316/100‑1 (MQ=60)
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGC                                                                                                   <  SRR3722030.1292926/100‑1 (MQ=60)
    GTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACC                                                                                               >  SRR3722030.37913/1‑100 (MQ=60)
      GGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCG                                                                                             >  SRR3722030.579885/1‑100 (MQ=60)
       GTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGC                                                                                            >  SRR3722030.2030389/1‑100 (MQ=60)
            ACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGC                                                                                       >  SRR3722030.56475/1‑100 (MQ=60)
            ACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGC                                                                                       >  SRR3722030.223586/1‑100 (MQ=60)
              CCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTC                                                                                     <  SRR3722030.2514416/100‑1 (MQ=60)
                   CCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACC                                                                                >  SRR3722030.418370/1‑100 (MQ=60)
                   CCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACC                                                                                >  SRR3722030.2347269/1‑100 (MQ=60)
                    agGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCG                                                                               <  SRR3722030.2649597/98‑1 (MQ=60)
                     AGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGC                                                                              >  SRR3722030.1802202/1‑100 (MQ=60)
                     AGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGC                                                                              >  SRR3722030.2064074/1‑100 (MQ=60)
                     AGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGC                                                                              >  SRR3722030.2050205/1‑100 (MQ=60)
                      GTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCT                                                                             >  SRR3722030.1333651/1‑100 (MQ=60)
                      GTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCT                                                                             >  SRR3722030.1784398/1‑100 (MQ=60)
                      GTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCT                                                                             >  SRR3722030.356315/1‑100 (MQ=60)
                      GTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCT                                                                             >  SRR3722030.609515/1‑100 (MQ=60)
                          gtcggcagcgtcagatgtgtataagagacaGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAG                                                                         <  SRR3722030.2768592/70‑1 (MQ=60)
                          GTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAG                                                                         >  SRR3722030.2030870/1‑100 (MQ=60)
                          GTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAG                                                                         >  SRR3722030.271434/1‑100 (MQ=60)
                          GTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAG                                                                         >  SRR3722030.426244/1‑100 (MQ=60)
                          GTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTAGGGCTACCCGCGTCGCTCAAACCGCTTCAG                                                                         >  SRR3722030.845594/1‑100 (MQ=60)
                                 GCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAAT                                                                  <  SRR3722030.2506935/100‑1 (MQ=60)
                                 GCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAAT                                                                  <  SRR3722030.1205686/100‑1 (MQ=60)
                                 GCAATGCAGTTTGGTATCCCCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAAT                                                                  <  SRR3722030.3034491/100‑1 (MQ=60)
                                   AATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTG                                                                <  SRR3722030.2923661/100‑1 (MQ=60)
                                         GTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATctgtctct                                                          >  SRR3722030.2392739/1‑92 (MQ=60)
                                         GTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGC                                                          >  SRR3722030.2770068/1‑100 (MQ=60)
                                         GTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGC                                                          >  SRR3722030.2242494/1‑100 (MQ=60)
                                                ATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTA                                                   <  SRR3722030.634833/100‑1 (MQ=60)
                                                      GCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACC                                             <  SRR3722030.63280/100‑1 (MQ=60)
                                                       CGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCC                                            >  SRR3722030.3040656/1‑100 (MQ=60)
                                                         GCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGG                                          >  SRR3722030.666092/1‑100 (MQ=60)
                                                           AAACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAA                                        >  SRR3722030.45722/1‑100 (MQ=60)
                                                            AACTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCctgtctcttatacacatctgac                                       >  SRR3722030.1623943/1‑78 (MQ=60)
                                                              CTACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGT                                     >  SRR3722030.1922934/1‑100 (MQ=60)
                                                                ACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCT                                   >  SRR3722030.432015/1‑100 (MQ=60)
                                                                ACCCCTTTATTGCCGACGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGCCT                                   >  SRR3722030.2087202/1‑100 (MQ=60)
                                                                     ataagagacagCGATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGG                              <  SRR3722030.373444/89‑1 (MQ=60)
                                                                                 GATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGctgtctcttatacacatctgacgctgccgacgacatggccagtgta                  >  SRR3722030.1069419/1‑54 (MQ=60)
                                                                                 GATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGG                  >  SRR3722030.1348446/1‑100 (MQ=60)
                                                                                 GATATGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGG                  >  SRR3722030.1141992/1‑100 (MQ=60)
                                                                                     TGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACG              >  SRR3722030.2791987/1‑100 (MQ=60)
                                                                                     TGGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACG              >  SRR3722030.777076/1‑100 (MQ=60)
                                                                                      GGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTctgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctc             >  SRR3722030.1711841/1‑48 (MQ=60)
                                                                                      GGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTctgtctcttatacacatctgacgctgccgacgacatggc             >  SRR3722030.200876/1‑61 (MQ=60)
                                                                                      GGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGT             >  SRR3722030.158886/1‑100 (MQ=60)
                                                                                      GGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGT             >  SRR3722030.2863459/1‑100 (MQ=60)
                                                                                      GGATAATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGT             >  SRR3722030.66732/1‑100 (MQ=60)
                                                                                           ATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATctgtctcttatacacatctgacgctgccgacgacatggccagtgtagatctcggtggt        >  SRR3722030.1600585/1‑42 (MQ=60)
                                                                                           ATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA        >  SRR3722030.1656489/1‑100 (MQ=60)
                                                                                           ATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA        >  SRR3722030.169649/1‑100 (MQ=60)
                                                                                           ATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA        >  SRR3722030.371784/1‑100 (MQ=60)
                                                                                           ATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA        >  SRR3722030.2760448/1‑100 (MQ=60)
                                                                                           ATCTGGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGA        >  SRR3722030.22610/1‑100 (MQ=60)
                                                                                               GGGGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAAC    <  SRR3722030.715299/100‑1 (MQ=60)
                                                                                                 GGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAG  <  SRR3722030.2050324/100‑1 (MQ=60)
                                                                                                 |                                                                                                   
CGCTGTGGTAAAACCCCCACCAGTCTGTATCTGGCAATGCAGTTTGGTATCCGCGCGGCAAACTACCCCTTTATTGCCGACGATATGGATAATCTGGTGCTACCCGCGTCGCTCAAACCGCTTCAGCATAAATTGTTCGGCCTGACTATCGACCCGGAACGTCTGGCGGCGATTCGCGAGGAACGTCGGGAGAACAG  >  NZ_CP009273/1782179‑1782375

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