Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 335,804 T→G 6.8% noncoding (31/89 nt) PP_mr04 → C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947335,8040TG 6.8% 174.3 / 4.6 59noncoding (31/89 nt)PP_mr04C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base G (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

ATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACA  >  NC_002947/335664‑335945
                                                                                                                                            |                                                                                                                                             
ctcctccGCGCGCCGCCGAAGGACCGCTGGTGGCGGCGCAGGGGCGCACGGATACCGGTCTCGGGTCGCGCTGGCGACGGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCCTGGTGGTCGTACGTGGGCAGACCTCGGTCTgg                                                                                                                                       <  8:62978/148‑1 (MQ=255)
   cTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGcc                                                                                                                                     <  1:230022/148‑1 (MQ=255)
         gcgcCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGCCAGACCGAGGTCTGccca                                                                                                                                    >  2:227835/1‑142 (MQ=255)
         gcgcCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCCTTACCCGCCATGGTGGTCGTACGTGGCCAGACCGAGGTCTGccca                                                                                                                                    <  1:227835/143‑2 (MQ=255)
          cgccgcCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTcc                                                                                                                             <  4:222488/149‑1 (MQ=255)
          cgccgcCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTcc                                                                                                                             <  5:160302/149‑1 (MQ=255)
           gccgccGAAGTAACGCTGTTGTCGGGGAAGTGGCTCCCGGATACCGTTTTCGGATCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCt                                                                                                                            >  5:226667/1‑149 (MQ=255)
               ccgAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTAtt                                                                                                                        >  8:160070/1‑149 (MQ=255)
                                tCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCAcc                                                                                                        <  7:87398/148‑1 (MQ=255)
                                  ggggAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGCCAGACCGAGGTCTGccca                                                                                                                                    <  5:289497/118‑2 (MQ=255)
                                  ggggAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGCCAGACCGAGGTCTGccca                                                                                                                                    >  6:289497/1‑117 (MQ=255)
                                          ggCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGc                                                                                             >  7:29175/1‑149 (MQ=255)
                                           gCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCt                                                                                            >  1:89345/1‑149 (MQ=255)
                                           gCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCt                                                                                            >  4:301046/1‑149 (MQ=255)
                                             tCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCg                                                                                                                                             <  3:255871/98‑1 (MQ=255)
                                             tCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCg                                                                                                                                             >  4:255871/1‑98 (MQ=255)
                                               aCGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTgccg                                                                                        >  4:255362/1‑149 (MQ=255)
                                               aCGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTgccg                                                                                        <  6:253703/149‑1 (MQ=255)
                                                    tACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGc                                                                                                    >  7:288170/1‑132 (MQ=255)
                                                    tACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGc                                                                                                    <  8:288170/132‑1 (MQ=255)
                                                         ttttCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATcc                                                                              >  5:292402/1‑149 (MQ=255)
                                                           ttCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAACTCCGt                                                                            >  2:258773/1‑149 (MQ=255)
                                                              gggTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTgcc                                                                                         <  6:212374/133‑1 (MQ=255)
                                                              gggTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTgcc                                                                                         >  5:212374/1‑133 (MQ=255)
                                                              gggTCGCGCTGGCGACTGGGGGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTga                                                                          <  6:226667/148‑1 (MQ=255)
                                                                    ctctGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGAt                                                                   >  7:2918/1‑149 (MQ=255)
                                                                         gCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGAttttt                                                               <  5:67378/148‑1 (MQ=255)
                                                                                ggTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGAACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCAATCGGATTCTCACCCTGCGTACGGCTCACGCCCAAAATCCGTGACAAACAATTTTCGCAGtct                                                       >  8:227425/1‑149 (MQ=255)
                                                                                 gtgtGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGTCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGtctc                                                      <  2:283217/149‑1 (MQ=255)
                                                                                 gtgtGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGtctc                                                      >  8:225214/1‑149 (MQ=255)
                                                                                  tgtgAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGc                                                                                                    <  4:251706/102‑1 (MQ=255)
                                                                                  tgtgAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTAACTCGGTTTCTCACCCTGc                                                                                                    >  3:251706/1‑102 (MQ=255)
                                                                                    tgAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCtt                                                   >  3:308072/1‑149 (MQ=255)
                                                                                      agaACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTcc                                                 <  3:255362/149‑1 (MQ=255)
                                                                                       gaACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCg                                                <  5:288477/149‑1 (MQ=255)
                                                                                        aaCCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGAtt                                                                  >  2:293724/1‑130 (MQ=255)
                                                                                        aaCCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGAtt                                                                  <  1:293724/130‑1 (MQ=255)
                                                                                        aaCCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGt                                               <  4:308072/149‑1 (MQ=255)
                                                                                          cccGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTg                                              >  4:306709/1‑148 (MQ=255)
                                                                                           ccGCCTTCAATAAGCCTTACCCGCCATGGTGGTCGTAAGTGGGCAGAACTCGGTATCGCATGGCTTCCTTATTCACTAGGTTTCTCAACCTGCGTACGGCTTCCGCCACAAATCCGTGAGAAAGATATTTGTCAGTCTCCTTAAGTGaa                                            >  1:15535/1‑149 (MQ=255)
                                                                                           ccGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGaa                                            >  1:158562/1‑149 (MQ=255)
                                                                                             gACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCg                                                                             <  8:111172/114‑1 (MQ=255)
                                                                                              aCTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATaa                                         <  7:225214/149‑1 (MQ=255)
                                                                                                ttGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAg                                        >  6:291933/1‑148 (MQ=255)
                                                                                                ttGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAg                                        >  5:8594/1‑148 (MQ=255)
                                                                                                 tGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGa                                      >  7:262272/1‑149 (MQ=255)
                                                                                                     tAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGGGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTg                                  <  7:227425/149‑1 (MQ=255)
                                                                                                     tAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTgg                                                             <  4:163265/122‑1 (MQ=255)
                                                                                                     tAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTgg                                                             >  3:163265/1‑122 (MQ=255)
                                                                                                     tAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTg                                  <  8:2918/149‑1 (MQ=255)
                                                                                                       aGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGc                                 <  5:291933/148‑1 (MQ=255)
                                                                                                              cccGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTAACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGaaa                          >  5:212520/1‑148 (MQ=255)
                                                                                                                         ggtCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGt                                               >  6:53895/1‑116 (MQ=255)
                                                                                                                                 gTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAAt      <  2:89345/149‑1 (MQ=255)
                                                                                                                                  ttggCAGACCTCGGTCTTGCCGGGCTTCCTTATTCACTGGGTTTCTCAACCTGCGTACGGCTGCCGCCCCAAATCCGTGATAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATc     <  8:29175/147‑1 (MQ=255)
                                                                                                                                   gggCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATca    >  1:32278/1‑149 (MQ=255)
                                                                                                                                     gCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATcaca  <  7:160070/149‑1 (MQ=255)
                                                                                                                                       aGACCTCGGTCTGGCCGGGCTTCCTTATTGACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGa            <  5:280759/137‑1 (MQ=255)
                                                                                                                                       aGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGa            >  6:280759/1‑137 (MQ=255)
                                                                                                                                            |                                                                                                                                             
ATCCTACGCGCGCCGCCGAAGTAACGCTGTTGTCGGGGAAGTGGCTCACGGATACCGTTTTCGGGTCGCTCTGGCGACTGGGTGTGAGAACCCGACTTGAATAAGCGTTACCCGCCATGGTGGTCGTACGTGGGCAGACCTCGGTCTGGCCGGGCTTCCTTATTCACTCGGTTTCTCACCCTGCGTACGGCTGCCGCCCCAAATCCGTGAGAAAGATTTTTGGCAGTCTCCTTCCGTGAATAAGGATTTGCCATGAAAAAAGATGTTTCTGAGCCAATCACA  >  NC_002947/335664‑335945

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: