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,283,874 Δ1 bp intergenic (‑350/‑105) BW25113_RS11435 ← / → narP IS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,283,8710T.94.7% 138.7 / 3.9 38intergenic (‑347/‑108)BW25113_RS11435/narPIS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (8/28);  total (8/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

TTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTC  >  NZ_CP009273/2283783‑2283955
                                                                                        |                                                                                    
ttttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtt                                                                                     <  2:1484874/90‑1 (MQ=255)
 tttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAGCGAACGACAAGGTTGTAGtttc                                                                                   <  2:1504988/90‑2 (MQ=255)
  ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtttc                                                                                   >  2:537199/1‑89 (MQ=255)
  ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtttc                                                                                   <  1:1520300/90‑2 (MQ=255)
  ttttAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtttc                                                                                   <  2:1096663/90‑2 (MQ=255)
       gttgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTAc                                                                             <  1:763357/90‑1 (MQ=255)
         tgttTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACAt                                                                           >  1:33138/1‑90 (MQ=255)
              ttCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCca                                                                      <  1:661962/90‑1 (MQ=255)
              ttCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCca                                                                      <  2:247722/90‑1 (MQ=255)
                        gaAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAAt                                                            >  2:766572/1‑90 (MQ=255)
                        gaAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAAt                                                            <  1:755340/90‑1 (MQ=255)
                              gATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtt                                                                                     >  1:387876/1‑60 (MQ=255)
                              gATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGtt                                                                                     <  2:387876/60‑1 (MQ=255)
                              gATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTa                                                      <  1:755784/90‑1 (MQ=255)
                                    aaaTTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTc                                                >  1:174006/1‑90 (MQ=255)
                                       ttCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCAcg                                             <  1:381789/90‑1 (MQ=255)
                                         ccGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCAcgcg                                           <  1:818052/90‑1 (MQ=255)
                                          cGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCAcgcg                                           <  1:1356724/90‑1 (MQ=255)
                                          cGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACgcgc                                          <  1:631098/90‑1 (MQ=255)
                                              tttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTgg                                      >  2:862449/1‑90 (MQ=255)
                                               ttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGt                                     <  2:1261936/90‑1 (MQ=255)
                                               ttATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGt                                     <  2:165628/90‑1 (MQ=255)
                                                tataATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTa                                    <  2:655304/90‑1 (MQ=255)
                                                   aaTGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCg                                 <  2:170766/90‑1 (MQ=255)
                                                     tGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGt                                <  2:528458/89‑1 (MQ=255)
                                                     tGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGt                                >  1:528458/1‑89 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTAcc                            <  2:598579/90‑1 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCa                           <  1:914679/90‑1 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCa                           <  1:1539507/90‑1 (MQ=255)
                                                         aaTGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCa                           <  2:154034/90‑1 (MQ=255)
                                                          atgatgCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCaa                          <  2:174006/90‑1 (MQ=255)
                                                              tgCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAacgct                      >  2:130426/1‑90 (MQ=255)
                                                                ccAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAccgctac                    <  2:1201268/90‑1 (MQ=255)
                                                                ccAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAacgctac                    <  2:1271035/90‑1 (MQ=255)
                                                                  aaaGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAacgctacgc                  <  2:648419/90‑1 (MQ=255)
                                                                    aGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTc                <  1:1400911/90‑1 (MQ=255)
                                                                      cGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCaa              <  2:1581297/90‑1 (MQ=255)
                                                                            aCAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAAGATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATaa        <  1:245379/90‑1 (MQ=255)
                                                                            aCAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATaa        <  1:521192/90‑1 (MQ=255)
                                                                               aGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGa     >  1:894418/1‑90 (MQ=255)
                                                                                 gTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGAtt   <  1:581493/90‑1 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  <  2:1231874/90‑1 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  <  2:41921/90‑1 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  >  1:1260552/1‑90 (MQ=255)
                                                                                  ttGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTc  <  1:1256853/90‑1 (MQ=255)
                                                                                        |                                                                                    
TTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTC  >  NZ_CP009273/2283783‑2283955

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ 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
GATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAG  >  NZ_CP009273/2283777‑2283972
                                                                                                 |                                                                                                  
GATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCAC                                                                                                 <  SRR3722078.9698/100‑1 (MQ=60)
GATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCAC                                                                                                 <  SRR3722078.1464959/100‑1 (MQ=60)
 ATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACT                                                                                                <  SRR3722078.407935/100‑1 (MQ=60)
  TTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGCCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTA                                                                                               <  SRR3722078.1584102/100‑1 (MQ=60)
     ATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACAT                                                                                            >  SRR3722078.33548/1‑100 (MQ=60)
        TTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTC                                                                                         <  SRR3722078.1545092/100‑1 (MQ=60)
             GTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATAC                                                                                    <  SRR3722078.773179/100‑1 (MQ=60)
                    TTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAAT                                                                             <  SRR3722078.669542/100‑1 (MQ=60)
                          ATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGctgtctcttatacacatctg                                                                       >  SRR3722078.392603/1‑80 (MQ=60)
                              GAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACAT                                                                   <  SRR3722078.765009/100‑1 (MQ=60)
                                AGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTC                                                                 >  SRR3722078.176043/1‑100 (MQ=60)
                                    GATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGC                                                             <  SRR3722078.765462/100‑1 (MQ=60)
                                             TTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAG                                                    <  SRR3722078.386449/100‑1 (MQ=60)
                                               CCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCG                                                  <  SRR3722078.828869/100‑1 (MQ=60)
                                                CGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCG                                                  <  SRR3722078.1378506/100‑1 (MQ=60)
                                                CGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGT                                                 <  SRR3722078.638358/100‑1 (MQ=60)
                                                 GTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTT                                                >  SRR3722078.534702/1‑100 (MQ=60)
                                                               AATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCT                                  <  SRR3722078.927091/100‑1 (MQ=60)
                                                               AATGATGCCAAAGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCT                                  <  SRR3722078.1564606/100‑1 (MQ=60)
                                                                          agacagCGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATG                       <  SRR3722078.1443609/94‑1 (MQ=60)
                                                                          agacagCGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATG                       <  SRR3722078.620950/94‑1 (MQ=60)
                                                                          AGCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATG                       <  SRR3722078.1423601/100‑1 (MQ=60)
                                                                           GCGAACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGA                      >  SRR3722078.906440/1‑100 (MQ=60)
                                                                              AACGACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTC                   >  SRR3722078.1280535/1‑100 (MQ=60)
                                                                                  ACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAAGATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAAT               <  SRR3722078.248379/100‑1 (MQ=60)
                                                                                  ACAAGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAAT               <  SRR3722078.527355/100‑1 (MQ=60)
                                                                                     AGGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAA            >  SRR3722078.1374133/1‑100 (MQ=60)
                                                                                      GGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAA           >  SRR3722078.1024235/1‑100 (MQ=60)
                                                                                      GGTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAA           >  SRR3722078.669551/1‑100 (MQ=60)
                                                                                       GTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAAC          <  SRR3722078.588273/100‑1 (MQ=60)
                                                                                       GTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAAC          >  SRR3722078.947409/1‑100 (MQ=60)
                                                                                       GTTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAAC          >  SRR3722078.356253/1‑100 (MQ=60)
                                                                                        TTGTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACC         <  SRR3722078.1276762/100‑1 (MQ=60)
                                                                                          GTAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTC       <  SRR3722078.1018538/100‑1 (MQ=60)
                                                                                           TAGTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACGTAATGATTCTAATAAAACCTCA      >  SRR3722078.774721/1‑100 (MQ=60)
                                                                                                TTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAG  <  SRR3722078.131932/100‑1 (MQ=60)
                                                                                                 |                                                                                                  
GATTGATTTTTTAGTTGTTTTTCTTGATGAGAAGCTGATGCAAAATTCCGTCTTTATAATGAAAATGATGCCAAAGCGAACGACAAGGTTGTAGTTTTCACTACATGTCCATACATAAAATGGGGTAACATTCACGCGCCTGGTAGCGTTACCAACGCTACGCTCAAACATAATGATTCTAATAAAACCTCAGGAG  >  NZ_CP009273/2283777‑2283972

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