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 F5 I97 R1 275 44.0 1733530 96.5% 1672856 122.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 574,600 A→C intergenic (‑251/+40) bor ← / ← ybcV serum resistance lipoprotein Bor/DUF1398 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273574,6000AC75.4% 114.6 / 43.1 57intergenic (‑251/+40)bor/ybcVserum resistance lipoprotein Bor/DUF1398 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (6/8);  new base C (19/24);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATA  >  NZ_CP009273/574464‑574733
                                                                                                                                        |                                                                                                                                     
gTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGccacca                                                                                                                                     >  2:241865/1‑139 (MQ=255)
gTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGccacca                                                                                                                                     >  1:281565/1‑139 (MQ=255)
 tGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCcaccac                                                                                                                                    >  2:195863/1‑139 (MQ=255)
   cGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCAccc                                                                                                                                  <  2:666379/139‑1 (MQ=255)
    gATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCAcccc                                                                                                                                 >  2:606723/1‑139 (MQ=255)
                        ccGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAAc                                                                                                             >  1:601253/1‑139 (MQ=255)
                        ccGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAAcc                                                                                                            >  1:536906/1‑139 (MQ=255)
                         cGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCg                                                                                                           <  2:536906/139‑1 (MQ=255)
                           gATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCg                                                                                                           >  2:283419/1‑138 (MQ=255)
                            aTTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGa                                                                                                         <  2:498891/139‑1 (MQ=255)
                             ttAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGAc                                                                                                        >  2:567237/1‑139 (MQ=255)
                                      tCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGa                                                                                                         <  2:499790/129‑1 (MQ=255)
                                      tCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGa                                                                                                         >  1:499790/1‑129 (MQ=255)
                                              ccTGACGGGGAACTTCTCTGCGGGAGTGTCCGGTAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCg                                                                                                           <  2:165364/119‑1 (MQ=255)
                                              ccTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCg                                                                                                           >  1:165364/1‑119 (MQ=255)
                                               cTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTgg                                                                                       >  2:215724/1‑138 (MQ=255)
                                               cTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                      <  1:372527/139‑1 (MQ=255)
                                               cTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                      <  1:510545/139‑1 (MQ=255)
                                               cTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                      <  1:606723/139‑1 (MQ=255)
                                                tGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                      >  2:168223/1‑138 (MQ=255)
                                                tGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGtt                                                                                                     <  2:163047/123‑1 (MQ=255)
                                                tGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                      >  2:576293/1‑138 (MQ=255)
                                                tGACGGGGAACTTCTCTGCGGGAGAGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGtt                                                                                                     >  1:163047/1‑123 (MQ=255)
                                                  aCGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGaaa                                                                                    <  1:567237/138‑1 (MQ=255)
                                                  aCGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGaaa                                                                                    <  1:215724/138‑1 (MQ=255)
                                                    ggggAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGAt                                                                                 <  2:431106/139‑1 (MQ=255)
                                                       gAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTtg                                                                              >  1:376951/1‑139 (MQ=255)
                                                            tctcTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAg                                                                         >  1:352223/1‑139 (MQ=255)
                                                             ctctGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTa                                                                          <  1:586925/137‑1 (MQ=255)
                                                                    ggAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGa                                                                 <  1:168223/139‑1 (MQ=255)
                                                                    ggAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGa                                                                 <  1:287233/139‑1 (MQ=255)
                                                                     gAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGaa                                                                >  1:296237/1‑139 (MQ=255)
                                                                        tgtCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGc                                                             <  2:376951/139‑1 (MQ=255)
                                                                          tCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctc                                                           <  1:283419/139‑1 (MQ=255)
                                                                          tCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGctc                                                           <  1:395917/139‑1 (MQ=255)
                                                                          tCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGct                                                            <  1:239172/138‑1 (MQ=255)
                                                                                            cGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGaa                                         <  2:296237/139‑1 (MQ=255)
                                                                                                 cacacaGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATc                                    <  1:241865/139‑1 (MQ=255)
                                                                                                   cacaGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCaa                                  <  2:352223/139‑1 (MQ=255)
                                                                                                    acaGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCaaa                                 <  1:576293/139‑1 (MQ=255)
                                                                                                        ggTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGta                             <  1:629603/139‑1 (MQ=255)
                                                                                                        ggTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGa                                                                                                         <  1:333995/63‑1 (MQ=255)
                                                                                                        ggTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGa                                                                                                         >  2:333995/1‑63 (MQ=255)
                                                                                                             aGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGt                         <  1:834897/138‑1 (MQ=255)
                                                                                                                  gTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAat                      >  2:15327/1‑136 (MQ=255)
                                                                                                                  gTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAat                      <  1:15327/136‑1 (MQ=255)
                                                                                                                    aCATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCtt                 >  2:564208/1‑139 (MQ=255)
                                                                                                                              tattatGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGt       >  1:353477/1‑139 (MQ=255)
                                                                                                                                ttatGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                      <  2:484193/58‑1 (MQ=255)
                                                                                                                                ttatGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGa                                                                                      >  1:484193/1‑58 (MQ=255)
                                                                                                                                ttatGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCaa                                  <  2:470481/110‑1 (MQ=255)
                                                                                                                                ttatGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCaa                                  >  1:470481/1‑110 (MQ=255)
                                                                                                                                  atGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGa    >  2:21171/1‑138 (MQ=255)
                                                                                                                                   tGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATg                                                              <  2:319782/79‑1 (MQ=255)
                                                                                                                                   tGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATg                                                              >  1:319782/1‑79 (MQ=255)
                                                                                                                                   tGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGct                                                            <  1:10773/81‑1 (MQ=255)
                                                                                                                                   tGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGct                                                            >  2:10773/1‑81 (MQ=255)
                                                                                                                                   tGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGata  >  1:829621/1‑139 (MQ=255)
                                                                                                                                   tGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGata  >  1:803612/1‑139 (MQ=255)
                                                                                                                                   tGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGat   <  1:21171/138‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
GTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATA  >  NZ_CP009273/574464‑574733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGG  >  NZ_CP009273/574453‑574748
                                                                                                                                                   |                                                                                                                                                    
ATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCAC                                                                                                                                                   >  SRR3722120.285604/1‑151 (MQ=60)
                        TCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACC                                                                                                                           >  SRR3722120.609607/1‑151 (MQ=60)
                        TCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCG                                                                                                                          >  SRR3722120.544498/1‑151 (MQ=60)
                                      caagcagaagacggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGATGCACACAGGGTTTAGCGCGTACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTT                                                                                                             <  SRR3722120.338778/85‑1 (MQ=60)
                                      GATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTT                                                                                                             >  SRR3722120.506961/1‑151 (MQ=60)
                                              ATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATctgtctctt                                                                                                      >  SRR3722120.167631/1‑141 (MQ=60)
                                                GTCCACAGCCCTGACGGGGAACTTCTCTGCGGCAGAGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGctgtct                                                                                                   >  SRR3722120.165287/1‑145 (MQ=60)
                                                       GCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGT                                                                                            >  SRR3722120.382506/1‑151 (MQ=60)
                                                         CCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGT                                                                                          <  SRR3722120.517827/151‑1 (MQ=60)
                                                         CCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGT                                                                                          <  SRR3722120.615144/151‑1 (MQ=60)
                                                         CCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGT                                                                                          <  SRR3722120.378030/151‑1 (MQ=60)
                                                            GACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAG                                                                                        <  SRR3722120.218743/150‑1 (MQ=60)
                                                            GACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAG                                                                                        <  SRR3722120.575165/150‑1 (MQ=60)
                                                            GACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGT                                                                                       >  SRR3722120.357322/1‑151 (MQ=60)
                                                                     CTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAAT                                                                              >  SRR3722120.300470/1‑151 (MQ=60)
                                                                       TCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAAT                                                                              <  SRR3722120.595105/149‑1 (MQ=60)
                                                                              GGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGT                                                                     <  SRR3722120.291360/151‑1 (MQ=60)
                                                                              GGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGT                                                                     <  SRR3722120.170522/151‑1 (MQ=60)
                                                                                   cgagattcgccttagtctcgtgggctcggagatgtgtataagagacagGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATA                                                                <  SRR3722120.10904/103‑1 (MQ=60)
                                                                                    GTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAG                                                               <  SRR3722120.287487/151‑1 (MQ=60)
                                                                                    GTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAG                                                               <  SRR3722120.401711/151‑1 (MQ=60)
                                                                                    GTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATA                                                                <  SRR3722120.242577/150‑1 (MQ=60)
                                                                                                           GCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGT                                        <  SRR3722120.245316/151‑1 (MQ=60)
                                                                                                              CACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAAT                                     <  SRR3722120.584310/151‑1 (MQ=60)
                                                                                                                  GGGTTTAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCT                                 <  SRR3722120.638321/151‑1 (MQ=60)
                                                                                                                       TAGCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTA                             <  SRR3722120.846148/150‑1 (MQ=60)
                                                                                                                         gtgtataagagacagGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATG                           <  SRR3722120.514067/135‑1 (MQ=60)
                                                                                                                         GCGCGTACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGT                          <  SRR3722120.15506/151‑1 (MQ=60)
                                                                                                                              TACATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTG                     >  SRR3722120.358599/1‑151 (MQ=60)
                                                                                                                                CATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaa                    >  SRR3722120.491166/1‑80 (MQ=60)
                                                                                                                                CATGTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAActgtctcttatacacatct                   >  SRR3722120.477302/1‑132 (MQ=60)
                                                                                                                                   GTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAActgtctcttatacacatctccgagcccacgagactaaggcgaatctcgt                 >  SRR3722120.324396/1‑101 (MQ=60)
                                                                                                                                   GTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATAT                >  SRR3722120.814461/1‑151 (MQ=60)
                                                                                                                                   GTATTGTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATAT                >  SRR3722120.840774/1‑151 (MQ=60)
                                                                                                                                        GTATTATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttg            >  SRR3722120.83381/1‑84 (MQ=60)
                                                                                                                                             ATGCCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCC       <  SRR3722120.21409/150‑1 (MQ=60)
                                                                                                                                                CCACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaaaaaaa    >  SRR3722120.159613/1‑76 (MQ=60)
                                                                                                                                                 CACCACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAActgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtct  >  SRR3722120.791961/1‑93 (MQ=60)
                                                                                                                                                   |                                                                                                                                                    
ATGTATGCCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGGGGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAAACGATGCACACAGGGTTTAGCGCGTACATGTATTGTATTATGCCAACACCCCGGTGCTGACACGGAAGAAACCGGACGTTATGATTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTATCAAAGGTATAGTAATATCTTTTATGTTCGTGGATATTTGTAATCCATCGG  >  NZ_CP009273/574453‑574748

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