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% 241.1 / NA 74*271* (TAG→TAAyaeIphosphodiesterase with model substrate bis‑pNPP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/39);  total (35/39)

GTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCGACA  >  NC_000913/184125‑184392
                                                                                                                                    |                                                                                                                                       
gtgCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACAcc                                                                                                                                   <  2:574929/139‑1 (MQ=255)
        tttAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAg                                                                                                                           >  2:785073/1‑139 (MQ=255)
                ttGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTa                                                                                                                   >  2:1674017/1‑139 (MQ=255)
                 tggcctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCa                                                                                                                         <  1:579831/132‑1 (MQ=255)
                 tggcctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCa                                                                                                                         >  2:579831/1‑132 (MQ=255)
                  ggcctggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACt                                                                                                                 >  1:1601571/1‑139 (MQ=255)
                      tggccCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCg                                                                                                             <  1:640807/139‑1 (MQ=255)
                             cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCa                                                                                                      >  2:1815456/1‑139 (MQ=255)
                             cccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCa                                                                                                      >  2:787562/1‑139 (MQ=255)
                                ggCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGtt                                                                                                   <  1:52292/139‑1 (MQ=255)
                                     gACCGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGAc                                                                                              <  1:1639704/139‑1 (MQ=255)
                                            aTCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAAcc                                                                                       >  1:1526559/1‑139 (MQ=255)
                                                 aGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACa                                                                                  >  1:29041/1‑139 (MQ=255)
                                                   aaaGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACaaa                                                                                <  1:676978/139‑1 (MQ=255)
                                                        cTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGcc                                                                           <  1:1561707/139‑1 (MQ=255)
                                                         tttttCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGccc                                                                          <  1:787562/139‑1 (MQ=255)
                                                           tttCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACAcc                                                                                                                                   >  2:823328/1‑80 (MQ=255)
                                                           tttCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACAcc                                                                                                                                   <  1:823328/80‑1 (MQ=255)
                                                           tttCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccac                                                                        >  2:1607911/1‑139 (MQ=255)
                                                                tcCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccac                                                                   <  1:1607911/139‑1 (MQ=255)
                                                                tcCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccac                                                                   <  1:785073/139‑1 (MQ=255)
                                                                   aGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                <  2:835851/139‑1 (MQ=255)
                                                                   aGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                <  2:29041/139‑1 (MQ=255)
                                                                        tCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTtgtg                                                           <  2:1526559/139‑1 (MQ=255)
                                                                               tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                >  1:917203/1‑127 (MQ=255)
                                                                               tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                <  2:917203/127‑1 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTacac                                                                                                                                    >  1:1368561/1‑57 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTacac                                                                                                                                    <  2:1368561/57‑1 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                          >  1:593959/1‑99 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                          <  2:593959/99‑1 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAAc                                                                                        >  1:732034/1‑101 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAAc                                                                                        <  2:732034/101‑1 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGccc                                                                          <  2:187104/115‑1 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGccc                                                                          >  1:187104/1‑115 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTtgt                                                            <  1:538533/129‑1 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTtgt                                                            >  2:1855528/1‑129 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTtgt                                                            >  2:538533/1‑129 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTtgt                                                            <  1:1855528/129‑1 (MQ=255)
                                                                                 gCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATgtgt                                                  >  2:459404/1‑139 (MQ=255)
                                                                                         aTCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGcc                                          >  1:1509410/1‑139 (MQ=255)
                                                                                                         cTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccac                                                                   >  2:855105/1‑98 (MQ=255)
                                                                                                         cTTTATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccac                                                                   <  1:855105/98‑1 (MQ=255)
                                                                                                            tATTCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCg                       <  2:1601571/139‑1 (MQ=255)
                                                                                                               tCTGGTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACg                    <  1:377208/139‑1 (MQ=255)
                                                                                                                   gTGATAAAGTTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAtg                                                     >  2:1648288/1‑102 (MQ=255)
                                                                                                                   gTGATAAAGTTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAtg                                                     <  1:1648288/102‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGggc                                                                                                          >  2:1552343/1‑49 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGggc                                                                                                          <  1:1552343/49‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                          <  1:274175/65‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                          <  1:173012/65‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                          >  2:173012/1‑65 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCaa                                                                                          >  2:274175/1‑65 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCa                                                                                           >  1:635374/1‑64 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCa                                                                                           <  2:635374/64‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgcc                                                                     >  1:1055541/1‑86 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgcc                                                                     <  2:1055541/86‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccac                                                                        <  1:1200603/83‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccac                                                                        >  2:1200603/1‑83 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                <  2:871373/91‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                >  1:871373/1‑91 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCa                                   <  2:1876979/120‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCa                                   >  1:1876979/1‑120 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCgg                         >  2:287566/1‑130 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCgg                         <  1:287566/130‑1 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAAc                >  2:1189897/1‑139 (MQ=255)
                                                                                                                   gTGATAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAAc                >  1:635428/1‑139 (MQ=255)
                                                                                                                      aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                <  1:1657915/88‑1 (MQ=255)
                                                                                                                      aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACggg                                                                >  2:1657915/1‑88 (MQ=255)
                                                                                                                      aTAAAATTCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGtt             <  1:1189897/139‑1 (MQ=255)
                                                                                                                             tCACGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTc      >  1:1840730/1‑139 (MQ=255)
                                                                                                                                cGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACATACTGCccacgccac                                                                   <  2:329416/75‑1 (MQ=255)
                                                                                                                                cGATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccacgccac                                                                   >  1:329416/1‑75 (MQ=255)
                                                                                                                                 gATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCGACa  <  1:1674017/139‑1 (MQ=255)
                                                                                                                                 gATTTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCGACa  <  1:459404/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTTCGACA  >  NC_000913/184125‑184392

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