Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,882,428 C→G 7.0% Y75* (TAC→TAG PP_2537 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,882,4280CG 7.0% 180.6 / 4.3 57Y75* (TAC→TAGPP_2537hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base G (2/2);  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 = 5.95e-01

GATGCCGAACTGGTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGA  >  NC_002947/2882282‑2882575
                                                                                                                                                  |                                                                                                                                                      
gATGCCGAACTGGTGCAAATGAAACCGCAACTGCCAATACGCACTGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGtacta                                                                                                                                                      <  5:44113/149‑1 (MQ=255)
     cGAACTGGTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCAGGGGCTGCAAAACCGAGATAAACCAAAAGGCCATCGACACGCCTGACGACCGCCAGGTGATGAAGAAGGGCCCCTACTACCAgg                                                                                                                                                 >  1:314776/1‑149 (MQ=255)
         cTGGGGCAAAGGAACCCGCAACTGCCAATACGCACCGCACTGCTGGAAGCGTCGGGCATCATGGGCTGCAGCACCGAGATCACCCGCCCGGCCAGCGACAGGCCGGACGACCGCCAGGGGATGATGAAGGCCACGTACTACCTGGCCgg                                                                                                                                             <  3:31607/149‑1 (MQ=255)
            gTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCAGGTAGTACGTGGccttc                                                                                                                                            <  3:15928/147‑4 (MQ=255)
            gTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGccatccc                                                                                                                                          <  8:335130/149‑6 (MQ=255)
             tGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCTTCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATg                                                                                                                                         <  1:317721/149‑1 (MQ=255)
             tGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGccatcccg                                                                                                                                         >  3:109065/1‑143 (MQ=255)
             tGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATg                                                                                                                                         <  7:224170/149‑1 (MQ=255)
             agCAAATGAAACCGCAACAGCCAATACGCACCGCACTTATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATg                                                                                                                                         <  7:336969/148‑1 (MQ=255)
                aaaTGAAACCGCAACTGCCAATACGCACCGCAGTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCCTCGACAGGCCTGCCGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGc                                                                                                                                       >  2:123400/1‑148 (MQ=255)
                   tGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAg                                                                                                                                   >  1:337409/1‑149 (MQ=255)
                   tGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAg                                                                                                                                   >  4:189213/1‑149 (MQ=255)
                    gAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGt                                                                                                                                  >  5:234384/1‑149 (MQ=255)
                       aCCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCAc                                                                                                                               <  4:262294/149‑1 (MQ=255)
                         cGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTg                                                                                                                             <  6:234384/149‑1 (MQ=255)
                                 ccAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCg                                                                                                                                              <  3:156975/124‑1 (MQ=255)
                                 ccAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCg                                                                                                                                              >  4:156975/1‑124 (MQ=255)
                                   aaTACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGcc                                                                                                                   >  2:39435/1‑149 (MQ=255)
                                   aaTACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCGGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGcc                                                                                                                   <  2:49097/149‑1 (MQ=255)
                                          accgcacTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAgctgct                                                                                                            >  3:299875/1‑149 (MQ=255)
                                          accgcacTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAgctgct                                                                                                            >  1:353816/1‑149 (MQ=255)
                                               acTGATGGATGCTTCGGGCATCAGGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGa                                                                                                       <  2:314776/149‑1 (MQ=255)
                                               acTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGa                                                                                                       <  4:192786/149‑1 (MQ=255)
                                                   aTGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACg                                                                                                  >  3:333350/1‑149 (MQ=255)
                                                     ggAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGt                                                                                                 >  2:29574/1‑148 (MQ=255)
                                                      gAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTc                                                                                                <  6:152312/148‑1 (MQ=255)
                                                        aGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAgg                                                                                             <  2:353816/149‑1 (MQ=255)
                                                         gCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGAGGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGgc                                                                                            <  6:390879/149‑1 (MQ=255)
                                                           ttCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGgcgc                                                                                          >  5:142467/1‑149 (MQ=255)
                                                                    catGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCg                                                                                >  2:254710/1‑149 (MQ=255)
                                                                             cATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGaa                                                                        <  5:155351/148‑1 (MQ=255)
                                                                                  ccGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAg                                                                  >  8:57438/1‑149 (MQ=255)
                                                                                         cACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACg                                                                                                                      <  2:341083/92‑1 (MQ=255)
                                                                                         cACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACg                                                                                                                      >  1:341083/1‑92 (MQ=255)
                                                                                            ccGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGc                                                        >  6:36811/1‑149 (MQ=255)
                                                                                            ccGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGc                                                        <  3:189213/149‑1 (MQ=255)
                                                                                              gCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATATAAGCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGcc                                                       >  1:324956/1‑149 (MQ=255)
                                                                                                  ggCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGCAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTggg                                                   <  1:254710/149‑1 (MQ=255)
                                                                                                       tCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGgcgc                                            <  2:337409/149‑1 (MQ=255)
                                                                                                           cAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAAc                                        >  7:127762/1‑149 (MQ=255)
                                                                                                             ggCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACtt                                      <  8:127762/149‑1 (MQ=255)
                                                                                                              gCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTg                                     >  7:233448/1‑149 (MQ=255)
                                                                                                                            aGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCgg                       >  2:276821/1‑149 (MQ=255)
                                                                                                                              gTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACTGCCAGCTGCTGCAGATTCGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGtt                     <  1:123400/149‑1 (MQ=255)
                                                                                                                                gatgatgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGttt                    >  7:356681/1‑148 (MQ=255)
                                                                                                                                  tgatgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCt                              <  7:364222/136‑1 (MQ=255)
                                                                                                                                  tgatgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCt                              >  8:364222/1‑136 (MQ=255)
                                                                                                                                    atgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGTTACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCCGCGAACTTGGCACCCTTGGGCGGTTTTaa                 <  1:29574/148‑1 (MQ=255)
                                                                                                                                    atgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATc               <  4:333350/149‑1 (MQ=255)
                                                                                                                                     tgaAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGGGAAGCGTAGCACGGGTTGCAGGGGGGCCGCAACATTGGAACCGTTGGGCGGTTTAAAtgt              >  3:220704/1‑147 (MQ=255)
                                                                                                                                       ccggCCAGGTAGTACGTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTc            >  7:39963/3‑149 (MQ=255)
                                                                                                                                       aaGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGCTACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGACCTTGGCACCCTTGGGCGGTTTTAATCGTc            >  8:295136/1‑149 (MQ=255)
                                                                                                                                        cggCCAGGTAGTACGTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGt                      <  3:94535/137‑1 (MQ=255)
                                                                                                                                        cggCCAGGTAGTACGTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGt                      >  4:94535/2‑138 (MQ=255)
                                                                                                                                        aGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCt           <  5:301136/149‑1 (MQ=255)
                                                                                                                                         ggCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTg          >  4:237090/1‑149 (MQ=255)
                                                                                                                                           ccACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGt         >  4:15487/1‑148 (MQ=255)
                                                                                                                                              cGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAaca     <  5:36811/149‑1 (MQ=255)
                                                                                                                                                  ctacCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGCCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGa  <  1:276821/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GATGCCGAACTGGTGCAAATGAAACCGCAACTGCCAATACGCACCGCACTGATGGAAGCTTCGGGCATCATGGGCTGCATCACCGAGATCACCCGCAAGGCCATCGACAGGCCTGACGACCGCCAGGTGATGATGAAGGCCACGTACTACCTGGCCGGGATGGCCAAGTCACTGATCGACGGCCAGCTGCTGCAGATACGTCAGGCGCGCGAAGGCGGGGGCGAAGCCTAGCCAGGCTTGCCTGGGGGGCGCGAACTTGGCACCCTTGGGCGGTTTTAATCGTCTGTAACACGA  >  NC_002947/2882282‑2882575

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: