Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,977,675 T→A 7.1% Q227L (CAG→CTG)  PP_2604 ← major facilitator family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,977,6750TA 7.1% 169.9 / 2.3 56Q227L (CAG→CTG) PP_2604major facilitator family transporter
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base A (2/2);  total (26/30)
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.76e-01

GCTGGTTCAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAACGCAGCCTGCCAGCCCCATACAGCCAGGATCAG  >  NC_002947/2977529‑2977817
                                                                                                                                                   |                                                                                                                                              
gCTGGTTCAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACAGGCGCGGAAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGCATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCctgctg                                                                                                                                               >  3:127747/1‑148 (MQ=255)
 cTGGTTCAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGa                                                                                                                                             >  8:253598/1‑149 (MQ=255)
    gTTCAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCAGCAGGATga                                                                                                                                          >  1:67774/1‑148 (MQ=255)
       cAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACGGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGGTGAACGGCACGGGGGTCGGCCCGGGCTCGTCATCCTGCTGGATGTAc                                                                                                                                        <  1:5071/148‑1 (MQ=255)
       cAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGCAAGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTa                                                                                                                                         <  1:102261/148‑1 (MQ=255)
       cAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACt                                                                                                                                       >  4:23609/1‑149 (MQ=255)
        aGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTg                                                                                                                                      >  3:106955/1‑149 (MQ=255)
        aGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCAGCAGGATgacga                                                                                                                                       <  2:67774/148‑5 (MQ=255)
         ggAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCAGCAGGATgacga                                                                                                                                       >  4:284897/1‑143 (MQ=255)
         ggAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCAGCAGGATgacga                                                                                                                                       <  3:284897/147‑5 (MQ=255)
          gAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCa                                                                                                                                    >  5:221059/1‑149 (MQ=255)
          gAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCa                                                                                                                                    >  3:219482/1‑149 (MQ=255)
                       cAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGacg                                                                                                                       >  1:87101/1‑149 (MQ=255)
                        aGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGacgc                                                                                                                      >  7:189272/1‑149 (MQ=255)
                              ggTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCgg                                                                                                                <  3:23609/149‑1 (MQ=255)
                               gTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCggg                                                                                                               >  3:306904/1‑149 (MQ=255)
                               gTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCggg                                                                                                               >  4:319743/1‑149 (MQ=255)
                                    aCCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGggtg                                                                                                           <  4:219482/148‑1 (MQ=255)
                                        gAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGAACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCtt                                                                                                      >  4:66984/1‑149 (MQ=255)
                                             ccGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCggg                                                                                                 <  2:109587/149‑1 (MQ=255)
                                                gcgcGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTg                                                                                              <  7:200382/149‑1 (MQ=255)
                                                 cgcgGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTgt                                                                                             <  4:324953/149‑1 (MQ=255)
                                                  gcgGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTgta                                                                                            <  4:106955/149‑1 (MQ=255)
                                                     gTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTgtagta                                                                                         <  8:71485/149‑1 (MQ=255)
                                                           gCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCa                                                                                    <  2:45834/148‑1 (MQ=255)
                                                             aTAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCa                                                                                 >  5:146799/1‑149 (MQ=255)
                                                                 gCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCAccc                                                                              <  6:221059/148‑1 (MQ=255)
                                                                  ccACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCa                                                                             >  7:281810/1‑148 (MQ=255)
                                                                        cAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAaca                                                                        <  8:46612/147‑1 (MQ=255)
                                                                         aaaaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCaaaa                                                                          <  5:53293/144‑1 (MQ=255)
                                                                         aaaaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCaaaa                                                                          >  6:53293/1‑144 (MQ=255)
                                                                         aaaaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGa                                                                     <  3:326866/149‑1 (MQ=255)
                                                                         aaaaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGa                                                                     <  7:253598/149‑1 (MQ=255)
                                                                         aaaaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGa                                                                     >  8:244684/1‑149 (MQ=255)
                                                                         aaaaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGa                                                                     >  6:109026/1‑149 (MQ=255)
                                                                           aaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGAcc                                                                   <  5:109026/149‑1 (MQ=255)
                                                                           aaGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGAcc                                                                   <  6:146799/149‑1 (MQ=255)
                                                                              cccAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCcac                                                                <  2:87101/149‑1 (MQ=255)
                                                                               ccAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCcaca                                                               >  5:183874/1‑149 (MQ=255)
                                                                                cAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCacac                                                              <  8:189272/149‑1 (MQ=255)
                                                                                cAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCacac                                                              >  1:288743/1‑149 (MQ=255)
                                                                                   gTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCCCCCAAACCACGAGCCACAccag                                                           >  8:138442/1‑149 (MQ=255)
                                                                                      ccACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTcc                                                                                                    >  4:324009/1‑105 (MQ=255)
                                                                                      ccACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTcc                                                                                                    <  3:324009/105‑1 (MQ=255)
                                                                                         cgc‑aggcgcaggATGCGGCTGTACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACAccaggcccag                                                     <  7:157049/149‑1 (MQ=255)
                                                                                             aggcgcaggATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCt                                                  <  3:66984/149‑1 (MQ=255)
                                                                                               gcgcaggATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGc                                                >  4:235147/1‑149 (MQ=255)
                                                                                                         cGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGa                                       <  4:306904/148‑1 (MQ=255)
                                                                                                         cGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGAt                                      <  7:244684/149‑1 (MQ=255)
                                                                                                          ggCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATg                                     <  8:43442/149‑1 (MQ=255)
                                                                                                                  ggCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAAGGCAga                             >  1:123338/1‑148 (MQ=255)
                                                                                                                                ctcgggctcgTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAACGCAGCCTGCCAGCCCCATa               <  8:281810/149‑1 (MQ=255)
                                                                                                                                      ctcgTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGCACGCAGCCTGCCAGCCCCATACAGccc         >  4:187613/1‑148 (MQ=255)
                                                                                                                                      ctcgTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAACGCAGCCTGCCAGCCCCATACAGCCa         <  3:96878/149‑1 (MQ=255)
                                                                                                                                         gTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAACGCAGCCTGCCAGCCCCATACAGCCAGGa      >  7:47425/1‑149 (MQ=255)
                                                                                                                                             tCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAACGCAGCCTGCCAGCCCCATACAGCCAGGATCAg  <  3:235147/149‑1 (MQ=255)
                                                                                                                                             tCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTCGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAACGCAGCCTGCCAGCCCCATACAGCCAGGATCAg  >  4:190302/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GCTGGTTCAGGAACGGCGGCAGCCAGTACAGGTAGAACCAGAACACCGGCGCGGTAATCGCATAGGCCACTGCAAAAGCCCAGGTGCCACGC‑AGGCGCAGGATGCGGCTGAACGGCACGCGGGTCGGCTCGGGCTCGTCATCCTGCTGGATGTACTGCACCTCGCTCTGACGCACGCGGGGGTGGTCTTCCGGGTTGTAGTATTTCAGCACCCAAAACACGACCCACACCAGGCCCAGGCTGCCCATGGCGATGAACGCAGCCTGCCAGCCCCATACAGCCAGGATCAG  >  NC_002947/2977529‑2977817

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