Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,905 A→G 9.7% D241G (GAC→GGC)  dnaEB → error‑prone DNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,528,9050AG 9.7% 134.7 / 5.4 54D241G (GAC→GGC) dnaEBerror‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (3/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ATGCACGTGCGTGGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCCGCTGCGC  >  NC_002947/3528757‑3529053
                                                                                                                                                     |                                                                                                                                                    
aTGCACGTGCGTGGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTATGCCAATGGCGAACGGCACCTGAGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGa                                                                                                                                                      >  5:1001/1‑149 (MQ=255)
 tGCACGTGCGTGGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGAc                                                                                                                                                     <  5:101416/149‑1 (MQ=255)
    aCGTGCGTGGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACctg                                                                                                                                                  <  3:50391/149‑1 (MQ=255)
        gcgtgGGCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTgg                                                                                                                                              <  8:342578/149‑1 (MQ=255)
          gtgGTCGCCGTGCCTTTCAGGAATGAATGACTGACATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATTGCGACCGGCACCTGCGCAGCCAGGAGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGcc                                                                                                                                            >  8:390441/1‑149 (MQ=255)
          gtgGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGcc                                                                                                                                            <  8:94762/149‑1 (MQ=255)
                gCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGAcc                                                                                                                                      >  1:435493/1‑149 (MQ=255)
                 ccGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGAccc                                                                                                                                     >  2:234320/1‑149 (MQ=255)
                   gTGCCTTGCAGGACTGCTTGTCTGCCATCGGGCAGCATTGCAAGTTGAGCGAGGCCGGGCGCTTTATGTTTGCCAATGGGGAACGGCAACTGCGCAGCGAGGCGCAACTGGCCGAGCTGTTTCGGCTCTTCCTGCTGGGCgattcgctg                                                                                                                                   <  6:1001/149‑8 (MQ=255)
                   gTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTg                                                                                                                                   <  5:70998/149‑1 (MQ=255)
                            aGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGcc                                                                                                                          <  6:330947/149‑1 (MQ=255)
                                   cATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGGCCAGCAGGTCgagcgga                                                                                                                                     <  3:398681/131‑5 (MQ=255)
                                   cATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGGCCAGCAGGTCgagcgga                                                                                                                                     <  7:46025/131‑5 (MQ=255)
                                   cATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGGCCAGCAGGTCgagcgga                                                                                                                                     >  4:398681/1‑127 (MQ=255)
                                   cATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGGCCAGCAGGTCgagcgga                                                                                                                                     >  8:46025/1‑127 (MQ=255)
                                   cATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCt                                                                                                                   >  8:191890/1‑149 (MQ=255)
                                    aTGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTg                                                                                                                  >  4:71042/1‑149 (MQ=255)
                                     tGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGc                                                                                                                 <  4:345771/149‑1 (MQ=255)
                                        cTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAg                                                                                                              <  7:19699/149‑1 (MQ=255)
                                         tGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTCGTCATTGCCCGGCGCTGCCAGt                                                                                                             >  4:299231/1‑149 (MQ=255)
                                           ccATCCGCCAGCATTGCCAGGTGAGCGAGTCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGttt                                                                                                           <  6:354539/149‑1 (MQ=255)
                                            cATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGttt                                                                                                           >  8:233592/1‑148 (MQ=255)
                                                 gCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGGCCCCTGCGCAGCCAGGCGCAACTGGCCGGGCTGTTTCCGCTCGCCCTCCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCTGTTTGAcgt                                                                                                      >  1:337172/1‑147 (MQ=255)
                                                 gCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGTCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTg                                                                                                     >  7:52388/1‑149 (MQ=255)
                                                 gCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGgc                                                                                                                      <  1:23373/132‑1 (MQ=255)
                                                 gCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGgc                                                                                                                      >  2:23373/1‑132 (MQ=255)
                                                 gCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTg                                                                                                     <  4:286243/149‑1 (MQ=255)
                                                    aGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGc                                                                                                  <  7:323242/149‑1 (MQ=255)
                                                       aTTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCgaga                                                                                                                                        <  5:194966/108‑1 (MQ=255)
                                                       aTTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCgaga                                                                                                                                        >  6:194966/1‑108 (MQ=255)
                                                       aTTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCgaa                                                                                               >  3:79487/1‑149 (MQ=255)
                                                                      aGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCTGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTAc                                                                                >  6:155656/1‑149 (MQ=255)
                                                                                cTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAAc                                                                      <  1:117569/149‑1 (MQ=255)
                                                                                cTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAAc                                                                      >  3:153124/1‑149 (MQ=255)
                                                                                   tCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTAc                                                                                      <  8:402083/130‑1 (MQ=255)
                                                                                   tCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTAc                                                                                      >  7:402083/1‑130 (MQ=255)
                                                                                     tGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTg                                                                 >  7:217163/1‑149 (MQ=255)
                                                                                       tttGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCTAACTGAACTACCAGTACCCACGCGAACTGGt                                                                  <  1:234320/146‑1 (MQ=255)
                                                                                              aTGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCTGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAgggg                                                        <  8:223834/149‑1 (MQ=255)
                                                                                              aTGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAgggg                                                        <  5:289038/149‑1 (MQ=255)
                                                                                               tGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAg                                                                                   <  6:162312/121‑1 (MQ=255)
                                                                                               tGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAg                                                                                   >  5:162312/1‑121 (MQ=255)
                                                                                               tGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGc                                                       <  7:191890/149‑1 (MQ=255)
                                                                                                           aCCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGc                                             >  3:226614/1‑148 (MQ=255)
                                                                                                           aCCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGc                                             >  3:220550/1‑148 (MQ=255)
                                                                                                              tGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAgc                                         >  6:405707/1‑149 (MQ=255)
                                                                                                                     ccAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTgc                                  <  4:79487/149‑1 (MQ=255)
                                                                                                                          cgcAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGaa                             <  4:153124/149‑1 (MQ=255)
                                                                                                                            cAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCTGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAAtt                           <  5:155656/149‑1 (MQ=255)
                                                                                                                                  gtcGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGa                     <  8:52388/147‑1 (MQ=255)
                                                                                                                                  gCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGc                       >  5:216705/1‑147 (MQ=255)
                                                                                                                                          gTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTGCCAGTACCGCCGCGAAGTGGTGCCCGAGGGGCATCCCCCGGCAAGCTGGCTGCGCGAATTGTGGGAGCAAGGcc             >  8:406812/1‑148 (MQ=255)
                                                                                                                                           tATCCGCTCGACCTGCTGTCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGATCGAACTGAACTACCAGTACCCACGCTAACTGGTGCCCGAGGGGCATACCCCTGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCa             <  7:390441/148‑1 (MQ=255)
                                                                                                                                           cagggtcTCGGCCAGCAGGTCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCAt            >  3:95884/7‑149 (MQ=255)
                                                                                                                                             tCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGc          >  6:20122/1‑149 (MQ=255)
                                                                                                                                                     aCCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCcgctgcgc  >  5:47142/1‑149 (MQ=255)
                                                                                                                                                     aCCTGCTGGCCGAGACCCTGGGCATGGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGGTGGCTGCGCGAATTGTGCGAGCAAGGCAGGCcgctgcgc  <  3:299231/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGCACGTGCGTGGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCCGCTGCGC  >  NC_002947/3528757‑3529053

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: