Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 253,369 G→A 14.5% intergenic (‑154/+17) PP_0202 ← / ← PP_0203 CBS domain‑containing protein/dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947253,3690GA14.5% 131.3 / 13.2 55intergenic (‑154/+17)PP_0202/PP_0203CBS domain‑containing protein/dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base A (4/4);  total (26/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 = 4.65e-01

TGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATGGC  >  NC_002947/253225‑253511
                                                                                                                                                |                                                                                                                                              
tGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGTCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                            <  8:175090/149‑1 (MQ=255)
  cGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAgg                                                                                                                                             <  5:242243/146‑1 (MQ=255)
  cGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTCTTCGCGGTTTTAGCCGCGAAAAGGcc                                                                                                                                           <  6:292571/148‑1 (MQ=255)
   gggTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCgg                                                                                                                                         >  6:185588/1‑149 (MQ=255)
                aCTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaa                                                                                                                            <  2:26775/149‑1 (MQ=255)
                 cTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaaa                                                                                                                           <  5:185588/149‑1 (MQ=255)
                        ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCg                                                                                                                     <  4:274444/148‑1 (MQ=255)
                               cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAaga                                                                                                                                               <  6:308571/115‑1 (MQ=255)
                               cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAaga                                                                                                                                               >  5:308571/1‑115 (MQ=255)
                               cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGc                                                                                                             <  1:148905/149‑1 (MQ=255)
                               cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGc                                                                                                             >  5:207966/1‑149 (MQ=255)
                                tAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTACCGGCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGCTAc                                                                                                                                      <  6:282157/123‑1 (MQ=37)
                                tAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTACCGGCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGCTAc                                                                                                                                      >  5:282157/1‑123 (MQ=37)
                                 aaaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGccc                                                                                                           >  8:50255/1‑149 (MQ=255)
                                  aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAagag                                                                                                                                              >  7:288281/1‑113 (MQ=255)
                                  aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAagag                                                                                                                                              <  8:288281/113‑1 (MQ=255)
                                  aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAg                                                                                                                          >  4:218053/1‑133 (MQ=255)
                                  aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAg                                                                                                                          <  3:218053/133‑1 (MQ=255)
                                  aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGccc                                                                                                           <  4:39744/148‑1 (MQ=255)
                                  aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGcc                                                                                                            <  1:243566/147‑1 (MQ=255)
                                   aGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCg                                                                                                          >  3:261391/1‑148 (MQ=255)
                                          ggTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCAt                                                                                                                               <  5:60504/120‑1 (MQ=255)
                                          ggTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCAt                                                                                                                               >  6:60504/1‑120 (MQ=255)
                                                 cgcTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCAtt                                                                                           <  7:50255/149‑1 (MQ=255)
                                                       gCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCcgc                                                                                                    <  1:27846/134‑1 (MQ=255)
                                                       gCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCcgc                                                                                                    >  2:27846/1‑134 (MQ=255)
                                                           ggCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcacca                                                                                 >  4:116152/1‑149 (MQ=255)
                                                              aaTGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCt                                                                              >  4:37708/1‑149 (MQ=255)
                                                               aTGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTg                                                                             >  6:194541/1‑149 (MQ=255)
                                                                tGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGc                                                                            <  6:7729/149‑1 (MQ=255)
                                                                  cAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCaa                                                                          <  4:261391/149‑1 (MQ=255)
                                                                           gCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCAt                                                                 <  6:207966/149‑1 (MQ=255)
                                                                                   gCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCt                                                         <  3:116152/149‑1 (MQ=255)
                                                                                          aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATGCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACgg                                                   >  4:43039/1‑148 (MQ=255)
                                                                                          aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTc                                                        <  4:46821/143‑1 (MQ=255)
                                                                                          aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACAGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTc                                                        >  3:46821/1‑143 (MQ=255)
                                                                                             aaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                               <  4:9190/149‑1 (MQ=255)
                                                                                             aaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                               <  3:37708/149‑1 (MQ=255)
                                                                                                      gCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgtt                                      >  1:242262/1‑149 (MQ=255)
                                                                                                       cTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttg                                     <  5:194541/149‑1 (MQ=255)
                                                                                                               ttGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                               <  8:103095/131‑1 (MQ=255)
                                                                                                               ttGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                               >  7:103095/1‑131 (MQ=255)
                                                                                                                     cGGCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                               >  5:144901/1‑125 (MQ=38)
                                                                                                                     cGGCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                               <  6:144901/125‑1 (MQ=38)
                                                                                                                      ggCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGc                            >  6:147294/1‑143 (MQ=39)
                                                                                                                      ggCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGc                            <  5:147294/143‑1 (MQ=39)
                                                                                                                       gCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgtg                      >  4:45895/1‑148 (MQ=255)
                                                                                                                          ttttCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttgtcgtt                                <  7:240050/135‑1 (MQ=255)
                                                                                                                          ttttCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttgtcgtt                                >  8:240050/1‑135 (MQ=255)
                                                                                                                                    tAAAACCGCTAAGCGGCCGGTACTGCGCATATAAGCTTCGGCTTCTGCCCGATCCGCGCCACCATTTCTCTCACCACCTGCAAGTCCATCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTTCCCGgtgt                     <  8:189544/136‑1 (MQ=255)
                                                                                                                                    tAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCgg        >  2:202173/1‑149 (MQ=255)
                                                                                                                                    tAAAACCGCGAAGAGGCCGGTACTGCGAATCAAAGCTTCGGCATCTGCCCGATCCGAGCCAACATTTCGATCACCACCTGCAAGTCCAGAATGAACTGCTCCACGTTCTTGAACTCGTTGTCGTTGTGCCCGgtgt                     >  7:189544/1‑136 (MQ=255)
                                                                                                                                     aaaaCCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCgg        >  2:254295/1‑148 (MQ=255)
                                                                                                                                        aCCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGCCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATg    <  2:287918/149‑1 (MQ=255)
                                                                                                                                        aCCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGCCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATg    >  1:287918/1‑149 (MQ=255)
                                                                                                                                         ccGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATgg   <  2:242262/149‑1 (MQ=255)
                                                                                                                                          cgcgAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCTCTCACCACCTGCAAGTCCAGCATGAACTGCTCAACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATggc  >  2:152349/1‑149 (MQ=255)
                                                                                                                                          cgcgAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACgg                                                   >  7:257928/1‑100 (MQ=255)
                                                                                                                                          cgcgAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTACCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATggc  >  8:34097/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATGGC  >  NC_002947/253225‑253511

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