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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,350,381 G→C G457G (GGC→GGG cadC ← lysine decarboxylation/transport transcriptional activator CadC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,350,3810GC100.0% 180.6 / NA 52G457G (GGC→GGGcadClysine decarboxylation/transport transcriptional activator CadC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (24/28);  total (24/28)

GGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTT  >  NZ_CP009273/4350296‑4350465
                                                                                     |                                                                                    
gggTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCaa                                                                                  >  2:846836/1‑90 (MQ=255)
gggTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCaa                                                                                  >  1:1202880/1‑90 (MQ=255)
  gtgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                <  2:1202880/90‑1 (MQ=255)
  gtgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                <  2:1163712/90‑1 (MQ=255)
  gtgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                <  2:859682/90‑1 (MQ=255)
  gtgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                >  1:1163712/1‑90 (MQ=255)
   tgtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCa                                                                               >  2:773489/1‑90 (MQ=255)
    gtTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCaa                                                                              <  1:516499/90‑1 (MQ=255)
        gCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAAcaca                                                                          <  1:331105/90‑1 (MQ=255)
        gCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAAcaca                                                                          <  2:1224654/90‑1 (MQ=255)
            cTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAAt                                                                      >  2:1141/1‑90 (MQ=255)
             tGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAAtt                                                                     <  1:812912/90‑1 (MQ=255)
              gggCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAAttt                                                                    >  1:229895/1‑90 (MQ=255)
                  gTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGcc                                                                >  1:926948/1‑90 (MQ=255)
                       ttAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGAc                                                           <  2:1223229/90‑1 (MQ=255)
                        tAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACa                                                          <  1:906556/90‑1 (MQ=255)
                         aaaGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACAt                                                         >  1:974054/1‑90 (MQ=255)
                               ggTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAg                                                   >  1:268857/1‑90 (MQ=255)
                                gTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                <  2:904128/60‑1 (MQ=255)
                                gTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGc                                                                                >  1:904128/1‑60 (MQ=255)
                                gTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGa                                                  >  2:571933/1‑90 (MQ=255)
                                     atatGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAAt                                             <  1:795143/90‑1 (MQ=255)
                                          cATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAg                                        >  2:1167023/1‑90 (MQ=255)
                                          cATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAg                                        >  1:1044643/1‑90 (MQ=255)
                                          cATAAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAg                                        <  1:1059601/90‑1 (MQ=255)
                                           aTCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCTGGACGTTTCAAGATCAATGCCAGt                                       >  1:427662/1‑90 (MQ=255)
                                            tCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTa                                      >  1:825552/1‑90 (MQ=255)
                                            tCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTa                                      >  1:149006/1‑90 (MQ=255)
                                              aGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACTCATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTAtt                                    >  1:498756/1‑90 (MQ=255)
                                              aGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTAtt                                    >  2:644287/1‑90 (MQ=255)
                                              aGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTAtt                                    >  1:912523/1‑90 (MQ=255)
                                              aGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTAtt                                    >  2:311572/1‑90 (MQ=255)
                                                  gctTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATc                                <  2:388091/90‑1 (MQ=255)
                                                    ttCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGc                              <  2:936112/90‑1 (MQ=255)
                                                    ttCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGc                              <  2:366108/90‑1 (MQ=255)
                                                     tCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGcc                             <  2:1096361/90‑1 (MQ=255)
                                                            tCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAg                      >  2:415370/1‑90 (MQ=255)
                                                               tCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACt                   >  2:117334/1‑90 (MQ=255)
                                                                  ccTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCAt                <  1:981071/90‑1 (MQ=255)
                                                                   cTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATc               <  2:1035052/90‑1 (MQ=255)
                                                                   cTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATc               <  1:806320/90‑1 (MQ=255)
                                                                   cTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATc               <  1:543860/90‑1 (MQ=255)
                                                                    ttcatttcatAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCt              >  2:161494/1‑90 (MQ=255)
                                                                         ttcatAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGtttt         <  2:13699/90‑1 (MQ=255)
                                                                         ttcatAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGtttt         <  2:912523/90‑1 (MQ=255)
                                                                         ttcatAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGtttt         <  1:509152/90‑1 (MQ=255)
                                                                           catAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGttttac       <  1:644287/90‑1 (MQ=255)
                                                                             tAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTGCAAGATCAATGCCAGTATTTATGGCCTGGTAAGACTCATCTGttttacct     >  2:949732/1‑90 (MQ=255)
                                                                              aaaCCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGttttacctt    <  1:949732/90‑1 (MQ=255)
                                                                                aCCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGttttacctttt  <  2:52971/90‑1 (MQ=255)
                                                                                aCCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGttttacctttt  <  1:67286/90‑1 (MQ=255)
                                                                                aCCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGttttacctttt  <  2:861846/90‑1 (MQ=255)
                                                                                     |                                                                                    
GGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTT  >  NZ_CP009273/4350296‑4350465

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCAGAGCAC  >  NZ_CP009273/4350286‑4350475
                                                                                               |                                                                                              
ATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAA                                                                                            >  SRR3722089.1221174/1‑100 (MQ=60)
  ctacactaagatggtcgtcggcagcgtcagatgtgtataagagacagATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGC                                                                                          <  SRR3722089.513566/53‑1 (MQ=60)
  CCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGC                                                                                          >  SRR3722089.1181407/1‑100 (MQ=60)
       AAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACAC                                                                                     <  SRR3722089.784615/100‑1 (MQ=60)
              GTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTT                                                                              <  SRR3722089.522562/100‑1 (MQ=60)
              GTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTT                                                                              >  SRR3722089.232813/1‑100 (MQ=60)
                  GCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCC                                                                          >  SRR3722089.940479/1‑100 (MQ=60)
                  GCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCC                                                                          <  SRR3722089.335247/100‑1 (MQ=60)
                       TGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGAC                                                                     <  SRR3722089.824278/100‑1 (MQ=60)
                         GGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACAT                                                                   >  SRR3722089.988438/1‑100 (MQ=60)
                               AATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAG                                                             >  SRR3722089.272298/1‑100 (MQ=60)
                                ATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATctgtctcttatacacatctg                                                            >  SRR3722089.917260/1‑80 (MQ=60)
                                  TAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATC                                                          <  SRR3722089.919726/100‑1 (MQ=60)
                                          GTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAG                                                  >  SRR3722089.1060287/1‑100 (MQ=60)
                                           TGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCTGGACGTTTCAAGATCAATGCCAGT                                                 >  SRR3722089.432816/1‑100 (MQ=60)
                                            GAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTA                                                >  SRR3722089.837192/1‑100 (MQ=60)
                                            GAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTA                                                >  SRR3722089.150856/1‑100 (MQ=60)
                                              GATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACTCATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATT                                              >  SRR3722089.504661/1‑100 (MQ=60)
                                              GATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATT                                              >  SRR3722089.925808/1‑100 (MQ=60)
                                               ATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTT                                             <  SRR3722089.806099/100‑1 (MQ=60)
                                                    CATAAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGC                                        <  SRR3722089.1075499/100‑1 (MQ=60)
                                                                            CCTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACC                <  SRR3722089.995597/100‑1 (MQ=60)
                                                                             CTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCT               <  SRR3722089.817539/100‑1 (MQ=60)
                                                                             CTTCATTTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCT               <  SRR3722089.550191/100‑1 (MQ=60)
                                                                                   TTCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACC         <  SRR3722089.515150/100‑1 (MQ=60)
                                                                                    TCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCA        >  SRR3722089.246339/1‑100 (MQ=60)
                                                                                    TCATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCA        >  SRR3722089.221002/1‑100 (MQ=60)
                                                                                     CATAAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCAG       <  SRR3722089.652526/100‑1 (MQ=60)
                                                                                        AAACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCAGAGC    <  SRR3722089.963688/100‑1 (MQ=60)
                                                                                          ACCTTCCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCAGAGCAC  <  SRR3722089.68101/100‑1 (MQ=60)
                                                                                               |                                                                                              
ATCCAGTAAAGGGTGTTTGCCCCTGGGCGTAAATTAAAGGCGGTGAGATATGCATCAGCTGCTTCCCGGTTCATCCCCTTCATTTCATAAACCTTGCCAAGCAACACATAATTTAGCCAGGACATTTCAAGATCAATGCCAGTATTTATCGCCTGGTAAGACTCATCTGTTTTACCTTTTACCAGAGCAC  >  NZ_CP009273/4350286‑4350475

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