Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,841 A→C 9.1% intergenic (+22/‑548) PP_0159 → / → PP_16SA family III CoA‑transferase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947170,8410AC 9.1% 138.7 / 7.8 55intergenic (+22/‑548)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base C (3/2);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-02

GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCA  >  NC_002947/170694‑170988
                                                                                                                                                    |                                                                                                                                                   
gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGAGCTCTAGATGAAAAGCGGGGAGGCCGTCAg                                                                                                                                                    <  1:31590/149‑1 (MQ=255)
  gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                                  <  8:230049/149‑1 (MQ=255)
     gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGTCAGGCCTc                                                                                                                                               >  1:172278/1‑149 (MQ=255)
     gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGGGTTGGGGCTGGATGCTGGCGAGGTACAGCGGCTGCGTGAGGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGGCAGGCCTc                                                                                                                                               >  5:289559/1‑149 (MQ=255)
       cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                             <  4:215131/149‑1 (MQ=255)
        aGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGGGGTGCGGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCTTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTcccc                                                                                                                                            <  8:288625/149‑1 (MQ=255)
        aGTGCAGTATCCCCGGCCGCCACCACTGCGGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGGGTTGGGGCGGGATGCGGGCGCGGTACAGCGTCTGCGTGATGCTGGGGGGCTCTAGATGAAACGCGGGGAGGCCTGACGGCCTcccc                                                                                                                                            <  4:65354/149‑1 (MQ=255)
           gtaGTTTCGCGGGGCGCACCCAGTGGTGGGGGAGCATACCGATGGGGTGCTGGGGTATGTGTGGGGGCTGGATGCTGTCGAGGTCCAGCGTCGGGGTGTGGCTGGGGTGCGCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCt                                                                                                                                         <  8:197504/147‑1 (MQ=255)
           ggAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCt                                                                                                                                         >  2:126095/1‑149 (MQ=255)
                            caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTa                                                                                                                         >  3:103255/1‑148 (MQ=255)
                             accccggctgGGCGAGCATACCGCCGTGGGGCTGGGGGAAGTGGTGGGGCTGGAGGCGGGCGAGGTACAGCGTCTGCGTGATGGTGGGGGGCGCTAGCTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATc                                                                                                                       <  6:289559/143‑1 (MQ=37)
                                                aCCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCGATCAAAGGTAATTGAAAttt                                                                                                      >  3:240489/1‑146 (MQ=255)
                                                  ctaaGTGGTGGTGGGGTATGTGTGGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTTATGCTGGGTTGCTCTATATCAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAgg                                                                                                  <  5:203582/147‑1 (MQ=38)
                                                  cGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGttttt                                                                                                                                 >  3:229882/1‑118 (MQ=255)
                                                          tgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGcc                                                                                          >  1:51937/1‑149 (MQ=255)
                                                            cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGTCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGccc                                                                                         >  5:113374/1‑148 (MQ=255)
                                                               ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTAAGCGGTATGAAAGGTAATTGAAATTTAAGGTTGACGCCctcgc                                                                                     >  4:152885/1‑147 (MQ=255)
                                                               ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTAAGGCCAACCCGCTTTTTTTTTAAGCGCTATAAAAGGTAATTGAAATTAAAGGTTTACGCCCtgtc                                                                                     >  8:56011/1‑149 (MQ=255)
                                                                 ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGGCGGCAGGCCTCTCCGCGTTGTTTTTGAGCGCGATAAAAGGTTATTGGAATTAAAGGTTTAAGCCCtctct                                                                                    >  3:10383/1‑147 (MQ=255)
                                                                 ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                                     <  8:107704/147‑1 (MQ=255)
                                                                 ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCg                                                                                    >  7:107704/1‑148 (MQ=255)
                                                                  gAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                    <  2:259286/115‑1 (MQ=255)
                                                                  gAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                                    >  1:259286/1‑115 (MQ=255)
                                                                    agggtTGGGGCTGGAGGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGTGGCTCTATATTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctct                                                                                      <  4:327261/140‑1 (MQ=39)
                                                                    aGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctct                                                                                      >  3:327261/1‑143 (MQ=255)
                                                                      tgtTGGGGCTGTAGGCTGCCGAGGTACAGGGTCTGGGTGAGGCGGGGGTGCTCTAGATCAAACGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcc                                                                               <  4:240489/148‑1 (MQ=38)
                                                                        ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                                                            >  3:344152/1‑149 (MQ=255)
                                                                          ggggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCATTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTTAAGGTTGAAGCGCTCTCGAATCCCCTta                                                                          >  4:239520/1‑149 (MQ=255)
                                                                             gCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata                                                                        >  2:65402/1‑148 (MQ=255)
                                                                             gCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata                                                                        <  1:65402/148‑1 (MQ=255)
                                                                                  aTGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc                                                                  <  4:103255/149‑1 (MQ=255)
                                                                                   tGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcc                                                                 >  7:113316/1‑149 (MQ=255)
                                                                                     cTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCGTTGTTTTTGAGCGCTTTCAAAGGTAATTGAAATTAAAGGTTGACGCCCTATCGAATCCCCTTATAATGCGccc                                                                >  3:61032/1‑148 (MQ=38)
                                                                                      tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCACCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTGAAATTTAACGTTGAAGCCCGCTCGAATCCCCTTATAATGCGccccc                                                              >  8:325115/1‑148 (MQ=39)
                                                                                        gCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACt                                                            <  4:344152/149‑1 (MQ=255)
                                                                                               aCAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTTTAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGAGCCCCACTTCCAGc                                                      >  3:62192/1‑148 (MQ=255)
                                                                                                   cGTCTTCGTGATGGGGGGGTGCCTATTAAGAAAAGCGGTGAGGGCGTGAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAt                                                 <  7:79794/149‑1 (MQ=21)
                                                                                                          gttaTGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                                                      <  1:341015/134‑1 (MQ=255)
                                                                                                          gTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTAAGGGCTATAAAAGGTAATTTAAATTAAAGGTTGACGCCCGCTCGAATGCCCTTATAATGAGCCCCACATCCAga                                                      >  2:341015/1‑136 (MQ=38)
                                                                                                           tGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTATCGAATCCCCGTATAATGCGCCCCACTTCCAGCGACATCGGAACg                                          >  5:7547/1‑148 (MQ=255)
                                                                                                             aTGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGGCGTCAGGCCTACCCGCTTTGTTTTTCAGCGCTTTCAAAGGTAATTGTAATTAAAGGTTGACGGCCTGGGGCATGCCCGTATAATGCGGCCCACTTCGCGCGCCATCGGAACgta                                        >  6:124868/1‑146 (MQ=255)
                                                                                                                      tGCTCTAGATTAAAAGCGGGCAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg                              <  3:152885/149‑1 (MQ=255)
                                                                                                                        ctctAGATGAAAAGCGGGGAGGCCGTCATGCCGACCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGCATCCCCGTAAAATGAGCCCCCCTTACAGCGACATCGGAACGAAAAACTCCtt                               >  3:324003/1‑146 (MQ=25)
                                                                                                                           tAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGAt                          >  1:250525/1‑148 (MQ=255)
                                                                                                                              gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCa                       <  2:30319/149‑1 (MQ=255)
                                                                                                                              gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTGTGTTTTTCACCGCAATCAAAGGTAATTGAAAATAAAGGTTTACGCCCACTCGAATCCCCATATAATGCGCCCCACTTCCAGCGACATAGGAACGAAAAACAACTTGAGAttg                        >  8:183597/1‑147 (MQ=21)
                                                                                                                               aaacaaaGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGCCGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                      >  6:265872/5‑149 (MQ=18)
                                                                                                                               aTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                      <  2:250525/149‑1 (MQ=255)
                                                                                                                                ttaaaaGGGGGCAGGCCGTCAGGCCGCCCCGCTTTTTTCATCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                      <  1:126095/146‑1 (MQ=18)
                                                                                                                                tGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt                     <  7:11399/149‑1 (MQ=255)
                                                                                                                                   aaaGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt                 >  3:66772/1‑149 (MQ=18)
                                                                                                                                   aaaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa                   <  1:114361/148‑1 (MQ=255)
                                                                                                                                    aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt                 <  3:239520/149‑1 (MQ=255)
                                                                                                                                    aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAAGCCCTCTCGAATCCCCTTATAATGCACCCAACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt                 >  2:246827/1‑149 (MQ=255)
                                                                                                                                        ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag             >  4:76201/1‑149 (MQ=255)
                                                                                                                                                    aGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCa  >  8:271753/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGGCA  >  NC_002947/170694‑170988

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