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,717,498 G→C G374A (GGC→GCC)  ydhK → FUSC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,717,4980GC100.0% 322.1 / NA 85G374A (GGC→GCC) ydhKFUSC family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (47/38);  total (47/38)

CACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGT  >  NZ_CP009273/1717411‑1717582
                                                                                       |                                                                                    
cACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGgcg                                                                                    <  2:1951595/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  1:1206149/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  2:592720/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  2:503078/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  1:1394775/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  2:2556393/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  2:1334947/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  2:1871773/90‑1 (MQ=255)
   ggACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCAt                                                                                 <  1:2330435/90‑1 (MQ=255)
      cATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGa                                                                              >  2:2598769/1‑90 (MQ=255)
        tCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGt                                                                            >  2:1053633/1‑90 (MQ=255)
        tCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGt                                                                            >  2:115090/1‑90 (MQ=255)
                ggCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTc                                                                    >  1:910904/1‑90 (MQ=255)
                ggCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTc                                                                    >  1:995025/1‑90 (MQ=255)
                ggCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTc                                                                    >  2:2948357/1‑90 (MQ=255)
                     gTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAAt                                                               >  1:2820275/1‑90 (MQ=255)
                     gTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATCATGATTGCCGCATGGAGTATTGCTTCGCAAt                                                               <  2:2259675/90‑1 (MQ=255)
                          aCCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGAt                                                          >  2:2210826/1‑90 (MQ=255)
                          aCCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGAt                                                          >  1:1264469/1‑90 (MQ=255)
                              aTAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCg                                                      >  2:1179873/1‑90 (MQ=255)
                              aTAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCg                                                      >  1:1445852/1‑90 (MQ=255)
                               tAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCgg                                                     >  1:2932040/1‑90 (MQ=255)
                               tAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCgg                                                     >  2:654447/1‑90 (MQ=255)
                                     cGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCaa                                               <  2:1759641/90‑1 (MQ=255)
                                     cGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCaa                                               >  1:535072/1‑90 (MQ=255)
                                     cGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCaa                                               >  1:633273/1‑90 (MQ=255)
                                      gAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAAt                                              <  1:2537238/90‑1 (MQ=255)
                                      gAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAAt                                              <  2:2794152/90‑1 (MQ=255)
                                            aTGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGcc                                                       <  1:2651605/75‑1 (MQ=255)
                                            aTGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGcc                                                       >  2:2651605/1‑75 (MQ=255)
                                            aTGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTa                                        >  2:1095216/1‑90 (MQ=255)
                                                 gAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATgg                                                                               >  2:1291922/1‑46 (MQ=255)
                                                 gAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATgg                                                                               <  1:1291922/46‑1 (MQ=255)
                                                 gAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCt                                   <  2:2751913/90‑1 (MQ=255)
                                                  aGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTg                                  >  1:95506/1‑90 (MQ=255)
                                                   gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATgg                                                                               >  2:608313/1‑44 (MQ=255)
                                                   gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATgg                                                                               <  1:608313/44‑1 (MQ=255)
                                                   gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                                 >  2:1987414/1‑90 (MQ=255)
                                                   gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                                 >  2:2635971/1‑90 (MQ=255)
                                                   gCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTgg                                 >  2:1334813/1‑90 (MQ=255)
                                                    cGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                                <  1:1715355/90‑1 (MQ=255)
                                                    cGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                                <  2:1972688/90‑1 (MQ=255)
                                                    cGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                                >  2:54440/1‑90 (MQ=255)
                                                    cGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                                >  2:523084/1‑90 (MQ=255)
                                                    cGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                                <  1:454244/90‑1 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCAt                                          <  2:1736264/78‑1 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCAt                                          >  1:1736264/1‑78 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgcag                              <  2:1482005/90‑1 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgcag                              >  1:388035/1‑90 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgcag                              >  1:461020/1‑90 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgcag                              >  1:603297/1‑90 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                                >  1:2434470/1‑88 (MQ=255)
                                                      ggCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgc                                <  2:2434470/88‑1 (MQ=255)
                                                        cTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAAt                                              <  2:490332/72‑1 (MQ=255)
                                                        cTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAAt                                              >  1:490332/1‑72 (MQ=255)
                                                        cTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGgcagca                            >  2:2555367/1‑90 (MQ=255)
                                                          gCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAt                          >  1:1141945/1‑90 (MQ=255)
                                                          gCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAt                          >  1:1269081/1‑90 (MQ=255)
                                                          gCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAt                          >  1:612838/1‑90 (MQ=255)
                                                           cGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAtt                         >  2:555786/1‑90 (MQ=255)
                                                            gTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTa                        <  1:367094/90‑1 (MQ=255)
                                                              aCATTTTGTGCGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGAATTAACGCTGGCAGCAATTAGc                      <  1:3003828/90‑1 (MQ=255)
                                                              aCATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGc                      <  1:2442508/90‑1 (MQ=255)
                                                              aCATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGc                      <  1:2938235/90‑1 (MQ=255)
                                                                aTTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTg                    <  1:702053/90‑1 (MQ=255)
                                                                aTTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTg                    >  2:1206820/1‑90 (MQ=255)
                                                                aTTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTg                    <  1:1741450/90‑1 (MQ=255)
                                                                 ttttGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGc                   >  1:895143/1‑90 (MQ=255)
                                                                  tttGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCg                                                      >  2:2140454/1‑54 (MQ=255)
                                                                  tttGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCg                                                      <  1:2140454/54‑1 (MQ=255)
                                                                  tttGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCg                  <  2:645890/90‑1 (MQ=255)
                                                                  tttGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCg                  <  2:388035/90‑1 (MQ=255)
                                                                   ttGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGt                 <  2:632208/90‑1 (MQ=255)
                                                                    tGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGta                >  1:2879675/1‑90 (MQ=255)
                                                                    tGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGta                >  2:464740/1‑90 (MQ=255)
                                                                    tGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGta                >  1:1213092/1‑90 (MQ=255)
                                                                     gTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtac               >  1:1797294/1‑90 (MQ=255)
                                                                        cGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactct            <  2:535072/90‑1 (MQ=255)
                                                                         gTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactcta           >  2:2647588/1‑90 (MQ=255)
                                                                         gTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactcta           >  2:1574594/1‑90 (MQ=255)
                                                                          ttGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGtactctac          >  2:2632534/1‑90 (MQ=255)
                                                                             atgatgatTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTcc       <  1:735011/90‑1 (MQ=255)
                                                                                atgatTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTccgcc    <  2:3003577/90‑1 (MQ=255)
                                                                                atgatTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTccgcc    <  1:2920696/90‑1 (MQ=255)
                                                                                  gatTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGt  >  2:772061/1‑90 (MQ=255)
                                                                                    ttGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCt                                   <  2:2343684/55‑1 (MQ=255)
                                                                                       |                                                                                    
CACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGT  >  NZ_CP009273/1717411‑1717582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ 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
TCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTT  >  NZ_CP009273/1717408‑1717595
                                                                                          |                                                                                                 
TCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAG                                                                                          <  SRR3722030.550577/100‑1 (MQ=60)
      nnnnnnnnnnnnnnnnnnnnCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722030.2564512/80‑1 (MQ=60)
      GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722030.2372749/100‑1 (MQ=60)
      GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722030.1223976/100‑1 (MQ=60)
      GGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGC                                                                                    <  SRR3722030.1416711/100‑1 (MQ=60)
        gtcgtcggcagcgtcagatgtgtataagagacaGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTT                                                                                  <  SRR3722030.1311486/67‑1 (MQ=60)
        atgtcgtcggcagcgtcagatgtgtataagagacagGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTT                                                                                  <  SRR3722030.617743/64‑1 (MQ=60)
        accaccgagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagCATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTT                                                                                  <  SRR3722030.371457/42‑1 (MQ=60)
         CATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTC                                                                                 >  SRR3722030.924871/1‑100 (MQ=60)
         CATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTC                                                                                 >  SRR3722030.1010105/1‑100 (MQ=60)
              GGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAAT                                                                            >  SRR3722030.2872334/1‑100 (MQ=60)
                  gcgaccaccgagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGA                                                                        <  SRR3722030.1173225/39‑1 (MQ=60)
                   GGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGAT                                                                       >  SRR3722030.1283372/1‑100 (MQ=60)
                       CGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCG                                                                   >  SRR3722030.1468821/1‑100 (MQ=60)
                        GTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGG                                                                  >  SRR3722030.2986106/1‑100 (MQ=60)
                              CCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAA                                                            >  SRR3722030.543312/1‑100 (MQ=60)
                              CCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAA                                                            >  SRR3722030.643102/1‑100 (MQ=60)
                                gacagCGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATG                                                          <  SRR3722030.2700427/95‑1 (MQ=60)
                                 cagcgtcagatgtgtataagagacaGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGC                                                         <  SRR3722030.2178733/75‑1 (MQ=60)
                                         GAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGC                                                 <  SRR3722030.2583700/100‑1 (MQ=60)
                                           AGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTG                                               >  SRR3722030.96925/1‑100 (MQ=60)
                                               ATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCct                                           >  SRR3722030.1764910/1‑98 (MQ=60)
                                               ATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAG                                           >  SRR3722030.2478874/1‑100 (MQ=60)
                                               ATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAG                                           >  SRR3722030.612649/1‑100 (MQ=60)
                                               ATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAG                                           >  SRR3722030.468021/1‑100 (MQ=60)
                                               ATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAG                                           >  SRR3722030.393848/1‑100 (MQ=60)
                                                 GTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCctgtctct                                         >  SRR3722030.497839/1‑92 (MQ=60)
                                                   GGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAT                                       >  SRR3722030.1158951/1‑100 (MQ=60)
                                                   GGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAT                                       >  SRR3722030.622332/1‑100 (MQ=60)
                                                   GGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAAT                                       >  SRR3722030.1288093/1‑100 (MQ=60)
                                                       CGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGC                                   <  SRR3722030.461148/100‑1 (MQ=60)
                                                       CGGGCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGC                                   <  SRR3722030.1743421/100‑1 (MQ=60)
                                                          GCTGCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGC                                >  SRR3722030.908882/1‑100 (MQ=60)
                                                             GCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCnnnnnnnnnnnnnnnnnnnn                             >  SRR3722030.1466285/1‑80 (MQ=60)
                                                             GCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTA                             >  SRR3722030.2932813/1‑100 (MQ=60)
                                                             GCGTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTA                             >  SRR3722030.1231024/1‑100 (MQ=60)
                                                              CGTACATTCTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTAC                            >  SRR3722030.1827526/1‑100 (MQ=60)
                                                               GTACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACT                           <  SRR3722030.372565/100‑1 (MQ=60)
                                                                 ACATTTTGTGCGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGAATTAACGCTGGCAGCAATTAGCTGCGTACTCT                         <  SRR3722030.3059171/100‑1 (MQ=60)
                                                                 ACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCT                         <  SRR3722030.2487081/100‑1 (MQ=60)
                                                                 ACATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCT                         <  SRR3722030.2992387/100‑1 (MQ=60)
                                                                   ATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTAC                       <  SRR3722030.1770231/100‑1 (MQ=60)
                                                                   ATTTTGTACGTTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTAC                       <  SRR3722030.712903/100‑1 (MQ=60)
                                                                             TTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTctgtctcttatacacatctgacgct             >  SRR3722030.2386281/1‑75 (MQ=60)
                                                                             TTGATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCG             >  SRR3722030.1081448/1‑100 (MQ=60)
                                                                                ATGATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAG          <  SRR3722030.746349/100‑1 (MQ=60)
                                                                                   ATGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCAC       <  SRR3722030.2974553/100‑1 (MQ=60)
                                                                                    TGATTGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACC      >  SRR3722030.998230/1‑100 (MQ=60)
                                                                                        TGCCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTT  <  SRR3722030.2730330/100‑1 (MQ=60)
                                                                                          |                                                                                                 
TCCCACGGACATCCGGCCTGGCTCGTCATACCGATAACGCCGAAGCTATGTGGAGCGGGCTGCGTACATTTTGTACGTTGATGATGATTGGCGCATGGAGTATTGCTTCGCAATGGGATGCCGGTGCCAATGCATTAACGCTGGCAGCAATTAGCTGCGTACTCTACTCCGCCGTCGCAGCACCGTTT  >  NZ_CP009273/1717408‑1717595

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