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 557,194 A→G Q681R (CAG→CGG)  sfmD → fimbrial biogenesis usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273557,1940AG97.6% 133.9 / ‑3.8 41Q681R (CAG→CGG) sfmDfimbrial biogenesis usher protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/17);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

CACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCA  >  NZ_CP009273/557067‑557308
                                                                                                                               |                                                                                                                  
cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCg                                                                                                         >  2:691071/1‑139 (MQ=255)
       gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGt                                                                                                  <  2:528103/139‑1 (MQ=255)
           gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCggg                                                                                              <  2:502550/139‑1 (MQ=255)
                aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATg                                                                                          >  1:537553/1‑138 (MQ=255)
                 tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGa                                                                                         >  1:415682/1‑138 (MQ=255)
                  aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTc                                                                                                 >  1:503487/1‑129 (MQ=255)
                  aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTc                                                                                                 <  2:503487/129‑1 (MQ=255)
                                    acGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagcc                                                                     <  2:138919/139‑1 (MQ=255)
                                    acGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagcc                                                                     <  2:537553/139‑1 (MQ=255)
                                         cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCt                                                                >  2:481921/1‑139 (MQ=255)
                                           gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGaaaa                                                                                      <  2:549209/115‑1 (MQ=255)
                                           gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGaaaa                                                                                      >  1:549209/1‑115 (MQ=255)
                                           gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagcc                                                                     >  1:141452/1‑132 (MQ=255)
                                           gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagcc                                                                     <  2:141452/132‑1 (MQ=255)
                                           gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCtt                                                               >  1:374433/1‑138 (MQ=255)
                                              ccACGGCAAACTGGCAGGCCGCAAACGGTACGATGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGATAACTGGCGGCTGTCTGGCGGGGCGGTCGGGCATGAAAAAGGCATAACGCTGAGCCAGCCTTTAgg                                                           >  1:664076/1‑139 (MQ=255)
                                                 cGGCAAACTGGCAGGCCGCATACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCGGGCTGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAgggta                                                        <  2:150327/139‑3 (MQ=255)
                                                           gCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa                                                                                 >  2:699562/1‑104 (MQ=255)
                                                           gCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa                                                                                 <  1:699562/104‑1 (MQ=255)
                                                               gCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGa                                          <  1:149927/139‑1 (MQ=255)
                                                                   cTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGt                                               >  1:608276/1‑130 (MQ=255)
                                                                   cTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGt                                               <  2:608276/130‑1 (MQ=255)
                                                                      aCGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTaaa                                     <  2:397612/137‑1 (MQ=255)
                                                                           aCGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTg                              >  1:257500/1‑139 (MQ=255)
                                                                             gCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGt                                                                                                  >  2:261175/1‑69 (MQ=255)
                                                                             gCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGt                                                                                                  <  1:261175/69‑1 (MQ=255)
                                                                              cTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGACGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATGAAAGCGCCTGgcg                           <  2:415682/139‑1 (MQ=255)
                                                                                 gggggCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGgcgc                          <  1:481921/137‑1 (MQ=255)
                                                                                 gggggCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGgcgc                          <  1:290086/137‑1 (MQ=255)
                                                                                    ggCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCg                     >  1:642169/1‑139 (MQ=255)
                                                                                        ggTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGgtg                  >  1:594705/1‑138 (MQ=255)
                                                                                                  cGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGa        >  1:419759/1‑138 (MQ=255)
                                                                                                    aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa     >  2:746115/1‑139 (MQ=255)
                                                                                                    aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa     >  1:244238/1‑139 (MQ=255)
                                                                                                    aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa     >  2:626179/1‑139 (MQ=255)
                                                                                                    aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa     >  2:509865/1‑139 (MQ=255)
                                                                                                    aCCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGaaaa     >  2:504848/1‑139 (MQ=255)
                                                                                                        cgATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCa  >  2:615555/1‑138 (MQ=255)
                                                                                                        cgATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCa  >  2:400672/1‑138 (MQ=255)
                                                                                                                   aCGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGgcgc                          >  1:202583/1‑103 (MQ=255)
                                                                                                                   aCGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGgcgc                          <  2:202583/103‑1 (MQ=255)
                                                                                                                               |                                                                                                                  
CACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCA  >  NZ_CP009273/557067‑557308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ 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
GAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTA  >  NZ_CP009273/557046‑557321
                                                                                                                                                    |                                                                                                                               
gcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacaGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGC                                                                                                                               <  SRR3722133.567072/99‑1 (MQ=60)
                          AAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGA                                                                                                      >  SRR3722133.544733/1‑150 (MQ=60)
                           gaagacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagacagGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAA                                                                                                    <  SRR3722133.264744/91‑1 (MQ=60)
                           AGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAA                                                                                                     >  SRR3722133.421426/1‑150 (MQ=60)
                            GCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAA                                                                                                    >  SRR3722133.510271/1‑150 (MQ=60)
                                            gctcggagatgtgtataagagacagACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC                                                                                   <  SRR3722133.708618/126‑1 (MQ=60)
                                                     GTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCctgtctcttatac                                                                           >  SRR3722133.556554/1‑137 (MQ=60)
                                                     GTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTT                                                                           >  SRR3722133.143313/1‑150 (MQ=60)
                                                     GTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTT                                                                           >  SRR3722133.379666/1‑150 (MQ=60)
                                                        CACGGGCAGCGCCACGGCAAACTGGCAGGCCGCAAACGGTACGATGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGATAACTGGCGGCTGTCTGGCGGGGCGGTCGGGCATGAAAAAGGCATAACGCTGAGCCAGCCTTTAGGG                                                                       >  SRR3722133.672734/1‑151 (MQ=60)
                                                                   ttttttgttttttttttcaagcagaagacggcatacgagatctagtacggtctcgtgggctcggagatgtgtataagagaCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATG                                                             <  SRR3722133.296405/70‑1 (MQ=60)
                                                                             CTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAA                                                   >  SRR3722133.616311/1‑150 (MQ=60)
                                                                                   GGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCT                                            <  SRR3722133.151893/151‑1 (MQ=60)
                                                                                     CCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGG                                          >  SRR3722133.261011/1‑151 (MQ=60)
                                                                                              GTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGG                                 >  SRR3722133.650565/1‑151 (MQ=60)
                                                                                                  GCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGT                              >  SRR3722133.602557/1‑150 (MQ=60)
                                                                                                     GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTAC                            <  SRR3722133.488440/149‑1 (MQ=60)
                                                                                                     GGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTAC                            <  SRR3722133.294030/149‑1 (MQ=60)
                                                                                                            GGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAA                    >  SRR3722133.425547/1‑150 (MQ=60)
                                                                                                              GTATAACTACGACCGCGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAAT                 >  SRR3722133.247543/1‑151 (MQ=60)
                                                                                                                             CGATCAACATGACGTTAACTGGCGGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACctgtctcttatacacatctccgagcc  >  SRR3722133.205281/1‑125 (MQ=60)
                                                                                                                                                    |                                                                                                                               
GAAGGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGCCAGCCTTTAGGGGATACCAATGTTTTGATTAAAGCGCCTGGCGCAGGCGGTGTACGCATTGAAAATCAAACTGGCATTTTA  >  NZ_CP009273/557046‑557321

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