Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,794,290 T→G 7.4% C45G (TGC→GGC) ‡ PP_2447 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,794,2900TG 7.4% 163.3 / 4.6 54C45G (TGC→GGC) ‡PP_2447hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base G (2/2);  total (28/26)
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.99e-01

CGCCGCTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGCCTAT  >  NC_002947/2794142‑2794429
                                                                                                                                                     |                                                                                                                                           
cgccgcTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCt                                                                                                                                             <  2:18219/149‑1 (MQ=255)
  ccgcTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTgc                                                                                                                                           <  6:301683/149‑1 (MQ=255)
    gcTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTgcgc                                                                                                                                         >  2:297227/1‑149 (MQ=255)
    gcTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTgcgc                                                                                                                                         <  1:59399/149‑1 (MQ=255)
     cTGTCGGACCTATGATTGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGcc                                                                                                                                        >  2:44111/1‑149 (MQ=255)
        tCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGgc                                                                                                                                     >  3:102129/1‑149 (MQ=255)
        tCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGgc                                                                                                                                     >  5:360883/1‑149 (MQ=255)
                aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGcc                                                                                                                             >  2:253680/1‑149 (MQ=255)
                aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGcc                                                                                                                             >  3:117890/1‑149 (MQ=255)
                aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgccaggg                                                                                                                                <  3:259481/146‑6 (MQ=255)
                aTGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgccaggg                                                                                                                                >  4:259481/1‑141 (MQ=255)
                          gCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgcca                                                                                                                                   <  7:143258/133‑3 (MQ=255)
                          gCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCGCAGCCGgcca                                                                                                                                   >  8:143258/1‑131 (MQ=255)
                            cAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTc                                                                                                                   >  3:80270/1‑147 (MQ=255)
                             agagCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAg                                                                                                                >  8:92809/1‑149 (MQ=255)
                              gagCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAgg                                                                                                               <  2:145400/149‑1 (MQ=255)
                               agCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAgg                                                                                                               <  1:256777/149‑1 (MQ=255)
                                gCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGc                                                                                                             >  8:177631/1‑149 (MQ=255)
                                  cGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCt                                                                                                           <  5:370311/149‑1 (MQ=255)
                                    ggCGCTGCCTGCCTTGCGTTGGATCAATCGTGTGTCCATGGATCAAAATTCCTTCAAGGGTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGc                                                                                                         <  7:285307/149‑1 (MQ=255)
                                    ggCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGc                                                                                                         >  8:28943/1‑149 (MQ=255)
                                        ctgcctgcCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCa                                                                                                     <  5:173273/149‑1 (MQ=255)
                                           cctgcctgGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGc                                                                                                   <  1:253680/148‑1 (MQ=255)
                                                    cGTTTGTACATGGGATTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCGACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGGAACCTCAAGGGCCTGCCGCAGCGAGTCGAGAt                                                                                         <  5:109485/149‑1 (MQ=255)
                                                      ttGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATca                                                                                       <  1:297227/149‑1 (MQ=255)
                                                          aaCAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATcagcag                                                                                   >  4:18057/1‑149 (MQ=255)
                                                           aCAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCagcagc                                                                                  <  4:80270/149‑1 (MQ=255)
                                                             aaGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTTCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGt                                                                                >  3:232183/1‑149 (MQ=255)
                                                                         aTGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCg                                                                    <  4:102129/149‑1 (MQ=255)
                                                                           ggAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTg                                                                  >  7:282046/1‑149 (MQ=255)
                                                                            gAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCTGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTAGCGTg                                                                  <  6:360883/148‑1 (MQ=255)
                                                                                       cTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAAcc                                                      >  1:52358/1‑149 (MQ=255)
                                                                                       cTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAAc                                                       <  1:173610/148‑1 (MQ=255)
                                                                                        ttCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAAc                                                             <  7:93932/141‑1 (MQ=255)
                                                                                        ttCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAAc                                                             >  8:93932/1‑141 (MQ=255)
                                                                                        ttCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACca                                                     >  1:333324/1‑149 (MQ=255)
                                                                                            aGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGc                                                 >  2:311981/1‑149 (MQ=255)
                                                                                             ggCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCg                                                >  2:217646/1‑149 (MQ=255)
                                                                                                tGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCgggg                                             >  5:70158/1‑149 (MQ=255)
                                                                                                 gCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGc                                            <  3:18057/149‑1 (MQ=255)
                                                                                                    ccccGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCAt                                         <  4:117890/149‑1 (MQ=255)
                                                                                                      ccGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATgg                                       >  8:86084/1‑149 (MQ=255)
                                                                                                      ccGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATgg                                       >  2:107257/1‑149 (MQ=255)
                                                                                                       cggcggCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGgc                                      <  7:92809/149‑1 (MQ=255)
                                                                                                                cTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTggcgg                              <  1:311981/148‑1 (MQ=255)
                                                                                                                    ctgcctgcctgcTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTg                         <  7:177631/149‑1 (MQ=255)
                                                                                                                        ctgcctgcTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAg                           <  3:45193/143‑1 (MQ=255)
                                                                                                                        ctgcctgcTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAg                           >  4:45193/1‑143 (MQ=255)
                                                                                                                                  gctgCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAggg           >  2:300939/1‑149 (MQ=255)
                                                                                                                                     gCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAgg            >  3:262207/1‑145 (MQ=255)
                                                                                                                                     gCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAgg            <  4:262207/145‑1 (MQ=255)
                                                                                                                                       ccACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTTCCGCAGCGCGGCGAGCTCAGCAGCGTGCCGGTCTACCGGGGCAACGCCCACCCCCGCGGGGCCAAGGCCCTGGCGGCCGTGCGTTCGCAGCAGGGGgcg       >  8:145904/1‑149 (MQ=255)
                                                                                                                                        cACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGgcgc      <  3:243138/149‑1 (MQ=255)
                                                                                                                                         aCCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGcc     >  5:218182/1‑149 (MQ=255)
                                                                                                                                            cTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGCCTAt  <  4:232183/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
CGCCGCTGTCGGACCTATGATCGGTGGCCAGAGCCGGGCGCTGCCTGCCTGGCGTTGGAACAAGCGAGTGTCCATGGAGCAACATTCCTTCAAGGCTGCCCCCCGGCGGCTACTGGCTGCCTGCCTGCTGGCTGCCACCCTGGCCGGCTGCGCCGGCGCACCGCCCGGCAACCTCAAGGGCCTGCCGCAGCGAGTCGAGATCAGCAGCGTGCCGTTCTACCGTGGCAACGCCAACCACAGCGGGGCCATGGCGCTGGCGGCAGTGCTGTCGCAGCAGGGGGCGCCTAT  >  NC_002947/2794142‑2794429

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: