Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 62,075 A→G 6.8% Y107C (TAT→TGT)  PP_0053 → sulfide:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294762,0750AG 6.8% 167.5 / 2.8 59Y107C (TAT→TGT) PP_0053sulfide:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base G (2/2);  total (30/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 = 7.42e-01

GGTGCCGGTGTGTTCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTGGTGCAAAAGCTCA  >  NC_002947/61930‑62218
                                                                                                                                                 |                                                                                                                                                
ggTGCCGGTGTGTTCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTAt                                                                                                                                                 <  5:207079/147‑1 (MQ=255)
  tGCCGGTGTGTTCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGaa                                                                                                                                              >  1:172909/1‑148 (MQ=255)
  tGCCGGTGTGTTCAAGGCACCAAGCACTGCCCGCACCAAGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAAc                                                                                                                                             <  8:192184/149‑1 (MQ=255)
          tgtTCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATc                                                                                                                                    <  2:207134/149‑1 (MQ=255)
            ttCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACTCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGt                                                                                                                                  <  1:292998/149‑1 (MQ=255)
            ttCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTGTTCATAGCTGAtggc                                                                                                                                  <  7:259414/149‑4 (MQ=255)
            ttCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTGTTCATAGCTGAtg                                                                                                                                    <  5:199370/147‑2 (MQ=255)
            ttCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGTACAGCGGAACGt                                                                                                                                  >  3:211958/1‑149 (MQ=255)
            ttCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCg                                                                                                                                   <  1:353601/148‑1 (MQ=255)
            ttCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCg                                                                                                                                   >  4:27157/1‑148 (MQ=255)
            ttCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCg                                                                                                                                   >  5:18602/1‑148 (MQ=255)
             tCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCAGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTc                                                                                                                                 >  3:49905/1‑149 (MQ=255)
               aaGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTg                                                                                                                               >  7:94664/1‑149 (MQ=255)
                    aCCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCAATGGCCCAGTTGGTCATCCTCCATGATGGTCGCGCCATCAGCAATGAACAGCTGATCGTCTGACCCg                                                                                                                          >  6:85645/1‑149 (MQ=255)
                     cgaaGCACTGCCCGAACCTCGGCGTGCACCATTCCTCGGGGTGGGCGCTGGGGCAAGTCTCGGGTAAAGGGTTTCGACTCATTGGCCCAGTTGTTCATCCTCGATGATGGTCGCGCCATCAGCTATTAACAGCTGATCGTCTGCCCCgg                                                                                                                         <  8:76396/147‑1 (MQ=255)
                               cccGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTg                                                                                                               <  3:207736/149‑1 (MQ=255)
                                 cGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGa                                                                                                             >  6:251463/1‑149 (MQ=255)
                                 cGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGa                                                                                                             <  8:254671/149‑1 (MQ=255)
                                  gCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGAc                                                                                                            <  3:27157/149‑1 (MQ=255)
                                      cATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTggg                                                                                                        >  8:216514/1‑149 (MQ=255)
                                       aTGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGggc                                                                                                       <  7:90244/149‑1 (MQ=255)
                                           cGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGcc                                                                                                   <  4:269620/149‑1 (MQ=255)
                                                aCCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGa                                                                                              <  5:85645/149‑1 (MQ=255)
                                                 ccATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGAt                                                                                             <  6:253819/149‑1 (MQ=255)
                                                    ttCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATgg                                                                                           >  6:290640/1‑148 (MQ=255)
                                                    ttCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATgg                                                                                           <  5:290640/148‑1 (MQ=255)
                                                       cTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACACGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCt                                                                                        <  1:92090/149‑1 (MQ=255)
                                                                  cgcTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCt                                                                            <  8:236551/149‑1 (MQ=255)
                                                                           aaGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCaa                                                                   <  3:166955/149‑1 (MQ=255)
                                                                           aaGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGAAGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCACGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACACTGGGCCGCaa                                                                   >  1:358652/1‑149 (MQ=255)
                                                                                gcgGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGgtg                                                              >  8:299725/1‑149 (MQ=255)
                                                                                    gTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCAc                                                          >  7:63326/1‑149 (MQ=255)
                                                                                          ggATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCaa                                                    >  8:274108/1‑149 (MQ=255)
                                                                                                  cccATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGcc                                                                       >  4:253909/1‑122 (MQ=255)
                                                                                                  cccATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGcc                                                                       <  3:253909/122‑1 (MQ=255)
                                                                                                   ccATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGggc                                                                                                       >  1:283039/1‑89 (MQ=255)
                                                                                                   ccATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGggc                                                                                                       <  2:283039/89‑1 (MQ=255)
                                                                                                   ccATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATc                                               <  5:207427/145‑1 (MQ=255)
                                                                                                   ccATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATc                                               >  6:207427/1‑145 (MQ=255)
                                                                                                     aTTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTAc                                          >  5:173146/1‑148 (MQ=255)
                                                                                                         gCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCa                                                           >  1:147836/1‑127 (MQ=255)
                                                                                                                      tCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATAcc                        >  4:143190/1‑149 (MQ=255)
                                                                                                                        cTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTg                      <  5:111660/149‑1 (MQ=255)
                                                                                                                         tCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTgg                     >  6:319550/1‑149 (MQ=255)
                                                                                                                         tCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTgg                     >  1:225014/1‑149 (MQ=255)
                                                                                                                         tCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTgg                     >  3:136982/1‑149 (MQ=255)
                                                                                                                         tCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTgg                     >  2:375708/1‑149 (MQ=255)
                                                                                                                           gatgatGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTgg                     >  6:144053/1‑147 (MQ=255)
                                                                                                                           gatgatGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTgg                     <  7:216514/147‑1 (MQ=255)
                                                                                                                             tgatgGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGc                 >  2:252606/1‑149 (MQ=255)
                                                                                                                               atgGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTg               >  8:338902/1‑149 (MQ=255)
                                                                                                                                tgGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTg               <  5:251463/148‑1 (MQ=255)
                                                                                                                                tgGGCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCGGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTg               <  8:94664/148‑1 (MQ=255)
                                                                                                                                       gcCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTGGTGCaaaa       <  3:143190/149‑1 (MQ=255)
                                                                                                                                         gaTCAGCTGTTCATAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTGGTGCAAAAGc     >  6:340630/2‑149 (MQ=255)
                                                                                                                                         gaTCAGCTGTTCATAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATAGCTACGACCTGGCACCCTATACCTGGGAGCTGGTGCAAAAGc     >  3:105356/2‑149 (MQ=255)
                                                                                                                                          aTCAGCTATGAACAGCTGATCGTCTGCCCCGCCCTCAAGCTGGACTGggcgg                                                                                                     <  2:348913/52‑1 (MQ=255)
                                                                                                                                          aTCAGCTATGAACAGCTGATCGTCTGCCCCGCCCTCAAGCTGGACTGggcgg                                                                                                     >  1:348913/1‑52 (MQ=255)
                                                                                                                                             aGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTGGTGCAAAAGCTca  <  1:252606/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGTGCCGGTGTGTTCAAGGCACCAAGCACTGCCCGCACCATGGCGTCGACCATTCCTCGCGGTGTGCGCTGGGTCAAGGCGCGGGTAAAGGGATTCGACCCATTGGCCCAGTTGGTCATCCTCGATGATGGTCGCGCCATCAGCTATGAACAGCTGATCGTCTGCCCCGGCCTCAAGCTGGACTGGGCGGCCATTGATGGGCTCAGCGAGACCCTGGGCCGCAACGGTGTCACCTCCAACTATCGCTACGACCTGGCACCCTATACCTGGGAGCTGGTGCAAAAGCTCA  >  NC_002947/61930‑62218

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: