Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,547,591	T→C	F307L (TTC→CTC)	murP →	N‑acetylmuramic acid permease, EIIBC component, PTS system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,547,591	0	T	C	100.0%	253.0 / NA	66	F307L (TTC→CTC)	murP	N‑acetylmuramic acid permease, EIIBC component, PTS system
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (33/33);  total (33/33)

TCACCCTGCTGATCACCGCTACGCTTGCGTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGC  >  NC_000913/2547457‑2547723                                                                                                                                       |                                                                                                                                     tCACCCTGCTGATCACCGCTACGCTTGCGTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCt                                                                                                                                  >  2:199298/1‑139 (MQ=255)            aTCACCGCTACGCTTGCGTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTgg                                                                                                                       >  1:790869/1‑139 (MQ=255)                           gCGTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGcc                                                                                                                           >  2:911741/1‑120 (MQ=255)                           gCGTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGcc                                                                                                                           <  1:911741/120‑1 (MQ=255)                           gCGTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATc                                                                                                        >  1:1750483/1‑139 (MQ=255)                             gTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAg                                                                                                      <  2:1272506/139‑1 (MQ=255)                                  tGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCtt                                                                                                 >  1:879635/1‑139 (MQ=255)                                   gATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCtt                                                                                                 >  2:658366/1‑138 (MQ=255)                                   gATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCtt                                                                                                 <  1:658366/138‑1 (MQ=255)                                    aTCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTa                                                                                               <  2:1006075/139‑1 (MQ=255)                                     tCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCt                                                                                                                                  <  2:1328830/103‑1 (MQ=255)                                     tCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCt                                                                                                                                  <  2:268608/103‑1 (MQ=255)                                     tCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCt                                                                                                                                  >  1:490845/1‑103 (MQ=255)                                     tCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCt                                                                                                                                  <  2:490845/103‑1 (MQ=255)                                     tCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCt                                                                                                                                  >  1:268608/1‑103 (MQ=255)                                     tCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCAGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCt                                                                                                                                  >  1:1328830/1‑103 (MQ=255)                                         tATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGTGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTa                                                                                               >  1:724217/1‑134 (MQ=255)                                         tATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGTGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTa                                                                                               <  2:724217/134‑1 (MQ=255)                                           tGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGt                                                                                        >  1:1765314/1‑139 (MQ=255)                                                cTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTAtt                                                                                                                   <  2:1187030/107‑1 (MQ=255)                                                cTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTAtt                                                                                                                   >  1:1187030/1‑107 (MQ=255)                                                cTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTAcc                                                                                   >  1:431651/1‑139 (MQ=255)                                                cTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTAcc                                                                                   >  1:1588333/1‑139 (MQ=255)                                                cTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTAcc                                                                                   >  1:1480118/1‑139 (MQ=255)                                                 tGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCt                                                                                  <  2:879635/139‑1 (MQ=255)                                                             tATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGa                                                                      <  2:1480118/139‑1 (MQ=255)                                                                 cGAAGGTATGTCATGGCTGTTTATGCTCCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAg                                                                                                      <  1:1282907/103‑1 (MQ=255)                                                                 cGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAg                                                                                                      >  2:1282907/1‑103 (MQ=255)                                                                       tATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGAGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGa                                                            <  2:1539560/139‑1 (MQ=255)                                                                         tGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCg                                                                                                              >  1:2151/1‑87 (MQ=255)                                                                         tGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCg                                                                                                              <  2:2151/87‑1 (MQ=255)                                                                         tGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGAtt                                                          >  1:1507299/1‑139 (MQ=255)                                                                          gTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGc                                                                                                   <  1:630623/97‑1 (MQ=255)                                                                          gTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGc                                                                                                   >  2:630623/1‑97 (MQ=255)                                                                          gTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGcc                                                                 <  2:216673/131‑1 (MQ=255)                                                                          gTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGcc                                                                 >  1:216673/1‑131 (MQ=255)                                                                            cATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCaa                                                       <  2:1765314/139‑1 (MQ=255)                                                                            cATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCaa                                                       <  2:1750483/139‑1 (MQ=255)                                                                               ggCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTAt                                                                                              >  1:919171/1‑97 (MQ=255)                                                                               ggCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTAt                                                                                              <  2:919171/97‑1 (MQ=255)                                                                                gCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCg                                                                                                              <  2:1022032/80‑1 (MQ=255)                                                                                gCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCg                                                                                                              >  1:1022032/1‑80 (MQ=255)                                                                                 cTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGa                                                                      >  1:1653366/1‑119 (MQ=255)                                                                                 cTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGa                                                                      <  2:1653366/119‑1 (MQ=255)                                                                                 cTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAAc                                                      <  1:325337/135‑1 (MQ=255)                                                                                 cTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAAc                                                      >  2:325337/1‑135 (MQ=255)                                                                                 cTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCt                                                  >  2:1681843/1‑139 (MQ=255)                                                                                  tGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCtt                                                 <  2:431651/139‑1 (MQ=255)                                                                                         gCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTccc                                          <  2:300147/139‑1 (MQ=255)                                                                                              tGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCt                                     >  1:444583/1‑139 (MQ=255)                                                                                               gAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCtt                                    <  1:1252067/139‑1 (MQ=255)                                                                                                aaCAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTcg                                                                                <  1:671414/94‑1 (MQ=255)                                                                                                aaCAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTcg                                                                                >  2:671414/1‑94 (MQ=255)                                                                                                aaCAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCt                                                  <  2:7652/124‑1 (MQ=255)                                                                                                aaCAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCt                                                  >  1:7652/1‑124 (MQ=255)                                                                                                     tAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATg                              >  1:1241387/1‑139 (MQ=255)                                                                                                          cGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCg                                                                                                                             >  1:1406585/1‑39 (MQ=255)                                                                                                          cGTTCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCg                                                                                                                             <  2:1406585/39‑1 (MQ=255)                                                                                                            ttCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGgcg                       <  2:444583/139‑1 (MQ=255)                                                                                                            ttCGGTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGgcg                       <  2:1588333/139‑1 (MQ=255)                                                                                                               ggTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGgcg                       <  1:910016/136‑1 (MQ=255)                                                                                                               ggTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGgcg                       >  2:910016/1‑136 (MQ=255)                                                                                                               ggTTGTGCGGTTTTAGCCGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCgg                    >  1:1705779/1‑139 (MQ=255)                                                                                                                               ccGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCAt                                        >  1:927848/1‑103 (MQ=255)                                                                                                                               ccGGGCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCAt                                        <  2:927848/103‑1 (MQ=255)                                                                                                                                 gggCTGCTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGggcgcggcgc  <  2:1507299/139‑1 (MQ=255)                                                                                                                                   gctgctCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGa                                                                      <  1:1734329/69‑1 (MQ=255)                                                                                                                                   gctgctCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGa                                                                      >  2:1734329/1‑69 (MQ=255)                                                                                                                                       |                                                                                                                                     TCACCCTGCTGATCACCGCTACGCTTGCGTACCTGATCATTATGCCGCTGGGTGGCTGGTTATTCGAAGGTATGTCATGGCTGTTTATGCACCTGAACAGTAATCCGTTCGGTTGTGCGGTTTTAGCCGGGCTGTTCCTGATCGCCGTGGTATTTGGCGTGCATCAGGGCTTTATTCCTGTTTACCTCGCGTTAATGGACAGCCAGGGATTCAACAGCTTATTTCCCATCCTTTCAATGGCAGGCGCGGGCCAGGTGGGCGCGGCGC  >  NC_000913/2547457‑2547723

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