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
A2 F5 I179 R1 182 40.3 1547118 96.8% 1497610 120.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,540,434 A→C V227G (GTT→GGT)  gntP ← gluconate permease GntP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,540,4340AC86.2% 169.2 / 18.3 58V227G (GTT→GGT) gntPgluconate permease GntP
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base C (21/29);  total (24/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTA  >  NZ_CP009273/4540300‑4540569
                                                                                                                                      |                                                                                                                                       
cAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGcc                                                                                                                                     >  1:692673/1‑139 (MQ=255)
   tAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGa                                                                                                                                   >  2:581343/1‑138 (MQ=255)
   tAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGAGGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGaa                                                                                                                                  <  2:208507/139‑1 (MQ=255)
    aaaCATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGaaa                                                                                                                                 <  2:692673/139‑1 (MQ=255)
         ttGCAATCGGCGAGGACCGGATAAAGTTAACCACTCCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAggg                                                                                                                            <  2:563260/139‑1 (MQ=255)
           gCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAggg                                                                                                                            >  2:758935/1‑137 (MQ=255)
            cAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAg                                                                                                                         <  1:677745/139‑1 (MQ=255)
                cGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGGAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGatta                                                                                                                     <  2:423741/139‑1 (MQ=255)
                           cGAAAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGCCCCCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTAc                                                                                                          <  2:100214/139‑1 (MQ=255)
                            gATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTAcc                                                                                                         >  2:165839/1‑139 (MQ=255)
                                             cccAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGatta                                                                                                                     >  1:67142/1‑110 (MQ=255)
                                             cccAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGatta                                                                                                                     <  2:67142/110‑1 (MQ=255)
                                             cccAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTc                                                                                          <  2:463932/137‑1 (MQ=255)
                                                     ggggTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAatgatg                                                                                >  1:311971/1‑139 (MQ=255)
                                                               ttACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTgg                                                                         <  2:278417/136‑1 (MQ=255)
                                                               ttACCAGCCAGATATTGGCGATGGTGGTAGAGATCAAAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTgg                                                                         <  2:328258/136‑1 (MQ=255)
                                                                tACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCTCCGGGATCAGCGGCCCCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGgcg                                                                      <  1:218203/138‑1 (MQ=255)
                                                                   cagccagATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTc                                                                   >  1:147373/1‑138 (MQ=255)
                                                                          atatTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTggg                                                                        <  1:543746/126‑1 (MQ=255)
                                                                          atatTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTggg                                                                        >  2:543746/1‑126 (MQ=255)
                                                                                gCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAgg                                                      <  1:165839/138‑1 (MQ=255)
                                                                                   aTGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACt                                                  <  1:499208/139‑1 (MQ=255)
                                                                                    tggtggTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTc                                                                                          >  1:65929/1‑98 (MQ=255)
                                                                                    tggtggTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTc                                                                                          <  2:65929/98‑1 (MQ=255)
                                                                                    tggtggTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACtt                                                 >  1:81736/1‑139 (MQ=255)
                                                                                       tggtAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCgg                                              <  1:309764/139‑1 (MQ=255)
                                                                                       tggtAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCgg                                              <  2:81736/139‑1 (MQ=255)
                                                                                        ggtAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGc                                             <  1:339269/139‑1 (MQ=255)
                                                                                         gtAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGc                                             <  1:758935/138‑1 (MQ=255)
                                                                                             agaTCATAATGATCGCCGGGATCAGCGCCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGAt                                        <  2:29801/139‑1 (MQ=255)
                                                                                             agaTCATAAGGATCGCCGGGATCAGCGGCACCAGAATGGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGAt                                        <  1:631508/139‑1 (MQ=255)
                                                                                                tCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGAt                                        >  2:328501/1‑136 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCCCCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTAcc                                                                                                         <  1:247096/70‑1 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATc                                                                                                             <  2:17813/66‑1 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATc                                                                                                             >  1:17813/1‑66 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTAcc                                                                                                         >  2:247096/1‑70 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCgg                                                                                                       >  2:25488/1‑72 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCgg                                                                                                       <  1:25488/72‑1 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGa                                    >  1:50859/1‑139 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGa                                    >  1:93343/1‑139 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGa                                    >  2:547815/1‑139 (MQ=255)
                                                                                                 cATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGa                                    >  1:569066/1‑139 (MQ=255)
                                                                                                  aTAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGAc                                   <  2:648817/139‑1 (MQ=255)
                                                                                                      ttaTCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTgc                               <  2:607265/137‑1 (MQ=255)
                                                                                                      tGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTgc                               <  2:593107/139‑1 (MQ=255)
                                                                                                        aTCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTggg                                                                        <  1:40197/96‑1 (MQ=255)
                                                                                                        aTCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTggg                                                                        >  2:40197/1‑96 (MQ=255)
                                                                                                         tCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGTGCGCTCAAGATTGCCGCGGAACTTCAGCAGGATCAGACCTgcg                              >  2:732659/1‑137 (MQ=255)
                                                                                                          cGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGa                                                                                                                        >  2:304248/1‑46 (MQ=255)
                                                                                                          cGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGa                                                                                                                        <  1:304248/46‑1 (MQ=255)
                                                                                                           gCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAgg                                                                                                                             <  2:466869/40‑1 (MQ=255)
                                                                                                           gCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAgg                                                                                                                             >  1:466869/1‑40 (MQ=255)
                                                                                                           gCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTgcg                              >  2:669030/1‑135 (MQ=255)
                                                                                                                     gCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGAt                                                              >  2:514800/1‑93 (MQ=255)
                                                                                                                     gCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGAt                                                              <  1:514800/93‑1 (MQ=255)
                                                                                                                      cGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTggg               <  1:669030/139‑1 (MQ=255)
                                                                                                                           ccAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGt          <  2:452813/139‑1 (MQ=255)
                                                                                                                                   gAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTa  <  2:4142/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                       
CAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGATCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTA  >  NZ_CP009273/4540300‑4540569

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGC  >  NZ_CP009273/4540289‑4540572
                                                                                                                                                    |                                                                                                                                          
AACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCG                                                                                                                                       >  SRR3722133.701660/1‑151 (MQ=60)
                      GCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATAT                                                                                                                 <  SRR3722133.686546/151‑1 (MQ=60)
                        tcttttttttttttttcaagcagaagacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATC                                                                                                                <  SRR3722133.308420/70‑1 (MQ=60)
                                      aagacggcatacgatatctagtacggtctcgtgggctcggagatgtgtataagagacaGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCCCCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGC                                                                                                 <  SRR3722133.250480/93‑1 (MQ=60)
                                        gacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacaGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTT                                                                                               <  SRR3722133.25812/95‑1 (MQ=60)
                                             GTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCtgtctcttatacacatct                                                                                          >  SRR3722133.67986/1‑133 (MQ=60)
                                                     CTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGG                                                                                  >  SRR3722133.316257/1‑151 (MQ=60)
                                                                   GTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCA                                                                     >  SRR3722133.149311/1‑150 (MQ=60)
                                                                       cagTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATT                                                                <  SRR3722133.550986/148‑1 (MQ=60)
                                                                          cgtgggctcggagatgtgtataagagacAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATT                                                                <  SRR3722133.40701/120‑1 (MQ=60)
                                                                          TTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCTCCGGGATCAGCGGCCCCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGC                                                              <  SRR3722133.221166/150‑1 (MQ=60)
                                                                                GCCAGATATTGGCGATGGTGGTAAGGAAGGA‑‑‑ACAGTTCCAATTAAAATTGGAATTCCGGCCCCGGGGGGAATTCCAAGGCCGGGGCCACCCCCAAGGAAAATTCCGGAAAACCCCGGCCCCGGAAAAAAGGAAGGGGGGCCAAGGAATTATTC                                                                                                                     <  SRR3722133.438095/64‑4294967272 (MQ=60)
                                                                                    cggtctcgtgggctcggagatgtgtataagagacagCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAAC                                                      <  SRR3722133.521716/115‑1 (MQ=60)
                                                                                    GATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCtgtctcttatacacatctccgagcccacg                                                    >  SRR3722133.66765/1‑121 (MQ=60)
                                                                                    GATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTC                                                   >  SRR3722133.82768/1‑151 (MQ=60)
                                                                                          GGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAG                                              <  SRR3722133.168048/150‑1 (MQ=60)
                                                                                             GATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATC                                          <  SRR3722133.505948/151‑1 (MQ=60)
                                                                                                 GTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaaaaaaaaaaaaaaaaaacaaat                                       >  SRR3722133.204232/1‑58 (MQ=60)
                                                                                                 GTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgc                                      >  SRR3722133.18027/1‑91 (MQ=60)
                                                                                                 GTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGAC                                      >  SRR3722133.576610/1‑151 (MQ=60)
                                                                                                 GTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGAC                                      >  SRR3722133.94544/1‑151 (MQ=60)
                                                                                                 GTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGAC                                      >  SRR3722133.51498/1‑151 (MQ=60)
                                                                                                 GTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGAC                                      <  SRR3722133.314020/151‑1 (MQ=60)
                                                                                                 GTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCAGAGGAACTTCGACAGGNN                                           >  SRR3722133.388290/1‑146 (MQ=60)
                                                                                                  GGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACC                                     <  SRR3722133.344037/151‑1 (MQ=60)
                                                                                                   GGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACC                                     <  SRR3722133.768752/150‑1 (MQ=60)
                                                                                                       GAGA‑‑‑TCATAAGGATCGCCGGGATCAGCGGCACCAGAATGGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGC                                <  SRR3722133.639779/151‑1 (MQ=60)
                                                                                                              TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCctgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtcttctgcttgaaaaaaaaaaataaacaaca                               >  SRR3722133.473221/1‑62 (MQ=60)
                                                                                                                                   GCGGCACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAG       <  SRR3722133.677732/151‑1 (MQ=60)
                                                                                                                                       CACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCtgtctcttatacacatctccgagcccacgagaccgtactagatctcgtatgccgtctt   >  SRR3722133.514371/1‑93 (MQ=60)
                                                                                                                                        ACCAGAATCGAACCGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCctgtctcttatacacatctccgagcccacgagacc  >  SRR3722133.261015/1‑116 (MQ=60)
                                                                                                                                                    |                                                                                                                                          
AACCACCATCGCAATAAACATTGCAATCGGCGAGGAACCGATAAAGTTAACCACTTCCCAGGCAGGGGTATCTTTTACCAGCCAGATATTGGCGATGGTGGTAGAGA‑‑‑TCATAATGATCGCCGGGATCAGCGGCACCAGAATCGAAACGCCGAAAGAGGGCAGATTATTCATATCTACCGGTTGATCTGCTTTCAGGAATGATGGCGTTGGGCGCTCAAGATTGCCGAGGAACTTCGGCAGGATCAGACCTGCGCAGATTACACTTGGGATCGTCACCAGTACGC  >  NZ_CP009273/4540289‑4540572

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