Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,390,205 G→T 60.3% G34V (GGG→GTG)  hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,390,2050GT60.3% 27.4 / 53.8 55G34V (GGG→GTG) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base G (13/9);  new base T (14/19);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

AATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTA  >  CP009273/4390059‑4390352
                                                                                                                                                  |                                                                                                                                                   
aaTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                                                                                                                                                   >  1:375793/1‑149 (MQ=255)
aaTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                                                                                                                                                   >  1:376142/1‑149 (MQ=255)
   cGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCaaa                                                                                                                                                >  2:572087/1‑149 (MQ=255)
               aGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTg                                                                                                                                    <  2:43624/149‑1 (MQ=255)
                       tAAACATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTCAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTc                                                                                                                            <  2:375731/149‑1 (MQ=255)
                            atatAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCATTGCGTCGGGAACGTGTTCCAGTTTATATTTATTTGGTGAATGGTATTAAGCTGCAAGTTCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                                                                                       <  2:563934/149‑1 (MQ=255)
                             tataAGGAACAGACAGAATGGCGAAGGGGCAATCTTTATAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATc                                                                                                                      <  2:862418/149‑1 (MQ=255)
                             tataAGGAAACGGGAGAATTGCCAAGGGGCAAGCTTTACAAGGTCCTTTCCTGAACGTACTGGGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATc                                                                                                                      <  1:476107/149‑1 (MQ=255)
                             tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATc                                                                                                                      <  1:93598/149‑1 (MQ=255)
                                     aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGGACGCACTGCGTCGGGAACGTGTCCCCGTCTCTGTTTATTTGGTGCATGGTATTAAACTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtttgt                                                                                                              >  2:735771/1‑146 (MQ=255)
                                     aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGAGCCGTTCCTGAACGCATTACGTCGGGAACGTGGTCCAGTCTCTAATTATTTGGTGCATGGGATTAAGCTGCAAGTGGAAATCGAGTCTTTTGATCAGTTCGTgcaccggtttaa                                                                                                              >  1:122052/1‑138 (MQ=255)
                                      aaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGaaa                                                                                                             >  2:560329/1‑149 (MQ=255)
                                            gaATGGCTAAGGGGCAATCTGTACAAGATCCGTTCCTGAACGCACGGCGTCGGGAACGTGTTACAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACg                                                                                                       <  2:589395/149‑1 (MQ=255)
                                             aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAGTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                                                                                      >  2:619012/1‑149 (MQ=255)
                                             aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGGTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAACCACgg                                                                                                      >  2:815652/1‑149 (MQ=255)
                                                             tCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                                                                                      >  2:176730/1‑149 (MQ=255)
                                                             tCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                                                                                      >  2:363746/1‑149 (MQ=255)
                                                              cTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                                                                                     <  1:842396/149‑1 (MQ=255)
                                                               tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                                                                                    >  2:353197/1‑149 (MQ=255)
                                                                ttACAAGATCCGTTCCTGGACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACaa                                                                                   <  1:127425/149‑1 (MQ=255)
                                                                ttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAACACACGGCCAGCCAGATGGTTTACaa                                                                                   >  2:826338/1‑149 (MQ=255)
                                                                      gATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgc                                                                             >  2:787909/1‑149 (MQ=255)
                                                                           gTTCCTGAACTCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAttt                                                                        <  1:497465/149‑1 (MQ=255)
                                                                           gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAttt                                                                        <  2:153192/149‑1 (MQ=255)
                                                                           gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAttt                                                                        <  2:376142/149‑1 (MQ=255)
                                                                           gTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAttt                                                                        <  2:375793/149‑1 (MQ=255)
                                                                               cTGAACGCACTGCGTCGGGAACGTGGTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACGCGGCCAGCCAGATGGTTTACAAGCACGCGATTTCTAc                                                                    <  2:553802/149‑1 (MQ=255)
                                                                                   aCGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtgtt                                                                <  1:572087/149‑1 (MQ=255)
                                                                                          gCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGt                                                         <  2:39773/149‑1 (MQ=255)
                                                                                                 gAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGccc                                                  >  2:860980/1‑149 (MQ=255)
                                                                                                   acctgtTCCAGTTTCCATTTATTTGGTCAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTCGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCgg                                                <  2:338618/146‑1 (MQ=255)
                                                                                                        ttCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTtct                                           >  2:804293/1‑149 (MQ=255)
                                                                                                           cactttCTAGTGCTTTGGCTAATGGTATTAAGCTGCAAGGGCAACTCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACCAGCACGCGATCTATACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                        <  2:13572/146‑1 (MQ=255)
                                                                                                           cAGTTTCTATTTATTTGGTGAACGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcat                                        <  2:627632/149‑1 (MQ=255)
                                                                                                             gTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCtcatca                                      >  2:313066/1‑149 (MQ=255)
                                                                                                               ttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                    <  1:342702/149‑1 (MQ=255)
                                                                                                               ttCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                    <  1:749054/149‑1 (MQ=255)
                                                                                                               tgcTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATcaca                                    <  1:775811/147‑1 (MQ=255)
                                                                                                                tCTATTTATTTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAg                                   <  1:860980/149‑1 (MQ=255)
                                                                                                                  tatttattTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCATCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCCTCACAGTa                                 >  1:207724/1‑149 (MQ=255)
                                                                                                                  tatttattTGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTTTCCCGTCTCGCCCGGTTTCtcatcaacgta                                 >  1:501330/1‑144 (MQ=255)
                                                                                                                    tttatttGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCACGGTTTCTCATCACAGTaac                               >  1:285569/1‑149 (MQ=255)
                                                                                                                       atttGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTaacaac                            <  2:695958/149‑1 (MQ=255)
                                                                                                                         ttGGTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCCGTTCGTGATCCTGTTGACAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGGTGCGCATCACAATaacaacac                          >  1:462417/1‑147 (MQ=255)
                                                                                                                            gTGAATGGTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCgg                       <  1:796700/149‑1 (MQ=255)
                                                                                                                             tGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGt                      >  2:660040/1‑149 (MQ=255)
                                                                                                                               aaTGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTgg                    >  1:87310/1‑149 (MQ=255)
                                                                                                                                  ggTATTAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCgg                 >  2:578141/1‑149 (MQ=255)
                                                                                                                                   gTATTAAGCTGCAAGGGCAAATCGAGTCTGTTGATCAGTTCGTGATCCTGTTGGAAAACACGGTCAGACAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCCGTTTCTCGTCACAGTAAATACGCCGGTGGCGGt                >  2:530450/1‑149 (MQ=255)
                                                                                                                                   gTATTAAGCAGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGCCAGACAGATGGGTTACAAGCACGCGATTTCTACCGTCGTCCCATCCCGCCCGCATTCTCACCACAGTAACAACGCCGGTGGccga                >  1:444943/1‑146 (MQ=255)
                                                                                                                                       tAAGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACca            >  1:410405/1‑149 (MQ=255)
                                                                                                                                        aaGCTGCAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcag           <  1:176730/149‑1 (MQ=255)
                                                                                                                                          gCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACcagca         <  1:560329/149‑1 (MQ=255)
                                                                                                                                           ctcgaaGTGCAAATGGAGTCTTTTGATCAGTTCGTGATCCTGCTGAACCGTTGGGTCGGCCAGATGGTTTACAAGCACGCGATGTCTACGGTTGTGCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCTGTGGCGGTACcagcag        <  1:613233/145‑1 (MQ=255)
                                                                                                                                              cAAGTGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAGAAGCACGCGATTTCTACTGTTGTCCCGACTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTaa     >  2:643194/1‑149 (MQ=255)
                                                                                                                                                 gggCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCCCGCGATGTCTCCGGTTGTCCCGTCTCGCCCGGTTTCCCATCACAGTAACACCGCCGGTGGCGGTAcccgacctaactc  >  1:368968/1‑138 (MQ=255)
                                                                                                                                                 gggCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTGCCAGCAGTAACTa  >  2:120219/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCCGTCTCGCCCGGTTTCTCATCACAGTAACAACGCCGGTGGCGGTACCAGCAGTAACTA  >  CP009273/4390059‑4390352

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