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,613,753 (A)6→5 coding (377/435 nt) marR → multiple antibiotic resistance transcriptional regulator MarR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,613,7480A.98.5% 297.6 / ‑3.9 68coding (372/435 nt)marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (36/31);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATG  >  NZ_CP009273/1613659‑1613833
                                                                                         |                                                                                     
ggcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACa                                                                                       >  1:2067413/1‑90 (MQ=255)
 gcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaa                                                                                      >  2:987891/1‑90 (MQ=255)
   gTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaa                                                                                    >  1:2170807/1‑90 (MQ=255)
           aaaaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGAcg                                                                           >  2:2556891/1‑90 (MQ=255)
           aaaaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGAcg                                                                           >  2:1344277/1‑90 (MQ=255)
                 tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGAc                                                                     >  1:1110249/1‑90 (MQ=255)
                 tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGAc                                                                     <  2:616566/90‑1 (MQ=255)
                 tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGAc                                                                     <  1:1984791/90‑1 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACg                                                                    >  2:1138777/1‑90 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACg                                                                    >  2:1676933/1‑90 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACg                                                                    >  2:1169475/1‑90 (MQ=255)
                    caccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGaa                                                                  <  2:403794/90‑1 (MQ=255)
                    caccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGaa                                                                  <  2:2005858/90‑1 (MQ=255)
                     accGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAg                                                                 <  1:1865420/90‑1 (MQ=255)
                      ccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGt                                                                >  1:1159411/1‑90 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  2:1505615/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  1:1912721/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  2:1857119/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  2:1005308/90‑1 (MQ=255)
                           ggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCaa                                                           <  1:987891/90‑1 (MQ=255)
                                ggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACtt                                                      <  1:588006/90‑1 (MQ=255)
                                 gcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                     <  2:296879/90‑1 (MQ=255)
                                 gcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                     <  1:1073418/90‑1 (MQ=255)
                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                               <  2:219933/90‑1 (MQ=255)
                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                               <  1:25733/90‑1 (MQ=255)
                                        gtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAttt                                              <  1:1918674/90‑1 (MQ=255)
                                         tgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  2:729541/75‑1 (MQ=255)
                                         tgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            >  1:729541/1‑75 (MQ=255)
                                           aaCAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCt                                           >  1:1221441/1‑90 (MQ=255)
                                           aaCAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCt                                           <  1:500031/90‑1 (MQ=255)
                                           aaCAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCt                                           >  2:1307155/1‑90 (MQ=255)
                                            aCAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCtt                                          <  2:1191001/90‑1 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTg                                                               <  2:2624858/66‑1 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTg                                                               >  1:2624858/1‑66 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  1:977559/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  1:2567125/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  1:918906/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  1:803613/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  1:2110856/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  2:1529101/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  1:534846/1‑90 (MQ=255)
                                                 gCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGaa                                     >  2:1215857/1‑90 (MQ=255)
                                                      cAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTcc                                <  1:972613/90‑1 (MQ=255)
                                                      cAATAAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTcc                                <  1:129700/90‑1 (MQ=255)
                                                           aGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCg                           >  2:2220301/1‑90 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                     >  1:1496975/1‑61 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                     <  2:311674/61‑1 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                     <  2:1496975/61‑1 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                     >  1:311674/1‑61 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTaa                        >  2:1403206/1‑90 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTaa                        >  1:236126/1‑90 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTaa                        <  2:2067413/90‑1 (MQ=255)
                                                                gCCAGGACCTGCACCAAGAATTAACAAAACCCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAAc                      <  2:1408069/90‑1 (MQ=255)
                                                                gCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAAc                      >  2:629055/1‑90 (MQ=255)
                                                                  cAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACaa                    >  2:191948/1‑90 (MQ=255)
                                                                    ggACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACaaaa                  >  1:1261315/1‑90 (MQ=255)
                                                                       ccTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAag               <  2:2166655/90‑1 (MQ=255)
                                                                        cTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAaga              <  1:1307155/90‑1 (MQ=255)
                                                                        cTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAaga              <  2:243556/90‑1 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag             >  2:2332536/1‑90 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag             >  1:538088/1‑90 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag             >  2:1828523/1‑90 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAAGATTTACTTAAGAAATGCCTGCGGGAAACAAAAAagag             >  1:1708715/1‑90 (MQ=255)
                                                                          gCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAgg            >  1:2652125/1‑90 (MQ=255)
                                                                          gCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAgg            <  2:1324820/90‑1 (MQ=255)
                                                                          gCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAgg            >  2:1855677/1‑90 (MQ=255)
                                                                              cAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATg        >  2:2338661/1‑90 (MQ=255)
                                                                                  aaTTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa    <  2:2817082/90‑1 (MQ=255)
                                                                                  aaTTAAAAAAAACCTGACGGCGGACGAAGGGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa    <  1:2663799/90‑1 (MQ=255)
                                                                                    ttAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATg  >  2:608380/1‑90 (MQ=255)
                                                                                    ttAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATg  >  1:110203/1‑90 (MQ=255)
                                                                                         |                                                                                     
GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATG  >  NZ_CP009273/1613659‑1613833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 30 ≤ 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
TGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGACGCAATACTGACG  >  NZ_CP009273/1613666‑1613852
                                                                                       |                                                                                                   
TGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGAC                                                                                        >  SRR3722030.1126823/1‑100 (MQ=60)
     AAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGT                                                                                   >  SRR3722030.1176646/1‑100 (MQ=60)
          TACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAA                                                                              <  SRR3722030.2019614/100‑1 (MQ=60)
              ACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACT                                                                          <  SRR3722030.1897362/100‑1 (MQ=60)
                   CGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGT                                                                     <  SRR3722030.1945852/100‑1 (MQ=60)
                    GGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTA                                                                    <  SRR3722030.1002892/100‑1 (MQ=60)
                        CGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTctgtc                                                                >  SRR3722030.740796/1‑95 (MQ=60)
                         GGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGC                                                               <  SRR3722030.597120/100‑1 (MQ=60)
                          GCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCT                                                              >  SRR3722030.1239467/1‑100 (MQ=60)
                          GCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCT                                                              <  SRR3722030.1089491/100‑1 (MQ=60)
                              TATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGctgtctcttataca                                                          >  SRR3722030.2673144/1‑86 (MQ=60)
                              TATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAG                                                          >  SRR3722030.815973/1‑100 (MQ=60)
                              TATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAG                                                          >  SRR3722030.992441/1‑100 (MQ=60)
                              TATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAG                                                          >  SRR3722030.543080/1‑100 (MQ=60)
                              TATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAG                                                          >  SRR3722030.2614204/1‑100 (MQ=60)
                              TATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAG                                                          >  SRR3722030.2148510/1‑100 (MQ=60)
                              TATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAG                                                          >  SRR3722030.932987/1‑100 (MQ=60)
                                TGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAA                                                        <  SRR3722030.26137/100‑1 (MQ=60)
                                 GTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAA                                                       <  SRR3722030.1951950/100‑1 (MQ=60)
                                    AACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTC                                                    <  SRR3722030.507686/100‑1 (MQ=60)
                                             ATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTctgtctcttatacacatct                                           >  SRR3722030.1520982/1‑81 (MQ=60)
                                             ATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTctgtctcttatacacatct                                           >  SRR3722030.316215/1‑81 (MQ=60)
                                             ATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAA                                           >  SRR3722030.239510/1‑100 (MQ=60)
                                               CAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAAC                                         <  SRR3722030.987429/100‑1 (MQ=60)
                                               CAATAAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAAC                                         <  SRR3722030.131571/100‑1 (MQ=60)
                                                   TAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAA                                     >  SRR3722030.1280140/1‑100 (MQ=60)
                                                    AGTTGGCCAGGACCTGCACCAAAAATAAAAAAAAACCTGACGGCGGACAAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGCCCTGCCGTAAACAAAAA                                    <  SRR3722030.1423203/100‑1 (MQ=60)
                                                        GGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAG                                >  SRR3722030.546394/1‑100 (MQ=60)
                                                        GGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAAGATTTACTTAAGAAATGCCTGCGGGAAACAAAAAAGAG                                >  SRR3722030.1736598/1‑100 (MQ=60)
                                                         GCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGG                               >  SRR3722030.2700955/1‑100 (MQ=60)
                                                                 CTCCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGA                       <  SRR3722030.1327061/100‑1 (MQ=60)
                                                                   GCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATG                     >  SRR3722030.111821/1‑100 (MQ=60)
                                                                          GAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGAC              >  SRR3722030.2751614/1‑100 (MQ=60)
                                                                           AATTAAAAAAAACCTGACGGCGGACGAAGGGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGACG             <  SRR3722030.2712882/100‑1 (MQ=60)
                                                                                  CAAAAACCTGACGGCGGACAAAGGGGCAACATTTGAGTATTTGCTAAAAAAAGCCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGACGCAATAC       <  SRR3722030.1252535/100‑1 (MQ=60)
                                                                                     AAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGACGCAATACTGA    <  SRR3722030.1425312/100‑1 (MQ=60)
                                                                                     AAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTCGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGACGCAATACTGA    >  SRR3722030.2957990/1‑100 (MQ=60)
                                                                                       ACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGACGCAATACTGACG  <  SRR3722030.1233817/100‑1 (MQ=60)
                                                                                       ACCTGACGGCGGACGAAGGGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGGATGACGATGTCCAGACGCAATACTGACG  >  SRR3722030.2926138/1‑100 (MQ=60)
                                                                                       |                                                                                                   
TGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATGTCCAGACGCAATACTGACG  >  NZ_CP009273/1613666‑1613852

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