Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,361,790 C→G 7.3% intergenic (+25/‑80) PP_t45 → / → yphA tRNA‑Leu/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,361,7900CG 7.3% 168.8 / 4.2 55intergenic (+25/‑80)PP_t45/yphAtRNA‑Leu/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base G (2/2);  total (28/27)
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.20e-01

AGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGATGGTC  >  NC_002947/4361653‑4361938
                                                                                                                                         |                                                                                                                                                    
aGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCttttt                                                                                                                                           <  4:175264/149‑1 (MQ=255)
    ggTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTtgctg                                                                                                                                       >  6:72475/1‑149 (MQ=255)
     gTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTtgctgc                                                                                                                                      >  5:111683/1‑149 (MQ=255)
       cGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCt                                                                                                                                    >  2:178524/1‑149 (MQ=255)
         ttGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGa                                                                                                                                  <  5:41069/149‑1 (MQ=255)
            cGACCCCTTCAAGCGTGCGAGTGTGGTGTAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAg                                                                                                                               <  6:228191/149‑1 (MQ=255)
              aCCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                         >  5:321173/1‑134 (MQ=255)
              aCCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                         <  6:321173/137‑4 (MQ=255)
                cccTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAAt                                                                                                                           >  2:92890/1‑149 (MQ=255)
                  cTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                         >  5:84427/1‑130 (MQ=255)
                  cTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                         <  6:84427/133‑4 (MQ=255)
                   ttCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTc                                                                                                                        >  7:106199/1‑149 (MQ=255)
                        gCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGgcg                                                                                                                   >  4:4151/1‑149 (MQ=255)
                          gTGCGAGTGTGTTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGgcgc                                                                                                                  <  5:249084/148‑1 (MQ=255)
                              gAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAgct                                                                                                              >  1:131090/1‑148 (MQ=255)
                               aGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAgctgc                                                                                                            <  4:234591/149‑1 (MQ=255)
                               aGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAgctgc                                                                                                            <  5:276972/149‑1 (MQ=255)
                                    ggtggATTTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAg                                                                                                       >  4:15029/1‑149 (MQ=255)
                                     gtggAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATCTCTGCCGCAGCTCCTAAAGt                                                                                                      >  4:162431/1‑149 (MQ=255)
                                      tggATTTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTg                                                                                                     <  3:15029/149‑1 (MQ=255)
                                        gAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGa                                                                                                   <  5:72475/149‑1 (MQ=255)
                                        gAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGa                                                                                                   <  8:106199/149‑1 (MQ=255)
                                          aTTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGc                                                                                                 >  2:78720/1‑149 (MQ=255)
                                                ataCACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTATTGCTGTCAGGCAGCaaaaa                                                                                                                         <  5:190093/119‑4 (MQ=255)
                                                ataCACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTATTGCTGTCAGGCAGCaaaaa                                                                                                                         >  6:190093/1‑116 (MQ=255)
                                                  acacAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTc                                                                                            >  3:258211/1‑146 (MQ=255)
                                                  acacAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTc                                                                                            <  4:258211/146‑1 (MQ=255)
                                                   cacaGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCg                                                                                        >  2:317841/1‑149 (MQ=255)
                                                     cagcagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAg                                                                                                                               <  6:101900/108‑1 (MQ=255)
                                                     cagcagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAg                                                                                                                               >  5:101900/1‑108 (MQ=255)
                                                      agcagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTc                                                                                      >  5:71537/1‑148 (MQ=255)
                                                        cagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCAc                                                                                   <  5:302973/149‑1 (MQ=255)
                                                          gACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTg                                                                                 <  3:4151/149‑1 (MQ=255)
                                                           aCTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTg                                                                                 >  7:245375/1‑148 (MQ=255)
                                                                          gCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAAcgc                                                                           >  1:141741/1‑139 (MQ=255)
                                                                          gCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAAcgc                                                                           <  2:141741/139‑1 (MQ=255)
                                                                                gCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTcc                                                           <  6:111683/149‑1 (MQ=255)
                                                                                     gCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGtt                                                      >  7:262341/1‑149 (MQ=255)
                                                                                      cGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGtt                                                      <  8:245375/148‑1 (MQ=255)
                                                                                             gAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGt                                               <  6:71537/148‑1 (MQ=255)
                                                                                               gTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTa                                                               >  5:15881/1‑130 (MQ=255)
                                                                                               gTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTa                                                               <  6:15881/130‑1 (MQ=255)
                                                                                                  ccGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAAcgc                                         >  1:4508/1‑149 (MQ=255)
                                                                                                         tCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGAcatca                                  >  5:198861/1‑149 (MQ=255)
                                                                                                           gCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGAcatca                                  <  1:92890/147‑1 (MQ=255)
                                                                                                              ccAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATc                             >  3:350528/1‑149 (MQ=255)
                                                                                                              ccAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATc                             <  1:178524/149‑1 (MQ=255)
                                                                                                              ccAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATc                             <  1:317841/149‑1 (MQ=255)
                                                                                                               cAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATcc                            >  2:323701/1‑149 (MQ=255)
                                                                                                                      cGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATc                       <  1:323701/147‑1 (MQ=255)
                                                                                                                            aaaGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCCTGCGGTACTCCCTGATCGATGGTCAACGCGACATCATTATCCTGATCGCCcgcg                >  3:257715/1‑148 (MQ=255)
                                                                                                                            aaaGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCcgcg                >  3:124300/1‑148 (MQ=255)
                                                                                                                                gcgcCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTg           <  2:4508/149‑1 (MQ=255)
                                                                                                                                gcgcCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTg           >  4:172500/1‑149 (MQ=255)
                                                                                                                                    ccTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGa       <  8:262341/149‑1 (MQ=255)
                                                                                                                                        gcggcgCCTTTTTGCTGCCTGACAGCTATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGATg     <  6:198861/147‑1 (MQ=255)
                                                                                                                                         cggcgCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGATGGTc  <  3:172500/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
AGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGATGGTC  >  NC_002947/4361653‑4361938

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