Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,087,682 C→T 8.7% intergenic (+133/‑97) oruR → / → aphA ornithine utilization transcriptional regulator/acetylpolyamine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,087,6820CT 8.7% 161.9 / 4.3 58intergenic (+133/‑97)oruR/aphAornithine utilization transcriptional regulator/acetylpolyamine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base T (2/3);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

CCATTTTCGTGGTTGACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGCACCACGGCCGCTGC  >  NC_002947/6087534‑6087826
                                                                                                                                                    |                                                                                                                                                 
ccATTTTCGTGGTTGACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGAcc                                                                                                                                                   >  6:322742/1‑149 (MQ=255)
        gTGGTTGACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAGGTCTTTGGccgg                                                                                                                                           <  2:201957/149‑3 (MQ=255)
         tGGTTGACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCata                                                                                                                                          <  4:177457/149‑1 (MQ=255)
          ggTTGACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCatat                                                                                                                                         >  8:242256/1‑149 (MQ=255)
           gTTGACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAGGTCTTTGGccggat                                                                                                                                         <  3:294709/148‑5 (MQ=255)
              gACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAGGTCTTTGGccggatggg                                                                                                                                      <  5:314861/148‑8 (MQ=38)
                ccTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATc                                                                                                                                    <  8:284148/148‑1 (MQ=255)
                 cTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTGACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATAAATgcc                                                                                                                                  >  6:252712/1‑149 (MQ=255)
                      cTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCttt                                                                                                                              <  2:56481/148‑1 (MQ=255)
                         gCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCtttt                                                                                                                             <  7:38688/146‑1 (MQ=255)
                                    ggTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCCTCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGcgtcg                                                                                                               >  7:14113/1‑149 (MQ=255)
                                                    cAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGAt                                                                                                <  3:234726/148‑1 (MQ=255)
                                                     aGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCa                                                                                             <  8:194715/149‑1 (MQ=255)
                                                       gggCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGc                                                                                           >  5:128303/1‑149 (MQ=255)
                                                        ggCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGc                                                                                           <  7:107714/148‑1 (MQ=255)
                                                              acaACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCa                                                                                        <  2:285325/145‑1 (MQ=255)
                                                              acaACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCAACCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCa                                                                                        >  1:285325/1‑145 (MQ=255)
                                                                aaCCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACaa                                                                                  <  7:71854/149‑1 (MQ=255)
                                                                aaCCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACaa                                                                                  >  3:98252/1‑149 (MQ=255)
                                                                 aCCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAAt                                                                                 <  6:128303/149‑1 (MQ=255)
                                                                   cccTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCttt                                                                                                                              >  6:14566/1‑103 (MQ=255)
                                                                   cccTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCttt                                                                                                                              <  5:14566/103‑1 (MQ=255)
                                                                    ccTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGa                                                                              <  8:155202/149‑1 (MQ=255)
                                                                      tGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTgg                                                                                                                           >  8:94495/1‑103 (MQ=255)
                                                                       gACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGAcc                                                                            >  1:176187/1‑148 (MQ=255)
                                                                        aCGCCCTTGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCTGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGAcc                                                                            <  5:252712/147‑1 (MQ=255)
                                                                           cccTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGc                                                                       <  7:242256/149‑1 (MQ=255)
                                                                            ccTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGc                                                                       >  7:213690/1‑148 (MQ=255)
                                                                                gtggCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTc                                                                  <  2:67771/149‑1 (MQ=255)
                                                                                         gTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGga                                                          <  8:137795/148‑1 (MQ=255)
                                                                                              aGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGaacaa                                                    <  8:213690/149‑1 (MQ=255)
                                                                                                  ggTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAAt                                                >  2:4088/1‑149 (MQ=255)
                                                                                                  ggTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAAt                                                <  5:322742/149‑1 (MQ=255)
                                                                                                  ggTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGACATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAAt                                                <  8:149734/149‑1 (MQ=255)
                                                                                                   gTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAAt                                                >  2:284490/1‑148 (MQ=255)
                                                                                                        accGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCTTGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGa                                           <  1:80292/148‑1 (MQ=255)
                                                                                                          cgcgAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGa                                        >  4:307359/1‑149 (MQ=255)
                                                                                                           gcgAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGa                                        >  6:230992/1‑148 (MQ=255)
                                                                                                           gcgAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGAt                                       <  2:276870/149‑1 (MQ=255)
                                                                                                             gAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATTACAATGGACCCAAGCCAGTCCAACCGGATAACAACAAAGCTGAcagaac                                     >  6:149328/1‑144 (MQ=255)
                                                                                                             gAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGa                                           <  5:149328/143‑1 (MQ=255)
                                                                                                                  ggCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTcc                                >  4:49622/1‑149 (MQ=255)
                                                                                                                    ccGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATGTATTCCGa                              >  2:26061/1‑149 (MQ=255)
                                                                                                                    ccGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGa                              >  5:116414/1‑149 (MQ=255)
                                                                                                                     cGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGAt                             <  8:14113/149‑1 (MQ=255)
                                                                                                                     cGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGAt                             <  5:16225/149‑1 (MQ=255)
                                                                                                                      ggATAGGTCACCACCCATCCGGCCAAAGACATTGGCCAAATAAATCCCCATTTTCCAGTTTGCGTCGTTCTCCCATGCTGATCAGCCCATGACAATGGACCCACTCCAGTCCCACCGAAGAACAACAATGATCACGATCTATTACGAt                             >  8:171991/1‑148 (MQ=255)
                                                                                                                        aTAGGTCACCACCCATCCGGCCAAAGACCTTGGCTATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGac                          >  2:34026/1‑149 (MQ=255)
                                                                                                                        aTAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGac                          <  2:176187/149‑1 (MQ=255)
                                                                                                                               accaccCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGc                    <  6:116414/148‑1 (MQ=255)
                                                                                                                                     cATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTc                                                                  <  7:188502/96‑1 (MQ=255)
                                                                                                                                     cATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTc                                                                  >  8:188502/1‑96 (MQ=255)
                                                                                                                                       tatggCCAAGGTCTTTGGACATATCAATCACCTTTTGGCCGTTTGCGTCGTTGTCCCATGCTGATCAGCCCATGACAATGGACACACGACAGTCCCACCGGAGAACAAAAATGATGACGATCTATTCCGATGACAACCGCCTGcaccac           >  8:20694/4‑149 (MQ=255)
                                                                                                                                       tCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGcaccac           >  1:224207/1‑149 (MQ=255)
                                                                                                                                        atggCCAAGGTCTTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGCACCACg          >  1:294861/3‑149 (MQ=255)
                                                                                                                                            ccAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGCACCACGGCCg      <  2:224207/149‑1 (MQ=255)
                                                                                                                                             cAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGCACCACGGCCGc     <  1:284490/149‑1 (MQ=255)
                                                                                                                                               aaGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGCACCACGGCCGCTg   >  3:38929/1‑149 (MQ=255)
                                                                                                                                                 gACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGCACCACGGCCGCTGc  >  2:306711/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCATTTTCGTGGTTGACCTGTTCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGGGCGTCACAACCCCTGACGCCCTGGTGGCCATTGTAGGAGCGGGTTCACCCGCGAAAAGGCCGGATAGGTCACCACCCATCCGGCCAAAGACCTTGGCCATATCAATCCCCTTTTGGCCGTTTGCGTCGTTCTCCCCTGCTGATCAGCCCATGACAATGGACCCACGCCAGTCCCACCGGAGAACAACAATGCTGACGATCTATTCCGATGACCACCGCCTGCACCACGGCCGCTGC  >  NC_002947/6087534‑6087826

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: