Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,855,540 C→A 8.6% H107N (CAT→AAT)  PP_4270 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,855,5400CA 8.6% 174.4 / 6.8 58H107N (CAT→AAT) PP_4270hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base A (2/3);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

GGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGA  >  NC_002947/4855396‑4855685
                                                                                                                                                |                                                                                                                                                    
ggCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATtcatca                                                                                                                                                  <  2:254083/149‑1 (MQ=255)
  cgcTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAgc                                                                                                                                                <  5:213393/149‑1 (MQ=255)
    cTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAgcgg                                                                                                                                              <  8:160091/149‑2 (MQ=255)
      cGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAGCGGGa                                                                                                                                            <  6:189111/149‑1 (MQ=255)
       gCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATGAATCAgcggg                                                                                                                                             <  7:203755/147‑3 (MQ=255)
         cccGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACg                                                                                                                                          >  7:308514/1‑148 (MQ=255)
                              cgcctgcgcctgGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCcgacga                                                                                                                    <  1:229387/149‑1 (MQ=255)
                                      cctgGAGGATGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCTAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCAccc                                                                                                            <  1:14898/149‑1 (MQ=255)
                                          gaggagGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAt                                                                                                        >  1:170603/1‑149 (MQ=255)
                                                   aTGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGa                                                                                               <  8:308514/149‑1 (MQ=255)
                                                          gcggcgATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTg                                                                                        >  6:176463/1‑149 (MQ=255)
                                                              cgATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCAc                                                                                    >  6:52138/1‑149 (MQ=255)
                                                              cgATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCAc                                                                                    >  6:52069/1‑149 (MQ=255)
                                                                  atTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCAcccc                                                                                 <  5:52069/148‑1 (MQ=255)
                                                                       tgtgCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGa                                                                        <  5:52138/149‑1 (MQ=255)
                                                                         tgCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATc                                                                      <  2:198346/149‑1 (MQ=255)
                                                                         tgCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATc                                                                      <  5:176463/149‑1 (MQ=255)
                                                                          gCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCg                                                                     >  7:204741/1‑149 (MQ=255)
                                                                               gCTAAGTCAGTACCCAGTGGCTCAACGCACGACCAAGGGATTTTTTGTGGAGTATACCGCTGATTCATAAGCGCGACGTACTGGTTGGCCCCCCCCGCGACAGCCACCCGCTGCCCCCCGAGCAACTGCCACACCAA‑‑‑GCGATCGGCa                                                                  >  4:294712/1‑147 (MQ=255)
                                                                                    gTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGc                                                                   <  6:4527/141‑1 (MQ=255)
                                                                                    gTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGc                                                                   >  5:4527/1‑141 (MQ=255)
                                                                                      cAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTCCTgg                                                         >  8:266530/1‑149 (MQ=255)
                                                                                       aGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTggg                                                        <  2:170603/149‑1 (MQ=255)
                                                                                            ccAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTaca                                                   >  1:138183/1‑149 (MQ=255)
                                                                                              aGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACAcc                                                 >  4:160094/1‑149 (MQ=255)
                                                                                                  gCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGCGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTAtt                                             >  2:320530/1‑149 (MQ=255)
                                                                                                   cTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTc                                            <  4:157333/149‑1 (MQ=255)
                                                                                                    tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCg                                                                                                                        >  5:255429/1‑75 (MQ=255)
                                                                                                    tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCg                                                                                                                        <  6:255429/75‑1 (MQ=255)
                                                                                                    tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTc                                            >  5:59719/1‑148 (MQ=255)
                                                                                                    tCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTc                                            >  2:107327/1‑148 (MQ=255)
                                                                                                         gACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGc                                                       >  2:353561/1‑132 (MQ=255)
                                                                                                         gACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGc                                                       <  1:353561/132‑1 (MQ=255)
                                                                                                           ccGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTg                                     >  2:225666/1‑148 (MQ=255)
                                                                                                            cGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGc                                    >  6:324501/1‑148 (MQ=255)
                                                                                                                cAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCa                               <  1:225666/149‑1 (MQ=255)
                                                                                                                cAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCa                               <  2:138183/149‑1 (MQ=255)
                                                                                                                 aaGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGACGCCCCGGCGACAGCACCCCGATGAACCCGGAGCACCTGCCACCCCAG‑‑‑GAGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCAc                              >  5:158447/1‑149 (MQ=255)
                                                                                                                  aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAg                                                                                                                                                 >  4:130670/1‑36 (MQ=255)
                                                                                                                  aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAg                                                                                                                                                 <  3:130670/36‑1 (MQ=255)
                                                                                                                  aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAt                                                                                                        >  7:87783/1‑77 (MQ=255)
                                                                                                                  aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAt                                                                                                        <  8:87783/77‑1 (MQ=255)
                                                                                                                  aGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCAct                             <  8:204741/149‑1 (MQ=255)
                                                                                                                  aGGGAATTTGTGTGGATTATCCCGCTGATTCATCAGATCGACGTACTGGTTGTCCGACCCGGCGACAGCACCCCGATGCCCCCGGACCAACTGCCACCCCAT‑‑‑GCCAACGGCACGGTGCAGGGCTACACCTATTCAACCCTGCAACCAct                             >  8:94588/1‑149 (MQ=255)
                                                                                                                   gggATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActg                            >  7:183332/1‑149 (MQ=255)
                                                                                                                   gggATTTTGTGTGGAGAATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActg                            >  4:180967/1‑149 (MQ=255)
                                                                                                                    ggATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActgc                           <  6:59719/149‑1 (MQ=255)
                                                                                                                      aTTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCActgctg                         <  8:183332/149‑1 (MQ=255)
                                                                                                                              tgGAGTATCGCGCTGATTCATCATCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCTATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTTCAACCACTGCTGGACCACg                  <  6:158447/148‑1 (MQ=255)
                                                                                                                                gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGc                >  7:95384/1‑148 (MQ=255)
                                                                                                                                gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCACCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGc                >  6:28073/1‑148 (MQ=255)
                                                                                                                                gAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGAAGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCAATTCAACCCTGCAACCACTGCTCGACCACGGc                >  6:22833/1‑148 (MQ=255)
                                                                                                                                   tATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGt            <  7:266530/149‑1 (MQ=255)
                                                                                                                                    gtcgcGCTGATGAATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGtt           >  5:328770/5‑149 (MQ=255)
                                                                                                                                       gcGCTGATGAATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCa        >  4:226823/2‑149 (MQ=255)
                                                                                                                                       ccGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCa        >  8:157835/1‑149 (MQ=255)
                                                                                                                                         gCTGATTCATCAGCGCGACGTACAGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAGCTACCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGc         <  1:107327/149‑1 (MQ=255)
                                                                                                                                           tGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGt    >  5:350231/1‑149 (MQ=255)
                                                                                                                                             aTTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGa  <  3:160094/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGCGCTCGCCCCGAGTTCCACGTCATGGCCCGCCTGCGCCTGGAGGAGGCCATGAGCAGCGGCGATATTGATGTGCGTTGCTACGTCAGTACCCAGTGGCTCAACGACCGTCCAAGGGATTTTGTGTGGAGTATCCCGCTGATTCATCAGCGCGACGTACTGGTTGGCCGCCCCGGCGACAGCACCCCGATGCCCCCGGAGCAACTGCCACCCCAG‑‑‑GCGATCGGCACGGTGCTGGGCTACACCTATTCAACCCTGCAACCACTGCTGGACCACGGCCAGTTGCACCGTGA  >  NC_002947/4855396‑4855685

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