Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,659,909 T→C 9.0% intergenic (‑44/+31) PP_3223 ← / ← PP_3224 ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,659,9090TC 9.0% 150.2 / 4.0 55intergenic (‑44/+31)PP_3223/PP_3224ABC transporter substrate‑binding protein/class II aldolase/adducin domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (2/3);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659761‑3660055
                                                                                                                                                    |                                                                                                                                                  
ggCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGt                                                                                                                                                    >  4:17871/1‑149 (MQ=255)
          aGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAAc                                                                                                                                          <  7:254726/149‑1 (MQ=255)
          aGGGGCTTGCTGGTTGAGGGGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCGAACTGGCCGGaa                                                                                                                                           <  2:119972/148‑1 (MQ=255)
            ggTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGGGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGGGCAGTGGGGGCATGGGATTGCTCCTTCAATGAAGGGGGCTCTCCTGTTCCGGCGACCTGGCCGGAACAg                                                                                                                                        <  1:83163/149‑1 (MQ=255)
             gTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAgg                                                                                                                                       <  2:155057/149‑1 (MQ=255)
              tCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGc                                                                                                                                      <  7:102481/149‑1 (MQ=255)
                      gTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTg                                                                                                                              >  2:91876/1‑149 (MQ=255)
                       ttGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGc                                                                                                                             <  7:261670/149‑1 (MQ=255)
                        tGAGGTCTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGACTGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGcc                                                                                                                            >  8:5214/1‑149 (MQ=255)
                           ggTGTTGTCGGCCGGCGAGGAAGCGGTGAGGCGCAGGGCCATGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCtt                                                                                                                          <  5:18289/148‑1 (MQ=255)
                            gtgtTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAATGTGCCTTg                                                                                                                         <  1:367464/148‑1 (MQ=255)
                                 gccggccggcGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGagt                                                                                                                     >  3:288418/3‑147 (MQ=255)
                                  tCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGagtcag                                                                                                                  >  6:98622/1‑149 (MQ=255)
                                  tCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGagtca                                                                                                                   <  2:103026/148‑1 (MQ=255)
                                                aGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    >  1:48518/1‑149 (MQ=255)
                                                aGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAg                                                                                                    <  3:26437/149‑1 (MQ=255)
                                                 gCGGTGCGGCCCAGGGCCAGGGGGCCGATCGCGAGCAGTTGGTGCAGTCGGGGCATGGGATTGATCCTTCAATGAATCGGGCTCTCCTGTTCCGGCCAGCTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAgg                                                                                                   <  4:288418/149‑1 (MQ=255)
                                                          cccAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGa                                                                                          >  2:74831/1‑149 (MQ=255)
                                                           ccagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGAt                                                                                         >  1:12382/1‑149 (MQ=255)
                                                            cagggccagggTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATc                                                                                        <  4:306205/149‑1 (MQ=255)
                                                                          cGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAAt                                                                            >  5:63273/1‑147 (MQ=255)
                                                                          cGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAAt                                                                            >  4:361889/1‑147 (MQ=255)
                                                                            aTCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAgc                                                                        >  5:191565/1‑149 (MQ=255)
                                                                                                   gggCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATAcc                                                 <  1:91876/149‑1 (MQ=255)
                                                                                                     gCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCtt                                               <  2:48518/149‑1 (MQ=255)
                                                                                                           taTTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTACTGGTAATCGAGTCCCTCTCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCt                                         >  5:308888/2‑149 (MQ=255)
                                                                                                           gATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTc                                                                                              >  5:377972/1‑96 (MQ=255)
                                                                                                           gATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTc                                                                                              <  6:377972/96‑1 (MQ=255)
                                                                                                            aTTGATCATTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTg                                        >  3:131313/1‑149 (MQ=255)
                                                                                                                aTCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAAt                                                              <  1:335189/123‑1 (MQ=255)
                                                                                                                aTCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAAt                                                              >  2:335189/1‑123 (MQ=255)
                                                                                                                aTCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTg                                                       <  2:54459/130‑1 (MQ=255)
                                                                                                                aTCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTg                                                       >  1:54459/1‑130 (MQ=255)
                                                                                                                 tCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCAt                                   >  3:358007/1‑149 (MQ=255)
                                                                                                                   cTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCg                                 <  3:17871/149‑1 (MQ=255)
                                                                                                                   cTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCg                                 <  1:74831/149‑1 (MQ=255)
                                                                                                                    ttCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGACGTTGGATACATTCCCCCTGGCCATCg                                 >  5:242536/1‑148 (MQ=255)
                                                                                                                      cAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGa                               >  6:302442/1‑148 (MQ=255)
                                                                                                                       aatgaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGaa                              >  8:116745/1‑148 (MQ=255)
                                                                                                                       aatgaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGaa                              >  3:46877/1‑148 (MQ=255)
                                                                                                                          gaatgGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAAcc                          >  8:278857/1‑149 (MQ=255)
                                                                                                                           aatgGGGCTCTCCTGTTCCGGACAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAAcc                          <  4:131313/148‑1 (MQ=255)
                                                                                                                               gggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTTCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCTCTTg                     <  4:358007/149‑1 (MQ=255)
                                                                                                                               gggCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCtt                      <  5:98622/148‑1 (MQ=255)
                                                                                                                                 gCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGa                   <  2:12382/149‑1 (MQ=255)
                                                                                                                                   tctcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATc                                  <  3:173284/132‑1 (MQ=255)
                                                                                                                                   tctcCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATc                                  >  4:173284/1‑132 (MQ=255)
                                                                                                                                      ccTGTTCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGc               >  6:358864/1‑148 (MQ=255)
                                                                                                                                      ccTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCt              >  5:46203/1‑149 (MQ=255)
                                                                                                                                       cTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTc             <  7:116745/149‑1 (MQ=255)
                                                                                                                                         gTTCCGGCGAACTGGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCAt           >  1:349060/1‑149 (MQ=255)
                                                                                                                                              ggCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCtt                      >  2:366684/1‑133 (MQ=255)
                                                                                                                                              ggCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCtt                      <  1:366684/133‑1 (MQ=255)
                                                                                                                                               gCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCt     >  4:305181/1‑149 (MQ=255)
                                                                                                                                                  aGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCGTGCACCTTGa  <  6:191565/149‑1 (MQ=255)
                                                                                                                                                  aGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGa  <  7:278857/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCGATGTTGAGGGTCTTGCTGGTTGAGGTGTTGTCGGCCGGCGAGCAAGCGGTGAGGCCCAGGGCCAGGGTGCCGATCGCGAGCAGTTGGTGCAGTGGGGGCATGGGATTGATCCTTCAATGAATGGGGCTCTCCTGTTCCGGCCAGTTCGCCGGAACAGGCAAAAGGTGCCTTGGAGTCAGTCCTGGTAATCGAGGCCCTCGCGATCCAACTGCCGAATCAGCGCATCCCAATGCCGTTGGATACCTTCCCCCTGGCCATCGGAAAACCGCTTGGACTGCTCATGCACCTTGA  >  NC_002947/3659761‑3660055

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: ‑