Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,954,599 C→T 6.3% A280V (GCC→GTC)  lysA_2 → Diaminopimelate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,954,5990CT 6.3% 205.0 / 3.2 63A280V (GCC→GTC) lysA_2Diaminopimelate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (17/42);  new base T (2/2);  total (19/44)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

TGGGTGGCGGTGTCGGCGTGCGCTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGG  >  ACB122/5954467‑5954733
                                                                                                                                     |                                                                                                                                      
tggGTGGCGGTGTCGGCGTGCGCTACCGCGCCGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCGGGTGTTCGAACCAGGCCGCTACATCGTCGCCac                                                                                                                                     >  7:383382/1‑135 (MQ=255)
   gggggggTGGCGGCGGGCGCTACCGGGACGATCCGCCACCACGGGTGGCCGACTACATCCAGGCCA‑TCCGCGAGCGCGGGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCTCCAACg                                                                                                                                   <  6:86934/130‑1 (MQ=255)
     ggCGGTGTCGGCGTGCGCGACCGCGCCGAACCGCCACCACGGGTGGCCGCCGACATCAAGGCCC‑TCCCCGCGCGCGTGGGCGAGCGCGCCCTGGCCCAGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGcc                                                                                                                                 >  6:163566/1‑135 (MQ=255)
      gCGGGGTCGGCGTGCGCTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCg                                                                                                                                <  2:150109/135‑1 (MQ=255)
      gCGGGGGCGGCGTGCGCTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCg                                                                                                                                <  3:69099/135‑1 (MQ=255)
                  tGCGCTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGa                                                                                                                    <  4:363504/135‑1 (MQ=255)
                    cgcTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGAccc                                                                                                                 <  7:38869/136‑1 (MQ=255)
                     gcTCCCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGAccc                                                                                                                 <  5:87300/135‑1 (MQ=255)
                     gcTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGAccc                                                                                                                 <  5:152688/135‑1 (MQ=255)
                     gcTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGAccc                                                                                                                 <  1:44150/135‑1 (MQ=255)
                       gaCCGCGCCGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCAGGCCCTGGGGTTGGAACCCGGCCGCGACATCGTCGCCACCGCCGGCGTGCTGCTGCCCcgcg                                                                                                              >  8:361699/2‑136 (MQ=255)
                         ccgctccgaACAGCCACCACGGGGGGCGGACTCCATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGt                                                                                                             <  5:171632/129‑1 (MQ=255)
                         ccGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGt                                                                                                             <  8:243488/135‑1 (MQ=255)
                                 aaCAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCt                                                                                                    >  7:219876/1‑136 (MQ=255)
                                  aCAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCt                                                                                                    <  1:35122/135‑1 (MQ=255)
                                    aGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCa                                                                                                  <  1:256977/135‑1 (MQ=255)
                                    aGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTCGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGACTACCTc                                                                                                   >  7:286630/1‑134 (MQ=255)
                                       caccacGGGGGGCCGATTACATCAAGGCCA‑GCCGCGCGCGCGGGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCCTCTCCAACGCCGGCGTGCTGCTGACCCGCGTGGAAAACCTCAAGc                                                                                               <  3:56350/135‑1 (MQ=255)
                                         ccacGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGcaca                                                                                            <  8:105105/136‑1 (MQ=255)
                                          cacGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGcaca                                                                                            <  3:301382/135‑1 (MQ=255)
                                          cacGGGGGCCGGACTACATCAAGCCCC‑TCCGGTCGCGCGTGGGCTTCCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCACCGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGcaca                                                                                            <  1:377223/135‑1 (MQ=255)
                                               gTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAac                                                                                      <  7:396509/136‑1 (MQ=255)
                                                tGGCCGACACCATCAAGGCAA‑TCCGCGAGGGGGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAac                                                                                      <  5:285373/135‑1 (MQ=255)
                                                 ggCCGACTTCATCAAGGCCACTCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAac                                                                                      <  1:162957/135‑1 (MQ=255)
                                                 ggCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTg                                                                                                            <  8:175085/112‑1 (MQ=255)
                                                 ggCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTg                                                                                                            >  7:175085/1‑112 (MQ=255)
                                                     gACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGGTCGAACCAGGCCGCGACATCGTCGCCAACGCCGGCGTGCTGCGGACCCGCGTGGAATACCTCAAGCCCACCGAACACAAg                                                                                 >  5:304751/1‑135 (MQ=255)
                                                          cATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACtt                                                                            >  4:231716/1‑135 (MQ=255)
                                                             cacggCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGc                                                                         <  5:309202/132‑1 (MQ=255)
                                                                ggCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCcat                                                                      >  8:41563/1‑135 (MQ=255)
                                                                  ccA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACAGAACACAAGGACTTCGCcatct                                                                    >  1:162154/1‑134 (MQ=255)
                                                                      tCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCg                                                                 <  3:293239/135‑1 (MQ=255)
                                                                      tCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCg                                                                 <  6:348678/135‑1 (MQ=255)
                                                                      tCCGCGAGCGCGGGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCg                                                                 <  5:33269/135‑1 (MQ=255)
                                                                         gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGTCAACGCCGGCGTGCTGCTGACCcgcg                                                                                                              <  6:358358/87‑1 (MQ=255)
                                                                         gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGTCAACGCCGGCGTGCTGCTGACCcgcg                                                                                                              <  6:364240/87‑1 (MQ=255)
                                                                         gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGTCAACGCCGGCGTGCTGCTGACCcgcg                                                                                                              >  5:364240/1‑87 (MQ=255)
                                                                         gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGTCAACGCCGGCGTGCTGCTGACCcgcg                                                                                                              >  5:358358/1‑87 (MQ=255)
                                                                         gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGAt                                                               >  5:182347/1‑134 (MQ=255)
                                                                         gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGAt                                                               >  3:391588/1‑134 (MQ=255)
                                                                         gcgAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCCCGTGGTCTACCTCAAGCACCCCGAACACAAGGACTTCGCCATCATCGATGc                                                             >  7:337164/1‑136 (MQ=255)
                                                                               gcgTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGa                                                        <  4:11269/135‑1 (MQ=255)
                                                                                cgTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGAt                                                       <  4:166093/135‑1 (MQ=255)
                                                                                        aGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGGGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCt                                              <  7:325762/136‑1 (MQ=255)
                                                                                        aGCGCGACCTGGCCCTGGTGGTCGAACCAGGCCGCTACATCGTCGCCACCGCCGGCGTGCTGCTGACCCGCGTGGACTCCCTCCCGCACACCGAACACAAGGCCTTCGCCCTCAGCGATGCCGCGATGAACGAcc                                               >  4:120956/1‑135 (MQ=255)
                                                                                              aCCTGGCCCTGGTGTTCGACCCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATcc                                         <  1:226563/135‑1 (MQ=255)
                                                                                              aCCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCACCGCCGGCGTGCTGCTGACCCGCGTGTAATACCTCACGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATcc                                         <  1:338251/135‑1 (MQ=255)
                                                                                                   gCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCgg                                   <  8:219876/136‑1 (MQ=255)
                                                                                                    cccTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCgg                                   <  4:391588/135‑1 (MQ=255)
                                                                                                         gtgtTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTg                              <  7:361699/135‑1 (MQ=255)
                                                                                                         gtgtTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTg                              <  7:87538/135‑1 (MQ=255)
                                                                                                          tgtTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGt                             <  5:163566/135‑1 (MQ=255)
                                                                                                               gAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTTGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCa                         <  4:153773/134‑1 (MQ=255)
                                                                                                               gAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCTATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCa                         <  3:314074/134‑1 (MQ=255)
                                                                                                                 aCCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAac                                                                                      >  7:317234/1‑71 (MQ=255)
                                                                                                                 aCCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAac                                                                                      <  8:317234/71‑1 (MQ=255)
                                                                                                                 aCCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGc                      >  2:205543/1‑135 (MQ=255)
                                                                                                                   cAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTg                   <  8:383382/136‑1 (MQ=255)
                                                                                                                             cATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGgtg           >  2:313773/1‑134 (MQ=255)
                                                                                                                              aTCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGgtgt          <  3:188286/134‑1 (MQ=255)
                                                                                                                              aTCGTCGCCAACGCCGGCGTGCGGCGGACCCGCGTGGAATACCTCAAGCACACCGCACACACGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGgtgt          <  1:124496/134‑1 (MQ=255)
                                                                                                                                   cgCCAACGCCGGCGGGCTGCGGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAgcgc    <  5:4703/135‑1 (MQ=255)
                                                                                                                                    gCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGcgcg   <  3:183179/135‑1 (MQ=255)
                                                                                                                                     ccAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCgg  >  3:200519/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TGGGTGGCGGTGTCGGCGTGCGCTACCGCGACGAACAGCCACCACGGGTGGCCGACTACATCAAGGCCA‑TCCGCGAGCGCGTGGGCGAGCGCGACCTGGCCCTGGTGTTCGAACCAGGCCGCTACATCGTCGCCAACGCCGGCGTGCTGCTGACCCGCGTGGAATACCTCAAGCACACCGAACACAAGGACTTCGCCATCATCGATGCAGCGATGAACGACCTGATCCGCCCGGCCCTGTACCAGGCCTGGATGGGTGTCAGCGCGG  >  ACB122/5954467‑5954733

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