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
A5 F13 I0 R1 213 60.4 2554362 98.6% 2518600 147.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 2,958,523 C→T 100% A428V (GCA→GTA)  PP_2589 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC2,958,5230CT100.0% 101.0 / NA 30A428V (GCA→GTA) PP_2589aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/16);  total (14/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTG  >  NC_002947_CJ‑RC/2958377‑2958660
                                                                                                                                                  |                                                                                                                                         
cTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAg                                                                                                                                         >  5:207395/1‑149 (MQ=255)
   cATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGc                                                                                                                                        <  8:102608/147‑1 (MQ=255)
         gCCGAGCATTTTTGCCGGTGTCGAGGGGCGGATGATCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCTGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGgtgt                                                                                                                                <  1:209137/149‑1 (MQ=255)
               cATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGcg                                                                                                                          <  5:278747/149‑1 (MQ=255)
                 tttttGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGcgac                                                                                                                        >  5:244804/1‑149 (MQ=255)
                        cGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTc                                                                                                                  >  6:249492/1‑148 (MQ=255)
                        cGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTCGCCAACGACAGCATCTATGGCCTGGTAGCCTCTGAGTGGAGCGACGACTTc                                                                                                                  >  1:132219/1‑148 (MQ=255)
                               gAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGt                                                                                                          >  6:276812/1‑149 (MQ=255)
                                             ccTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGACGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAgcgc                                                                                            <  3:87117/149‑1 (MQ=255)
                                                        aggaggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGaaa                                                                                 <  2:44444/149‑1 (MQ=255)
                                                         ggaggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAg                                                                                >  8:200201/1‑149 (MQ=255)
                                                            ggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCg                                                                             <  6:86951/149‑1 (MQ=255)
                                                            ggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCg                                                                             <  2:284731/149‑1 (MQ=255)
                                                                     cccGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTAt                                                                    <  5:91136/149‑1 (MQ=255)
                                                                         gTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGt                                                                <  5:239148/149‑1 (MQ=255)
                                                                            cTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGaa                                                             <  3:230644/149‑1 (MQ=255)
                                                                            cTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGa                                                              <  3:222646/148‑1 (MQ=255)
                                                                               gCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAaca                                                           >  1:103876/1‑148 (MQ=255)
                                                                                      gcgcCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGAt                                                   >  3:242731/1‑149 (MQ=255)
                                                                                         ccTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATgcg                                                >  3:178014/1‑149 (MQ=255)
                                                                                               acacGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGAt                                          >  8:279368/1‑149 (MQ=255)
                                                                                                   ggAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCTTCTTTGGCCTGGTAGCCTCGGTGTGGAGCGAAGACTTCAACCAGGTGCAGGGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGa                                      <  4:263813/149‑1 (MQ=255)
                                                                                                                   cgcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTg                       <  7:200201/148‑1 (MQ=255)
                                                                                                                   cgcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTg                       <  3:264109/148‑1 (MQ=255)
                                                                                                                    gcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAATGCCGGTACGGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTggg                     <  5:229840/149‑1 (MQ=255)
                                                                                                                    gcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTggg                     >  3:76308/1‑149 (MQ=255)
                                                                                                                    gcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTggg                     >  3:71823/1‑149 (MQ=255)
                                                                                                                    gcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTggg                     >  3:58163/1‑149 (MQ=255)
                                                                                                                         gCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCaa                <  6:207395/149‑1 (MQ=255)
                                                                                                                                       cTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTg  >  5:239864/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                         
CTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTG  >  NC_002947_CJ‑RC/2958377‑2958660

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGG  >  NC_002947_CJ‑RC/2958376‑2958661
                                                                                                                                                   |                                                                                                                                          
ACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGC                                                                                                                                         >  NB501016:282:H2V7GAFX2:3:21411:1633:7827/1‑151 (MQ=60)
         AGCCGAGCATTTTTGCCGGTGTCGAGGGGCGGATGATCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCTGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTG                                                                                                                                <  NB501016:282:H2V7GAFX2:1:21110:10152:14862/151‑1 (MQ=60)
               CCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGA                                                                                                                          <  NB501016:282:H2V7GAFX2:3:21603:25650:12140/151‑1 (MQ=60)
                 TTTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACG                                                                                                                        >  NB501016:282:H2V7GAFX2:3:21508:24288:2642/1‑151 (MQ=60)
                        CCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTCGCCAACGACAGCATCTATGGCCTGGTAGCCTCTGAGTGGAGCGACGACTTCA                                                                                                                  >  NB501016:282:H2V7GAFX2:1:11305:4132:9360/1‑150 (MQ=60)
                                             GCCTGGCGCGCGAGGAGGGGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGACGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCC                                                                                            <  NB501016:282:H2V7GAFX2:2:11208:14579:17123/151‑1 (MQ=60)
                                                                     GCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATC                                                                    <  NB501016:282:H2V7GAFX2:3:11509:16232:8025/151‑1 (MQ=60)
                                                                         GGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTG                                                                <  NB501016:282:H2V7GAFX2:3:21506:23255:16193/151‑1 (MQ=60)
                                                                            ACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAAC                                                             <  NB501016:282:H2V7GAFX2:2:21207:5488:2260/151‑1 (MQ=60)
                                                                            ACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAA                                                              <  NB501016:282:H2V7GAFX2:2:21205:20340:6121/150‑1 (MQ=60)
                                                                               GGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACAC                                                           >  NB501016:282:H2V7GAFX2:1:11211:16558:2095/1‑150 (MQ=60)
                                                                                      AGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATG                                                   >  NB501016:282:H2V7GAFX2:2:21209:4896:16381/1‑151 (MQ=60)
                                                                                       ccgctgtatcacgaacgttttctGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGC                                                  <  NB501016:282:H2V7GAFX2:4:21402:3522:6306/128‑1 (MQ=60)
                                                                                         GCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGC                                                >  NB501016:282:H2V7GAFX2:2:21104:17682:15136/1‑151 (MQ=60)
                                                                                                    ggcagcagcagccgtgcgcctgggctgcgagagggTCTATGGGCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGCTGCACCGGGTAGCGCGGGGGTGGAAAGCCGGTACCGAATCGGTGACCACAGTCGAGGCGCTGGATGTGACG                                     >  NB501016:282:H2V7GAFX2:1:11211:25937:4120/36‑151 (MQ=60)
                                                                                                                   GCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGG                       <  NB501016:282:H2V7GAFX2:2:21303:7219:5599/150‑1 (MQ=60)
                                                                                                                   GCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGG                       <  NB501016:282:H2V7GAFX2:4:21408:22525:17835/150‑1 (MQ=60)
                                                                                                                    GGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAATGCCGGTACGGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGG                     <  NB501016:282:H2V7GAFX2:3:21504:23502:12566/151‑1 (MQ=60)
                                                                                                                    CGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGG                     >  NB501016:282:H2V7GAFX2:2:11206:26694:9512/1‑151 (MQ=60)
                                                                                                                    CGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGG                     >  NB501016:282:H2V7GAFX2:2:11205:1277:9518/1‑151 (MQ=60)
                                                                                                                    CGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGG                     >  NB501016:282:H2V7GAFX2:2:11202:2012:7993/1‑151 (MQ=60)
                                                                                                                                       TCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGG  >  NB501016:282:H2V7GAFX2:3:21506:16612:19031/1‑151 (MQ=60)
                                                                                                                                                   |                                                                                                                                          
ACTACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGG  >  NC_002947_CJ‑RC/2958376‑2958661

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