Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 184,257 C→T *271* (TAG→TAA yaeI ← phosphodiesterase with model substrate bis‑pNPP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913184,2570CT100.0% 236.2 / NA 72*271* (TAG→TAAyaeIphosphodiesterase with model substrate bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (39/33);  total (39/33)

TGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTA  >  NC_000913/184124‑184383
                                                                                                                                     |                                                                                                                              
tgtgCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTacac                                                                                                                           <  1:76132/139‑1 (MQ=255)
tgtgCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTacac                                                                                                                           <  2:1662555/139‑1 (MQ=255)
     ttATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGtt                                                                                                                      >  2:1550476/1‑139 (MQ=255)
     ttATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTAATTCTGGTGATAAAATTCACGATTAACACCAGtt                                                                                                                      >  1:1110650/1‑139 (MQ=255)
                   ggcctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                                 <  1:616133/130‑1 (MQ=255)
                   ggcctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                                 >  2:616133/1‑130 (MQ=255)
                      ctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTcc                                                                                                     >  1:1277208/1‑139 (MQ=255)
                              cccgGCTCGACTGTGATCGCAGAGAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAGATCAAAAAATTACCTGCTTTCTTCTGGTGATAAAATTCACGATTTacac                                                                                                                           <  2:1122014/109‑1 (MQ=255)
                              cccgGCTCGACTGTGATCGCAGAGAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTacac                                                                                                                           >  1:1122014/1‑109 (MQ=255)
                              cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                                 <  2:260429/119‑1 (MQ=255)
                              cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGc                                                                                                                 >  1:260429/1‑119 (MQ=255)
                              cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCg                                                                                                    <  1:1123398/132‑1 (MQ=255)
                              cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCg                                                                                                    >  2:1123398/1‑132 (MQ=255)
                               ccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACAcc                                                                                                                          <  2:115307/109‑1 (MQ=255)
                               ccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACAcc                                                                                                                          >  1:115307/1‑109 (MQ=255)
                                   cTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCa                                                                                        >  2:494606/1‑139 (MQ=255)
                                        cTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGc                                                                                   <  2:1110650/139‑1 (MQ=255)
                                         tgtgATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCa                                                                                  <  2:1277208/139‑1 (MQ=255)
                                          gtgATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                 <  1:494606/139‑1 (MQ=255)
                                                gCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATa                                                                           >  2:573045/1‑139 (MQ=255)
                                                     aaGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAg                                                                                            >  1:1324602/1‑117 (MQ=255)
                                                     aaGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAg                                                                                            <  2:1324602/117‑1 (MQ=255)
                                                     aaGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAAc                                                                      >  1:1213642/1‑139 (MQ=255)
                                                     aaGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAAc                                                                      >  2:675579/1‑139 (MQ=255)
                                                         cTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGcc                                                                  <  1:573045/139‑1 (MQ=255)
                                                            tttCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccac                                                               >  2:409413/1‑139 (MQ=255)
                                                               cttcCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgcc                                                            <  1:1550476/139‑1 (MQ=255)
                                                                 tcCAGACTACTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACGGCccacgccac                                                          >  1:1125812/1‑139 (MQ=255)
                                                                    aGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                       <  1:409413/139‑1 (MQ=255)
                                                                     gACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGt                                                      >  2:891999/1‑139 (MQ=255)
                                                                           tATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTaca                                                                                                                            <  2:830245/63‑1 (MQ=255)
                                                                           tATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTaca                                                                                                                            >  1:830245/1‑63 (MQ=255)
                                                                                  gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGCCGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATgtgt                                         >  1:1372630/1‑139 (MQ=255)
                                                                                  gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGccc                                                                 >  2:1388186/1‑115 (MQ=255)
                                                                                  gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGccc                                                                 <  1:1388186/115‑1 (MQ=255)
                                                                                  gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAtgtg                                          >  2:1610581/1‑138 (MQ=255)
                                                                                  gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAtgtg                                          <  1:1610581/138‑1 (MQ=255)
                                                                                  gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATgtgt                                         >  1:834045/1‑139 (MQ=255)
                                                                                     tAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCtt                                      >  1:1099903/1‑139 (MQ=255)
                                                                                             aaaaaaTTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                       <  1:264925/114‑1 (MQ=255)
                                                                                             aaaaaaTTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                       >  2:264925/1‑114 (MQ=255)
                                                                                                             tATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCg              >  1:1343931/1‑139 (MQ=255)
                                                                                                                tCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACg           <  1:675579/139‑1 (MQ=255)
                                                                                                                  tGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTa         <  2:1343931/139‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGtt                                                                                          <  2:31023/56‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGtt                                                                                          >  1:31023/1‑56 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                 <  2:609841/65‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                 >  1:609841/1‑65 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaaa                                                                       <  2:1069744/75‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaaa                                                                       >  1:1069744/1‑75 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccacg                                                         <  2:1070960/89‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccacg                                                         >  1:1070960/1‑89 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTa                                                >  1:436125/1‑98 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAt                                             <  2:424404/101‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAt                                             >  1:424404/1‑101 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCAtt                        >  1:26724/1‑122 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCAtt                        <  2:26724/122‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCg                 >  2:1011825/1‑129 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCg                 <  1:1011825/129‑1 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAAc       >  2:108578/1‑139 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAAc       >  2:923868/1‑139 (MQ=255)
                                                                                                                    gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACGGCCCACGCCACGGGTTGTGTa                                                <  2:436125/98‑1 (MQ=255)
                                                                                                                       aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAAc                                                                               >  1:734125/1‑64 (MQ=255)
                                                                                                                       aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAAc                                                                               <  2:734125/64‑1 (MQ=255)
                                                                                                                       aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAg                             >  2:495204/1‑114 (MQ=255)
                                                                                                                       aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAg                             <  1:495204/114‑1 (MQ=255)
                                                                                                                       aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGACCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGtt    <  2:1125812/139‑1 (MQ=255)
                                                                                                                         aaaaTTCACGATTTACACCAGTTCCAGCCTCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTa         >  2:1530994/1‑132 (MQ=255)
                                                                                                                         aaaaTTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTa         <  1:1530994/132‑1 (MQ=255)
                                                                                                                         aaaaTTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTa  >  2:218293/1‑139 (MQ=255)
                                                                                                                                 cGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaaa                                                                       <  1:847810/62‑1 (MQ=255)
                                                                                                                                 cGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaaa                                                                       >  2:847810/1‑62 (MQ=255)
                                                                                                                                     |                                                                                                                              
TGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTA  >  NC_000913/184124‑184383

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