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 I199 R1 220 57.5 3163680 97.0% 3068769 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,814,267 A→G A368A (GCT→GCC alaS ← alanine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,814,2670AG93.3% 126.7 / 2.0 45A368A (GCT→GCCalaSalanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (17/25);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAC  >  NZ_CP009273/2814185‑2814350
                                                                                  |                                                                                   
aaGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACt                                                                              <  2:608360/90‑1 (MQ=255)
  gCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGc                                                                            <  1:1378740/90‑1 (MQ=255)
    aGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctc                                                                          <  2:103204/90‑1 (MQ=255)
     gTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctct                                                                         <  1:1320645/90‑1 (MQ=255)
     gTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctct                                                                         <  1:1566868/90‑1 (MQ=255)
            cATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcag                                                                  <  2:737827/90‑1 (MQ=255)
             aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt                                                                 >  2:843219/1‑90 (MQ=255)
             aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt                                                                 <  2:403632/90‑1 (MQ=255)
             aTCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagt                                                                 <  1:424101/90‑1 (MQ=255)
                  gCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTtcttcagtcttca                                                            >  1:545159/1‑90 (MQ=255)
                      aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  1:1228485/90‑1 (MQ=255)
                      aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  2:124545/90‑1 (MQ=255)
                      aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  1:684936/90‑1 (MQ=255)
                      aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                        <  1:327893/90‑1 (MQ=255)
                       tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCttcttc                                                                    >  2:600847/1‑77 (MQ=255)
                       tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCttcttc                                                                    <  1:600847/77‑1 (MQ=255)
                       tCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAcc                                                       >  2:943356/1‑90 (MQ=255)
                           cAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCt                                                   <  2:375624/90‑1 (MQ=255)
                            aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  >  2:1319464/1‑90 (MQ=255)
                            aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  >  2:566668/1‑90 (MQ=255)
                            aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  <  2:1544249/90‑1 (MQ=255)
                            aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                                  <  2:276061/90‑1 (MQ=255)
                            aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTtcttcag                                                                  >  1:1292335/1‑74 (MQ=255)
                            aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTtcttcag                                                                  <  2:1292335/74‑1 (MQ=255)
                                          gCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCctgct                                    >  2:1300267/1‑90 (MQ=255)
                                           ctctTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCctgctg                                   <  1:398287/90‑1 (MQ=255)
                                           ctctTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCctgctg                                   <  2:308236/90‑1 (MQ=255)
                                                  tCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTc                            <  1:778365/90‑1 (MQ=255)
                                                  tCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTc                            <  1:798355/90‑1 (MQ=255)
                                                    cAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAg                          >  1:688005/1‑90 (MQ=255)
                                                     aGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAgg                         <  2:951247/90‑1 (MQ=255)
                                                      gCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGt                        >  1:1367359/1‑90 (MQ=255)
                                                         aCGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCtt                     >  1:1056070/1‑90 (MQ=255)
                                                              aGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCAccc                >  2:81426/1‑90 (MQ=255)
                                                              aGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCAccc                >  1:276076/1‑90 (MQ=255)
                                                                   gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCa                                                         >  2:795378/1‑44 (MQ=255)
                                                                   gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCa                                                         <  1:795378/44‑1 (MQ=255)
                                                                   gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  2:1450185/1‑90 (MQ=255)
                                                                   gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  1:206883/1‑90 (MQ=255)
                                                                   gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  1:269384/1‑90 (MQ=255)
                                                                   gcgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  2:752289/1‑90 (MQ=255)
                                                                   gcgcTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCaga           >  1:1289922/1‑90 (MQ=255)
                                                                         cAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCAt     <  2:1428839/90‑1 (MQ=255)
                                                                            agTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAc  <  2:269384/90‑1 (MQ=255)
                                                                            agTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAc  <  2:206883/90‑1 (MQ=255)
                                                                                  |                                                                                   
AAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAAC  >  NZ_CP009273/2814185‑2814350

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGATCAGCGGACCA  >  NZ_CP009273/2814169‑2814366
                                                                                                  |                                                                                                   
GTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGG                                                                                                    >  SRR3722078.964609/1‑100 (MQ=60)
 TGTCGTACAGACGGAAAGCAGTTTCACCACCCAGCGTATCACCAGAAAGTTTTGCCAGCTCCTCATCCAGCAACGCCAGACCGCGCTCAAGAGTACGGGC                                                                                                   <  SRR3722078.1369203/100‑1 (MQ=60)
  GTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCA                                                                                                  <  SRR3722078.332366/100‑1 (MQ=60)
       ACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTG                                                                                             <  SRR3722078.258870/100‑1 (MQ=60)
        CAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGC                                                                                            <  SRR3722078.131617/100‑1 (MQ=60)
                  GCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAG                                                                                  <  SRR3722078.1400978/100‑1 (MQ=60)
                     GTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCT                                                                               <  SRR3722078.1341754/100‑1 (MQ=60)
                     GTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCT                                                                               <  SRR3722078.1592378/100‑1 (MQ=60)
                        TCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCA                                                                            >  SRR3722078.551597/1‑100 (MQ=60)
                          cagATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGC                                                                          <  SRR3722078.607795/97‑1 (MQ=60)
                             ATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACC                                                                       <  SRR3722078.429211/100‑1 (MQ=60)
                                  GCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGtct                                                                  >  SRR3722078.1312906/1‑97 (MQ=60)
                                     ccatcgtcggcagcgtcagatgtgtataagagacagACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAAC                                                               <  SRR3722078.805830/64‑1 (MQ=60)
                                      ATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGTTCAACC                                                              <  SRR3722078.692973/100‑1 (MQ=60)
                                      ATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACC                                                              <  SRR3722078.331914/100‑1 (MQ=60)
                                      ATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACC                                                              <  SRR3722078.1247838/100‑1 (MQ=60)
                                                          GCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAG                                          >  SRR3722078.696124/1‑100 (MQ=60)
                                                           CTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGG                                         <  SRR3722078.403119/100‑1 (MQ=60)
                                                            TCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGT                                        >  SRR3722078.1389362/1‑100 (MQ=60)
                                                               TCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTT                                     >  SRR3722078.1071930/1‑100 (MQ=60)
                                                                  TCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCAC                                  <  SRR3722078.788460/100‑1 (MQ=60)
                                                                  TCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCAC                                  <  SRR3722078.808857/100‑1 (MQ=60)
                                                                    CAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCC                                >  SRR3722078.279443/1‑100 (MQ=60)
                                                                         ACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGA                           >  SRR3722078.209382/1‑100 (MQ=60)
                                                                         ACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGA                           >  SRR3722078.272677/1‑100 (MQ=60)
                                                                         ACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGA                           >  SRR3722078.1310451/1‑100 (MQ=60)
                                                                                      CTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCTTGCTGGCGTTTCAGGTCTTCACCCGCGGAGCCCATAACGTCG              >  SRR3722078.832113/1‑100 (MQ=60)
                                                                                      CTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCG              >  SRR3722078.1385179/1‑100 (MQ=60)
                                                                                        CCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGtctcttatacacatctgacgctgccgacgatggcggtcgtgtagatctcggtggtcg            >  SRR3722078.288650/1‑43 (MQ=60)
                                                                                             GTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGATCAGCG       >  SRR3722078.869062/1‑100 (MQ=60)
                                                                                               ACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGATCAGCGGA     >  SRR3722078.674972/1‑100 (MQ=60)
                                                                                                  GGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGATCAGCGGACCA  >  SRR3722078.718063/1‑100 (MQ=60)
                                                                                                  |                                                                                                   
GTGTCGTACAGACGGAAAGCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCATAACGTCGATCAGCGGACCA  >  NZ_CP009273/2814169‑2814366

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: