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
A4 F13 I1 R1 59 72.5 3655346 98.4% 3596860 137.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 1,419,549 +C intergenic (‑38/+98) PP_1242 ← / ← PP_1244 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC1,419,5481.C96.3% 92.3 / ‑3.4 27intergenic (‑37/+99)PP_1242/PP_1244hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (10/16);  total (10/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 = 1.00e+00

CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGC  >  NC_002947_CJ‑RC/1419401‑1419622
                                                                                                                                                     |                                                                          
cTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggc                                                                           <  5:198317/149‑1 (MQ=255)
   gCAGGGGATGCGGGGTTAGATGCGGTACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGTAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGTTTGTGGGTTGCCTGggctggc                                                                           <  2:258740/146‑1 (MQ=255)
   gCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGATTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGAGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctcgcctc                                                                        >  5:240156/1‑146 (MQ=255)
      ggggATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTc                                                                      >  5:256764/1‑149 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCgg                                                                  <  6:256764/149‑1 (MQ=255)
          aTGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCgg                                                                  <  3:371880/149‑1 (MQ=255)
                  ttCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGccc                                                          <  8:36342/149‑1 (MQ=255)
                           gttctgCGGGTGAGTGTTGTTTTGGGCTTGGTAGAAGCAGGTCTTGCGGACTGTGGTGGGTTTTGTTGTGTTGTGGGTGGGCTTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGGTGGGCTCTTCGCGGGCACGCCCGCTCCcag                                                  <  6:240156/143‑2 (MQ=255)
                                 cGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                           >  6:103893/1‑149 (MQ=255)
                                  gggTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACt                                           <  2:57289/148‑1 (MQ=255)
                                    gTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggc                                                                           <  8:214052/113‑1 (MQ=255)
                                    gTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggc                                                                           >  7:214052/1‑113 (MQ=255)
                                        gtgtTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGt                                    <  7:306951/149‑1 (MQ=255)
                                                   gCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCt                         >  4:313794/1‑149 (MQ=255)
                                                         gTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTcgc                                                                    >  4:23463/1‑100 (MQ=255)
                                                         gTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTcgc                                                                    <  3:23463/100‑1 (MQ=255)
                                                          tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                       >  1:379101/1‑144 (MQ=255)
                                                          tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgt                       <  2:379101/144‑1 (MQ=255)
                                                          tAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGa                  <  3:313794/149‑1 (MQ=255)
                                                            gAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAc                 <  5:103893/148‑1 (MQ=255)
                                                                cAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtgt            <  4:341855/149‑1 (MQ=255)
                                                                        gCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAc                 <  3:231596/136‑1 (MQ=255)
                                                                        gCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAc                 >  4:231596/1‑136 (MQ=255)
                                                                           gACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGt                                >  3:320497/1‑118 (MQ=255)
                                                                           gACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGt                                <  4:320497/118‑1 (MQ=255)
                                                                           gACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGATGGCTgtggtg                   <  1:132603/131‑1 (MQ=255)
                                                                           gACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGATGGCTgtggtg                   >  2:132603/1‑131 (MQ=255)
                                                                                       ttCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgtg                      >  7:153499/1‑116 (MQ=255)
                                                                                       ttCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCtgtg                      <  8:153499/116‑1 (MQ=255)
                                                                                       ttCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtggt                    <  4:191129/118‑1 (MQ=255)
                                                                                       ttCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTgtggt                    >  3:191129/1‑118 (MQ=255)
                                                                                                  tgtCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGc  <  3:335329/125‑1 (MQ=255)
                                                                                                                                                     |                                                                          
CTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGC  >  NC_002947_CJ‑RC/1419401‑1419622

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
CCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAG  >  NC_002947_CJ‑RC/1419400‑1419663
                                                                                                                                                      |                                                                                                                  
CCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCC                                                                                                                   <  NB501016:282:H2V7GAFX2:3:11607:3352:16373/151‑1 (MQ=60)
      AGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCG                                                                                                              >  NB501016:282:H2V7GAFX2:3:21404:23489:18935/1‑151 (MQ=60)
          GATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGATGCAGGTGTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGG                                                                                                          <  NB501016:282:H2V7GAFX2:2:21302:7727:16668/151‑1 (MQ=60)
                                    GGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCC                                                                                                                   >  NB501016:282:H2V7GAFX2:4:11609:1610:19725/1‑115 (MQ=60)
                                        AGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTG                                                                            <  NB501016:282:H2V7GAFX2:4:21411:12887:9563/151‑1 (MQ=60)
                                                         GGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACTGTGCAGTGCGTGACGGCTGTGG                                                              >  NB501016:282:H2V7GAFX2:1:21212:1317:3111/1‑148 (MQ=60)
                                                         GGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCG                                                                                                            <  NB501016:282:H2V7GAFX2:2:11104:23618:16721/102‑1 (MQ=60)
                                                          GTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTGGGTGGGCTTTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAC                                                          <  NB501016:282:H2V7GAFX2:2:21205:18738:4614/151‑1 (MQ=60)
                                                            AGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCccaaccccacaaactgtggggggcggggcgctaggaagagca                                                         >  NB501016:282:H2V7GAFX2:3:11506:21837:14425/1‑108 (MQ=60)
                                                            AGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACC                                                         <  NB501016:282:H2V7GAFX2:3:11505:10326:1684/150‑1 (MQ=60)
                                                                    GTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCcagatacgccggggcccccccccccacccaccccgatgatggagggttggcggcgtgggtggccggggg                                               >  NB501016:282:H2V7GAFX2:3:11507:18367:12012/1‑82 (MQ=60)
                                                                        TGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgtcaccacagccgtcacgcactgcacagtacctgtgggagcg                                              <  NB501016:282:H2V7GAFX2:1:21311:5269:17949/149‑43 (MQ=60)
                                                                        TGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACC                                                         <  NB501016:282:H2V7GAFX2:2:21101:16936:17660/138‑1 (MQ=60)
                                                                          CCGTCTGTGGGTGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGAATCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGgtcaccacagccgtcacgcactgcacagtacctgtgggagcgggcg                                          <  NB501016:282:H2V7GAFX2:1:11310:2402:16518/151‑47 (MQ=60)
                                                                          CCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGATGGCTGTGGTGAC                                                          <  NB501016:282:H2V7GAFX2:1:11209:6148:12506/135‑1 (MQ=60)
                                                                           CGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTG                                                                        >  NB501016:282:H2V7GAFX2:2:21206:8981:5946/1‑120 (MQ=60)
                                                                                 TGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCggggtaagggggggggccccccccccccccccacaaggaggggggggggtgggggggggggggggggggggggggggggggg                                  >  NB501016:282:H2V7GAFX2:4:11604:11376:2200/1‑69 (MQ=60)
                                                                                       GTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTG                                                            >  NB501016:282:H2V7GAFX2:2:11307:23453:8558/1‑120 (MQ=60)
                                                                                       GTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGG                                                              >  NB501016:282:H2V7GAFX2:4:11512:17238:8186/1‑118 (MQ=60)
                                                                                                  TTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCG                                          <  NB501016:282:H2V7GAFX2:2:21208:17557:12451/127‑1 (MQ=60)
                                                                                                      CGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCccaaagtcgcgagaccaccccaccccaccaatgtgtagatttgggatgggtctgggttggccaagccggccagcgggacggcgccagacgcggaggagctaa              >  NB501016:282:H2V7GAFX2:4:11609:5899:17465/1‑48 (MQ=60)
                                                                                                             CATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCcacttagagggcgccgcccccaccacaaaaattatgaaacaatgatgcccagactggttagaaagggggtaggccgtggagaccgaagcggcagccaaaaaagagaaa        >  NB501016:282:H2V7GAFX2:1:21109:18327:16075/1‑41 (MQ=60)
                                                                                                                 GTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCcacattgcgggggcggccccccccccaaagccgattgtgtggggggtggggggtttgtgggcgcgtggggggggggggggcccccgcatggcggggcgcggaataaaaagtggt  >  NB501016:282:H2V7GAFX2:3:21501:13671:2984/1‑37 (MQ=60)
                                                                                                                                                      |                                                                                                                  
CCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAG  >  NC_002947_CJ‑RC/1419400‑1419663

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