Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 287,977 A→G M397V (ATG→GTG)  yagH → CP4‑6 prophage, putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913287,9770AG100.0% 268.1 / NA 74M397V (ATG→GTG) yagHCP4‑6 prophage, putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/39);  total (35/39)

ACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCAT  >  NC_000913/287846‑288089
                                                                                                                                   |                                                                                                                
acacCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGtt                                                                                                           <  2:597344/139‑1 (MQ=255)
 cacCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTtc                                                                                                          >  1:1273009/1‑139 (MQ=255)
    ccTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCg                                                                                                       >  2:1446424/1‑139 (MQ=255)
              ctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACt                                                                                             >  2:889843/1‑139 (MQ=255)
                  ccGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCca                                                                                         >  2:1190406/1‑139 (MQ=255)
                          ccGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                 <  1:1190406/139‑1 (MQ=255)
                          ccGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                 <  1:1446424/139‑1 (MQ=255)
                                   ttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        <  1:626890/139‑1 (MQ=255)
                                    tACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                       >  2:553392/1‑139 (MQ=255)
                                    tACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                       <  2:1273009/139‑1 (MQ=255)
                                        gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGc                                                                   <  1:889843/139‑1 (MQ=255)
                                           ctATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctac                                                                >  2:859228/1‑139 (MQ=255)
                                            tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                <  1:142337/106‑1 (MQ=255)
                                            tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                >  2:142337/1‑106 (MQ=255)
                                            tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        >  2:152445/1‑130 (MQ=255)
                                            tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        <  1:152445/130‑1 (MQ=255)
                                                   aCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                     >  1:1317809/1‑126 (MQ=255)
                                                   aCGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                     <  2:1317809/126‑1 (MQ=255)
                                                      aCTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        >  2:1068974/1‑120 (MQ=255)
                                                      aCTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        <  1:1068974/120‑1 (MQ=255)
                                                          gCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        <  2:907969/116‑1 (MQ=255)
                                                          gCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        >  1:907969/1‑116 (MQ=255)
                                                           cTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACt                                                <  2:1044780/139‑1 (MQ=255)
                                                              aaTTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagca                                                                                      <  1:699512/98‑1 (MQ=255)
                                                              aaTTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCcagca                                                                                      >  2:699512/1‑98 (MQ=255)
                                                                  cGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTa                                         <  1:553392/139‑1 (MQ=255)
                                                                        tCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCtt                                   >  2:1475739/1‑139 (MQ=255)
                                                                        tCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCtt                                   >  2:383553/1‑139 (MQ=255)
                                                                        tCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCtt                                   >  1:715260/1‑139 (MQ=255)
                                                                        tCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCtt                                   >  2:796863/1‑139 (MQ=255)
                                                                          aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGt                                                                                                  <  1:10439/74‑1 (MQ=255)
                                                                          aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGt                                                                                                  >  2:10439/1‑74 (MQ=255)
                                                                            ccAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCGTtgtg                               <  2:913130/139‑1 (MQ=255)
                                                                               aTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGt                                                                                                  >  1:61538/1‑69 (MQ=255)
                                                                               aTCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGt                                                                                                  <  2:61538/69‑1 (MQ=255)
                                                                                 cGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTa                          <  1:1475739/139‑1 (MQ=255)
                                                                                   aCCGTGGCGCGCCGCTGGCAGCACTTCTCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                >  1:943109/1‑67 (MQ=255)
                                                                                   aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                <  2:943109/67‑1 (MQ=255)
                                                                                   aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGa                             <  1:21609/134‑1 (MQ=255)
                                                                                   aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGa                             >  2:21609/1‑134 (MQ=255)
                                                                                     cGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgag                      >  2:442978/1‑139 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGcgcg                                                                               >  1:276246/1‑81 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGcgcg                                                                               <  2:276246/81‑1 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                     >  1:1451891/1‑91 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                     >  2:1476924/1‑91 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                     <  2:1451891/91‑1 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCt                                                                     <  1:1476924/91‑1 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgagg                     <  1:442978/139‑1 (MQ=255)
                                                                                      gTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgagg                     <  2:1236223/139‑1 (MQ=255)
                                                                                          cgcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                <  1:1650851/60‑1 (MQ=255)
                                                                                          cgcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGc                                                                                                >  2:1650851/1‑60 (MQ=255)
                                                                                           gcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa                <  1:383553/139‑1 (MQ=255)
                                                                                           gcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa                <  1:796863/139‑1 (MQ=255)
                                                                                                  tGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                 >  2:449296/1‑67 (MQ=255)
                                                                                                  tGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                 <  1:449296/67‑1 (MQ=255)
                                                                                                    gcagcaCTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                 >  2:184453/1‑65 (MQ=255)
                                                                                                    gcagcaCTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                 <  1:184453/65‑1 (MQ=255)
                                                                                                         aCTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCt                                                                          <  2:1237786/67‑1 (MQ=255)
                                                                                                         aCTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCt                                                                          >  1:1237786/1‑67 (MQ=255)
                                                                                                         aCTTCGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCAt  <  2:715260/139‑1 (MQ=255)
                                                                                                             cGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        <  2:320481/65‑1 (MQ=255)
                                                                                                             cGCCTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        >  1:320481/1‑65 (MQ=255)
                                                                                                                cTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                       >  1:1572489/1‑63 (MQ=255)
                                                                                                                cTTCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                       <  2:1572489/63‑1 (MQ=255)
                                                                                                                  tCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        >  1:869715/1‑60 (MQ=255)
                                                                                                                  tCCGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                        <  2:869715/60‑1 (MQ=255)
                                                                                                                  tCCGGGCAGAAACGCGGGT‑GAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTa                                         >  1:1484169/1‑90 (MQ=255)
                                                                                                                  tCCGGGCAGAAACGCGGGT‑GAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTa                                         <  2:1484169/90‑1 (MQ=255)
                                                                                                                    cGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGcgcg                                                                               >  2:1447892/1‑51 (MQ=255)
                                                                                                                    cGGGCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGcgcg                                                                               <  1:1447892/51‑1 (MQ=255)
                                                                                                                     gggCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCaaaa                                                   >  2:1353693/1‑78 (MQ=255)
                                                                                                                     gggCAGAAACGCGGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCaaaa                                                   <  1:1353693/78‑1 (MQ=255)
                                                                                                                                cgGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCAt  >  1:441489/1‑116 (MQ=255)
                                                                                                                                cgGGTGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCAt  <  2:441489/116‑1 (MQ=255)
                                                                                                                                   |                                                                                                                
ACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCAT  >  NC_000913/287846‑288089

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