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 F18 I1 R1 59 77.7 3515218 98.5% 3462489 145.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 1,126,647 +C intergenic (‑101/+48) tdcG‑II ← / ← gcvP‑I L‑serine dehydratase/glycine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC1,126,6451.C95.7% 76.6 / ‑3.0 23intergenic (‑99/+50)tdcG‑II/gcvP‑IL‑serine dehydratase/glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (15/7);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947_CJ‑RC/1126497‑1126788
                                                                                                                                                       |                                                                                                                                               
cTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTaaa                                                                                                                                                  <  8:328900/149‑1 (MQ=255)
   gAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAAc                                                                                                                                                <  6:299096/147‑1 (MQ=255)
      gggCCTTTGCCGATCTTTAACAGGTCGAAGACGCTCAGTTACATGTGGACCTCCTGCCCGGTTGTTGGGTGGACATTGAGCCCTGTTGGCGCGGGTTTACCCGCTAATGCGGTGTTATAGCCACCATCGCATTCGGGGTTAAAcccgct                                                                                                                                           <  2:405799/149‑5 (MQ=255)
                             gTCCAACACGCTCAGTCAATTGTCGACCTCCGGCCTGGTTGTGGGGGGCACTTTGCGCCCTTTCGGGGCGGTTTAACCCGCTATTGGGCCGTTATCCCCACCACCGCATCCGCGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCaa                                                                                                                     <  2:44871/149‑1 (MQ=255)
                                                                      ttgGGGGGACTTTCCGCCCTGTCGGCGGGGTTTTCCCGGCTATGGCGGCGGTATCCCCACCATCCAATGCGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGctc                                                                             <  6:162378/148‑1 (MQ=255)
                                                                        ggggggACATTCCCCCCTGTCGGCGGGGGTTAACCCGCGAATGGGTAGTTATACCCACCATCGCATTCGGGGTAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCtcg                                                                          <  8:31971/149‑1 (MQ=255)
                                                                                      cccTGTCGGCGGGGGTTACCCGCCGACTGCGGCGTATTACCCACCACCCCATGCGCGGTAAACCCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGGGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGgc                                                            <  2:182/149‑1 (MQ=255)
                                                                                                    gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCAGGTTGc                                              >  3:95097/1‑149 (MQ=255)
                                                                                                    gTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGc                                              >  3:158830/1‑149 (MQ=255)
                                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa                                                                                 <  3:20189/108‑1 (MQ=255)
                                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATa                                                                                 >  4:20189/1‑108 (MQ=255)
                                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCGGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                        >  3:152804/1‑149 (MQ=255)
                                                                                                          ccGCGAATGCGGAGGTATAGCCACCATCGCATTCCCCGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGCCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGc                                        >  3:146519/1‑149 (MQ=255)
                                                                                                                     ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGtt                                                <  1:385874/131‑1 (MQ=255)
                                                                                                                     ggAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGtt                                                >  2:385874/1‑131 (MQ=255)
                                                                                                                            taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                       >  5:86466/1‑149 (MQ=255)
                                                                                                                            taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                       >  3:183224/1‑149 (MQ=255)
                                                                                                                            taGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAcc                       >  7:301691/1‑149 (MQ=255)
                                                                                                                              gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                     >  3:378902/1‑149 (MQ=255)
                                                                                                                              gCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCGCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCg                     >  4:371562/1‑149 (MQ=255)
                                                                                                                               ccaccaTCCCATTCGGGGTTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCgg                    <  2:152317/149‑1 (MQ=255)
                                                                                                                                 accaTCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGACGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGcc                  >  1:154736/1‑149 (MQ=255)
                                                                                                                                        cATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGCCcggc           >  4:58349/1‑149 (MQ=255)
                                                                                                                                         aTTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACcggcg          >  8:115429/1‑149 (MQ=255)
                                                                                                                                              cgGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAg     <  3:371562/149‑1 (MQ=255)
                                                                                                                                                 gTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAc  >  3:348963/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
CTGGAGGGGCCGATGCCGATCTTGAACAGGTCGAAGACGCTCAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTAC  >  NC_002947_CJ‑RC/1126497‑1126788

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCC  >  NC_002947_CJ‑RC/1126538‑1126793
                                                                                                              |                                                                                                                                                  
agctcacgggggcccccccgccggggtttggggggcaaactccccccgggtcgggggggttcccccgcccacgcggctgtaacccccaccccccactgggggtaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGG                                                                                                            <  NB501016:282:H2V7GAFX2:2:21306:24409:7452/42‑1 (MQ=60)
  gccactgcgccccccccgccgtgttgggggggcccaatcgcccccgggggggggggtttcccccggctttggggtttaaaccccccccccccctggggggtaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGG                                                                                                            <  NB501016:282:H2V7GAFX2:2:11209:8560:14804/42‑1 (MQ=60)
  cgccgggggccccccccgccttgtgggggggggaaccccgcccccgtgggggggggggtcccccggaccgggggctttacccccccccccccctgggggttaaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAG                                                                                                          <  NB501016:282:H2V7GAFX2:4:21409:25097:4136/44‑1 (MQ=60)
  ccctccgcccccccccccccgggggggggggcccccccccccccggggcgggggggtcccccccgccccgggccccacaccccccccccccctcgcgccccaaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAG                                                                                                          <  NB501016:282:H2V7GAFX2:4:21602:20735:7854/44‑1 (MQ=60)
            caccccctccctgttggtgggggccaattccccccgtggcgggggggtataccccccaatggcgctgtaacccccaccagcccttcgggtttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGG                                                                                                <  NB501016:282:H2V7GAFX2:4:21607:24230:9792/54‑1 (MQ=60)
                                    catcgccccccgtgtggcggggggccccccggctccgcgcgcgtacccccccccccccctgcgggggaaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACG                                                                        <  NB501016:282:H2V7GAFX2:4:11604:22569:18796/78‑1 (MQ=60)
                                          cccccgtgtcggggggggtcacccggctattgggtcttaacccccccaaccccttgggggttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGG                                                                  <  NB501016:282:H2V7GAFX2:3:21510:2076:12945/84‑1 (MQ=60)
                                            cccctgtcgggggggtttccccggctttgggttttataaccccccccccgattgggggtaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCG                                                                <  NB501016:282:H2V7GAFX2:4:11403:4404:13378/86‑1 (MQ=60)
                                                ggcgtggtgggttcaccctgtcatggggggtcaaacccacccctctctttgggttaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCAC                                                             <  NB501016:282:H2V7GAFX2:1:11305:26122:19335/89‑1 (MQ=60)
                                                    gcgggggcctccgcgtttgtgtctgtcattccccccgccgcttctgggttcaaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAG                                                        <  NB501016:282:H2V7GAFX2:4:21607:15496:1181/94‑1 (MQ=60)
                                                         ggtttccccggttattgcgtctttaacccacccacctattggggtttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGC                                                   <  NB501016:282:H2V7GAFX2:3:21605:5550:2434/99‑1 (MQ=60)
                                                          GGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGCCGGGCAGGCGCACACCAGGTTGCG                                                  >  NB501016:282:H2V7GAFX2:2:11204:12298:16901/1‑151 (MQ=60)
                                                          GGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCG                                                  >  NB501016:282:H2V7GAFX2:2:11303:15539:1732/1‑151 (MQ=60)
                                                             taccccgcattggggttttaaacccccccacccattggggtttaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGT                                                <  NB501016:282:H2V7GAFX2:2:11312:20026:2571/102‑1 (MQ=60)
                                                                CCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAG                                                                                     <  NB501016:282:H2V7GAFX2:2:11104:15230:9335/110‑1 (MQ=60)
                                                                CCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCGGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCC                                            >  NB501016:282:H2V7GAFX2:2:11302:26667:2891/1‑151 (MQ=60)
                                                                CCCGCGAATGCGGAGGTATAGCCACCATCGCATTCCCCGGAAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGCCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCC                                            >  NB501016:282:H2V7GAFX2:2:11301:1222:2908/1‑151 (MQ=60)
                                                                          gttgttaaacccacctccctttggggtgtaaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGTGGAGGGCTATCAGGCGTCCTTATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTG                                  <  NB501016:282:H2V7GAFX2:2:21307:24979:13088/116‑1 (MQ=60)
                                                                          CGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTG                                                    <  NB501016:282:H2V7GAFX2:1:21304:26362:5712/133‑1 (MQ=60)
                                                                               tttaagcaaccatcgtttgcgggttaaaccgCGCTCCTACATAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGA                              <  NB501016:282:H2V7GAFX2:1:11205:22883:17687/119‑1 (MQ=60)
                                                                                tacccccccatagccaggggtgttaacccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGAC                             <  NB501016:282:H2V7GAFX2:4:21410:12244:15923/121‑1 (MQ=60)
                                                                                aaccccccccccccctggggtttcaaaccCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACC                            <  NB501016:282:H2V7GAFX2:2:21203:16258:12025/122‑1 (MQ=60)
                                                                                 ATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCC                           >  NB501016:282:H2V7GAFX2:2:11306:1742:16938/1‑151 (MQ=60)
                                                                                 ATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCC                           >  NB501016:282:H2V7GAFX2:3:11503:2244:7539/1‑151 (MQ=60)
                                                                                 ATAGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCC                           >  NB501016:282:H2V7GAFX2:4:21412:12904:2508/1‑151 (MQ=60)
                                                                                   AGCCACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGG                         >  NB501016:282:H2V7GAFX2:2:21305:10536:13327/1‑151 (MQ=60)
                                                                                      AACCATCGCATTCGCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGACGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCC                      >  NB501016:282:H2V7GAFX2:1:11301:13783:7579/1‑151 (MQ=60)
                                                                                                   GCGGGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGT         <  NB501016:282:H2V7GAFX2:2:21304:19898:9245/151‑1 (MQ=60)
                                                                                                      GGTAAACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACT      >  NB501016:282:H2V7GAFX2:2:21212:18842:4528/1‑151 (MQ=60)
                                                                                                         AACCCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCC  >  NB501016:282:H2V7GAFX2:4:11509:15413:17679/1‑151 (MQ=60)
                                                                                                              |                                                                                                                                                  
CAGTGACATGTGGACCTCCTGCCAGGTTGTTGGGTGGACATTGAGCCCTGTAGGCGCGGGTTTACCCGCGAATGCGGAGGTATAGCCACCATCGCATTCGCGGGTAAACCGCTCCTACACAGAGCCCGCCACCCAAGCCAGGTGGCGGGGAGGGCTATCAGGCGTCCTGATAGCTCTCGATCGACGGGCAGGCGCACACCAGGTTGCGGTCGCCGAACACGTTGTCGACCCGGCCGACCGGCGGCCAGTACTTGCC  >  NC_002947_CJ‑RC/1126538‑1126793

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