Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,918 2 bp→CG coding (348‑349/1236 nt) aslB → putative arylsulfatase regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9180GC98.6% 380.2 / NA 148W116C (TGG→TGCaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base C (65/80);  total (66/82)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,980,9190CG98.7% 392.5 / NA 150P117A (CCG→GCG) aslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (66/82);  minor base T (0/2);  total (66/84)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 150 total aligned reads displayed.

ATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAAT  >  NC_000913/3980849‑3980988
                                                                     ||                                                                     
atatGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttt                                                                  >  2:1783343/1‑76 (MQ=255)
atatGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttt                                                                  >  2:4490554/1‑76 (MQ=255)
 tatGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtc                                                                 <  2:4024256/76‑1 (MQ=255)
  atGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtct                                                                <  1:983203/76‑1 (MQ=255)
   tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                               <  1:3544620/76‑1 (MQ=255)
   tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                               <  2:3548289/76‑1 (MQ=255)
   tGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                                               <  2:35584/76‑1 (MQ=255)
    ggTGCTGGCAGGAAGATAAGTAAGAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcg                                                              <  2:1067416/76‑1 (MQ=255)
     gTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgc                                                             <  2:1733622/76‑1 (MQ=255)
      tGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgcg                                                            <  2:3631192/76‑1 (MQ=255)
      gggtGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgcg                                                            <  2:1078746/73‑1 (MQ=255)
       gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                           >  2:4412672/1‑76 (MQ=255)
       gCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                                           >  2:4480327/1‑76 (MQ=255)
        cTGGCAGGAATATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                          <  1:4859129/76‑1 (MQ=255)
        cTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                                          <  1:3180563/76‑1 (MQ=255)
         tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                                         >  2:1962807/1‑76 (MQ=255)
         tGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                                         <  1:1160025/76‑1 (MQ=255)
          ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGc                                                        <  1:4100248/76‑1 (MQ=255)
          ggCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGc                                                        <  2:1714846/76‑1 (MQ=255)
           gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGca                                                       <  2:3486034/76‑1 (MQ=255)
           gCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGca                                                       <  1:4213237/76‑1 (MQ=255)
             aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatc                                                     <  2:1318877/76‑1 (MQ=255)
             aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatc                                                     <  1:54201/76‑1 (MQ=255)
             aGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatc                                                     >  2:3863605/1‑76 (MQ=255)
                aaGATAAGTAGCGGCTTCCAGACTAACGGCGTGTGACTGGATGACGAATTGTGCGCGTTTCTCGCGGAGCAGCAtt                                                  <  1:4296610/76‑1 (MQ=255)
                 aGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAttt                                                 <  2:263174/76‑1 (MQ=255)
                  gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCCTCAtttt                                                >  2:2006660/1‑76 (MQ=255)
                  gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtttt                                                <  2:1743973/76‑1 (MQ=255)
                  gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCAtttt                                                >  1:3103520/1‑76 (MQ=255)
                   aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTc                                               >  2:494686/1‑76 (MQ=255)
                   aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTc                                               >  1:3655132/1‑76 (MQ=255)
                     aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCtt                                             <  2:792308/76‑1 (MQ=255)
                        tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                          >  1:79661/1‑76 (MQ=255)
                        tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                          >  1:3463004/1‑76 (MQ=255)
                        tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCttgtt                                          >  1:4311876/1‑76 (MQ=255)
                         aaCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGGCCGCGTTTCCCGCGGACCACCATTCTCttgttg                                         >  2:310575/1‑76 (MQ=255)
                          aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTgg                                        >  1:3995354/1‑76 (MQ=255)
                          aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTgg                                        >  1:2478962/1‑76 (MQ=255)
                           cAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                       <  1:3804077/76‑1 (MQ=255)
                           cAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                                       >  2:2529262/1‑76 (MQ=255)
                            accTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGCGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                      >  2:2015321/3‑76 (MQ=255)
                            aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTggat                                      >  1:2158277/1‑74 (MQ=255)
                            aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                      >  2:2142091/1‑76 (MQ=255)
                            aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGtt                                      >  1:4540393/1‑76 (MQ=255)
                              ctgccAGACTAACGGCGGGGTGCTGGATGACGAATGGTGCGCGTGTCGCGCGGAGCATCATTGTCTTGTTGGTTTa                                    <  2:2864534/73‑1 (MQ=255)
                              cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTa                                    <  2:3422861/76‑1 (MQ=255)
                                tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATc                                  <  1:4579970/76‑1 (MQ=255)
                                tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATc                                  <  1:1955210/76‑1 (MQ=255)
                                tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATc                                  >  1:2640718/1‑76 (MQ=255)
                                 ctaGACGAGCGGCGTTCTGCTGGATTACGAATGGTGCGCGTTTCTCGCGGATCATCATTTTCTTGTTGGTTTATCg                                 <  2:3281223/74‑1 (MQ=255)
                                 ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCg                                 <  1:1356690/76‑1 (MQ=255)
                                 ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCg                                 >  2:3841881/1‑76 (MQ=255)
                                 ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCg                                 >  1:3635357/1‑76 (MQ=255)
                                  cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGc                                <  1:4151793/76‑1 (MQ=255)
                                  cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGc                                <  2:2026088/76‑1 (MQ=255)
                                   aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCt                               >  1:2915843/1‑76 (MQ=255)
                                   aGACTAACGGCGGGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCt                               <  2:2244585/76‑1 (MQ=255)
                                    gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTg                              <  2:3123815/76‑1 (MQ=255)
                                    gACTAACGGCGTGCTGCTGGATGACGAAGGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTg                              <  1:3101404/76‑1 (MQ=255)
                                     aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCggg                             >  1:2431159/1‑76 (MQ=255)
                                     aCTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTgg                             <  1:83744/76‑1 (MQ=255)
                                       tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGAt                           >  2:2772848/1‑76 (MQ=255)
                                        aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATg                          >  2:421122/1‑76 (MQ=255)
                                           ggCGTGCTGCTGGATGCCGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTGTGTTTATCGCTGGATggac                       >  1:3392518/1‑74 (MQ=255)
                                           ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                       >  1:4489493/1‑76 (MQ=255)
                                           ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                       >  1:4002955/1‑76 (MQ=255)
                                           ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                       >  1:3355790/1‑76 (MQ=255)
                                           ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                       <  1:1572011/76‑1 (MQ=255)
                                           cgggTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc                       <  2:3865615/73‑1 (MQ=255)
                                            gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGccc                      <  2:2252186/76‑1 (MQ=255)
                                            gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGccc                      >  1:3027784/1‑76 (MQ=255)
                                             gttgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCg                     <  1:3192017/74‑1 (MQ=255)
                                              gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGc                    >  2:3374319/1‑76 (MQ=255)
                                              gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGc                    >  1:2620337/1‑76 (MQ=255)
                                               tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGcc                   <  1:1565268/76‑1 (MQ=255)
                                               tgctgcTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGcc                   >  1:1150888/1‑76 (MQ=255)
                                                gctgctGGATGACGAATGATGCTCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCt                  <  1:4114615/76‑1 (MQ=255)
                                                 ctgctgGATGTCTGATGGTTGGCTTTTCTCGCGGATGATCGTTGTCTTGTTGGTTTATCTCTGGATGGCCCGCCTg                 <  1:2174068/76‑1 (MQ=255)
                                                  tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTga                <  1:2406104/76‑1 (MQ=255)
                                                   gctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTgag               <  1:1406876/76‑1 (MQ=255)
                                                   gctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTgag               >  1:4695349/1‑76 (MQ=255)
                                                    gtgGTTGAGGAAGGGTGCGCGTTTGTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTgaga              <  2:4749695/75‑1 (MQ=255)
                                                    ctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTgaga              <  1:3863605/76‑1 (MQ=255)
                                                     tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGAt             >  2:3187217/1‑76 (MQ=255)
                                                     tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGAt             <  2:3060640/76‑1 (MQ=255)
                                                      ggATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATc            >  2:1631907/1‑76 (MQ=255)
                                                       gaggACGAAGGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc           <  2:3400667/73‑1 (MQ=255)
                                                       gATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc           <  1:2478991/76‑1 (MQ=255)
                                                       gATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATcc           >  2:2902029/1‑76 (MQ=255)
                                                        aTGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCca          >  1:1709207/1‑76 (MQ=255)
                                                         tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcac         >  1:1864308/1‑76 (MQ=255)
                                                         tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcac         >  1:1622680/1‑76 (MQ=255)
                                                         tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcac         >  1:3768188/1‑76 (MQ=255)
                                                          gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGTCGGATGGCCCGCCTGAGTTCcaca        >  2:4137611/1‑76 (MQ=255)
                                                          gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcaca        >  2:3418200/1‑76 (MQ=255)
                                                          gACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCcaca        >  2:2735377/1‑76 (MQ=255)
                                                             gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc     <  2:1880216/76‑1 (MQ=255)
                                                             gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTCGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatc     >  1:4577957/1‑76 (MQ=255)
                                                              aaTGGTGCTCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatca    <  2:465433/76‑1 (MQ=255)
                                                              aaTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCAcaatca    >  1:3909163/1‑76 (MQ=255)
                                                                     ||                                                                     
ATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAAT  >  NC_000913/3980849‑3980988

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