Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,826,413 A→G K8E (AAA→GAA)  srlE → glucitol/sorbitol‑specific enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,826,4130AG97.2% 241.9 / ‑2.8 71K8E (AAA→GAA) srlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (34/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.56e-01

TCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGA  >  NC_000913/2826286‑2826546
                                                                                                                               |                                                                                                                                     
tcACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGt                                                                                                                            >  2:1673602/1‑139 (MQ=255)
  aCCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTgg                                                                                                                          >  1:1570843/1‑139 (MQ=255)
        ttcttcCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTgggg                                                                                                                    <  2:1649952/139‑1 (MQ=255)
           ttcCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGggcgg                                                                                                                 <  1:435508/139‑1 (MQ=255)
             ccGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGcc                                                                                                               <  2:1028895/139‑1 (MQ=255)
                       gTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCt                                                                                                     >  1:1503565/1‑139 (MQ=255)
                           ccGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAg                                                                                                       >  1:836380/1‑133 (MQ=255)
                           ccGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAg                                                                                                       <  2:836380/133‑1 (MQ=255)
                                  gACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGc                                                                                          >  1:1020558/1‑139 (MQ=255)
                                        cGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaa                                                                                    >  2:197189/1‑139 (MQ=255)
                                        cGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaa                                                                                    >  2:1066078/1‑139 (MQ=255)
                                             tttttGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaaaaaaa                                                                               <  2:1734310/139‑1 (MQ=255)
                                              ttttGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATACTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaaaaaaaa                                                                              <  2:1278997/139‑1 (MQ=255)
                                                 tGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCg                                                                           <  2:722190/139‑1 (MQ=255)
                                                  gagaAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATTACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGt                                                                          <  1:422151/139‑1 (MQ=255)
                                                  gagaAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGt                                                                          >  1:1303874/1‑139 (MQ=255)
                                                           aTGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCAc                                                                 <  1:197189/139‑1 (MQ=255)
                                                             gggCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCg                                                                                        <  2:1888665/114‑1 (MQ=255)
                                                             gggCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCg                                                                                        >  1:1888665/1‑114 (MQ=255)
                                                               gCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCAccgcc                                                             >  1:1290744/1‑139 (MQ=255)
                                                                  ttCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGt                                                          >  2:1626927/1‑139 (MQ=255)
                                                                    cAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTAc                                                        >  1:239704/1‑139 (MQ=255)
                                                                          gAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGcc                                                  >  1:730085/1‑139 (MQ=255)
                                                                          gAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCg                                                     >  2:1562756/1‑136 (MQ=255)
                                                                          gAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCg                                                     <  1:1562756/136‑1 (MQ=255)
                                                                           aaCAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCt                                                 >  2:1628661/1‑139 (MQ=255)
                                                                               aaaaGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGa                                             <  2:1570843/139‑1 (MQ=255)
                                                                               aaaaGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGa                                             <  2:1290744/139‑1 (MQ=255)
                                                                               aaaaGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGa                                             <  1:1626927/139‑1 (MQ=255)
                                                                                     tCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGaa                                                                                                                                 <  1:1317017/49‑1 (MQ=255)
                                                                                     tCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGaa                                                                                                                                 >  2:1317017/1‑49 (MQ=255)
                                                                                     tCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTg                                                <  2:108094/130‑1 (MQ=255)
                                                                                     tCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTg                                                >  1:108094/1‑130 (MQ=255)
                                                                                      cACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTatat                                                                    >  1:1663685/1‑109 (MQ=255)
                                                                                      cACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTatat                                                                    <  2:1663685/109‑1 (MQ=255)
                                                                                      cACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGa                                      <  1:1628661/139‑1 (MQ=255)
                                                                                       aCCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGAc                                     >  2:141511/1‑139 (MQ=255)
                                                                                            gCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAgg                                                                                                                                   <  2:1703333/40‑1 (MQ=39)
                                                                                            gCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAgg                                                                                                                                   >  1:1703333/1‑40 (MQ=255)
                                                                                            gCAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACt                                >  1:1031234/1‑139 (MQ=255)
                                                                                             cAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCtt                                                                                                          >  2:1710743/1‑64 (MQ=255)
                                                                                             cAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCtt                                                                                                          <  1:1710743/64‑1 (MQ=255)
                                                                                             cAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGcc                                                                                         >  2:2051554/1‑81 (MQ=255)
                                                                                             cAGGAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGcc                                                                                         <  1:2051554/81‑1 (MQ=255)
                                                                                               ggAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGc                             <  2:730085/139‑1 (MQ=255)
                                                                                                gAGCAACATCATGACGCATATTCGGATCGAAGAAGGa                                                                                                                                  >  1:33059/1‑37 (MQ=255)
                                                                                                gAGCAACATCATGACGCATATTCGGATCGAAGAAGGa                                                                                                                                  <  2:33059/37‑1 (MQ=255)
                                                                                                gAGCAACATCATGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGca                            <  1:1066078/139‑1 (MQ=255)
                                                                                                      catcatGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCa                                                                  <  2:1149907/95‑1 (MQ=255)
                                                                                                      catcatGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCa                                                                  >  1:1149907/1‑95 (MQ=255)
                                                                                                          atGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTActg                  >  1:1411684/1‑139 (MQ=255)
                                                                                                           tGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGAt                                            <  2:1387274/112‑1 (MQ=255)
                                                                                                           tGACGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGAt                                            >  1:1387274/1‑112 (MQ=255)
                                                                                                              cGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCAccgcc                                                             >  2:964468/1‑92 (MQ=255)
                                                                                                              cGCATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCAccgcc                                                             <  1:964468/92‑1 (MQ=255)
                                                                                                                cATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTg                                                >  2:1418328/1‑103 (MQ=255)
                                                                                                                cATATTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTg                                                <  1:1418328/103‑1 (MQ=255)
                                                                                                                 atatTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaaaaaaaa                                                                              <  2:1957817/72‑1 (MQ=255)
                                                                                                                 atatTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaaaaaaaa                                                                              >  1:1957817/1‑72 (MQ=255)
                                                                                                                 atatTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTc                                                                         >  1:468448/1‑77 (MQ=255)
                                                                                                                 atatTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTc                                                                         <  2:468448/77‑1 (MQ=255)
                                                                                                                   atTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGcac                           <  1:713943/121‑1 (MQ=255)
                                                                                                                   atTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGcac                           >  2:713943/1‑121 (MQ=255)
                                                                                                                   atTCGGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAgg         <  2:1031234/139‑1 (MQ=255)
                                                                                                                   atTCAGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaaaaaaaa                                                                              >  1:427233/1‑70 (MQ=255)
                                                                                                                   atTCAGATCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGaaaaaaaa                                                                              <  2:427233/70‑1 (MQ=255)
                                                                                                                         aTCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGca                            >  2:270164/1‑114 (MQ=255)
                                                                                                                         aTCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGca                            <  1:270164/114‑1 (MQ=255)
                                                                                                                          tCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAgg         >  2:299118/1‑132 (MQ=255)
                                                                                                                          tCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAgg         <  1:299118/132‑1 (MQ=255)
                                                                                                                          tCGAAGAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGa  <  1:141511/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                     
TCACCAATTTCTTCCGCGGCTGGGTGACCGATCTGACCACCGCCATTTTTGAGAAAAAGATGGGCATTCAACTTGAACAAAAAGTTCACCTGGCAGGAGCAACATCATGACGCATATTCGGATCGAAAAAGGAACGGGTGGCTGGGGCGGCCCGCTTGAGCTGAAAGCCACGCCGGGAAAAAAAATCGTCTATATCACCGCCGGTACCCGGCCTGCGATTGTTGACAAACTGGCACAGCTTACTGGCTGGCAGGCTATTGA  >  NC_000913/2826286‑2826546

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