Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,092 T→C 100% A245A (GCA→GCG U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,0920TC100.0% 241.5 / NA 66A245A (GCA→GCGU712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/37);  total (29/37)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGA  >  CP006881/3723957‑3724219
                                                                                                                                       |                                                                                                                               
aaaaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                              <  1:25913/139‑1 (MQ=255)
aaaaaaCTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                              <  6:314986/139‑1 (MQ=255)
        tCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTCGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCAAg                                                                                                                      >  7:126496/1‑139 (MQ=255)
         caagagATCACTCCCTGAAAGTCAAAGCGGTTATTCAGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAg                                                                                                                      <  8:126496/136‑1 (MQ=255)
           agagATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGAACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgcgc                                                                                                                   <  6:277314/139‑1 (MQ=255)
               aTCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACAt                                                                                                               >  3:213183/1‑139 (MQ=255)
                 cACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACAtt                                                                                                              >  8:88403/1‑138 (MQ=255)
                        tGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGc                                                                                                       >  6:231738/1‑138 (MQ=255)
                          cAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGcga                                                                                                     <  1:84445/138‑3 (MQ=255)
                              tCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGaa                                                                                                >  2:183851/1‑139 (MQ=255)
                                  aGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACt                                                                                            >  5:275839/1‑139 (MQ=255)
                                        tATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACttttttt                                                                                      >  2:305211/1‑139 (MQ=255)
                                            ccGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                                   <  4:213183/138‑1 (MQ=255)
                                            ccGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                                   <  4:26793/138‑1 (MQ=255)
                                             cGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGAt                                                                                  <  5:126011/138‑1 (MQ=255)
                                               agagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                                >  7:168191/1‑138 (MQ=255)
                                               agagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                                >  6:135130/1‑138 (MQ=255)
                                               agagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTATGGAtta                                                                                >  8:43034/1‑136 (MQ=255)
                                                 agaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAg                                                                                                                      >  2:180496/1‑98 (MQ=255)
                                                 agaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAg                                                                                                                      <  1:180496/98‑1 (MQ=255)
                                                  gaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATa                                                                            <  7:88403/139‑1 (MQ=255)
                                                  gaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATa                                                                            <  4:61630/139‑1 (MQ=255)
                                                   aTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATa                                                                            >  1:194123/1‑138 (MQ=255)
                                                      ttGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGc                                                                        <  1:120284/139‑1 (MQ=255)
                                                        gAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGc                                                                        <  7:2790/137‑1 (MQ=255)
                                                           gACGCTGGACATGCTGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACCTTGATAGGCGAATGAATACTTTTTTTGGATGAATTAAAGCAtctc                                                                   <  7:151176/139‑1 (MQ=255)
                                                           gACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtctc                                                                   <  4:130477/139‑1 (MQ=255)
                                                           gACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtctc                                                                   <  1:96389/139‑1 (MQ=255)
                                                           gACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtctc                                                                   <  8:168191/139‑1 (MQ=255)
                                                             cGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCATAAACTCCGCCAGCTCCAGCTCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcct                                                                 <  7:43034/139‑1 (MQ=255)
                                                             cGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcct                                                                 >  6:169526/1‑139 (MQ=255)
                                                              gCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTTATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcctc                                                                >  8:119642/1‑139 (MQ=255)
                                                              gCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcctc                                                                >  7:301084/1‑139 (MQ=255)
                                                              gCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcctc                                                                <  2:250559/139‑1 (MQ=255)
                                                              gCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTATTCAAAAAATCCGCCAGCTCCAGCGCACAATGATACGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcct                                                                 >  4:199682/1‑138 (MQ=255)
                                                                tGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                              >  2:255738/1‑139 (MQ=255)
                                                                tGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTGTTTTTGGATGAAATAAAGCATCTCCTCGt                                                              >  6:83300/1‑139 (MQ=255)
                                                                  gACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                              <  8:313661/137‑1 (MQ=255)
                                                                   aCATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCa                                                           <  8:265534/139‑1 (MQ=255)
                                                                   aCATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCa                                                           <  3:319426/139‑1 (MQ=255)
                                                                     aTGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATa                                                         <  1:110245/139‑1 (MQ=255)
                                                                          gCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATg                                                    >  1:312132/1‑139 (MQ=255)
                                                                            gTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGa                                                  <  5:292270/139‑1 (MQ=255)
                                                                             tGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                                 >  1:213961/1‑139 (MQ=255)
                                                                             tGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                                 <  8:301084/139‑1 (MQ=255)
                                                                             tGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                                 <  7:316416/139‑1 (MQ=255)
                                                                             tGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                                 <  8:33338/139‑1 (MQ=255)
                                                                                   aGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCtt                                            >  5:292490/1‑138 (MQ=255)
                                                                                     ccGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                                 >  6:24440/1‑131 (MQ=255)
                                                                                     ccGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                                 <  5:24440/131‑1 (MQ=255)
                                                                                     ccGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCAt                                         >  1:177356/1‑139 (MQ=255)
                                                                                       ggCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCAttt                                       >  8:153511/1‑139 (MQ=255)
                                                                                        gCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                              >  4:175654/1‑51 (MQ=255)
                                                                                        gCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                              <  3:175654/51‑1 (MQ=255)
                                                                                          cTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                                   >  6:141667/1‑92 (MQ=255)
                                                                                          cTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                                   <  5:141667/92‑1 (MQ=255)
                                                                                               attCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCaa                               <  6:275839/137‑1 (MQ=255)
                                                                                                          tGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCggctgaccg                    <  1:257/139‑3 (MQ=255)
                                                                                                            tAATTGATATTGATTTTTCAAAAACTCCTCCAGCTCCAGCGCACATTGATATGCGAATGAATACTTTTTTTGGATGATATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGAATCCAACCGGCTGACcgcg                  <  4:93681/139‑1 (MQ=255)
                                                                                                              aTTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                              >  7:80420/1‑93 (MQ=255)
                                                                                                              aTTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                              <  8:80420/93‑1 (MQ=255)
                                                                                                                  atatTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcgc            <  7:180444/139‑3 (MQ=255)
                                                                                                                  atatTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTAtcgc            <  6:155490/139‑3 (MQ=255)
                                                                                                                   tatTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCa           <  3:260374/139‑1 (MQ=255)
                                                                                                                      tCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAAt        >  6:100558/1‑139 (MQ=255)
                                                                                                                            ttCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGGCTGACCGCGTATCGCAAATGCGTGa  >  1:45930/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                               
AAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCAAATGCGTGA  >  CP006881/3723957‑3724219

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: