Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 3,065,508 G→A V53M (GTG→ATG)  rcsC_6 → Sensor histidine kinase RcsC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,065,5080GA98.4% 215.6 / NA 64V53M (GTG→ATG) rcsC_6Sensor histidine kinase RcsC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/37);  minor base C (1/0);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTT  >  ACB122/3065380‑3065638
                                                                                                                                |                                                                                                                                  
gggCGGTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACg                                                                                                                              <  7:439104/135‑1 (MQ=255)
        gCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTc                                                                                                                      <  8:112235/135‑1 (MQ=255)
          cGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACTGTATTTTCCa                                                                                                                    <  5:190316/135‑1 (MQ=255)
            ggCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                                  <  4:230169/135‑1 (MQ=255)
            ggCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                                  <  6:426657/135‑1 (MQ=255)
               tgctgGGGCTGGGTGACCACTCGGCGCGTAAAAGCCAATACCCGGAGCTGTCCGCCCGCCTGGACGAACTAGAAGCCGACCGCAACCGTTACAAATGGCTGTTCGAAAACGCCAAGCACGGTATTTTCCAGGccc                                                                                                               >  4:351766/1‑134 (MQ=255)
                gctggggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGCACGCAACCGTTACAAATGGCTGTTCGAAAACGCCCTGCACGGTATTTTCCAGGcccg                                                                                                              >  3:19161/1‑133 (MQ=255)
                gctggggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAg                                                                                                              >  1:369048/1‑135 (MQ=255)
                 ctggggcggggGGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGGTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGc                                                                                                             <  8:445271/135‑1 (MQ=255)
                  tggggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGcc                                                                                                            <  8:341560/135‑1 (MQ=255)
                     ggctgggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAAAGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTc                                                                                                         >  5:221712/1‑135 (MQ=255)
                         gggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCCGGa                                                                                                     >  2:397337/1‑135 (MQ=255)
                         gggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGa                                                                                                     >  1:291598/1‑135 (MQ=255)
                          ggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGAc                                                                                                    <  8:30256/135‑1 (MQ=255)
                          ggTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGAc                                                                                                    >  5:397481/1‑135 (MQ=255)
                                                       cccGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAg                                                                                                       <  8:187023/103‑1 (MQ=255)
                                                       cccGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAg                                                                                                       >  7:187023/1‑103 (MQ=255)
                                                       cccGGAGCTGTCCGCCCGCCGGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGGGCCGCCAACCCGGCGCTGgc                                                                       <  8:277372/135‑1 (MQ=255)
                                                        ccGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGGGCCGCCAACCCGGCGCTGgc                                                                       <  8:482505/134‑1 (MQ=255)
                                                             gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                                 <  6:191941/135‑1 (MQ=255)
                                                             gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                                 <  6:505973/135‑1 (MQ=255)
                                                             gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                                 <  2:369048/135‑1 (MQ=255)
                                                             gCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCAt                                                                 <  2:220126/135‑1 (MQ=255)
                                                               tGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGc                                                               >  3:33906/1‑135 (MQ=255)
                                                               tGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGc                                                               >  6:455581/1‑135 (MQ=255)
                                                               tGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGc                                                               >  5:542942/1‑135 (MQ=255)
                                                               tGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGc                                                               >  4:186332/1‑135 (MQ=255)
                                                                ggccgcccgcGTGGACGAAGTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCt                                                              <  8:340792/133‑1 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAACACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCGGGCGCGCATGCTCgg                                                           >  5:339035/1‑134 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCCCGGTATTTTCCAGGCCAGCCTTCAGGCCGGCATGCGTGCCGCCCGCCCGGCGCTGGCGCGCCTGCTCGGc                                                          >  6:55752/1‑135 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                                  >  7:423652/1‑79 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAgg                                                                                                                  <  8:423652/79‑1 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCgg                                                           >  4:33053/1‑134 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCgg                                                           >  2:435998/1‑134 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCgg                                                           >  6:497645/1‑134 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGATGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCACCCCGGCGCTGGCGCGCATGCTCgg                                                           >  1:103159/1‑134 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGCAACAGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGCAATGCGTGACGCCAACCCGGAGATGGCGCGAATGCTCgg                                                           >  7:29151/1‑134 (MQ=255)
                                                                    gcccgccTGGACGAACTGGAAGCCGAACGAAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACTGCAAGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCgg                                                           >  8:107720/1‑134 (MQ=255)
                                                                       cgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTAt                                                       <  7:353386/135‑1 (MQ=255)
                                                                       cgccTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTGCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGAGGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTAt                                                       <  4:224554/135‑1 (MQ=255)
                                                                              aCGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGAcc                                                <  8:367476/135‑1 (MQ=255)
                                                                                gAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGAccc                                               >  8:315488/1‑134 (MQ=255)
                                                                                 aaCTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGc                                             <  4:229650/135‑1 (MQ=255)
                                                                                 aaCTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGc                                             <  5:198851/135‑1 (MQ=255)
                                                                                 aaCTGCAACCCGAACGCAACCGTTACAAAGGGCTGTTCCAAAACGCCAGGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGGGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGc                                             <  5:11465/135‑1 (MQ=255)
                                                                                       aaGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCgg                                       >  7:268862/1‑135 (MQ=255)
                                                                                       aaGCCGAACGCAACCGTTACAAATGGCTGTGCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCgg                                       <  3:33053/135‑1 (MQ=255)
                                                                                           cGAACGCAACCGTTACAAATGGCGGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCt                                   <  6:135714/135‑1 (MQ=255)
                                                                                              aCGCAACCGTTACAAATGGCTGTGCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGtt                                <  3:308759/135‑1 (MQ=255)
                                                                                               cGCAACCGTTACAAATGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCGCGGCTATGACGACCCGCAAGCGGTGCTGttt                               <  5:497645/135‑1 (MQ=255)
                                                                                                   aCCGTTACTAATGGCGGTTCGAAAACGGCATGCACGGTATTTTCCAGGCCAGCCTGCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTccc                           <  6:504068/135‑1 (MQ=255)
                                                                                                           aactGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCctgactga                    <  2:98555/131‑1 (MQ=255)
                                                                                                           aaaTGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCctgactga                    <  8:347354/134‑1 (MQ=255)
                                                                                                              tGGCTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTg                <  4:115476/135‑1 (MQ=255)
                                                                                                                 cTGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgcc             <  7:315488/135‑1 (MQ=255)
                                                                                                                  tGTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGCCTGACCTGgccg            >  2:437291/1‑135 (MQ=255)
                                                                                                                   gTTCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgccgc           >  1:214860/1‑135 (MQ=255)
                                                                                                                    ttCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgccgcc          <  8:326153/135‑1 (MQ=255)
                                                                                                                    ttCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgccgc           >  4:503950/1‑134 (MQ=255)
                                                                                                                     tCGAAAACGCCATGCACGGTTTTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgccgcc          <  7:285313/134‑1 (MQ=255)
                                                                                                                     tCGAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGgccgcc          <  6:397481/134‑1 (MQ=255)
                                                                                                                       gAAAACGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAAc       >  3:260664/1‑135 (MQ=255)
                                                                                                                          aaCGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTg    <  6:542942/135‑1 (MQ=255)
                                                                                                                            cGCCATGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGtt  <  6:392852/135‑1 (MQ=255)
                                                                                                                                |                                                                                                                                  
GGGCGTTGGCCGGGCTGCTGGGGCTGGGTGACCACTCGGCGCGTAAAAGCCATTACCCGGAGCTGTCCGCCCGCCTGGACGAACTGGAAGCCGAACGCAACCGTTACAAATGGCTGTTCGAAAACGCCGTGCACGGTATTTTCCAGGCCAGCCTTCAGGACGGCATGCGTGCCGCCAACCCGGCGCTGGCGCGCATGCTCGGCTATGACGACCCGCAAGCGGTGCTGTTTTCCCTGACTGACCTGGCCGCCAACCTGTT  >  ACB122/3065380‑3065638

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