Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,025 G→A *118* (TAG→TAA ybaA → DUF1428 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,0250GA100.0% 220.6 / NA 68*118* (TAG→TAAybaADUF1428 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/31);  total (37/31)

TTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCTGGCGT  >  NC_000913/476893‑477160
                                                                                                                                    |                                                                                                                                       
ttCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCg                                                                                                                                   >  1:982976/1‑139 (MQ=255)
  cAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTg                                                                                                                                 <  2:73505/139‑1 (MQ=255)
   aaaaGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGa                                                                                                                                >  1:680541/1‑139 (MQ=255)
   aaaaGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGa                                                                                                                                >  2:1249120/1‑139 (MQ=255)
    aaaGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGAc                                                                                                                               <  1:1421424/139‑1 (MQ=255)
       gagGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACggg                                                                                                                            <  2:677765/139‑1 (MQ=255)
        agGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGc                                                                                                                           >  1:1593734/1‑139 (MQ=255)
          gTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCtg                                                                                                                         >  2:378228/1‑139 (MQ=255)
                       ctAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGt                                                                                                            >  1:964165/1‑139 (MQ=255)
                               aaGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGt                                                                                                            <  2:161680/131‑1 (MQ=255)
                               aaGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGt                                                                                                            >  1:161680/1‑131 (MQ=255)
                               aaGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaa                                                                                                    >  1:277475/1‑139 (MQ=255)
                                    gatgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                                      <  2:1685117/116‑1 (MQ=255)
                                    gatgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                                      >  1:1685117/1‑116 (MQ=255)
                                           gACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAAt                                                                                        >  2:126680/1‑139 (MQ=255)
                                                 cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                                  <  1:126680/139‑1 (MQ=255)
                                                  ggATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                                  <  1:74396/138‑1 (MQ=255)
                                                  ggATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                                  >  2:74396/1‑138 (MQ=255)
                                                    aTGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCgg                                                                               <  1:1249120/139‑1 (MQ=255)
                                                        aaGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCa                                                                           >  1:698867/1‑139 (MQ=255)
                                                         agagTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAg                                                                          <  2:964165/139‑1 (MQ=255)
                                                             ttCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTa                                                                      >  1:895331/1‑139 (MQ=255)
                                                                  cGAGTCCATGTCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGaaa                                                                 >  1:1599528/1‑139 (MQ=255)
                                                                  cGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGaaa                                                                 >  2:1255407/1‑139 (MQ=255)
                                                                   gAGTCCATGTCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAt                                                                <  2:1599528/139‑1 (MQ=255)
                                                                     gTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                                       >  1:1206578/1‑114 (MQ=255)
                                                                     gTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                                       <  2:1206578/114‑1 (MQ=255)
                                                                     gTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtct                                                              <  2:277475/139‑1 (MQ=255)
                                                                     gTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtct                                                              >  1:797022/1‑139 (MQ=255)
                                                                         aTGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGtctgtctcttatac                                                          <  2:680541/139‑14 (MQ=255)
                                                                                ttGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                             <  1:376653/113‑1 (MQ=255)
                                                                                ttGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                             >  2:376653/1‑113 (MQ=255)
                                                                                     ggCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                                      <  1:1469874/67‑1 (MQ=255)
                                                                                     ggCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                                      >  2:1469874/1‑67 (MQ=255)
                                                                                         aGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTa                                                                                                                                      >  1:1285467/1‑47 (MQ=255)
                                                                                         aGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTa                                                                                                                                      <  2:1285467/47‑1 (MQ=255)
                                                                                                gATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagc                                   <  2:1593734/139‑1 (MQ=255)
                                                                                                    tATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagccagc                               <  2:895331/139‑1 (MQ=255)
                                                                                                     aTGGCGGATTCTAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCGCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagccagcc                              <  1:1549710/139‑1 (MQ=255)
                                                                                                      tGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAt                                                                <  1:1157912/104‑1 (MQ=255)
                                                                                                      tGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAt                                                                >  2:1157912/1‑104 (MQ=255)
                                                                                                       ggCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCt                            >  2:174656/1‑139 (MQ=255)
                                                                                                        gCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTg                           >  1:1148325/1‑139 (MQ=255)
                                                                                                           gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                             >  1:1181600/1‑86 (MQ=255)
                                                                                                           gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                             <  2:1181600/86‑1 (MQ=255)
                                                                                                           gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGAcc                        >  2:1430295/1‑139 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaaa                                                                                                 <  1:172891/65‑1 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaaa                                                                                                 >  2:172891/1‑65 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGcc                                               >  1:783026/1‑115 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGcc                                               <  2:783026/115‑1 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagc                                   <  2:32409/127‑1 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagc                                   >  1:32409/1‑127 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCg                       <  2:1148325/139‑1 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCg                       >  2:823841/1‑139 (MQ=255)
                                                                                                            aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCg                       >  2:510502/1‑139 (MQ=255)
                                                                                                             ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                                  <  2:1534461/79‑1 (MQ=255)
                                                                                                             ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                                  >  1:1534461/1‑79 (MQ=255)
                                                                                                             ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGt                      >  1:1138912/1‑139 (MQ=255)
                                                                                                              tCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTa                     >  2:1265716/1‑139 (MQ=255)
                                                                                                                     aaTTATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                                       <  2:1089538/66‑1 (MQ=255)
                                                                                                                     aaTCATCGACGAATAAGTAGCGTGACGGGCTGAGCGGCGTGCCGTCAGCCCAAAAATTCACAAATg                                                                                       >  1:1089538/1‑66 (MQ=255)
                                                                                                                      atcatcGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtcttct                                                           >  2:468727/1‑93 (MQ=255)
                                                                                                                      atcatcGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtcttct                                                           <  1:468727/93‑1 (MQ=255)
                                                                                                                            gacgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCt                            <  2:619065/118‑1 (MQ=255)
                                                                                                                            gacgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCt                            >  1:619065/1‑118 (MQ=255)
                                                                                                                            gacgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCt       >  1:490150/1‑139 (MQ=255)
                                                                                                                              cgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCTgg     <  1:510502/139‑1 (MQ=255)
                                                                                                                                 aTAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCTGGCGt  <  1:378228/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
TTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCTGGCGT  >  NC_000913/476893‑477160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: