Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,862,518 Δ1 bp pseudogene (1086/1617 nt) glvC ← pseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,862,5160C.100.0% 245.6 / NA 60pseudogene (1088/1617 nt)glvCpseudogene, arbutin specific enzyme IIC component of PTS,enzyme, Transport of small molecules: Carbohydrates, organic acids, alcohols, PTS system, arbutin‑like IIB component, PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (28/32);  total (28/32)

AGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCA  >  NC_000913/3862378‑3862648
                                                                                                                                          |                                                                                                                                    
aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGc                                                                                                                                      >  2:1306931/1‑139 (MQ=255)
aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGc                                                                                                                                      >  2:727759/1‑139 (MQ=255)
aGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGc                                                                                                                                      >  1:1275999/1‑139 (MQ=255)
 gAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                     <  1:727759/139‑1 (MQ=255)
    tCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                  <  2:1575684/139‑4 (MQ=255)
        tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                                                                             <  1:613084‑M1/139‑3 (MQ=255)
        tgtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTca                                                                                                                             <  2:324560‑M1/139‑3 (MQ=255)
         gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                                            <  1:1008946‑M1/139‑4 (MQ=255)
         gtgCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                                            <  1:1526587‑M1/139‑4 (MQ=255)
                aaaCCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                                     >  1:1295062‑M1/1‑129 (MQ=255)
                   ccACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                                  <  2:265430‑M1/139‑14 (MQ=255)
                   ccACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaag                                                                                                                  <  1:1306931‑M1/139‑14 (MQ=255)
                      cGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                                               <  1:1112665‑M1/139‑17 (MQ=255)
                      cGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCACGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagaga                                                                                                               >  1:745337‑M1/1‑123 (MQ=255)
                                   ggTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                                                                  >  1:330600‑M1/1‑110 (MQ=255)
                                   ggTAAAGCACAACCCAATGGCGATATGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcg                                                                                                  >  1:1092386‑M1/1‑110 (MQ=255)
                                     tAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgac                                                                                                >  2:1013280‑M1/1‑108 (MQ=255)
                                                  aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                          >  2:262973‑M1/1‑95 (MQ=255)
                                                  aaTGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgag                                                                                          <  1:262973‑M1/132‑38 (MQ=255)
                                                       cGATTTGGGTCAGCATCATATTCGCGTGGTTGCTGAACATGGGATCCAGTTTTGCGGAAAAACTTAGTCAATCAGACCTCCGCCA‑‑‑TGCTCACCACACCAAAGCGATACATTAAGGTAGACATTGAGGCCGCCAGCACgg                                                                          >  1:260847/1‑139 (MQ=255)
                                                       cGATTTGGGTCAGCATCATATCCACGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                              <  2:745337‑M1/139‑50 (MQ=255)
                                                         aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggc                                                                            <  2:1195713‑M1/139‑52 (MQ=255)
                                                          tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                  <  2:1133892/85‑4 (MQ=255)
                                                          tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGccatg                                                                                                                                  >  1:1133892/1‑82 (MQ=255)
                                                          tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcg                                                                           >  1:1441083‑M1/1‑87 (MQ=255)
                                                          tttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcg                                                                           >  1:1126604‑M1/1‑87 (MQ=255)
                                                           ttGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                                          <  2:1275999‑M1/139‑54 (MQ=255)
                                                                tCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtacc                                                                     <  2:1295062‑M1/139‑59 (MQ=255)
                                                                    catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                     >  1:240603/1‑72 (MQ=255)
                                                                    catcatATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                     <  2:240603/72‑1 (MQ=255)
                                                                        atatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacag                                                             <  1:1013280‑M1/139‑67 (MQ=255)
                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                     >  1:703821/1‑67 (MQ=255)
                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                     <  2:703821/67‑1 (MQ=255)
                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcactacaccaaa                                                                                                                   >  2:1422425‑M1/1‑72 (MQ=255)
                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcactacaccaaa                                                                                                                   <  1:1422425‑M1/84‑13 (MQ=255)
                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                              <  2:644337‑M1/121‑50 (MQ=255)
                                                                         tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacg                                                                              >  1:644337‑M1/1‑72 (MQ=255)
                                                                          atCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagca                                                           <  2:1441083‑M1/139‑69 (MQ=255)
                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagca                                                                                >  2:206080‑M1/1‑67 (MQ=255)
                                                                              gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagca                                                                                <  1:206080‑M1/114‑48 (MQ=255)
                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                     >  2:1571902/1‑57 (MQ=255)
                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGcc                                                                                                                                     <  1:1571902/57‑1 (MQ=255)
                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgatattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaa                                                  >  1:465584‑M1/1‑62 (MQ=255)
                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaa                                                  >  1:1197039‑M1/1‑62 (MQ=255)
                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                                            <  2:733451‑M1/65‑4 (MQ=255)
                                                                                   gTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcac                                                                                                                            >  1:733451‑M1/1‑62 (MQ=255)
                                                                                     tGCTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatg                                                <  1:890199‑M1/139‑80 (MQ=255)
                                                                                       cTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaa                                              <  2:1476234‑M1/139‑82 (MQ=255)
                                                                                       cTGAACATCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacatggaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgat                                              <  1:1006127‑M1/139‑82 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcac                                                                                                                            <  2:1030776/54‑4 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGttcac                                                                                                                            >  1:1030776/1‑51 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaag                                       >  2:85088‑M1/1‑51 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaag                                       >  1:118481‑M1/1‑51 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                                          <  2:157104‑M1/104‑54 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgt                                                                          >  1:157104‑M1/1‑51 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCaaa                                                                                                                >  2:256295/1‑63 (MQ=255)
                                                                                              tCGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTT‑‑‑CACCACACCaaa                                                                                                                <  1:256295/63‑1 (MQ=255)
                                                                                               cGGGATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaagg                                      >  2:1697‑M1/1‑50 (MQ=255)
                                                                                                 ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagg                                          >  1:357883‑M1/1‑48 (MQ=255)
                                                                                                 ggATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacagg                                          <  2:357883‑M1/133‑86 (MQ=255)
                                                                                                  gATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                                                        >  1:1192617‑M1/1‑47 (MQ=255)
                                                                                                  gATCCAGTTTTGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacg                                                        <  2:1192617‑M1/118‑72 (MQ=255)
                                                                                                         ttttGCGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccag                            <  2:118481‑M1/139‑100 (MQ=255)
                                                                                                              cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggtt                       >  1:862544‑M1/1‑35 (MQ=255)
                                                                                                              cGGTAAAACTTAGTCAATCAGACCTCCGCCATGTTcaccacacca                                                                                                                     <  1:167410/45‑5 (MQ=255)
                                                                                                                   aaaCTTAGTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggta                  >  2:1109495‑M1/1‑30 (MQ=255)
                                                                                                                      cttagTCAATCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaatt               <  1:938323‑M1/137‑113 (MQ=255)
                                                                                                                               tCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                 >  2:1348887‑M1/1‑18 (MQ=255)
                                                                                                                               tCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                 <  1:1348887‑M1/64‑47 (MQ=255)
                                                                                                                               tCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                 >  2:577991‑M1/1‑18 (MQ=255)
                                                                                                                               tCAGACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagc                                                                                 <  1:577991‑M1/64‑47 (MQ=255)
                                                                                                                                  gACCTCCGCCATGTTcaccacaccaaagagatacattacggtcgatattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcccaccagcatg   <  2:465584‑M1/139‑125 (MQ=255)
                                                                                                                                   aCCTCCGCCATGTTcaccacaccaaagagatacattacggtcgacattgaggccgccagcacggcgtgtaccgcaaacagcaacggtgaaatgaacaggaaggtaaattccagcggttcggtaattcccaccagcatgg  >  2:1596142‑M1/1‑14 (MQ=255)
                                                                                                                                          |                                                                                                                                    
AGAATCAGTGTGCGGAAAACCACGAAGTACAGCAGGGTAAAGCACAACCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCA  >  NC_000913/3862378‑3862648

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