Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,430,392 A→G 8.6% intergenic (‑74/+87) PP_4769 ← / ← PP_4770 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,430,3920AG 8.6% 154.2 / 4.6 58intergenic (‑74/+87)PP_4769/PP_4770hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base G (2/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

AAATGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAA  >  NC_002947/5430245‑5430530
                                                                                                                                                   |                                                                                                                                           
aaaTGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAg                                                                                                                                            >  1:223652/1‑149 (MQ=255)
  aTGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCGGGAGTGCCGGGGGCATCGTCCCCACACCTACGCGGCCCCCCCCGCTCCCCCAgg                                                                                                                                           >  6:205017/1‑148 (MQ=255)
   tGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGGCACGCCCGCTCCCACGGGtg                                                                                                                                         >  1:62225/1‑148 (MQ=255)
       aGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATcgc                                                                                                                                    >  3:113716/1‑149 (MQ=255)
            tGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTa                                                                                                                               >  4:49775/1‑149 (MQ=255)
              ggCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTac                                                                                                                              <  3:49775/148‑1 (MQ=255)
                    gcgcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGt                                                                                                                        <  8:14314/148‑1 (MQ=255)
                    gcgcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                       <  6:7128/149‑1 (MQ=255)
                        aCGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGGCACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAg                                                                                                                   <  2:62225/149‑1 (MQ=37)
                          ggTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCg                                                                                                                 <  2:86456/149‑1 (MQ=255)
                           gtgtAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCgg                                                                                                                <  5:227077/149‑1 (MQ=255)
                           gtgtACTAGAGGTTCTGCACATCGGCGAACAGCGCGATCCTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCgg                                                                                                                <  2:60718/149‑1 (MQ=255)
                            tgtAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCggg                                                                                                               <  2:39360/149‑1 (MQ=37)
                            tgtAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCggg                                                                                                               >  4:318949/1‑149 (MQ=255)
                               agtagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGt                                                                                                            >  6:213153/1‑149 (MQ=255)
                                           gCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAg                                                                                                                                            <  1:159037/106‑1 (MQ=255)
                                           gCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAg                                                                                                                                            >  2:159037/1‑106 (MQ=255)
                                           gCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                >  1:13522/1‑149 (MQ=255)
                                            cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                       <  4:200687/125‑1 (MQ=255)
                                            cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                       >  3:200687/1‑125 (MQ=255)
                                             acaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgc                                                                                              <  2:13522/149‑1 (MQ=255)
                                             acaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgc                                                                                              <  3:318949/149‑1 (MQ=255)
                                              caTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCGCCCCCAGGTATCGGGCGTACCCGCGGGGGAGCGGGCGTGCCCGCGAAGGGgcg                                                                                             >  8:285430/1‑149 (MQ=255)
                                               aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAgggg                                                                                               >  8:316791/1‑146 (MQ=255)
                                               aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAgggg                                                                                               <  7:316791/146‑1 (MQ=255)
                                               aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgcgc                                                                                            <  5:213153/149‑1 (MQ=255)
                                               aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgcg                                                                                             >  7:46738/1‑148 (MQ=255)
                                                 cGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCa                                                                                           >  3:292466/1‑148 (MQ=255)
                                                      aaCAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCGGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAgcgc                                                                                     <  7:285430/149‑1 (MQ=255)
                                                             cgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGa                                                                              <  4:292466/149‑1 (MQ=255)
                                                                tCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                           >  3:143100/1‑149 (MQ=255)
                                                                 cTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATg                                                                          >  1:340394/1‑149 (MQ=255)
                                                                  tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGc                                                                         >  2:241729/1‑149 (MQ=255)
                                                                    tttCCCCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                           <  6:38826/145‑1 (MQ=255)
                                                                    tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAGGGGGCGCTTTGCGccccttcgcgg                                                                           <  3:324461/145‑8 (MQ=38)
                                                                    tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                           <  2:223652/145‑1 (MQ=255)
                                                                     ttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCca                                                                       <  3:261733/148‑1 (MQ=255)
                                                                      tCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                >  8:170164/1‑122 (MQ=255)
                                                                      tCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                <  7:170164/122‑1 (MQ=255)
                                                                        aCCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCGCCCACACGTAGCGCGCGTACACCCGTGGGAGCGGGCGGGCCCGCGAACGGGCGCAAAGCGCCCCCGGAAATGCcaccag                                                                   >  8:245740/1‑149 (MQ=255)
                                                                                  cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                         >  1:50372/1‑149 (MQ=255)
                                                                                      tttttGGCCTGGAGTGCCTGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACATGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGt                                                     <  1:247086/149‑1 (MQ=255)
                                                                                        tttGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTc                                                    <  4:113716/148‑1 (MQ=255)
                                                                                                  aGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCcag                                                                      <  4:259947/120‑1 (MQ=255)
                                                                                                  aGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCcag                                                                      >  3:259947/1‑120 (MQ=255)
                                                                                                   gTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGgc                                                                                              <  3:87067/95‑1 (MQ=255)
                                                                                                   gTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGGGGGAGCGGGCGTGCCCGCGAAGGGgc                                                                                              >  4:87067/1‑95 (MQ=255)
                                                                                                         gggggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATcgtcgt                                  <  8:46738/149‑1 (MQ=255)
                                                                                                              ctgctgCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCgag                             <  2:50372/149‑1 (MQ=255)
                                                                                                                    cgcAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAAc                       <  4:143100/149‑1 (MQ=255)
                                                                                                                         aCCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCAcc                  <  1:241729/149‑1 (MQ=255)
                                                                                                                           cTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGt                >  1:357525/1‑149 (MQ=255)
                                                                                                                               gcgGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATc            >  1:6369/1‑149 (MQ=255)
                                                                                                                                       gCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGTGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCa     >  3:107198/1‑148 (MQ=255)
                                                                                                                                       gCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCa     >  4:212487/1‑148 (MQ=255)
                                                                                                                                         ccGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaaaa  <  2:6369/149‑1 (MQ=255)
                                                                                                                                           gCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAgaagaa                        >  7:282035/1‑125 (MQ=25)
                                                                                                                                           gCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAgaagaa                        <  8:282035/125‑1 (MQ=25)
                                                                                                                                                   |                                                                                                                                           
AAATGGCAGCCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAA  >  NC_002947/5430245‑5430530

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: