Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 165,266 A→G 10.5% I64T (ATC→ACC)  pntB ← pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947165,2660AG10.5% 149.9 / 8.0 57I64T (ATC→ACC) pntBpyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (3/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CACGGCGGCGATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGCACAC  >  NC_002947/165118‑165414
                                                                                                                                                    |                                                                                                                                                    
cACGGCGGCGATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGa                                                                                                                                                      >  4:138169/1‑149 (MQ=255)
     cggcgATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCg                                                                                                                                                 >  2:111610/1‑149 (MQ=255)
         gATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGGTATCGGCta                                                                                                                                             <  5:187385/149‑2 (MQ=255)
             gCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGc                                                                                                                                         <  8:80471/149‑1 (MQ=255)
               aaTGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCaa                                                                                                                                       <  1:315834/149‑1 (MQ=255)
                     cacCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAg                                                                                                                                 >  6:4871/1‑149 (MQ=255)
                               aGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCtt                                                                                                                       <  3:138169/149‑1 (MQ=255)
                                   cGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGACCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTtata                                                                                                                    >  2:356842/1‑148 (MQ=255)
                                     aTCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGACCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAAt                                                                                                                 <  1:356842/149‑1 (MQ=255)
                                        aTGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACTAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGa                                                                                                               <  1:111610/148‑1 (MQ=255)
                                          gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGGTATCGGCTATGTCa                                                                                                                                        <  6:50401/121‑1 (MQ=255)
                                          gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGGTATCGGCTATGTCa                                                                                                                                        >  5:50401/1‑121 (MQ=255)
                                          gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGa                                                                                                               <  7:347642/146‑1 (MQ=255)
                                          gCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGa                                                                                                               >  8:347642/1‑146 (MQ=255)
                                            tgtgCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGcc                                                                                                           >  2:22002/1‑148 (MQ=255)
                                               gCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATa                                                                                                                                                    <  4:108807/104‑1 (MQ=255)
                                               gCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATa                                                                                                                                                    >  3:108807/1‑104 (MQ=255)
                                                cATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCg                                                                                                      >  2:369589/1‑149 (MQ=255)
                                                 aTGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCg                                                                                                      >  7:250583/1‑148 (MQ=255)
                                                  tGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCaa                                                                                                                           >  1:131423/1‑126 (MQ=255)
                                                  tGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCaa                                                                                                                           <  2:131423/126‑1 (MQ=255)
                                                   gAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAg                                                                                                              >  5:316380/1‑138 (MQ=255)
                                                   gAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAg                                                                                                              <  6:316380/138‑1 (MQ=255)
                                                    aaaGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAg                                                                                                   <  6:221486/148‑1 (MQ=255)
                                                      aGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGa                                                                                                <  4:272534/149‑1 (MQ=255)
                                                                  gggCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTg                                                                                                 <  4:214114/136‑1 (MQ=255)
                                                                  gggCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTg                                                                                                 >  3:214114/1‑136 (MQ=255)
                                                                  gggCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGcc                                                                                    >  2:197027/1‑149 (MQ=255)
                                                                      aTTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCAt                                                                                >  7:215393/1‑149 (MQ=255)
                                                                         ttGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGc                                                                              <  1:369589/148‑1 (MQ=255)
                                                                             tCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGAt                                                                          >  7:85071/1‑148 (MQ=255)
                                                                                tttCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATg                                                                      <  1:197027/149‑1 (MQ=255)
                                                                                   cGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGcc                                                                    <  5:4871/148‑1 (MQ=255)
                                                                                             ggCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTg                                                          <  8:250583/148‑1 (MQ=255)
                                                                                                    aTCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTgccgcgccg                                                  >  4:111058/1‑149 (MQ=255)
                                                                                                        aaCCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAg                                              >  4:217167/1‑149 (MQ=255)
                                                                                                        aaCCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAg                                              <  8:215393/149‑1 (MQ=255)
                                                                                                        aaCCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAg                                              >  5:288230/1‑149 (MQ=255)
                                                                                                         gtcggcggCACCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGAt                                                            <  5:9182/132‑1 (MQ=255)
                                                                                                         gtcggcggCACCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGAt                                                            >  6:9182/3‑134 (MQ=255)
                                                                                                             gcggTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAggtggt                                         <  3:111058/149‑1 (MQ=255)
                                                                                                             gcggTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAggtggt                                         >  8:210624/1‑149 (MQ=255)
                                                                                                                gTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAAc                                                                >  1:32402/1‑123 (MQ=255)
                                                                                                                gTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAAc                                                                <  2:32402/123‑1 (MQ=255)
                                                                                                                     gccgACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGc                                  >  8:153576/1‑148 (MQ=255)
                                                                                                                             gcagGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGACAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGccc                         >  5:218621/1‑149 (MQ=255)
                                                                                                                                 gCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGa                     <  8:33905/149‑1 (MQ=255)
                                                                                                                                        ttgccacaGCCGGTATCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTg              >  6:264987/8‑149 (MQ=255)
                                                                                                                                         tgccacaGCCGGTATCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTg              >  1:44561/7‑148 (MQ=255)
                                                                                                                                           aCATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGa            <  3:217167/148‑1 (MQ=255)
                                                                                                                                            cacaGCCGGTATCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATg          <  5:264987/146‑1 (MQ=255)
                                                                                                                                             aTAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCt               <  5:139032/143‑1 (MQ=255)
                                                                                                                                             aTAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCt               >  6:139032/1‑143 (MQ=255)
                                                                                                                                               aGCCGATACCGGCTGTGGCAAGTTCAGCGCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAg       >  1:211778/1‑149 (MQ=255)
                                                                                                                                               aGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAg       >  7:298971/1‑149 (MQ=255)
                                                                                                                                                gCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAg       >  3:17144/1‑148 (MQ=255)
                                                                                                                                                 ccGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGca     <  1:22002/149‑1 (MQ=255)
                                                                                                                                                 ccGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGca     <  4:17144/149‑1 (MQ=255)
                                                                                                                                                    aTACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGCacac  <  2:211778/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACGGCGGCGATGGCAATGAACACCGCGGCCAGGCCGATCATGCTGTGCATGAAAGCGACCAGTTCGGGCATTTTGGTCATTTCGACGCGCTTGGCCATGATCGAACCGGCGGTGCCGCCGACCAGCAGGCCGACGATGACATAGCCGATACCGGCTGTGGCAAGTTCAGCCCCAAGCTTATAAATGAGGCCAACCGTAGTGAGGATGGCGAGCCCCATGCCGATCATGCCGAACAGATTGCCGCGCCGCGAGGTGGTCGGGTGCGACAGGCCCTTGAGCGCCTGGATGAAGCACAC  >  NC_002947/165118‑165414

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