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 I176 R1 180 42.9 1710034 96.5% 1650182 119.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 918,237 T→A N14K (AAT→AAA yljB → protein YljB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273918,2370TA100.0% 143.1 / NA 37N14K (AAT→AAAyljBprotein YljB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (20/17);  total (20/17)

GGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGACAAT  >  NZ_CP009273/918106‑918370
                                                                                                                                   |                                                                                                                                     
ggggCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGa                                                                                                                                   >  2:681468/1‑136 (MQ=255)
ggggCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGa                                                                                                                                   >  2:617598/1‑136 (MQ=255)
       ttCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGc                                                                                                                         >  1:532594/1‑139 (MQ=255)
        tCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAagtg                                                                                                                                    >  2:50859/1‑127 (MQ=255)
        tCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAagtg                                                                                                                                    <  1:50859/127‑1 (MQ=255)
                 cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  2:92814/1‑139 (MQ=255)
                 cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  1:840704/1‑139 (MQ=255)
                 cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  1:478747/1‑139 (MQ=255)
                 cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  1:378567/1‑139 (MQ=255)
                 cGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGAGCAGATATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAAc                                                                                                               >  2:586166/1‑139 (MQ=255)
                            gAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAgatga                                                                                                      >  1:581301/1‑137 (MQ=255)
                               aTACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGtt                                                                                                 <  1:488846/139‑1 (MQ=255)
                               aTACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGtt                                                                                                 <  1:681468/139‑1 (MQ=255)
                                    aTTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatcta                                                                                            <  1:27791/139‑1 (MQ=255)
                                    aTTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTtatcta                                                                                            <  2:532594/139‑1 (MQ=255)
                                         gagaATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCAt                                                                                        <  1:586166/138‑1 (MQ=255)
                                             aTTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGa                                                                                   <  1:92814/139‑1 (MQ=255)
                                                    cGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGAtt                                                                            >  1:636630/1‑139 (MQ=255)
                                                      aGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCaa                                                                                                                >  2:655225/1‑101 (MQ=255)
                                                      aGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCaa                                                                                                                <  1:655225/101‑1 (MQ=255)
                                                           tCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAgg                                                                     <  2:770968/139‑1 (MQ=255)
                                                           tCAAGCATTTGACAGGGGACAAGAGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAgg                                                                     <  1:253722/139‑1 (MQ=255)
                                                                    tGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAg                                                             <  1:123475/138‑1 (MQ=255)
                                                                     gACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGct                                                           >  2:420132/1‑139 (MQ=255)
                                                                     gACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGct                                                           >  1:211310/1‑139 (MQ=255)
                                                                          gggACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGt                                                      >  2:640702/1‑139 (MQ=255)
                                                                          gggACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGt                                                      >  2:720423/1‑139 (MQ=255)
                                                                           ggACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGt                                                      >  1:233883/1‑138 (MQ=255)
                                                                                    gCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGa                                                                                   <  2:589395/100‑1 (MQ=255)
                                                                                    gCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGa                                                                                   >  1:589395/1‑100 (MQ=255)
                                                                                     cAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGa                                                                                                             >  2:116898/1‑73 (MQ=255)
                                                                                     cAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGa                                                                                                             <  1:116898/73‑1 (MQ=255)
                                                                                         tATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGtt                                       <  1:720423/139‑1 (MQ=255)
                                                                                         tATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGtt                                       <  1:617598/139‑1 (MQ=255)
                                                                                         tATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGtt                                       <  2:840704/139‑1 (MQ=255)
                                                                                             aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAg                                                                      >  2:532921/1‑104 (MQ=255)
                                                                                             aaTCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAg                                                                      <  1:532921/104‑1 (MQ=255)
                                                                                                                         ttGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAAc         <  2:233883/137‑1 (MQ=255)
                                                                                                                              aGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGACAAt  >  2:207435/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAACCGTAACCGAAGATGATAACTGACAAT  >  NZ_CP009273/918106‑918370

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAAC  >  NZ_CP009273/918101‑918344
                                                                                                                                        |                                                                                                           
agGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGC                                                                                               <  SRR3722130.51556/149‑1 (MQ=60)
 AGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCT                                                                                              >  SRR3722130.540728/1‑151 (MQ=60)
           GTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACG                                                                                    >  SRR3722130.852874/1‑151 (MQ=60)
           GTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACG                                                                                    >  SRR3722130.384579/1‑151 (MQ=60)
           GTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACG                                                                                    >  SRR3722130.486172/1‑151 (MQ=60)
                     gggctcggagatgtgtataagagacagGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGAT                                                                           <  SRR3722130.665025/123‑1 (MQ=60)
                      CGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGAT                                                                           >  SRR3722130.590092/1‑149 (MQ=60)
                       acggcatacgagattcgccttagtctcgtgggctcggagatgtgtataagaGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGT                                                                        <  SRR3722130.118556/100‑1 (MQ=60)
                              CTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATctgtctcttatacacatctccgagcccacgagact                                                                 >  SRR3722130.299021/1‑116 (MQ=60)
                                   AATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGT                                                            <  SRR3722130.691591/151‑1 (MQ=60)
                                   AATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGT                                                            <  SRR3722130.496419/151‑1 (MQ=60)
                                        AATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGT                                                       <  SRR3722130.28155/151‑1 (MQ=60)
                                             AGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGAT                                                   <  SRR3722130.595003/150‑1 (MQ=60)
                                              GAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTG                                                 >  SRR3722130.646154/1‑151 (MQ=60)
                                                 AATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTC                                              <  SRR3722130.94111/151‑1 (MQ=60)
                                                              gctcggagatgtgtataagagacAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCT                                 <  SRR3722130.541059/128‑1 (MQ=60)
                                                               GTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTC                                >  SRR3722130.214510/1‑151 (MQ=60)
                                                               GTCAAGCATTTGACAGGGGACAAGAGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTC                                <  SRR3722130.257694/151‑1 (MQ=60)
                                                                     CATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTC                           >  SRR3722130.237474/1‑150 (MQ=60)
                                                                        TTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGC                        <  SRR3722130.125246/150‑1 (MQ=60)
                                                                              GGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCctgtctcttatacacatctccgagcccac                 >  SRR3722130.598279/1‑122 (MQ=60)
                                                                                        ccttagtctcgtgggctcggagatgtgtataagagacagTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCG       <  SRR3722130.681629/112‑1 (MQ=60)
                                                                                             GTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAAC  <  SRR3722130.626876/151‑1 (MQ=60)
                                                                                             GTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAAAGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAAC  <  SRR3722130.731054/151‑1 (MQ=60)
                                                                                                                                        |                                                                                                           
GAGAGGGGGCTGTTCAAAACCTCGCCAACTCTAGAAATACAATTTAGAGAATTAGGGCGAGCCGTCAAGCATTTGACAGGGGACAAGGGGCAGGTATGAATCAAAAATTTGAAGCAGTTAACGCTATTGACAGGAATGTGACAGATGTCGCTGATGCCAACGATAGATGATAGTTATCTATCATGTGGAGTAGATTGGTCAGGCAAATAAGCTCTTGTCAGCGGCAGGGCGTTCTGCCGATAAC  >  NZ_CP009273/918101‑918344

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