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 F2 I194 R1 197 32.2 1845014 95.2% 1756453 85.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,356,378 (A)8→7 coding (213/2382 nt) yhcD → fimbrial biogenesis outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,356,3710A.96.5% 221.9 / 2.1 58coding (206/2382 nt)yhcDfimbrial biogenesis outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (46/10);  total (48/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.69e-01

GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA  >  NZ_CP009273/3356284‑3356455
                                                                                       |                                                                                    
gtgtAGATCCTAATATATCGCATTATTTTTCCCGTTCAACCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaa                                                                                    >  2:461924/1‑90 (MQ=255)
   tAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaa                                                                                 <  2:318370/90‑1 (MQ=255)
    aGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGaaaaaaa                                                                                >  2:61586/1‑90 (MQ=255)
         cTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCaa                                                                          >  2:79598/1‑90 (MQ=255)
          tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  2:915292/1‑90 (MQ=255)
          tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  1:459539/1‑90 (MQ=255)
          tAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  1:532268/1‑90 (MQ=255)
            atatatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACAt                                                                       <  2:829429/90‑1 (MQ=255)
            atatatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACAt                                                                       >  1:737554/1‑90 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         <  1:386286/85‑1 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAc                                                                         >  2:386286/1‑85 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGc                                                                    >  2:41731/1‑90 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGc                                                                    >  1:320980/1‑90 (MQ=255)
               tatCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGc                                                                    >  2:15394/1‑90 (MQ=255)
                atCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCt                                                                   >  1:252375/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:232007/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:297905/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:262485/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:917398/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  1:830070/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  <  1:517775/90‑1 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:559225/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  2:770527/1‑90 (MQ=255)
                 tCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTa                                                                  >  1:388672/1‑90 (MQ=255)
                   gCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTAcg                                                                >  1:746805/1‑90 (MQ=255)
                   gCATTATTTTTCCCGTCCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTAcg                                                                >  2:524255/1‑90 (MQ=255)
                        attTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCtt                                                           <  2:513683/90‑1 (MQ=255)
                        attTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCtt                                                           >  2:692403/1‑90 (MQ=255)
                        attTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCtt                                                           <  1:232865/90‑1 (MQ=255)
                             tCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTtgatg                                                      >  2:437212/1‑90 (MQ=255)
                                cGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGaaa                                                   >  1:682051/1‑90 (MQ=255)
                                  ttCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAAt                                                 >  2:112546/1‑90 (MQ=255)
                                      gCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCg                                             >  2:81728/1‑90 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        >  1:431691/1‑90 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        <  2:362845/90‑1 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        <  1:614410/90‑1 (MQ=255)
                                           ttttttGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACAtt                                        <  2:508286/90‑1 (MQ=255)
                                                    aGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAt                               >  1:773432/1‑90 (MQ=255)
                                                    aGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAt                               >  1:45114/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTc                     >  2:620000/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:594855/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:239391/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  2:8807/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:388260/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:436270/1‑90 (MQ=255)
                                                               cACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCt                    >  1:856747/1‑90 (MQ=255)
                                                                  tGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGca                 >  2:738893/1‑90 (MQ=255)
                                                                    aTAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGcaac               >  1:680691/1‑90 (MQ=255)
                                                                       gTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAg            >  2:399675/1‑90 (MQ=255)
                                                                       gTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAg            >  1:366880/1‑90 (MQ=255)
                                                                        tATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTc                                    <  2:424188/65‑1 (MQ=255)
                                                                        tATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTc                                    >  1:424188/1‑65 (MQ=255)
                                                                         aTCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGcc          <  2:12229/90‑1 (MQ=255)
                                                                           cAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCg         >  1:724258/1‑90 (MQ=255)
                                                                           cAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCgg        >  1:565378/1‑90 (MQ=255)
                                                                             gTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGtt      >  2:189686/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:504158/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:674431/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:717026/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  2:180715/1‑90 (MQ=255)
                                                                                 aTGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTaaa  >  1:157084/1‑90 (MQ=255)
                                                                                       |                                                                                    
GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA  >  NZ_CP009273/3356284‑3356455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ 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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAA  >  NZ_CP009273/3356284‑3356472
                                                                                              |                                                                                              
GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAC                                                                                          >  SRR3722072.541862/1‑100 (MQ=60)
GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAAC                                                                                          >  SRR3722072.467965/1‑100 (MQ=60)
  GTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACAT                                                                                        >  SRR3722072.750629/1‑100 (MQ=60)
     GATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGC                                                                                     >  SRR3722072.326592/1‑100 (MQ=60)
      ATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCT                                                                                    >  SRR3722072.256648/1‑100 (MQ=60)
       TCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTA                                                                                   >  SRR3722072.844476/1‑100 (MQ=60)
       TCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTA                                                                                   >  SRR3722072.395640/1‑100 (MQ=60)
         CTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACG                                                                                 >  SRR3722072.760000/1‑100 (MQ=60)
          TAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGC                                                                                <  SRR3722072.393196/100‑1 (MQ=60)
                 TCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGA                                                                         <  SRR3722072.527095/100‑1 (MQ=60)
                      TTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAA                                                                    >  SRR3722072.694253/1‑100 (MQ=60)
                        ATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAAT                                                                  <  SRR3722072.236756/100‑1 (MQ=60)
                                 GTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATT                                                         >  SRR3722072.439552/1‑100 (MQ=60)
                                          GTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAT                                                >  SRR3722072.45937/1‑100 (MQ=60)
                                          GTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGAT                                                >  SRR3722072.787018/1‑100 (MQ=60)
                                           TTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATC                                               <  SRR3722072.625365/100‑1 (MQ=60)
                                                     GGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCT                                     >  SRR3722072.871535/1‑100 (MQ=60)
                                                     GGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCT                                     >  SRR3722072.444219/1‑100 (MQ=60)
                                                     GGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCT                                     >  SRR3722072.243403/1‑100 (MQ=60)
                                                     GGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCT                                     >  SRR3722072.605466/1‑100 (MQ=60)
                                                     GGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCT                                     >  SRR3722072.395213/1‑100 (MQ=60)
                                                          ATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAAC                                >  SRR3722072.692869/1‑100 (MQ=60)
                                                             TTCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAG                             >  SRR3722072.373400/1‑100 (MQ=60)
                                                              TCACTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCctgtctcttatacac                            >  SRR3722072.431901/1‑85 (MQ=60)
                                                                 CTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCG                          >  SRR3722072.737122/1‑100 (MQ=60)
                                                                 CTGATAGTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGG                         >  SRR3722072.575510/1‑100 (MQ=60)
                                                                       GTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA                   >  SRR3722072.686508/1‑100 (MQ=60)
                                                                       GTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA                   >  SRR3722072.729785/1‑100 (MQ=60)
                                                                       GTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA                   >  SRR3722072.159834/1‑100 (MQ=60)
                                                                       GTATCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAA                   >  SRR3722072.513256/1‑100 (MQ=60)
                                                                          TCAGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAAT                >  SRR3722072.688426/1‑100 (MQ=60)
                                                                            AGTAAATGGCGAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTC              >  SRR3722072.567762/1‑100 (MQ=60)
                                                                                         AAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAA  >  SRR3722072.352951/1‑100 (MQ=60)
                                                                                              |                                                                                              
GTGTAGATCCTAATATATCGCATTATTTTTCCCGTTCAGCCCGTTTTTTGCCAGGTGAATATTCACTGATAGTATCAGTAAATGGCGAAAAAAAAGGCAACATTGCTACGCGCTTTGATGAAAATGGCGACATTTGTCTTGATCAGGCATTTCTGCAACAAGCCGGTTTAAAAATTCCTTCTGAAGAAA  >  NZ_CP009273/3356284‑3356472

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: