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
A1 F4 I88 R1 217 30.5 1742422 95.9% 1670982 87.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,545,596 T→C E181G (GAG→GGG)  yjiC ← YjfZ/YjiC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,545,5960TC100.0% 196.7 / NA 51E181G (GAG→GGG) yjiCYjfZ/YjiC family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (32/19);  total (32/19)

ATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAG  >  NZ_CP009273/4545513‑4545683
                                                                                   |                                                                                       
aTGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGt                                                                                  >  1:366214/1‑91 (MQ=255)
aTGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGt                                                                                  >  2:293646/1‑91 (MQ=255)
                  cGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTgaga                                                                <  1:831264/91‑1 (MQ=255)
                    gaagaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGaa                                                               >  1:109177/1‑90 (MQ=255)
                    gaagaaTTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGaa                                                               <  2:109177/90‑1 (MQ=255)
                          tttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGaa                                                        <  1:499907/91‑1 (MQ=255)
                          tttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGaa                                                        <  2:449238/91‑1 (MQ=255)
                           ttGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAg                                                       <  1:344690/91‑1 (MQ=255)
                              cTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGt                                                    <  1:207097/91‑1 (MQ=255)
                                 ccccTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTg                                                  <  1:387230/90‑1 (MQ=255)
                                 ccccTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTg                                                  >  2:387230/1‑90 (MQ=255)
                                        tACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAatatat                                          <  1:865921/91‑1 (MQ=255)
                                           cGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCt                                       <  1:585684/91‑1 (MQ=255)
                                            gTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTa                                      >  2:335496/1‑91 (MQ=255)
                                              cATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATc                                    >  2:351280/1‑91 (MQ=255)
                                                ttATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGt                                  >  2:546395/1‑91 (MQ=255)
                                                     tGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTc                             >  2:243240/1‑91 (MQ=255)
                                                     tGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTc                             <  2:659160/91‑1 (MQ=255)
                                                      gATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGaa                                                        >  1:466583/1‑63 (MQ=255)
                                                      gATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGaa                                                        <  2:466583/63‑1 (MQ=255)
                                                      gATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTcc                            >  2:248256/1‑91 (MQ=255)
                                                      gATCAATATAGATATCAATCCTTGGCTGCCCAGTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTcc                            >  1:861783/1‑91 (MQ=255)
                                                        tCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCgg                          >  1:601628/1‑91 (MQ=255)
                                                           atatAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCAt                       <  1:655089/91‑1 (MQ=255)
                                                           atatAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCAt                       <  1:358195/91‑1 (MQ=255)
                                                                gATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                           <  2:84622/66‑1 (MQ=255)
                                                                gATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAatata                                           >  1:84622/1‑66 (MQ=255)
                                                                gATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCa                  >  2:816320/1‑91 (MQ=255)
                                                                gATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCa                  >  2:297679/1‑91 (MQ=255)
                                                                gATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCa                  >  1:228812/1‑91 (MQ=255)
                                                                gATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCa                  >  1:646465/1‑91 (MQ=255)
                                                                gATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCa                  >  2:620486/1‑91 (MQ=255)
                                                                  tatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                >  1:802863/1‑91 (MQ=255)
                                                                  tatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                >  1:769806/1‑91 (MQ=255)
                                                                  tatCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAgg                >  1:697994/1‑91 (MQ=255)
                                                                   atCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGa               >  1:128073/1‑91 (MQ=255)
                                                                   atCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGa               >  2:678017/1‑91 (MQ=255)
                                                                   atCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGa               >  2:459029/1‑91 (MQ=255)
                                                                     cAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                         <  1:433771/47‑1 (MQ=255)
                                                                     cAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGa                                                         >  2:433771/1‑47 (MQ=255)
                                                                     cAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAc             >  2:617164/1‑91 (MQ=255)
                                                                     cAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAc             >  1:73629/1‑91 (MQ=255)
                                                                       aTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGaa                                                        >  1:501680/1‑46 (MQ=255)
                                                                       aTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGaa                                                        <  2:501680/46‑1 (MQ=255)
                                                                       aTCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACAt           >  2:653195/1‑91 (MQ=255)
                                                                         ccTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGaaa                                               <  1:727941/53‑1 (MQ=255)
                                                                         ccTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGaaa                                               >  2:727941/1‑53 (MQ=255)
                                                                         ccTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTg         <  2:746779/91‑1 (MQ=255)
                                                                          cTTCGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGc        <  1:839558/91‑1 (MQ=255)
                                                                            tGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCaa      >  2:382545/1‑91 (MQ=255)
                                                                            tGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCaa      >  2:805179/1‑91 (MQ=255)
                                                                                tgccCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAg  <  1:266720/91‑1 (MQ=255)
                                                                                   |                                                                                       
ATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAG  >  NZ_CP009273/4545513‑4545683

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGT  >  NZ_CP009273/4545503‑4545693
                                                                                             |                                                                                                 
GCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGT                                                                                            >  SRR3722118.373837/1‑101 (MQ=60)
 CATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTC                                                                                           <  SRR3722118.783767/101‑1 (MQ=60)
      AGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCT                                                                                      <  SRR3722118.89726/101‑1 (MQ=60)
                    GGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAG                                                                        >  SRR3722118.111201/1‑101 (MQ=60)
                            CGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGG                                                                <  SRR3722118.849579/101‑1 (MQ=60)
                                   gtctcgtgggctcggagatgtgtataagagacagATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAA                                                         <  SRR3722118.442924/67‑1 (MQ=60)
                                    TTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAAT                                                        <  SRR3722118.510512/101‑1 (MQ=60)
                                     TTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATA                                                       <  SRR3722118.351835/101‑1 (MQ=60)
                                        CTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATAT                                                    <  SRR3722118.211239/101‑1 (MQ=60)
                                          ACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCC                                                  <  SRR3722118.395338/101‑1 (MQ=60)
                                             tgggctcggagatgtgtataagagacaggGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTAT                                               <  SRR3722118.743926/72‑1 (MQ=60)
                                                  TACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGG                                          <  SRR3722118.885205/101‑1 (MQ=60)
                                                     CGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTC                                       <  SRR3722118.598112/101‑1 (MQ=60)
                                                      GTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATctgtctcttatacacatc                                      >  SRR3722118.476425/1‑83 (MQ=60)
                                                      GTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCAGTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCC                                      >  SRR3722118.880954/1‑101 (MQ=60)
                                                        CATTATCTGATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGG                                    >  SRR3722118.614462/1‑101 (MQ=60)
                                                                GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGctgtctcttatacac                            >  SRR3722118.86190/1‑86 (MQ=60)
                                                                GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCA                            >  SRR3722118.660415/1‑101 (MQ=60)
                                                                GATCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCA                            >  SRR3722118.233433/1‑101 (MQ=60)
                                                                  TCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGG                          >  SRR3722118.713227/1‑101 (MQ=60)
                                                                  TCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGG                          >  SRR3722118.786716/1‑101 (MQ=60)
                                                                  TCAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGG                          >  SRR3722118.820496/1‑101 (MQ=60)
                                                                   CAATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGA                         >  SRR3722118.130472/1‑101 (MQ=60)
                                                                     ATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAC                       >  SRR3722118.74995/1‑101 (MQ=60)
                                                                     ATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAC                       <  SRR3722118.669239/101‑1 (MQ=60)
                                                                     ATATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAAC                       <  SRR3722118.365637/101‑1 (MQ=60)
                                                                       ATAGATATCAATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATctgtctcttatacacatctccgagcccacgagacg                     >  SRR3722118.512321/1‑66 (MQ=60)
                                                                                 ATCCTTGGCTGCCCTTTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCActgtctcttataaacatctccga           >  SRR3722118.807812/1‑78 (MQ=60)
                                                                                 ATCCTTGGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGG           >  SRR3722118.87427/1‑101 (MQ=60)
                                                                                    CTTCGCTGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCC        <  SRR3722118.858105/101‑1 (MQ=60)
                                                                                          TGCCCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGT  <  SRR3722118.272148/101‑1 (MQ=60)
                                                                                             |                                                                                                 
GCATTAAGTAATGCGGCAAGGGCGAATCCGGAAGAATTTGCTACCCCTGTTACCGTCATTATCTGATCAATATAGATATCAATCCTTGGCTGCTCATTGGTCCATCTTTTTGCCTGAGAAGCTGGAAGGGTTGAAATATATCCTATCGTGGTTCCGGCATAACCAGGAACATTGCAACCAGGGCCTAATGT  >  NZ_CP009273/4545503‑4545693

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