Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,528,914 C→T 9.3% A244V (GCC→GTC)  dnaEB → error‑prone DNA polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,528,9140CT 9.3% 149.1 / 4.8 55A244V (GCC→GTC) dnaEBerror‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base T (3/2);  total (29/26)
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.41e-01

GTGGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCCGCTGCGCTGGCCAGA  >  NC_002947/3528767‑3529061
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                             >  4:398681/1‑127 (MQ=255)
                         cATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGGCCAGCAGGTCgagcgga                                                                                                                                             <  3:398681/131‑5 (MQ=255)
                         cATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGGCCAGCAGGTCgagcgga                                                                                                                                             <  7:46025/131‑5 (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                                                                                              >  7:402083/1‑130 (MQ=255)
                                                                         tCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTAc                                                                                              <  8:402083/130‑1 (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)
                                                                                                                                             ctgctgGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCcgctgcgctg        >  3:82164/1‑149 (MQ=255)
                                                                                                                                                ctgGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCCGCTGCGCTGGc      <  8:217163/148‑1 (MQ=255)
                                                                                                                                                   gCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGGCCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCCGCTGCGCTGGCCAGa  >  6:173682/1‑149 (MQ=255)
                                                                                                                                                   gCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCCGCTGCGCTGGCCAGa  >  3:259336/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTGGTCGCCGTGCCTTGCAGGACTGCATGACTGCCATCCGCCAGCATTGCCAGGTGAGCGAGGCCGGGCGCTTTCTGTTTGCCAATGGCGAACGGCACCTGCGCAGCCAGGCGCAACTGGCCGAGCTGTATCCGCTCGACCTGCTGGCCGAGACCCTGGTCATTGCCCGGCGCTGCCAGTTTGACCTGAGCGAACTGAACTACCAGTACCCACGCGAACTGGTGCCCGAGGGGCATACCCCGGCAAGCTGGCTGCGCGAATTGTGCGAGCAAGGCATGCCGCTGCGCTGGCCAGA  >  NC_002947/3528767‑3529061

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: