Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 738,853 G→A *254* (TAG→TAA ybfD → H repeat‑associated putative transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913738,8530GA100.0% 225.8 / NA 69*254* (TAG→TAAybfDH repeat‑associated putative transposase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (34/35);  total (34/35)

TCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCG  >  NC_000913/738721‑738987
                                                                                                                                    |                                                                                                                                      
tCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                  <  1:1505140/139‑1 (MQ=255)
tCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                  <  2:17382/139‑1 (MQ=255)
      gAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTa                                                                                                                            <  2:850734/139‑1 (MQ=255)
      gAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTa                                                                                                                            <  2:339197/139‑1 (MQ=255)
           aaTTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGtt                                                                                                                       <  2:227657/139‑1 (MQ=255)
                 ataatCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTactgtctcttat                                                                                                                 <  2:1495384/139‑12 (MQ=17)
                 ataatCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAAc                                                                                                                 <  2:1459962/139‑1 (MQ=255)
                       ccAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTg                                                                                                           >  2:794512/1‑139 (MQ=255)
                          aGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACg                                                                                                        >  1:138465/1‑139 (MQ=255)
                            cATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTc                                                                                                      <  1:109449/139‑1 (MQ=255)
                              tGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAg                                                                                                    <  1:1087692/139‑1 (MQ=255)
                              tGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAg                                                                                                    <  1:820170/139‑1 (MQ=255)
                               gACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACACGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGt                                                                                                   >  1:1138156/1‑139 (MQ=255)
                                       cGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCgg                                                                                           >  2:303145/1‑139 (MQ=255)
                                                      gagTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTa                                                                                                                            <  2:935704/91‑1 (MQ=255)
                                                      gagTAACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTa                                                                                                                            >  1:935704/1‑91 (MQ=255)
                                                              gCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                    <  1:794512/139‑1 (MQ=255)
                                                               cAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTc                                                                   <  1:882452/139‑1 (MQ=255)
                                                                 gagaAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                            >  2:367403/1‑112 (MQ=255)
                                                                 gagaAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                            <  1:367403/112‑1 (MQ=255)
                                                                 gagaAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCt                                                                 >  2:637828/1‑139 (MQ=255)
                                                                  agaagaAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTg                                                                >  1:34469/1‑139 (MQ=255)
                                                                  agaagaAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGAGTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTg                                                                >  1:992395/1‑139 (MQ=255)
                                                                   gaagaaACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGc                                                               <  1:637828/139‑1 (MQ=255)
                                                                     agaaACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTa                                                             >  1:561807/1‑139 (MQ=255)
                                                                        aaCCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTacaa                                                          <  1:303145/139‑1 (MQ=255)
                                                                             gTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGa                                                     >  1:696890/1‑139 (MQ=255)
                                                                              tCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                    >  2:902653/1‑123 (MQ=255)
                                                                              tCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                    <  1:902653/123‑1 (MQ=255)
                                                                                ttCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGC‑‑TCGTAAATTTGTCCTGCTACAACAGGATTAAc                                                >  2:1415672/1‑139 (MQ=255)
                                                                                  cATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAAc                                                <  2:138465/139‑1 (MQ=255)
                                                                                         tttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                  <  1:1321782/50‑1 (MQ=255)
                                                                                         tttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAggg                                                                                                                                  >  2:1321782/1‑50 (MQ=255)
                                                                                         tttGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaa                                         >  1:197865/1‑139 (MQ=255)
                                                                                           tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAat                                       >  1:235722/1‑139 (MQ=255)
                                                                                           tGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAat                                       >  2:1490599/1‑139 (MQ=255)
                                                                                            gCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                            >  2:530615/1‑69 (MQ=255)
                                                                                            gCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                            <  1:530615/69‑1 (MQ=255)
                                                                                                  tCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                    <  1:672921/103‑1 (MQ=255)
                                                                                                  tCCCTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                    >  2:672921/1‑103 (MQ=255)
                                                                                                   cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                                        >  2:581560/1‑130 (MQ=255)
                                                                                                   cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACaaa                                        <  1:581560/130‑1 (MQ=255)
                                                                                                   cccTGATGAGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtc                               >  1:399443/1‑139 (MQ=255)
                                                                                                         tgaGCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGC‑‑TCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTc                       <  1:1415672/139‑1 (MQ=255)
                                                                                                            gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                            >  2:648576/1‑53 (MQ=255)
                                                                                                            gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCt                                                                                                            <  1:648576/53‑1 (MQ=255)
                                                                                                            gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCt                      >  2:1121523/1‑139 (MQ=255)
                                                                                                            gCTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCt                      <  2:561807/139‑1 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                            >  1:532300/1‑68 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCg                                                                                            <  2:532300/68‑1 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCg                                                                             >  2:1139667/1‑83 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                             <  2:1237115/131‑1 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTtctc                             >  1:1237115/1‑131 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTa                     >  1:931723/1‑139 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTa                     >  2:153038/1‑139 (MQ=255)
                                                                                                             cTTATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGATTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCg                                                                             <  1:1139667/83‑1 (MQ=255)
                                                                                                               tATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTaca                   <  2:34469/139‑1 (MQ=255)
                                                                                                               tATTGATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTaca                   >  1:690036/1‑139 (MQ=255)
                                                                                                                   gATTTCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTa               <  2:1138156/139‑1 (MQ=255)
                                                                                                                     tttCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAAc                                                <  1:75861/104‑1 (MQ=255)
                                                                                                                     tttCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAAc                                                >  2:75861/1‑104 (MQ=255)
                                                                                                                       tCACGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTcc           >  1:871968/1‑139 (MQ=255)
                                                                                                                          cGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                    <  1:135043/79‑1 (MQ=255)
                                                                                                                          cGTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGt                                                                    >  2:135043/1‑79 (MQ=255)
                                                                                                                           gTTTGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTa       >  2:97414/1‑139 (MQ=255)
                                                                                                                              tGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATTTCATTTCTCAACGTCTACACTTACTCCTGTAAAc    <  2:781097/139‑1 (MQ=255)
                                                                                                                              tGAATAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAAc    <  2:690036/139‑1 (MQ=255)
                                                                                                                                aaTAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCg  >  1:380211/1‑139 (MQ=255)
                                                                                                                                aaTAAAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCg  <  1:1121523/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TCCGCTGAAAGAATTAAATAATCCCAAGCATGACAGTTACGCAATTAGTGAAAAGAGTCACGGCAGAGAAGAAACCCGTCTTCATATTGTTTGCGATGTCCCTGATGAGCTTATTGATTTCACGTTTGAATAGAAAGGGCAGCTATGGTTAGAAACTACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCGTAAATTTGTCCTGCTACAACAGGATTAACTTCACAAATATCATTTCTCAACGTCTACACTTACTCCTGTAAACCG  >  NC_000913/738721‑738987

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