Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 99,667 T→A 7.4% intergenic (+16/+91) PP_0093 → / ← PP_0094 zinc protease/5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294799,6670TA 7.4% 154.2 / 2.6 54intergenic (+16/+91)PP_0093/PP_0094zinc protease/5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base A (2/2);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

CGGCGCTGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGCGAAGTGGAAAGCGG  >  NC_002947/99530‑99813
                                                                                                                                            |                                                                                                                                                  
cggcgCTGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCAAAGCGGcccca                                                                                                                                          <  1:274815/149‑3 (MQ=255)
  gcgcTGGTTCGTGAGGGGGGTGGGGATTGGGTTCGGGTTTACGCGCGGGGCGAGCATTTGGGGCGGTTGCCGAGTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATcgc                                                                                                                                       <  7:341129/149‑1 (MQ=39)
    gcTGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGa                                                                                                                                     <  1:2674/149‑1 (MQ=255)
    gcTGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGa                                                                                                                                     >  2:118865/1‑149 (MQ=255)
     cTGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCAAAGCGGcccc                                                                                                                                           <  1:168570/142‑2 (MQ=255)
     cTGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCAAAGCGGcccc                                                                                                                                           >  2:168570/1‑141 (MQ=255)
      tGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGaca                                                                                                                                   >  2:7074/1‑149 (MQ=255)
          tCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGACCATCGCGACACAAg                                                                                                                               >  8:44616/1‑149 (MQ=255)
            gTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGc                                                                                                                             >  7:48591/1‑149 (MQ=255)
                  tgggtgggGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCTGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCcac                                                                                                                       <  5:174163/149‑1 (MQ=255)
                       ggggATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTTTATCGGGCCCCGTTTCGGGCCAAAGCGACACAACGGCCCCCCaa                                                                                                                     >  1:312300/1‑148 (MQ=255)
                          gATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTAAGATTT‑‑TGAGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGccc                                                                                                               >  7:340031/1‑149 (MQ=255)
                              gggTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGc                                                                                                            <  8:48591/148‑1 (MQ=255)
                                gTTCGGGTTTACGCGCTGGTCGATCATTTGGGGAGGTTGCCGTTTGAAGAAAGGATAGCGGCGTTGCAGCGTCTGATGTTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGAGACAAGGGGCCGTCAAAGGCCCTGCGTa                                                                                                         <  8:340031/149‑1 (MQ=37)
                                   cGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCATAAAGGCCCTGCGTACg                                                                                                       <  7:44616/148‑1 (MQ=255)
                                      ttttACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCg                                                                                                    >  8:309478/1‑148 (MQ=255)
                                        ttACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGc                                                                                                      >  6:351239/1‑144 (MQ=255)
                                        ttACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGc                                                                                                      <  5:351239/144‑1 (MQ=255)
                                          aCGCGCGGGTCGAGCAGTTGGGGCGGTTGCCGATCGACGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCg                                                                                                <  8:63609/148‑1 (MQ=255)
                                           cgcgcgGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAg                                                                                              <  3:170502/149‑1 (MQ=255)
                                               cgGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTATGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTgg                                                                                          <  1:118865/149‑1 (MQ=255)
                                                     gAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑AGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGc                                                                                    <  5:76229/149‑1 (MQ=255)
                                                          tttGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGtt                                                                               <  7:309478/149‑1 (MQ=255)
                                                                gCGGTGGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGGTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAag                                                                         <  3:207295/149‑1 (MQ=255)
                                                                  ggTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAagca                                                                       <  1:200680/149‑1 (MQ=255)
                                                                       ccGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCct                                                                  >  2:138577/1‑149 (MQ=255)
                                                                                     gaaagCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGc                                                    >  8:137940/1‑149 (MQ=255)
                                                                                     gaaagCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGc                                                    <  7:137940/149‑1 (MQ=255)
                                                                                       aagCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCt                                                  >  2:175391/1‑149 (MQ=255)
                                                                                                 gCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAg                                         <  1:138577/148‑1 (MQ=255)
                                                                                                   aGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCt                                      >  2:48053/1‑149 (MQ=255)
                                                                                                     gCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCt                                      >  8:123346/1‑148 (MQ=255)
                                                                                                      cGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGa                                    <  7:123346/149‑1 (MQ=255)
                                                                                                       gTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGAt                                   <  1:7074/149‑1 (MQ=255)
                                                                                                          tGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGc                                 >  5:72982/1‑148 (MQ=255)
                                                                                                           gATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCa                                >  2:324167/1‑148 (MQ=255)
                                                                                                           gATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCa                                >  8:38390/1‑148 (MQ=255)
                                                                                                                ttAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCgg                          <  2:312300/149‑1 (MQ=255)
                                                                                                                   gCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTa                       >  5:89423/1‑149 (MQ=255)
                                                                                                                     cTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCtt                                                 >  7:260748/1‑121 (MQ=255)
                                                                                                                     cTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCtt                                                 <  8:260748/121‑1 (MQ=255)
                                                                                                                      tttGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCg                    >  7:320923/1‑149 (MQ=255)
                                                                                                                      tttGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCg                    >  1:53026/1‑149 (MQ=255)
                                                                                                                      tttGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCg                    >  2:54261/1‑149 (MQ=255)
                                                                                                                      tttGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCg                    >  5:246253/1‑149 (MQ=255)
                                                                                                                        tGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGt                   <  7:38390/148‑1 (MQ=255)
                                                                                                                           gaTTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGCTCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGc                >  7:180070/1‑148 (MQ=255)
                                                                                                                            aTTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGAGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGCg               >  3:284793/1‑148 (MQ=255)
                                                                                                                                  tGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGCGAAGt           <  8:320923/148‑1 (MQ=255)
                                                                                                                                  atgggCCGCAAAGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGCGAAGTg          >  8:130742/3‑149 (MQ=255)
                                                                                                                                      gCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGt                                                                                            <  3:107789/63‑1 (MQ=255)
                                                                                                                                      gCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGt                                                                                            >  4:107789/1‑63 (MQ=255)
                                                                                                                                      gCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGCGAAGTGGa        >  6:188824/1‑147 (MQ=255)
                                                                                                                                          cTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGCGAAGTGGAAAGCgg  >  5:282994/1‑149 (MQ=255)
                                                                                                                                           tttGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCaaa                                                                          <  8:120506/76‑1 (MQ=255)
                                                                                                                                           tttGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCaaa                                                                          >  7:120506/1‑76 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CGGCGCTGTTTCGTGAGGTGGGTGGGGATTGGGTTCGGTTTTACGCGCGGGTCGAGCATTTGGGGCGGTTGCCGATTGAAGAAAGGAAAGCGGCGTTGCAGCGTCTGATGGTTAGCCTTTGAGATTT‑‑TGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCCCCACAAAGGCCCTGCGTACGCCGATCGAGGTGGGAGGGCCGGTTGCAAAGCAGCCCTGCTGGCGCTCAAGCCTTGATAAAAGCCTGATGCAAATCGGCTACCGTGGCGAAGTGGAAAGCGG  >  NC_002947/99530‑99813

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