Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,678,128 T→C 6.8% E78E (GAA→GAG yphE ← putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,1280TC 6.8% 177.3 / 3.2 59E78E (GAA→GAGyphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base C (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCTTTATTGAGCGTGAAGT  >  NC_000913/2677984‑2678273
                                                                                                                                                |                                                                                                                                                 
cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGt                                                                                                                                               >  6:170254/1‑149 (MQ=255)
cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGt                                                                                                                                               <  8:312225/149‑1 (MQ=255)
 aTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                                                                                                              <  5:348474/149‑1 (MQ=255)
         gTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGTCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCa                                                                                                                                      <  6:43631/149‑1 (MQ=255)
               attcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgc                                                                                                                                 <  5:241298/146‑1 (MQ=255)
               aatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgcg                                                                                                                                >  5:124797/1‑149 (MQ=255)
                atcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGCCTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGt                                                                                                                               >  6:167435/1‑149 (MQ=255)
                atcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAGCTCCTGATAAACCGCGCGCCCCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCCCCTTCCAGTCGCGt                                                                                                                               >  1:241466/1‑149 (MQ=255)
                   aTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                                                                                                            <  5:87479/149‑1 (MQ=255)
                           gcggcgGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATc                                                                                                                     <  5:170254/148‑1 (MQ=255)
                               cggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                                                                                                                <  6:296138/149‑1 (MQ=255)
                                      aCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgc                                                                                                                                 >  5:194460/1‑125 (MQ=255)
                                      aCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgc                                                                                                                                 <  6:194460/125‑1 (MQ=255)
                                                                    tCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCAtt                                                                           >  8:376892/1‑149 (MQ=255)
                                                                     cAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTc                                                                          <  5:167435/149‑1 (MQ=255)
                                                                      aGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTc                                                                          <  6:341468/148‑1 (MQ=255)
                                                                       gCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCg                                                                         >  2:48070/1‑148 (MQ=255)
                                                                        cccTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTc                                                                          <  1:139657/146‑1 (MQ=255)
                                                                        cccTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTc                                                                          >  2:139657/1‑146 (MQ=255)
                                                                         ccTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAAt                                                                      <  6:124797/149‑1 (MQ=255)
                                                                         ccTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAAt                                                                      <  8:351519/149‑1 (MQ=255)
                                                                           ttCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCAtt                                                                           <  6:148214/142‑1 (MQ=255)
                                                                           ttCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCAtt                                                                           >  5:148214/1‑142 (MQ=255)
                                                                           ttCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAAt                                                                      >  1:263611/1‑147 (MQ=255)
                                                                             ccaccaGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATga                                                                    <  7:376892/147‑1 (MQ=255)
                                                                               accaGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATga                                                                    <  4:78593/145‑1 (MQ=255)
                                                                               accaGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATga                                                                    >  3:78593/1‑145 (MQ=255)
                                                                               accaGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCCTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATga                                                                    <  7:340005/145‑1 (MQ=255)
                                                                               accaGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCCTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATga                                                                    >  4:49010/1‑145 (MQ=255)
                                                                               accaGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCCTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATga                                                                    >  8:340005/1‑145 (MQ=255)
                                                                               accaGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCCTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATga                                                                    <  3:49010/145‑1 (MQ=255)
                                                                                 caGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTc                                                               >  3:193711/1‑148 (MQ=255)
                                                                                     cTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGAttt                                                          <  2:263611/149‑1 (MQ=255)
                                                                                      ttAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATgaga                                                                  >  1:55829/1‑140 (MQ=255)
                                                                                      ttAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATgaga                                                                  <  2:55829/140‑1 (MQ=255)
                                                                                      ttAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTg                                                         >  7:89204/1‑149 (MQ=255)
                                                                                      ttAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCACCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTg                                                         >  1:230197/1‑149 (MQ=255)
                                                                                           tcgtGATAAACCGCGCGCACCCCCAGGTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATgaga                                                                  <  2:344505/132‑1 (MQ=255)
                                                                                           tCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTg                                                                                       <  5:281467/114‑1 (MQ=255)
                                                                                           tCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTg                                                                                       >  6:281467/1‑114 (MQ=255)
                                                                                              tGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCa                                                                             <  1:143123/121‑1 (MQ=255)
                                                                                              tGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCa                                                                             >  2:143123/1‑121 (MQ=255)
                                                                                              tGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGcc                                                 >  8:80756/1‑149 (MQ=255)
                                                                                                    aCCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGa                                                                                                                 <  2:281853/79‑1 (MQ=255)
                                                                                                    aCCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGa                                                                                                                 >  1:281853/1‑79 (MQ=255)
                                                                                                     ccGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTa                                          >  4:327071/1‑149 (MQ=255)
                                                                                                     ccGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTa                                          >  7:22414/1‑149 (MQ=255)
                                                                                                               ccccAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAgcg                                >  4:156993/1‑149 (MQ=255)
                                                                                                                 ccAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCg                                                              <  2:163722/117‑1 (MQ=255)
                                                                                                                 ccAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCg                                                              >  1:163722/1‑117 (MQ=255)
                                                                                                                        aGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCACCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACt                        <  2:230197/148‑1 (MQ=255)
                                                                                                                             cACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCt                  <  3:327071/149‑1 (MQ=255)
                                                                                                                                  cGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCTTTAtt             <  2:241466/149‑1 (MQ=255)
                                                                                                                                      tCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCTTTATTGAGc         <  8:22414/149‑1 (MQ=255)
                                                                                                                                          cGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCTTTATTGAGCGTg      >  6:183967/1‑148 (MQ=255)
                                                                                                                                           gTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCTTTATTGAGCGTGaa    <  1:48070/149‑1 (MQ=255)
                                                                                                                                             aGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCt                                       >  7:98988/1‑112 (MQ=255)
                                                                                                                                             aGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCt                                       <  8:98988/112‑1 (MQ=255)
                                                                                                                                              gCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCTTTATTGAGCGTGAAGt  <  3:156993/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGCCGGTAAGCATTCGAATGAGAGTCGATTTGCCCGCGCCGTTTTTACCTAACAGCGCACGAACTTCGCCTTTATTGAGCGTGAAGT  >  NC_000913/2677984‑2678273

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