Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,095,484 C→T intergenic (‑47/‑21) insE1 ← / → ymdE IS3 transposase A/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,095,4840CT100.0% 214.3 / NA 67intergenic (‑47/‑21)insE1/ymdEIS3 transposase A/pseudogene
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/29);  total (38/29)

AGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAA  >  NC_000913/1095367‑1095615
                                                                                                                     |                                                                                                                                   
aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTg                                                                                                                <  2:449315/139‑1 (MQ=14)
aGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTg                                                                                                                >  2:913606/1‑139 (MQ=14)
     ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTgg                                                                                                           <  2:472866/139‑1 (MQ=25)
     ttCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTgg                                                                                                           >  2:613456/1‑139 (MQ=25)
      tCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGc                                                                                                          >  2:993060/1‑139 (MQ=37)
       cACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGca                                                                                                         >  1:550214/1‑139 (MQ=37)
            cGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACAcgc                                                                                                    >  1:541254/1‑139 (MQ=255)
                    aGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATg                                                                                            <  2:541254/139‑1 (MQ=255)
                           tGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACg                                                                                     <  1:913606/139‑1 (MQ=255)
                                 ttACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTc                                                                               <  1:440331/139‑1 (MQ=255)
                                       ggTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCaa                                                                         <  1:613456/139‑1 (MQ=255)
                                        gTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAAt                                                                        <  2:1500844/139‑1 (MQ=255)
                                         ttttttACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATg                                                                       <  1:1778708/139‑1 (MQ=255)
                                                 tGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTa                                                               >  2:1520174/1‑139 (MQ=255)
                                                  ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGtt                                                                                                 <  1:1711428/104‑1 (MQ=255)
                                                  ggTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGtt                                                                                                 >  2:1711428/1‑104 (MQ=255)
                                                  ggTTGATACTGTTTTTGGCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTAt                                                              >  2:1123941/1‑139 (MQ=255)
                                                   gTTGATCCTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGcc                                                                                                              <  1:1386675/90‑1 (MQ=255)
                                                   gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGcc                                                                                                              >  2:1386675/1‑90 (MQ=25)
                                                   gTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                             >  1:507026/1‑139 (MQ=255)
                                                      gATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTg                                                          <  2:1435832/139‑1 (MQ=255)
                                                           tGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGa                                                     <  2:1331819/139‑1 (MQ=255)
                                                             tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGcaca                                                                                                       >  1:280776/1‑87 (MQ=255)
                                                             tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGcaca                                                                                                       <  2:280776/87‑1 (MQ=255)
                                                             tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                                   >  2:944591/1‑139 (MQ=255)
                                                             tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                                   >  2:91174/1‑139 (MQ=255)
                                                             tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                                   >  1:1576369/1‑139 (MQ=255)
                                                             tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                                   >  2:1776803/1‑139 (MQ=255)
                                                             tttttGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                                   >  1:1737407/1‑139 (MQ=255)
                                                                 tGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGat                                               <  1:944591/139‑1 (MQ=255)
                                                                  gTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGata                                              >  1:855508/1‑139 (MQ=255)
                                                                  gTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGata                                              >  1:291980/1‑139 (MQ=255)
                                                                      tgtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                          <  1:91174/139‑1 (MQ=255)
                                                                       gtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGcac                                                                                                        >  2:1714742/1‑76 (MQ=255)
                                                                       gtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGcac                                                                                                        <  1:1714742/76‑1 (MQ=255)
                                                                       gtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGtt                                                                                <  2:389241/100‑1 (MQ=255)
                                                                       gtgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGtt                                                                                >  1:389241/1‑100 (MQ=255)
                                                                        tgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATg                                                                       >  1:1325419/1‑108 (MQ=255)
                                                                        tgAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATg                                                                       <  2:1325419/108‑1 (MQ=255)
                                                                            tCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGcc                                    <  1:1123941/139‑1 (MQ=255)
                                                                             cACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCa                                   <  2:1576369/139‑1 (MQ=255)
                                                                             cACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCa                                   >  1:833103/1‑139 (MQ=255)
                                                                               ccTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCaga                                 <  2:507026/139‑1 (MQ=255)
                                                                                   tgactgaGAGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCAt                             <  2:291980/139‑1 (MQ=255)
                                                                                        gagaGTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGt                        <  1:993060/139‑1 (MQ=255)
                                                                                         agagTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                       >  1:1370336/1‑139 (MQ=255)
                                                                                         agagTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                       >  2:710051/1‑139 (MQ=255)
                                                                                         agagTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                       >  1:802359/1‑139 (MQ=255)
                                                                                         agagTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                       >  2:413283/1‑139 (MQ=255)
                                                                                          gagTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTc                      >  1:466360/1‑139 (MQ=255)
                                                                                          gagTTTACTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTc                      >  1:503053/1‑139 (MQ=255)
                                                                                                 cTCACTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGaa               <  2:833103/139‑1 (MQ=255)
                                                                                                     cTTAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTa           >  2:1657741/1‑139 (MQ=255)
                                                                                                       tAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTc                                                                               >  1:236691/1‑69 (MQ=255)
                                                                                                       tAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTc                                                                               <  2:236691/69‑1 (MQ=255)
                                                                                                       tAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGat                                               <  1:1108726/101‑1 (MQ=255)
                                                                                                       tAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                       <  2:17859/125‑1 (MQ=255)
                                                                                                       tAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGtt                       >  1:17859/1‑125 (MQ=255)
                                                                                                       tAGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGAGCGATTATCTTGCCGGACTGGat                                               >  2:1108726/1‑101 (MQ=255)
                                                                                                        aGCCGCGTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACa        >  2:1490802/1‑139 (MQ=255)
                                                                                                            gcgTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCaa    >  1:828172/1‑139 (MQ=255)
                                                                                                            gcgTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCaa    >  2:1189263/1‑139 (MQ=255)
                                                                                                            gcgTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCa     <  2:1640707/138‑1 (MQ=255)
                                                                                                            gcgTGTCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCa     >  1:1640707/1‑138 (MQ=255)
                                                                                                              gtgtCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCaaaa  >  1:1463349/1‑139 (MQ=255)
                                                                                                               tgtCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACTGCGTTCAGGCAATGCGCGATTATCTTGCCTGACTGGATATCGCTTCGCCAGAGCATc                            >  1:1183316/1‑112 (MQ=255)
                                                                                                               tgtCCATTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATc                            <  2:1183316/112‑1 (MQ=255)
                                                                                                                     |                                                                                                                                   
AGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAA  >  NC_000913/1095367‑1095615

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: