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 F33 I1 R1 23 35.1 1378992 95.0% 1310042 124.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,184,543 C→A P1100Q (CCG→CAG)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,184,5430CA100.0% 154.0 / NA 44P1100Q (CCG→CAG) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/23);  total (21/23)

GGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGC  >  NC_000913/4184397‑4184687
                                                                                                                                                  |                                                                                                                                                
ggTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGc                                                                                                                                                >  8:78387/1‑149 (MQ=255)
             gTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                                   <  5:109669/149‑1 (MQ=255)
                    gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                                 <  5:169440/144‑1 (MQ=255)
                    gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                                 >  6:169440/1‑144 (MQ=255)
                               ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                                   >  6:39678/1‑147 (MQ=255)
                                cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                                   >  5:135397/1‑146 (MQ=255)
                                      aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                          <  7:78387/149‑1 (MQ=255)
                                        aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                                   >  8:128424/1‑122 (MQ=255)
                                        aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                                   <  7:128424/122‑1 (MQ=255)
                                          gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                          >  8:79189/1‑145 (MQ=255)
                                          gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                          <  7:79189/145‑1 (MQ=255)
                                                        aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAt                                                                                                                       >  7:43454/1‑118 (MQ=255)
                                                        aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAt                                                                                                                       <  8:43454/118‑1 (MQ=255)
                                                                aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGa                                                                                 >  2:79770/1‑148 (MQ=255)
                                                                  gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTc                                                                                                          <  7:112383/121‑1 (MQ=255)
                                                                  gggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTc                                                                                                          >  8:112383/1‑121 (MQ=255)
                                                                       aaaTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                             <  6:39217/80‑1 (MQ=255)
                                                                          tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                             <  3:7632/126‑1 (MQ=255)
                                                                          tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                             >  4:7632/1‑126 (MQ=255)
                                                                          tttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGc                                                                      >  8:13059/1‑149 (MQ=255)
                                                                                  aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                 <  5:135787/114‑1 (MQ=255)
                                                                                  aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAc                                                                                                 >  6:135787/1‑114 (MQ=255)
                                                                                  aTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                              <  5:39678/149‑1 (MQ=255)
                                                                                     aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCaa                                                           >  8:96807/1‑149 (MQ=255)
                                                                                      aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                          <  4:2613/149‑1 (MQ=255)
                                                                                      aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                          >  3:2613/1‑149 (MQ=255)
                                                                                          gATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                                      <  7:13059/149‑1 (MQ=255)
                                                                                                        cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACATGATCAACGCCATGCTGAAAcagcag                                        <  1:79770/149‑1 (MQ=255)
                                                                                                        cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAAc                                                          <  8:105702/131‑1 (MQ=255)
                                                                                                            cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                                   >  1:30067/1‑70 (MQ=255)
                                                                                                            cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                                   <  2:30067/70‑1 (MQ=255)
                                                                                                             gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                                    >  6:166919/1‑148 (MQ=255)
                                                                                                                               gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                                >  3:11003/1‑86 (MQ=255)
                                                                                                                               gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                                <  4:11003/86‑1 (MQ=255)
                                                                                                                                 aGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                                    <  2:110374/128‑1 (MQ=255)
                                                                                                                                 aGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                                    >  1:110374/1‑128 (MQ=255)
                                                                                                                                      tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATAGGTCAGATCCTCGAAACCAACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                       >  5:150824/1‑120 (MQ=255)
                                                                                                                                           cTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgc      <  5:166919/148‑1 (MQ=255)
                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt               <  2:52301/136‑1 (MQ=255)
                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt               >  1:52301/1‑136 (MQ=255)
                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcgtgc  <  6:135397/149‑1 (MQ=255)
                                                                                                                                              aaCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAgcgtgc  <  7:96807/149‑1 (MQ=255)
                                                                                                                                               aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                          <  4:133304/108‑1 (MQ=255)
                                                                                                                                               aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCAATGCTGAAAcagc                                          >  3:133304/1‑108 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GGTATATCTGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGC  >  NC_000913/4184397‑4184687

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ 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
GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCG  >  NC_000913/4184409‑4184688
                                                                                                                                      |                                                                                                                                                 
GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCT                                                                                                                                   <  NB501016:228:HFK2FAFXY:3:21408:18939:3482/151‑1 (MQ=60)
       CGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCG                                                                                                                                 <  NB501016:228:HFK2FAFXY:3:21608:10011:10702/146‑1 (MQ=60)
                   CCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTC                                                                                                                   >  NB501016:228:HFK2FAFXY:3:21506:9210:14164/1‑148 (MQ=60)
                         GACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCG                                                                                                          <  NB501016:228:HFK2FAFXY:4:11612:4854:17124/151‑1 (MQ=60)
                           CAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCT                                                                                                                                   <  NB501016:228:HFK2FAFXY:4:21511:23757:18066/124‑1 (MQ=60)
                             AGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCG                                                                                                          <  NB501016:228:HFK2FAFXY:4:21401:12394:5570/147‑1 (MQ=60)
                                           CACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATC                                                                                                                       >  NB501016:228:HFK2FAFXY:4:11508:1831:3934/1‑120 (MQ=60)
                                                     AGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTCG                                                                                                          <  NB501016:228:HFK2FAFXY:4:21503:21959:19028/123‑1 (MQ=60)
                                                          GAAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGagatcggaagagcacacgtctgaacccagtcactaaagcgcatctcgtaagccgtcttctgcttgaa                                                                          >  NB501016:228:HFK2FAFXY:3:11504:26377:9345/1‑83 (MQ=60)
                                                             ATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGG                                                                                             <  NB501016:228:HFK2FAFXY:2:11104:20488:11358/128‑1 (MQ=60)
                                                                     GATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATC                                                              <  NB501016:228:HFK2FAFXY:3:11504:17253:12890/151‑1 (MQ=60)
                                                                     GATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACC                                                                                                 <  NB501016:228:HFK2FAFXY:3:21506:5182:17350/116‑1 (MQ=60)
                                                                         AACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACG                                                          >  NB501016:228:HFK2FAFXY:2:11102:19170:5347/1‑151 (MQ=60)
                                                                             CGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCAT                                                      <  NB501016:228:HFK2FAFXY:4:11406:24272:14690/151‑1 (MQ=60)
                                                                                           CCTTACGATGATAACGGTACGCCGGTAGACATCGTACTGAACCAGGTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACATGATCAACGCCATGCTGAAACAGCAGC                                        <  NB501016:228:HFK2FAFXY:1:11305:2429:20205/151‑1 (MQ=60)
                                                                                           CCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAACGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGagatcggaaggagcacac                                        >  NB501016:228:HFK2FAFXY:4:21412:6990:14879/1‑133 (MQ=60)
                                                                                               ACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTC                                                                                                                   >  NB501016:228:HFK2FAFXY:1:11111:10706:13931/1‑72 (MQ=60)
                                                                                                       AACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGT                                                                                                                                     <  NB501016:228:HFK2FAFXY:3:11501:6139:12744/46‑1 (MQ=60)
                                                                                                                  GGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAA                                                                                >  NB501016:228:HFK2FAFXY:2:11106:2032:4324/1‑88 (MQ=60)
                                                                                                                    TAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGA                                    >  NB501016:228:HFK2FAFXY:1:21105:12975:2643/1‑130 (MQ=60)
                                                                                                                         ATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATAGGTCAGATCCTCGAAACCAACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCA                                       >  NB501016:228:HFK2FAFXY:3:21601:24533:2078/1‑122 (MQ=60)
                                                                                                                              ACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCG      <  NB501016:228:HFK2FAFXY:3:21607:24154:9930/150‑1 (MQ=60)
                                                                                                                                 GAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCG  <  NB501016:228:HFK2FAFXY:4:21409:7190:2986/151‑1 (MQ=60)
                                                                                                                                 GAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATT               >  NB501016:228:HFK2FAFXY:1:11207:15333:17839/1‑138 (MQ=60)
                                                                                                                                  AACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCAATGCTGAAACAGCA                                          >  NB501016:228:HFK2FAFXY:2:21302:25954:18411/1‑110 (MQ=60)
                                                                                                                                      AGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGA                                                 <  NB501016:228:HFK2FAFXY:1:11311:23750:19004/99‑1 (MQ=60)
                                                                                                                                      |                                                                                                                                                 
GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAGCGTGCG  >  NC_000913/4184409‑4184688

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