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 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,6870GC100.0% 580.9 / NA 156G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (80/76);  total (80/76)

Only 100 of 156 total aligned reads displayed.

TCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785601‑3785772
                                                                                      |                                                                                     
tCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCg                                                                                    <  1:352309/90‑1 (MQ=255)
  ttCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGaa                                                                                  >  1:951750/1‑90 (MQ=255)
    ccATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAc                                                                                <  1:1652165/90‑1 (MQ=255)
          ccAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgag                                                                          <  2:2841464/90‑1 (MQ=255)
           cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAc                                                                                >  1:2826585/1‑83 (MQ=255)
           cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                         >  1:2081333/1‑90 (MQ=255)
             ggAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAACCCGTGAGAAt                                                                       <  1:2573197/90‑1 (MQ=255)
             ggAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAAt                                                                       >  2:1372632/1‑90 (MQ=255)
             ggAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAAt                                                                       >  1:882267/1‑90 (MQ=255)
              gAATGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGACGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      <  1:1771116/90‑1 (MQ=255)
                  gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:1285520/1‑90 (MQ=255)
                  gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:2492921/1‑90 (MQ=255)
                  gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:229125/1‑90 (MQ=255)
                  gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:2174944/1‑90 (MQ=255)
                  gccgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCAtt                                                                  >  1:2137642/1‑90 (MQ=255)
                    cgccAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTcc                                                                <  2:1518633/90‑1 (MQ=255)
                           ttttGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCttt                                                         >  1:1416715/1‑90 (MQ=255)
                             ttGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgc                                                       <  1:1940257/90‑1 (MQ=255)
                                  cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  2:2064555/90‑1 (MQ=255)
                                  cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  2:1900794/90‑1 (MQ=255)
                                  cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  2:407522/90‑1 (MQ=255)
                                  cAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCACCCCGAACCCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  <  1:2650585/90‑1 (MQ=255)
                                   aaGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCcgccgc                                                 >  1:1042475/1‑90 (MQ=255)
                                      ctctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAat                                              <  2:1285520/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCTGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:2815645/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:2169815/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:806414/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:1133409/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:2137642/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:2379047/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:2653066/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAAAAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:1367553/90‑1 (MQ=255)
                                       tctTTATGGCTGTCCTGAGTGCCGCAAAAAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:1133408/90‑1 (MQ=255)
                                         tttATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTcc                                                                <  2:2529328/69‑1 (MQ=255)
                                         tttATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcag                                           >  2:404417/1‑90 (MQ=255)
                                         tttATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:2581340/88‑1 (MQ=255)
                                               gCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGccgccg                                                  >  2:2483337/1‑77 (MQ=255)
                                                cTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGga                                    >  2:2731883/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATATAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1371870/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTga                                                                           <  2:468430/49‑1 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTga                                                                           >  1:468430/1‑49 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTg                                                                            >  2:635499/1‑48 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTgccgc                                                    >  1:1283217/1‑72 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatcagatc                                        <  1:2769425/84‑1 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  1:2607760/79‑1 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             >  2:2607760/1‑79 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:2753939/79‑1 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:205588/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:762994/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:440114/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:1804330/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:16421/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:2784531/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:2619578/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  2:1305877/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:809287/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:372181/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1859730/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1878014/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1486735/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:578353/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:145169/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:1442070/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:961954/1‑90 (MQ=255)
                                                  gTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatg                                  >  1:987458/1‑90 (MQ=255)
                                                       gAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCg                             >  1:1590242/1‑90 (MQ=255)
                                                         gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAACCCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           <  2:1456729/90‑1 (MQ=255)
                                                         gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                           <  1:2404712/90‑1 (MQ=255)
                                                           gCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAg                         >  2:2193044/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      <  2:2255951/44‑1 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATc                                                                      >  1:2255951/1‑44 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCTGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  2:1173342/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             >  1:1803634/1‑69 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                             <  2:1803634/69‑1 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:1065991/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:2962358/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:267634/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:1686439/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:2500948/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:2420388/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  1:698431/1‑90 (MQ=255)
                                                            ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGc                        >  2:2666466/1‑90 (MQ=255)
                                                              gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:160929/90‑1 (MQ=255)
                                                              gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  2:1242816/1‑90 (MQ=255)
                                                              gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      <  2:2377136/90‑1 (MQ=255)
                                                              gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  1:592234/1‑90 (MQ=255)
                                                              gCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGccc                      >  2:511250/1‑90 (MQ=255)
                                                                 aaTAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCgag                   <  2:1207681/90‑1 (MQ=255)
                                                                   tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                 <  1:205588/90‑1 (MQ=255)
                                                                         gCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgt           >  1:177869/1‑90 (MQ=255)
                                                                          cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          <  1:511250/90‑1 (MQ=255)
                                                                          cACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAgttgtt          <  2:1368866/90‑1 (MQ=255)
                                                                               aaGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGcac     <  2:1590242/90‑1 (MQ=255)
                                                                                 gCCATCCGGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  1:1917833/90‑1 (MQ=255)
                                                                                 gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  2:2896569/90‑1 (MQ=255)
                                                                                 gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  1:2674170/90‑1 (MQ=255)
                                                                                 gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  2:1821445/90‑1 (MQ=255)
                                                                                 gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  2:590909/90‑1 (MQ=255)
                                                                                 gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  2:692487/90‑1 (MQ=255)
                                                                                 gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa   <  2:1065991/90‑1 (MQ=255)
                                                                                      |                                                                                     
TCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAA  >  NZ_CP009273/3785601‑3785772

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTGG  >  NZ_CP009273/3785588‑3785786
                                                                                                   |                                                                                                   
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                     >  SRR3722030.2893927/1‑100 (MQ=60)
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                     >  SRR3722030.1715225/1‑100 (MQ=60)
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                     >  SRR3722030.2707017/1‑100 (MQ=60)
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATC                                                                                                     >  SRR3722030.2079210/1‑100 (MQ=60)
     GAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAA                                                                                                >  SRR3722030.966274/1‑100 (MQ=60)
       nnnnnnnnnnnnnnnnnnnnGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                              <  SRR3722030.1902481/80‑1 (MQ=60)
       ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                              <  SRR3722030.3070520/100‑1 (MQ=60)
       ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                              <  SRR3722030.269248/100‑1 (MQ=60)
       ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                              <  SRR3722030.581741/100‑1 (MQ=60)
       ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                              <  SRR3722030.1526277/100‑1 (MQ=60)
       ATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAAC                                                                                              <  SRR3722030.917385/100‑1 (MQ=60)
             TCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAG                                                                                        <  SRR3722030.357527/100‑1 (MQ=60)
             TCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAG                                                                                        <  SRR3722030.2192871/100‑1 (MQ=60)
              CTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGA                                                                                       >  SRR3722030.2118327/1‑100 (MQ=60)
              CTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGA                                                                                       >  SRR3722030.2878759/1‑100 (MQ=60)
               actggccatgtcgtcggcagcgtcagatgtgtataagagacagTGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAA                                                                                      <  SRR3722030.1113670/57‑1 (MQ=60)
               TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAA                                                                                      <  SRR3722030.2627849/100‑1 (MQ=60)
               TTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAA                                                                                      >  SRR3722030.1985253/1‑100 (MQ=60)
                TCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAAT                                                                                     >  SRR3722030.895799/1‑100 (MQ=60)
                 CCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATC                                                                                    <  SRR3722030.1678988/100‑1 (MQ=60)
                     cgtcggcagcgtcagatgtgtataagagacagCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                                <  SRR3722030.645372/68‑1 (MQ=60)
                     CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                                >  SRR3722030.2175857/1‑100 (MQ=60)
                     CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                                >  SRR3722030.2213953/1‑100 (MQ=60)
                     CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                                >  SRR3722030.1304929/1‑100 (MQ=60)
                     CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                                >  SRR3722030.232401/1‑100 (MQ=60)
                     CCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATT                                                                                >  SRR3722030.2538518/1‑100 (MQ=60)
                      CCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTC                                                                               >  SRR3722030.818815/1‑100 (MQ=60)
                          GGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAACCCGTGAGAATCCATTCCCGC                                                                           <  SRR3722030.2620401/100‑1 (MQ=60)
                           GAATGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGACGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCC                                                                          <  SRR3722030.1800648/100‑1 (MQ=60)
                              GGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTT                                                                       >  SRR3722030.1439078/1‑100 (MQ=60)
                                      GTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGC                                                               >  SRR3722030.2422313/1‑100 (MQ=60)
                                      GTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGC                                                               >  SRR3722030.1058174/1‑100 (MQ=60)
                                      GTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGC                                                               >  SRR3722030.2949993/1‑100 (MQ=60)
                                          TTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATC                                                           <  SRR3722030.1974064/100‑1 (MQ=60)
                                            GTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCtgtctcttatacacatctgacgctgccgacgacatg                                                         >  SRR3722030.2933518/1‑64 (MQ=60)
                                            GTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCctgtctcttat                                                         >  SRR3722030.2575632/1‑89 (MQ=60)
                                               agacagCTTTATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                      <  SRR3722030.1836700/94‑1 (MQ=60)
                                               CAGGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                      <  SRR3722030.2528733/100‑1 (MQ=60)
                                               CAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCACCCCGAACCCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATC                                                      <  SRR3722030.2699392/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.2171349/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCACCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.2701912/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.2628691/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.2208724/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.1150308/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAAAAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.1388911/100‑1 (MQ=60)
                                                    TCTTTATGGCTGTCCTGAGTGCCGCAAAAAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.1150307/100‑1 (MQ=60)
                                                    GCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGAT                                                 <  SRR3722030.2655648/100‑1 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATATAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1393315/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCtgtctcttatacacatctgacgctgccgac                                                >  SRR3722030.475558/1‑70 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTTCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1494968/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCtgtctcttatacaca                                                >  SRR3722030.1234173/1‑85 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGActgtctct                                                >  SRR3722030.1302567/1‑92 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.265998/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.976606/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.377729/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.587314/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1002456/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.978297/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.821724/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1891532/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1910284/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1464963/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.147257/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1510536/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.192294/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATG                                                >  SRR3722030.1782180/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATC                                                >  SRR3722030.1826483/1‑100 (MQ=60)
                                                     CTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATC                                                >  SRR3722030.2804802/1‑100 (MQ=60)
                                                         ATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTC                                            <  SRR3722030.2820575/100‑1 (MQ=60)
                                                          TGGCTGTCCTGAGTGCCGCAGATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCG                                           >  SRR3722030.2702284/1‑100 (MQ=60)
                                                          TGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCG                                           >  SRR3722030.1616020/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCctgtctcttatacacatctgacgctgccgacgacat                                      >  SRR3722030.2296713/1‑64 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATctgtctcttat                                      >  SRR3722030.1834024/1‑89 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.3063231/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.1713820/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.709232/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.930255/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.1081985/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.3016971/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.2546718/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.271516/1‑100 (MQ=60)
                                                               GTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGC                                      >  SRR3722030.2464509/1‑100 (MQ=60)
                                                                 CCTGAGTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCC                                    >  SRR3722030.601419/1‑100 (MQ=60)
                                                                      GTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAT                               <  SRR3722030.2448523/100‑1 (MQ=60)
                                                                      GTGCCGCAAATAAAACCCACCGAAGCCACCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAT                               <  SRR3722030.1445012/100‑1 (MQ=60)
                                                                      GGGCCGAAAATAAAACGCACCGAACCCATCCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAT                               <  SRR3722030.1652473/100‑1 (MQ=60)
                                                                           GCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTG                          <  SRR3722030.413819/100‑1 (MQ=60)
                                                                            CAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGT                         >  SRR3722030.180431/1‑100 (MQ=60)
                                                                                TAAAACGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGC                     <  SRR3722030.208516/100‑1 (MQ=60)
                                                                                       CACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG              >  SRR3722030.2993898/1‑100 (MQ=60)
                                                                                       CACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG              <  SRR3722030.519088/100‑1 (MQ=60)
                                                                                       CACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG              >  SRR3722030.665220/1‑100 (MQ=60)
                                                                                       CACCGAAGCCATCCCAAACCCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAG              <  SRR3722030.2218749/100‑1 (MQ=60)
                                                                                              GCCATCCGGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA       <  SRR3722030.1951086/100‑1 (MQ=60)
                                                                                              GCCATCCCGAAACGGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA       <  SRR3722030.2723469/100‑1 (MQ=60)
                                                                                              GCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA       <  SRR3722030.525989/100‑1 (MQ=60)
                                                                                              GCCATCCCAAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGA       <  SRR3722030.1394088/100‑1 (MQ=60)
                                                                                               CCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGAC      <  SRR3722030.1637847/100‑1 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGctgtc   >  SRR3722030.2547738/1‑95 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACtgt   >  SRR3722030.2991733/1‑97 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG   >  SRR3722030.2698037/1‑100 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG   >  SRR3722030.1567722/1‑100 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG   >  SRR3722030.1803055/1‑100 (MQ=60)
                                                                                                  TCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTG   >  SRR3722030.2804665/1‑100 (MQ=60)
                                                                                                   CCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTGG  >  SRR3722030.2812754/1‑100 (MQ=60)
                                                                                                   |                                                                                                   
GTAGAGAATCGCATCTTCCATCCCCAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAAAAGTTAATGACGTGG  >  NZ_CP009273/3785588‑3785786

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